我正在使用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和另外两个相关文件。
答案 0 :(得分:0)
为了存档完整性,此问题也在GDAL邮件列表中被询问为线程save raster as point-vector file。 Chaitanya似乎为它提供了解决方案。