Запись массива в растровый файл растерио создает прозрачный TIFF - PullRequest
0 голосов
/ 17 сентября 2018

Я использую 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)
...