python - 使用 rasterio 将数组保存为 Geotiff

Question

我有以下 numpy 数组:

supervised.shape
(1270, 1847)

我正在尝试使用以下代码使用rasterio将其保存到GeoTIFF:

with rasterio.open('/my/path/ReferenceRaster.tif') as src:
    ras_meta = src.profile

with rasterio.open('/my/output/path/output_supervised.tif', 'w', **ras_meta) as dst:
    dst.write(supervised)

其中 ras_meta 是:

{'driver': 'GTiff', 'dtype': 'float32', 'nodata': None, 'width': 1847, 'height': 1270, 'count': 1, 'crs': CRS.from_epsg(32736), 'transform': Affine(10.0, 0.0, 653847.1979372115,
       0.0, -10.0, 7807064.5603836905), 'tiled': False, 'interleave': 'band'}

我面临以下错误，我无法理解，因为引用栅格和我的监督数组具有相同的形状

ValueError: Source shape (1270, 1847) is inconsistent with given indexes 1

知道这里出了什么问题吗？我不完全理解该错误的含义。

Answer 1

write 需要一个形状为 (band, row, col) 的数组。您可以 reshape 数组，也可以使用write(supervised,indexes=1)。