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)
"""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()