来自XYZ csv窗体的光栅

Question

我有一个包含3列xy和z值的csv文件。我想做以下事情： 1如何将csv转换为矢量形状文件（点）。 2如何将步骤中获得的2个点转换为光栅。 三。假设一组“shp”格式的独立点文件如何将步骤2中获得的光栅值提取到新的csv文件中，以便对其执行so统计。在

问题是我不熟悉使用python进行地理处理，例如，对于步骤2，我习惯于在R中使用“graster”包轻松地完成这项工作，对于第3步，我还使用了“raster”包中的“extract”函数。然而，在python中，我可以在geopandas中执行步骤2，但是步骤2和3没有简单的答案。在

Answer 1

这个确切的例子存在于GDAL documentation。在

假设您的数据dem.csv，如下所示：

Easting,Northing,Elevation
86943.4,891957,139.13
87124.3,892075,135.01
86962.4,892321,182.04
87077.6,891995,135.01
...

创建一个VRT来描述它：

^{pr2}$

或者，如果CSV没有列，您可以如下指定GeometryField：

<GeometryField encoding="PointFromColumns" x="field_1" y="field_2" z="field_3"/>

将其另存为dem.vrt以供以后的命令使用。在

然后使用gdal_grid插值并生成GeoTIFF输出：

gdal_grid -a invdist:power=2.0:smoothing=1.0 -txe 85000 89000 -tye 894000 890000 -outsize 400 400 -of GTiff -ot Float64 -l dem dem.vrt dem.tiff