Convierta vértices de entidades poligonales en puntos usando Virtual Layer en QGIS

Dec 06 2020

Tengo una capa de polígono y quiero usar una capa virtual que crea puntos a partir de los vértices de la capa de polígono. La pregunta es ¿cómo escribir una consulta que haga esto?

Soy consciente de la herramienta "Extract vertices", pero eso no es lo que estoy buscando.

Lo que probé:

 select fid, nodes_to_points(p.geometry) as vertices
 from polygon as p

nodes_to_points()es una expresión QGIS, pero la documentación dice :

Las funciones de las expresiones QGIS también se pueden utilizar en una consulta de capa virtual.

Sin embargo, no dice nada sobre cómo incluirlo en la sintaxis de consulta (SQL). Actualización: parece que no todas las expresiones de QGIS se pueden usar en la capa virtual, consulte aquí una lista de funciones que se pueden usar: Lista de todas las funciones disponibles en la capa virtual de QGIS

También busqué en la lista de referencia de funciones SQL de SpatiaLite 4.2.0 , pero no pude encontrar ninguna función para crear puntos individuales para cada vértice. Lo más cerca que vine fue usando ST_DissolvePoints():
```
 select fid, ST_DissolvePoints (p.geometry) as vertices
 from polygon as p
```

Esto crea el efecto visual de vértices, pero no como puntos individuales, sino como una única entidad multipunto por polígono:

Respuestas

2 J.Monticolo Dec 18 2020 at 06:01

Aquí mi solución:

Requirements

QGIS 3. x
Una capa (multi) poligonal con un idcampo con valores distintos

Features

Una consulta
Funciona con multipolígonos
Configure el nombre de la capa en un solo lugar
Muestra la capa id, el número de pieza, el número de anillo y el número de vértice
Suelta el último vértice duplicado (en un polígono, inicio del vértice = final del vértice)

QGIS abierta, vaya a Database > Database > DB Manager... > Virtual Layers
Abra una nueva SQL Window
Copie el script SQL a continuación y reemplace, en la lyrparte, polygonscon el nombre de su capa real (o cambie el nombre en QGIS a su capa polygons):

-- list parts
WITH RECURSIVE gs_part(id, part) AS (
    SELECT conf_p.id, conf_p.start
    FROM conf_p
    UNION ALL
    SELECT conf_p.id, part + 1
    FROM gs_part gs, conf_p
    WHERE
        gs.id = conf_p.id
        AND gs.part + 1 <= conf_p.stop
),
-- list interior rings
gs_ring(id, part, ring) AS (
    SELECT conf_i.id, conf_i.part, conf_i.start
    FROM conf_i
    UNION ALL
    SELECT conf_i.id, conf_i.part, ring + 1
    FROM gs_ring gs, conf_i
    WHERE
        gs.id = conf_i.id
        AND gs.part = conf_i.part
        AND gs.ring + 1 <= conf_i.stop
),
-- list vertices
gs_vert(id, part, ring, vert) AS (
    SELECT conf_v.id, conf_v.part, conf_v.ring, conf_v.start
    FROM conf_v
    UNION ALL
    SELECT conf_v.id, conf_v.part, conf_v.ring, vert + 1
    FROM gs_vert gs, conf_v
    WHERE
        gs.id = conf_v.id
        AND gs.part = conf_v.part
        AND gs.ring = conf_v.ring
        AND gs.vert + 1 < conf_v.stop
),
--
parts AS (
    SELECT lyr.id, gs_part.part, st_geometryn(lyr.geometry, gs_part.part) AS geometry
    FROM gs_part, lyr
    WHERE lyr.id = gs_part.id
),
--
rings AS (
    SELECT
        'interior' AS info,
        parts.id,
        parts.part,
        gs_ring.ring,
        interior_ring_n(parts.geometry, ring) AS geometry
    FROM gs_ring, parts
    WHERE
        gs_ring.ring > 0
        AND gs_ring.id = parts.id
        AND gs_ring.part = parts.part
    UNION ALL
    SELECT
        'exterior',
        parts.id,
        parts.part,
        0,
        exterior_ring(parts.geometry)
    FROM parts
),
-- configuration
-- for parts
conf_p AS (
    SELECT id, 1 AS start, st_numgeometries(geometry) AS stop
    FROM lyr
),
-- for interior rings
conf_i AS (
    SELECT id, part, 0 AS start, num_interior_rings(geometry) AS stop
    FROM parts
),
-- for vertices
conf_v AS (
    SELECT id, part, ring, 1 AS start, st_npoints(geometry) AS stop
    FROM rings
),
-- for layer
lyr AS (
    SELECT
        'polygons' AS lyr_name, -- ## replace here with the 'layer name' ##
        p.*
    FROM polygons p --  ## replace here with the layer name ##
),
-- get layer crs
crs AS (
    SELECT CAST(
        SUBSTR(layer_property(lyr_name, 'crs'), 6, 15)
        AS integer
    ) AS id
    FROM lyr
    LIMIT 1
)

SELECT
    rings.info,
    rings.id,
    rings.part,
    rings.ring,
    gs_vert.vert AS vertex,
    SetSRID(st_pointn(rings.geometry, gs_vert.vert), crs.id) AS geometry
FROM gs_vert, rings, crs
WHERE
    gs_vert.id = rings.id
    AND gs_vert.part = rings.part
    AND gs_vert.ring = rings.ring

7 Taras Dec 07 2020 at 04:16

Desafortunadamente, las funciones ST_Dump()y ST_DumpPoints()aún no están implementadas en Virtual Layer. Por lo tanto, ofrezco una solución alternativa sobre cómo extraer vértices de entidades poligonales.

Esta solución incorpora varias funciones, a saber NumInteriorRings(), ST_InteriorRingN(), ST_ExteriorRing(), ST_PointN(), y ST_NPoints().

Alguien puede encontrar que esta solución no es perfecta en términos de rendimiento y complejidad, debido a la gran cantidad de resultados de la CROSS JOINcláusula que produce índices redundantes (tanto para anillos interiores como para vértices) que son realmente necesarios. Si alguien sabe cómo desentrañarlo, no dude en mejorar la consulta o dar una pista.