Converta vértices de feições de polígono em pontos usando a Camada Virtual no QGIS
Eu tenho uma camada poligonal e quero usar uma camada virtual que cria pontos a partir dos vértices da camada poligonal. A questão é como escrever uma consulta que faça isso?
Conheço a ferramenta "Extract vertices"
, mas não é isso que procuro.
O que eu tentei:
select fid, nodes_to_points(p.geometry) as vertices from polygon as p
nodes_to_points()
é uma expressão QGIS, mas a documentação diz :
Funções de expressões QGIS também podem ser usadas em uma consulta de camada virtual.
No entanto, não diz nada como incluí-lo na sintaxe de consulta (SQL). Atualização: parece que nem todas as expressões QGIS podem ser usadas na camada virtual, veja aqui uma lista de funções que podem ser usadas: Listando todas as funções disponíveis na Camada Virtual QGIS
Também procurei na lista de referência de funções SQL do SpatiaLite 4.2.0 , mas não consegui encontrar nenhuma função para criar pontos individuais para cada vértice. O mais próximo que cheguei foi usando
ST_DissolvePoints()
:select fid, ST_DissolvePoints (p.geometry) as vertices from polygon as p
Isso cria o efeito visual de vértices, mas não como pontos individuais, mas um único recurso multiponto por polígono:

Respostas
Aqui está minha solução:
Requirements
- QGIS 3. x
- Uma camada de (multi) polígonos com um
id
campo com valores distintos
Features
- Uma consulta
- Funciona com multipolígonos
- Configure o nome da camada em apenas um lugar
- Exibe a camada
id
, número da peça, número do anel e número do vértice - Solte o último vértice duplicado (em um polígono, início do vértice = final do vértice)
- Abrir QGIS, ir para a Database > Database > DB Manager... > Virtual Layers
- Abra um novo SQL Window
- Copie o script SQL abaixo e substitua, na lyrparte,
polygons
pelo nome real da sua camada (ou renomeie sob QGIS sua camadapolygons
):
-- 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
Infelizmente, as funções ST_Dump()e ST_DumpPoints()ainda não foram implementadas na Camada Virtual. Portanto, estou oferecendo uma solução alternativa sobre como extrair vértices de recursos de polígono.
Esta solução incorpora várias funções, nomeadamente NumInteriorRings(), ST_InteriorRingN(), ST_ExteriorRing(), ST_PointN(), e ST_NPoints().
Alguém pode achar que esta solução não é perfeita em termos de desempenho e complexidade, por causa da vasta saída da CROSS JOINcláusula que produz índices redundantes (para anéis internos, bem como para vértices) que são realmente necessários. Se alguém souber como desvendar, sinta-se à vontade para melhorar a consulta ou dar uma dica.
Vamos supor que existe uma camada de polígono chamada 'polygons'
, veja a imagem abaixo.

Com a seguinte consulta
-- generating the second series required to extract all vertices from polygons' contours
WITH RECURSIVE generate_series2(vert) AS (
SELECT 1
UNION ALL
SELECT vert + 1
FROM generate_series2, max_num_points_per_feature
WHERE vert + 1 <= stop2
),
-- finding max number of vertices in all features
max_num_points_per_feature AS (
SELECT MAX(st_npoints(geom1)) AS stop2
FROM inter_outer_contours
),
-- union of geometries of both outer and inner contours, represented as polylines
inter_outer_contours AS (
-- interior rings from all polygons
WITH interior AS (
-- generating the first series required to extract all interior rings
WITH RECURSIVE generate_series1(ring) AS (
SELECT 1
UNION ALL
SELECT ring + 1
FROM generate_series1, max_num_rings_per_polys
WHERE ring + 1 <= stop1
),
-- finding max number of interior rings within all polygons
max_num_rings_per_polys AS (
SELECT MAX(numrings) AS stop1
FROM num_rings_per_poly
),
-- finding how many interior rings each polygon has
num_rings_per_poly AS (
SELECT id, NumInteriorRings(geometry) AS numrings
FROM "polygons"
GROUP BY id
)
-- query to extract all interior rings from all polygons
SELECT p.id AS origin, -- a field represents the original polygon id
'interior' AS info, -- a text field represents the interior attribute
interior_ring_n(geometry, ring) AS geom1 -- setting geometry for a polyline
FROM "polygons" AS p
CROSS JOIN generate_series1
WHERE geom1 IS NOT NULL -- no null geometries accepted, needed because of the cross join
),
-- exterior rings from all polygons
exterior AS (
SELECT p.id AS origin, -- a field represents the original polygon id
'exterior' AS info, -- a text field represents the exterior attribute
st_exteriorring(geometry) AS geom1 -- setting geometry for a polyline
FROM "polygons" AS p
)
-- a union between interior and exterior rings from all polygons
SELECT *
FROM interior
UNION
SELECT *
FROM exterior
)
-- query to achieve all vertices from all polygons
SELECT ioc.origin,
ioc.info,
gs2.vert, -- a field represents the vertice index
st_pointn(ioc.geom1, gs2.vert) AS geom2 -- setting geometry for a point
FROM inter_outer_contours AS ioc
JOIN generate_series2 AS gs2
WHERE geom2 IS NOT NULL -- no null geometries accepted, needed because of the cross join
é possível obter todos os vértices

O resultado corresponde à saída do geoalgoritmo "Extrair vértices" .
Explicações:
- Na consulta acima , não se esqueça de alterar a
"polygons"
string de quatro vezes para o nome da camada poligonal real
FAÇAM:
- estenda a seção de explicação
- tente com multipolígonos
- adicione um novo campo
"id"
Referências:
- Usando RECURSIVE na Camada Virtual