Мой вопрос об открытии файлов TIFF с растерио.Я пытаюсь выполнить расчеты NDVI с помощью снимков Landsat 8 из USGS Earth Explorer.Я выполнил расчет с использованием растерио, а также с помощью программы Qgis.
После того, как я выполнил расчет (мин. Ndvi: -0.03, макс. Ndvi: 0.60), сохранил результат на диск.
Затем я читаю tiff с растерио:
ndvi_dir = "/Users/username/Downloads/ndvi_qgis.tif"
with rasterio.open(ndvi_dir) as src:
ndvi_qgis = src.read()
plt.imshow(ndvi_qgis[0])
Полученное изображение здесь:
tiff, открытое растерио
ЕслиЯ хочу посмотреть на профиль растра:
src.profile
{'driver': 'GTiff', 'dtype': 'float32', 'nodata': -3.4028234663852886e+38, 'width': 7801, 'height': 7901, 'count': 1, 'crs': CRS({'init': 'epsg:32630'}), 'transform': Affine(30.0, 0.0, 512385.0,
0.0, -30.0, 5846415.0), 'tiled': False}
Когда я проверил минимальное и максимальное числа этого изображения, min: -3.4028234663852886e + 38, max: 1.0
фактическое изображение, которое я хочу увидеть:
nvdi - результаты из qgis
Я думаю, что мне не хватает чего-то связанного с типами данных при чтении результата qgis tiff с библиотекой растерио.Любая помощь приветствуется.