Selezionare vertici duplicati↵

Ogni volta che estraiamo i vertici in una geometria (lineare o poligonale) otteniamo quasi sempre dei vertici duplicati, in questo esempio vedremo come selezionarli usando le espressioni.

Per selezionare i vertici duplicati vedremo due soluzioni diverse:

prima soluzione usando una funzione di aggregazione count

count(
    expression:=geom_to_wkt( $geometry,2),
    group_by:=geom_to_wkt( $geometry,2)
      ) > 1

Ogni espressione viene valutata riga per riga nella tabella attributi, questa espressione sarà vera solo in quelle righe dove il conteggio delle geometrie sono maggiori di uno e quindi duplicate: prima riga, conteggio con aggregazione, eistono più geometrie uguali?? se si il count è >1 e seleziona la riga.

seconda soluzione usando gli array

array_length(
        array_agg( 
            expression:=$geometry,
            group_by:=geom_to_wkt($geometry,2))
            ) > 1

Il concetto è simile a quello precedente, ma usando gli array: quindi quando un array contiene più di un elemento (geometria=vertice), la lunghezza dell'array sarà > 1 e quindi duplicato.

per entrambe stesso risultato:

Selezionare solo uno dei due punti doppi↵

espressione da utilizzare:

seleziona primo punto

if($id = 
  array_first( -- primo punto
    array_agg(
    expression:= $id,
    filter:=
      array_length(
      array_agg( 
      expression:=$geometry,
      group_by:=geom_to_wkt($geometry)))>1,
    group_by:=geom_to_wkt($geometry)
    )
  )
, $id, '')

selezionare ultimo punto

if($id = 
  array_last( -- ultimo punto
    array_agg(
    expression:= $id,
    filter:=
      array_length(
      array_agg( 
      expression:=$geometry,
      group_by:=geom_to_wkt($geometry)))>1,
    group_by:=geom_to_wkt($geometry)
    )
  )
, $id, '')

nel caso di molti punti sovrapposti si potrebbe utilizzare la notazione degli array, cioè le parentesi quadre:

if($id = 
    array_agg(
    expression:= $id,
    filter:=
      array_length(
      array_agg( 
      expression:=$geometry,
      group_by:=geom_to_wkt($geometry)))>1,
    group_by:=geom_to_wkt($geometry)
    )[0] -- primo vertice doppio; [-1] per ultimo; [x] per xesimo posto
, $id, '')

NB1: La selezione potrebbe non essere visibile a schermo a causa della sovrapposizione dei punti, per ovviare alternare le funzioni array_first con array_last.
NB2: I calcoli da fare sono molti e per dataset pesanti potrebbe impiegarci molto tempo (dipende anche dalle risorse hardware)

In collaborazione con Valerio Pinna che ringrazio per la disponibilità.

Prova tu↵

il geopackage è scaricabile qui

Funzioni e variabili utilizzate: