def test_no_squeeze_band(shared_datadir):
f = Path(shared_datadir) / TEST_FILE
g = GeoTiff(f)
assert g.da.band.size == 3
"""An example of using the GeoTiff data component.
The example is adapted from http://xarray.pydata.org/en/stable/examples/visualization_gallery.html#imshow()-and-rasterio-map-projections,
and it uses the sample GeoTIFF image "RGB.byte.tif" from the rasterio project.
"""
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
from bmi_geotiff import GeoTiff
tif_file = "RGB.byte.tif"
g = GeoTiff(tif_file)
print(g.da)
crs = ccrs.UTM("18N")
ax = plt.subplot(projection=crs)
g.da.plot.imshow(ax=ax, rgb="band", transform=crs)
plt.show()
def test_open(shared_datadir):
g = GeoTiff()
f = Path(shared_datadir) / TEST_FILE
g.open(f)
assert g.filename == f
assert isinstance(g.da, DataArray)
def test_squeeze_band():
g = GeoTiff(TEST_URL)
assert g.da.ndim != 3
def test_with_filename_as_url():
g = GeoTiff(TEST_URL)
assert isinstance(g, GeoTiff)
assert g.filename == TEST_URL
assert isinstance(g.da, DataArray)
def test_with_filename_as_keyword(shared_datadir):
f = Path(shared_datadir) / TEST_FILE
g = GeoTiff(filename=f)
assert isinstance(g, GeoTiff)
assert g.filename == f
assert isinstance(g.da, DataArray)
def test_bad_filename(tmpdir, shared_datadir):
with pytest.raises(RasterioIOError):
GeoTiff("foo.tif")
def test_instantiate_without_filename():
g = GeoTiff()
assert isinstance(g, GeoTiff)
assert g.filename is None
assert g.da is None
"""An example of using the GeoTiff data component. The data used are band 1 of Landsat 7 ETM scene LE71030642021102ASA00. """ import matplotlib.pyplot as plt from bmi_geotiff import GeoTiff url = "https://csdms.colorado.edu/data/LE07_L1GT_103064_20210412_20210412_02_RT_B1.TIF" g = GeoTiff() g.open(url) print(g.da) g.da.plot.imshow() plt.show()