python - 如何使用光谱 python 处理多光谱栅格文件？

Question

我有兴趣使用Spectral Python (SPy) 来可视化和分类多波段栅格 GeoTIFF（不是高光谱数据）。目前看来，只有.lan,.gis文件格式是可读的。

我尝试将文件转换为.lanwithgdal_translate但不支持图像格式（IOError: Unable to determine file type or type not supported）。

知道如何将此库用于非超光谱数据集吗？

Answer 1

将 GeoTIFF 文件转换为兼容格式（例如 LAN）。这可以通过以下两种方式之一来完成。从系统外壳，使用gdal_translate：

gdal_translate -of LAN file.tif file.lan

或在 Python 中类似：

from osgeo import gdal

src_fname = 'file.tif'
dst_fname = 'file.lan'
driver = gdal.GetDriverByName('LAN')

sds = gdal.Open(src_fname)
dst = driver.CreateCopy(dst_fname, sds)
dst = None  # close dataset; the file can now be used by other processes

请注意，第一种方法实际上更好，因为它还传输其他元数据，例如空间参考系统和可能的其他数据。要在 Python 中正确执行相同的操作，需要添加更多代码行。