Convertir les sommets d'entités surfaciques en points à l'aide de Virtual Layer dans QGIS

Dec 06 2020

J'ai une couche de polygones et je souhaite utiliser une couche virtuelle qui crée des points à partir des sommets de la couche de polygones. La question est de savoir comment écrire une requête qui fait cela?

Je connais l'outil "Extract vertices", mais ce n'est pas ce que je recherche.

Ce que j'ai essayé:

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

nodes_to_points()est une expression QGIS, mais la documentation dit :

Les fonctions des expressions QGIS peuvent également être utilisées dans une requête de couche virtuelle.

Cependant, il ne dit rien comment l'inclure dans la syntaxe de requête (SQL). Mise à jour: il semble que toutes les expressions QGIS ne puissent pas être utilisées dans la couche virtuelle, voir ici pour une liste des fonctions pouvant être utilisées: Liste de toutes les fonctions disponibles dans la couche virtuelle de QGIS

J'ai également regardé dans la liste de référence des fonctions SQL de SpatiaLite 4.2.0 , mais je n'ai trouvé aucune fonction pour créer des points individuels pour chaque sommet. Le plus proche que je suis venu utilisait ST_DissolvePoints():
```
 select fid, ST_DissolvePoints (p.geometry) as vertices
 from polygon as p
```

Cela crée l'effet visuel des sommets, mais pas en tant que points individuels, mais en une seule entité multipoint par polygone:

Réponses

2 J.Monticolo Dec 18 2020 at 06:01

Voici ma solution:

Requirements

QGIS 3. x
Une couche (multi) -polygone avec un idchamp avec des valeurs distinctes

Features

Une requête
Fonctionne avec des multipolygones
Configurer le nom de la couche à un seul endroit
Sortie du calque id, du numéro de pièce, du numéro de l'anneau et du numéro de sommet
Déposer le dernier sommet dupliqué (dans un polygone, début du sommet = fin du sommet)

Ouvrez QGIS, allez à Database > Database > DB Manager... > Virtual Layers
Ouvrez un nouveau SQL Window
Copiez le script SQL ci-dessous et remplacez, dans la lyrpièce, polygonspar votre vrai nom de couche (ou renommez sous QGIS votre couche 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

Malheureusement, les fonctions ST_Dump()et ST_DumpPoints()ne sont pas encore implémentées dans Virtual Layer. Par conséquent, je propose une solution de contournement sur la façon d'extraire des sommets d'entités polygonales.

Cette solution intègre plusieurs fonctions, à savoir NumInteriorRings(), ST_InteriorRingN(), ST_ExteriorRing(), ST_PointN(), et ST_NPoints().

Quelqu'un peut trouver que cette solution n'est pas parfaite en termes de performances et de complexité, en raison de la vaste sortie de la CROSS JOINclause qui produit des index redondants (pour les anneaux intérieurs ainsi que pour les sommets) qui sont réellement nécessaires. Si quelqu'un sait comment le démêler, n'hésitez pas à améliorer la requête ou à donner un indice.