1

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

*gdal_polygonize.py:*

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

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

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

是的 - 对不起，我忘了添加这些信息。

一般来说，在准备好 CSV 文件后（使用带有 -csv 选项的 gdal2xyz.py），我需要在它的开头添加一行：“Longitude,Latitude,Value”（不带引号）

然后我需要创建一个 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 和另外两个相关文件。

4

1 回答 1

0

为了存档的完整性，在 GDAL 邮件列表中也提出了这个问题，作为线程保存光栅作为点向量文件。柴坦亚似乎为它提供了解决方案。

于 2011-08-26T09:38:07.270 回答