【发布时间】:2011-10-23 01:46:19
【问题描述】:
我正在使用 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 文件后(使用带有 -csv 选项的 gdal2xyz.py), 我需要在它的开头添加一行: “经度、纬度、值”(不带引号)
然后我需要创建一个 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 和另外两个相关文件。
【问题讨论】: