主页 / coding / 问题 / 78567751
Accepted
Gery
Asked: 2024-06-03 05:35:54 +0800 CST

对相交点的值求和并将结果保存在多边形列中

  • 772

使用此链接我可以计算与多边形相交的所有点：

ogrinfo -sql "select count(*) from tmppol, tmppoi where st_intersects(tmppol.geometry,tmppoi.geometry) group by tmppol.fid" -dialect sqlite input.vrt

我的疑问是如何修改它以获得所有点值的总和并将结果放入多边形形状文件的列中？也许作为更新调用，例如：

ogrinfo tmppol.shp -dialect sqlite -sql "update tmppol set sum = sum(tmpoi.tn) from tmppol where st_intersects(tmppol.geom, tmppoi.geom) ???"

非常感谢任何提示，

MWE 文件可以从这里下载。

sql

1 个回答

  1. Best Answer

    计算每个多边形的点值总和

    ogrinfo -sql "SELECT tmppol.fid, SUM(tmppoi.value) as point_sum FROM tmppol, tmppoi WHERE ST_Intersects(tmppol.geometry, tmppoi.geometry) GROUP BY tmppol.fid" -dialect sqlite input.vrt

    创建一个名为 input.vrt 的文件

    <OGRVRTDataSource>
    <OGRVRTLayer name="tmppol">
        <SrcDataSource>tmppol.shp</SrcDataSource>
        <GeometryType>wkbPolygon</GeometryType>
        <LayerSRS>WGS84</LayerSRS>
    </OGRVRTLayer>
    <OGRVRTLayer name="tmppoi">
        <SrcDataSource>tmppoi.shp</SrcDataSource>
        <GeometryType>wkbPoint</GeometryType>
        <LayerSRS>WGS84</LayerSRS>
    </OGRVRTLayer>
</OGRVRTDataSource>

    使用：

    ogr2ogr -f "ESRI Shapefile" tmppol.shp -dialect sqlite -sql "UPDATE tmppol SET sum = (SELECT SUM(tmppoi.value) FROM tmppoi WHERE ST_Intersects(tmppol.geometry, tmppoi.geometry))" input.vrt
    • 1

相关问题