我正在使用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 和其他两个相关文件。
I am using QGIS software. I would like to show value of each raster cell as label.
My idea (I don't know any plugin or any functionality from QGIS which allow to it easier) is to export raster using gdal2xyz.py into coordinates-value format and then save it as vector (GML or shapefile). For this second task, I try to use
*gdal_polygonize.py:*
gdal_polygonize.py rainfXYZ.txt rainf.shp Creating output rainf.shp of
format GML.
0...10...20...30...40...50...60...70...80...90...100 - done.
unfortunately I am unable to load created file (even if I change the extension to .gml)
ogr2ogr tool don't even recognize this format.
yes - sorry I forgot to add such information.
In general after preparing CSV file (using gdal2xyz.py with -csv option),
I need to add one line at begining of it:
"Longitude,Latitude,Value" (without the quotes)
Then I need to create a VRT file which contain
*> <OGRVRTDataSource>
> <OGRVRTLayer name="Shapefile_name">
> <SrcDataSource>Shapefile_name.csv</SrcDataSource>
> <GeometryType>wkbPoint</GeometryType>
>
> <GeometryField encoding="PointFromColumns" x="Longitude"
> y="Latitude"/>
> </OGRVRTLayer> </OGRVRTDataSource>*
Run the command "ogr2ogr -select Value Shapefile_name.shp Shapefile_name.vrt". I got the file evap_OBC.shp and two other associated files.
发布评论
评论（1）
为了存档的完整性，这个问题也在 GDAL 邮件列表上作为线程 将光栅保存为点矢量文件。看来Chaitanya为此提供了解决方案。
For the sake of archive completeness, this question has also been asked on GDAL mailing list as thread save raster as point-vector file. It seems Chaitanya provided solution for it.