Question

我正在使用QGIS软件。我想将每个栅格单元的值显示为标签。我的想法（我不知道任何插件或QGIS的任何功能允许更容易）是使用gdal2xyz.py将栅格导出为坐标值格式，然后将其保存为矢量（GML或shapefile）。对于第二项任务，我尝试使用

* gdal_polygonize.py：*

gdal_polygonize.py rainfXYZ.txt rainf.shp创建输出rainf.shp 格式化GML。 0 ... 10 ... 20 ... 30 ... 40 ... 50 ... 60 ... 70 ... 80 ... 90 ... 100 - 完成。

不幸的是我无法加载创建的文件（即使我将扩展名更改为.gml）

ogr2ogr 工具甚至无法识别此格式。

是的 - 抱歉，我忘记添加此类信息。

一般在准备好CSV文件后（使用gdal2xyz.py和-csv选项），我需要在开头添加一行： “经度，纬度，价值”（不含引号）

然后我需要创建一个包含

的VRT文件

*> <OGRVRTDataSource>
>     <OGRVRTLayer name="Shapefile_name">
>         <SrcDataSource>Shapefile_name.csv</SrcDataSource>
>         <GeometryType>wkbPoint</GeometryType>
> 
>         <GeometryField encoding="PointFromColumns" x="Longitude"
> y="Latitude"/>
>     </OGRVRTLayer> </OGRVRTDataSource>*

运行命令“ogr2ogr -select Value Shapefile_name.shp Shapefile_name.vrt”。我得到了文件evap_OBC.shp和另外两个相关文件。

Answer 1

为了存档完整性，此问题也在GDAL邮件列表中被询问为线程save raster as point-vector file。 Chaitanya似乎为它提供了解决方案。

将文件以XYZ格式保存为矢量（GML或shp）

1 个答案: