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

Question

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

我尝试用.lan将文件转换为gdal_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中正确执行此操作，需要添加更多代码行。