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
|
"""
Tests of band mask creation, both .msk sidecar and internal.
"""
import numpy as np
import pytest
import rasterio
from rasterio.enums import MaskFlags
from rasterio.errors import RasterioIOError
from rasterio.windows import Window
from .conftest import gdal_version
def test_create_internal_mask(data):
"""Write an internal mask to the fixture's RGB.byte.tif."""
with rasterio.Env(GDAL_TIFF_INTERNAL_MASK=True):
with rasterio.open(str(data.join('RGB.byte.tif')), 'r+') as dst:
blue = dst.read(1, masked=False)
mask = 255 * (blue == 0).astype('uint8')
dst.write_mask(mask)
# There should be no .msk file
assert data.join('RGB.byte.tif').exists()
assert not data.join('RGB.byte.tif.msk').exists()
# Check that the mask was saved correctly.
with rasterio.open(str(data.join('RGB.byte.tif'))) as src:
assert (mask == src.read_masks()).all()
for flags in src.mask_flag_enums:
assert MaskFlags.per_dataset in flags
assert MaskFlags.alpha not in flags
assert MaskFlags.nodata not in flags
@pytest.mark.xfail(gdal_version.at_least("3.9"), reason="Internal mask are the default since 3.9.0.")
def test_create_sidecar_mask(data):
"""Write a .msk sidecar mask."""
with rasterio.open(str(data.join('RGB.byte.tif')), 'r+') as dst:
blue = dst.read(1, masked=False)
mask = 255 * (blue == 0).astype('uint8')
dst.write_mask(mask)
# There should be a .msk file in this case.
assert data.join('RGB.byte.tif').exists()
assert data.join('RGB.byte.tif.msk').exists()
# Check that the mask was saved correctly.
with rasterio.open(str(data.join('RGB.byte.tif'))) as src:
assert (mask == src.read_masks()).all()
for flags in src.mask_flag_enums:
assert MaskFlags.per_dataset in flags
assert MaskFlags.alpha not in flags
assert MaskFlags.nodata not in flags
# Check the .msk file, too.
with rasterio.open(str(data.join('RGB.byte.tif.msk'))) as msk:
assert (mask == msk.read(1, masked=False)).all()
def test_create_mask_windowed_sidecar(data):
"""Writing masks by window succeeds with sidecar mask
"""
with rasterio.Env(GDAL_TIFF_INTERNAL_MASK=False):
with rasterio.open(str(data.join('RGB.byte.tif')), 'r+') as dst:
for ij, window in dst.block_windows():
blue = dst.read(1, window=window, masked=False)
mask = 255 * (blue == 0).astype('uint8')
dst.write_mask(mask, window=window)
def test_create_mask_windowed_internal(data):
"""Writing masks by window with internal mask
"""
with rasterio.Env(GDAL_TIFF_INTERNAL_MASK=True):
with rasterio.open(str(data.join('RGB.byte.tif')), 'r+') as dst:
for ij, window in dst.block_windows():
blue = dst.read(1, window=window, masked=False)
mask = (blue != 0)
dst.write_mask(mask, window=window)
def test_create_mask_windowed_internal_spillover(data):
"""Writing mask to a window that spills over dataset raises an error
"""
with rasterio.Env(GDAL_TIFF_INTERNAL_MASK=True):
with rasterio.open(str(data.join('RGB.byte.tif')), 'r+') as dst:
mask = np.ones((1024, 1024), dtype="bool")
with pytest.raises(RasterioIOError):
dst.write_mask(mask, window=Window(800, 800, 1024, 1024))
|
