Я использую Rasterio для чтения GeoTIFF из Sentinel 2 и вычисления NDVI, но когда я пытаюсь сохранить его как новый TIF-файл, файл становится прозрачным изображением.Чего мне не хватает?
Это код:
class MapService:
def __init__(self):
dsRed = rasterio.open('data/bands/B04.TIF')
bandRed = dsRed.read(1)
dsNIR = rasterio.open('data/bands/B08.TIF')
bandNIR = dsNIR.read(1)
ndvi = np.zeros(dsRed.shape, dtype=rasterio.float32)
ndvi = (bandNIR.astype(float)-bandRed.astype(float))/(bandNIR+bandRed)
kwargs = dsRed.meta
kwargs.update(
dtype=rasterio.float32,
count=1,
compress='lzw')
self.ndvi = ndvi
with rasterio.open('result/example-total.tif', 'w', **kwargs) as dst:
dst.write(ndvi.astype(rasterio.float32),1)