python - 将大彩色图像保存为 `GTiff` 和 `gdal`

Question

我正在尝试保存尺寸为 (15000, 80000, 3) 的大图像。这个数组是一个 numpy 数组，我将其初始化为 im_final = np.zeros((15000,80000,,3))。为了保存，我使用 gdal 如下:

dst_ds = gdal.GetDriverByName('GTiff').Create('val.tif', 80000, 15000, 3, gdal.GDT_Byte)
dst_ds.GetRasterBand(1).WriteArray(im_final[:,:,0])   # write r-band to the    raster
dst_ds.GetRasterBand(2).WriteArray(im_final[:,:,1])   # write g-band to the raster
dst_ds.GetRasterBand(3).WriteArray(im_final[:,:,2])   # write b-band to the raster
dst_ds.FlushCache()                     # write to disk
dst_ds = None

当我保存它时，生成的图像是黑白的。但是，我需要图像为 RGB，有人知道问题出在哪里吗？此外，im_final 中的值为 uint16。

Answer 1

问题是您正在尝试将 uint16 写入 uint8 (gdal.GDT_Byte) 图像。如果您确实需要 8 位图像(例如，如果您想在非 GIS 程序中查看此图像)，最佳做法是将 im_final 缩放到 0-255 之间。这可以是从 0-65535 到 0-255 的映射，或者每个波段的最小值/最大值到 0-255 或任何其他数量的方式。

如果 im_final 中的值很重要，则在 driver.Create() 中使用 gdal.GDT_UInt16。