1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
|
import subprocess
import uuid
import numpy as np
import pytest
import rasterio
@pytest.fixture(scope='function')
def tempfile():
"""A temporary filename in the GDAL '/vsimem' filesystem"""
return f"/vsimem/{uuid.uuid4()}"
def complex_image(height, width, dtype):
"""An array with sequential elements"""
return np.array(
[complex(x, x) for x in range(height * width)],
dtype=dtype).reshape(height, width)
@pytest.mark.parametrize("dtype", ["complex", "complex64", "complex128"])
@pytest.mark.parametrize("height,width", [(20, 20)])
def test_read_array(tempfile, dtype, height, width):
"""_io functions read and write arrays correctly"""
in_img = complex_image(height, width, dtype)
with rasterio.open(tempfile, 'w+', driver='GTiff', dtype=dtype,
height=height, width=width, count=1) as dataset:
dataset.write(in_img, 1)
out_img = dataset.read(1)
assert (in_img == out_img).all()
def test_complex_nodata(tmpdir):
"""A complex dataset can be created with a real nodata value"""
import numpy as np
import rasterio
from rasterio.transform import Affine
x = np.linspace(-4.0, 4.0, 240)
y = np.linspace(-3.0, 3.0, 180)
X, Y = np.meshgrid(x, y)
Z1 = np.ones_like(X) + 1j
res = (x[-1] - x[0]) / 240.0
transform1 = Affine.translation(x[0] - res / 2, y[-1] - res / 2) * Affine.scale(res, -res)
tempfile = str(tmpdir.join("test.tif"))
with rasterio.open(tempfile, 'w', driver='GTiff', height=Z1.shape[0], width=Z1.shape[1], nodata=0, count=1, dtype=Z1.dtype, crs='+proj=latlong', transform=transform1) as dst:
dst.write(Z1, 1)
with rasterio.open(tempfile) as dst:
assert dst.nodata == 0
@pytest.mark.gdalbin
def test_complex_int16(tmpdir):
"""A cint16 dataset can be created"""
import numpy as np
import rasterio
from rasterio.transform import Affine
x = np.linspace(-4.0, 4.0, 240)
y = np.linspace(-3.0, 3.0, 180)
X, Y = np.meshgrid(x, y)
Z1 = np.ones_like(X) + 1j
res = (x[-1] - x[0]) / 240.0
transform1 = Affine.translation(x[0] - res / 2, y[-1] - res / 2) * Affine.scale(
res, -res
)
tempfile = str(tmpdir.join("test.tif"))
with rasterio.open(
tempfile,
"w",
driver="GTiff",
height=Z1.shape[0],
width=Z1.shape[1],
nodata=None,
count=1,
dtype="complex_int16",
crs="+proj=latlong",
transform=transform1,
) as dst:
dst.write(Z1, 1)
assert "Type=CInt16" in subprocess.check_output(["gdalinfo", tempfile]).decode(
"utf-8"
)
with rasterio.open(tempfile) as dst:
assert dst.nodatavals == (None,)
data = dst.read()
assert data.dtype == np.complex64
|
