Início / coding / Perguntas / 78567751
Accepted
Gery
Asked: 2024-06-03 05:35:54 +0800 CST

Somar valores de pontos de interseção e salvar o resultado em uma coluna poligonal

  • 772

Usando este link eu poderia contar todos os pontos que cruzam polígonos:

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

Minha dúvida é como ele pode ser modificado para obter a soma de todos os valores dos pontos e colocar o resultado em uma coluna do shapefile do polígono? Talvez como uma chamada de atualização, como:

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

Qualquer sugestão é muito apreciada,

Os arquivos MWE podem ser baixados aqui .

sql

1 respostas

  1. Best Answer

    calcular a soma dos valores dos pontos para cada polígono

    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

    crie um arquivo chamado 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>

    usar:

    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

relate perguntas