使用python进行光栅到点转换

Question

我想使用 python 脚本根据光栅中的值字段将多个光栅文件转换为点文件。但是，我想排除具有“0”值的单元格。有没有办法一步完成，还是我需要将原始栅格复制到新栅格并将 no_data 值设置为 0？有什么建议吗？谢谢。

Answer 1

免责声明：我是以下库的作者。

你可以试试lidario.Translator：它允许快速将多个 .tif 文件转换为任何类型的点云（np.array、pd.dataframe、csv 文件，...）。它在引擎盖下与 GDAL 一起工作。

Geotiff 到 np.array：

import lidario as lio
    
# Create a Translator object which will take a tif file and return a np.array
translator = lio.Translator("geotiff", "np")
    
# Translate the tif file and get the np.array
point_cloud = translator.translate("/path/to/file.tif", no_data=0, band=1)

直接保存到 CSV：

import lidario as lio

translator = lio.Translator("geotiff", "csv")
translator.translate("/path/to/file.tif", out_file='output.csv', no_data=0, band=1)