Pseudoknoten in freier GIS-Software finden?

Die Software gvSIG OA Digital Edition 2010 verfügt über Tools zur Topologie zum Auffinden von Pseudoknoten in linearer Geometrie. Ich habe die Clustertoleranz auf 0,00002 und die maximale Anzahl von Fehlern auf -10000 für die lineare Geometrie mit 20000 Verbindungszahlen festgelegt. Aber erfolgloses Ergebnis.

Gibt es Lösungen, die Pseudoknoten in freier GIS-Software finden?

Ich muss Pseudoknoten überlagern (eine Lösung für dieses Problem ist die Verwendung der Tools-Topologie von ArcInfo, aber für mich ist die Verwendung von freier Software die Priorität). Mit der linearen Geometrie wurden mehrere Benutzer in QGIS 1.8.0 in der PostGIS-Datenbank (Version 2.0.1) erstellt.

Neues Bild hinzufügen: 12 lineare Features mit drei Pseudoknoten in A (Zeile 4/5), B (Zeile 6/7), C (Zeile 9/10). Pseudoknoten sollten stattdessen Punkte sein - zwei lineare Merkmale mit Schnittpunkt in einem Punkt (Knoten) sollten ein lineares Merkmal sein (Linie 4/5 - Linie 4, ...).

Ist es möglich, eine Anfrage in PostGIS zu stellen, die zu einer Schicht von Pseudoknoten führt?

Neues Bild von Beispiel-Pseudoknoten hinzufügen: Wenn ich für die lineare Ebene Punktebene-Pseudoknoten (blaue Rechtecke) erhalte, habe ich folgende Fehler in der linearen Ebene korrigiert: A - fehlende Geometrie hinzufügen, B - Schnittlinie einrasten lassen, C - Pseudoknoten entfernen.

Hier eine generische Lösung, die Sie mit PostGIS oder einer anderen OGC-konformen Software umsetzen können.

HINWEIS: wie ich sage vor , ein Schlüsselbegriff in FOSS und GIS ist Standardisierung : Die besten Lösungen verabschieden Standards wie OGC diejenigen.

Ihr Problem ist, "Pseudoknoten zu finden" ... Aber ich denke, dass es ein bisschen mehr ist, "Nicht-Pseudoknoten zu finden und Linien von Pseudoknoten zu verbinden". Meine Lösung kann für beide verwendet werden.

OGC-Standards bieten:

• ST_Boundary (geom) : um die Knoten der Linien zu erkennen

• ST_Dump (geom) : um jeden einzelnen Knoten in einen SQL-Tabellendatensatz einzufügen .

• ST_DWithin, ST_Equals, ST_SnapToGrid, ST_Snap können für die Änderungstoleranz verwendet werden. Ich benutze ST_DWithin.

Wir können davon ausgehen, dass Ihr Hauptproblem mit diesen Objekten und Eigenschaften spezifiziert werden kann.

• nur Liniensegmente (eine Tabelle sind Liniensegment ), durch eine Linienfolge Geometrie dargestellt ... ich nicht mit MULTILNE getestet, wenn Sie GeometryType haben = MEHRFACH, können Sie teilen und Gussmultilinien mit ST_Dump und ST_LineMerge;

• Jedes Liniensegment hat eine (Geometrie-ID) -ID und eine (Farb-ID) -ID .

Der erste Schritt besteht also darin, die Knoten zu ermitteln, die von Verbindungslinien stammen.

CREATE TABLE cache_bounds AS
  SELECT gid as gid_seg, (ST_Dump(ST_Boundary(the_geom))).geom AS the_geom,
         gid as color 
         -- if you not have something for "color label" of lines, use gid.
  FROM linesegment;
ALTER TABLE cache_bounds ADD column gid serial PRIMARY KEY;

CREATE TABLE cache_joinnodes AS
  -- Use your TOLERANCE instead "1" at ST_DWithin and ST_Buffer.
  SELECT *, array_length(colors,1) as ncolors FROM (
   SELECT gid, array_distinct(array_cat(a_colors,b_colors)) as colors, the_geom FROM (
    SELECT 
      a.gid, array_agg(a.color) as a_colors, array_agg(b.color) as b_colors
      , st_buffer(a.the_geom,1) as the_geom -- any one to represent the join point.
    FROM cache_bounds a, cache_bounds b 
    WHERE a.gid>b.gid AND ST_DWithin(a.the_geom,b.the_geom,1)
    -- use ST_equals(a.the_geom,b.the_geom) if no tolerance.
    GROUP BY a.gid, a.the_geom
   ) as t
  ) as t2;

HINWEIS: Verwenden von Caches, da diese schneller sind als Ansichten. Verwenden Sie "EXPLAIN SELECT ...", um die CPU-Zeit zu überprüfen. Dies kann sehr lange dauern.

Hier werden Zyklen und durchgehende (gleichfarbige) Linien als ncolors=1Punkte und die Pseudoknoten als Punkte erkannt. ncolors=2Sie haben also eine Ebene mit diesen Punkten.

Ihre Tabelle mit "guten Knoten" enthält die ursprünglichen "Begrenzungspunkte" und keine "Pseudoknoten".

CREATE VIEW vw_joinnodes_full AS
  SELECT b.*, j.ncolors
  FROM cache_joinnodes j INNER JOIN cache_bounds b 
       ON j.gid=b.gid;

CREATE TABLE cache_good_nodes AS
  SELECT *  
  FROM vw_joinnodes_full 
  WHERE ncolors=1 OR ncolors>2;

-- IF NEED ... CREATE VIEW vw_correct_linesegment AS ... 

Refractions Research hat ein Linienreinigungs-Tool entwickelt , das zu tun scheint, was Sie wollen.

Line Cleaner bereinigt Netzwerke durch Vereinfachung komplexer, zyklischer, sehr kurzer und längenloser Geometrien und durch Entfernen von Pseudoknoten und unbedeutenden Scheitelpunkten. Am wichtigsten ist, dass in der Bereinigungsphase sichergestellt werden kann, dass Feature-Übereinstimmungen automatisch berücksichtigt werden

Der Quellcode ist bei GitHub zu finden.

Non-Free-Lösung: FME + MRF + SmartCleaner-Transformator

Kostenlose Lösung GRASS v.clean (das neueste QGIS 1.8.0 mit GRASS-Tools ist die einfachste Möglichkeit, es zu verwenden) und andere Tools zur Reinigung der Topologie

Hier sind Schritte, um Ihre Pseudoknoten mithilfe von OpenJump, einem kostenlosen GIS, zu finden.
QGIS und gvSIG haben das Sextante-Plugin, daher sollten die gleichen Schritte auch funktionieren, da
der Spatial-Join möglicherweise etwas anders ist.
Ich habe zum Testen Version 1.2 verwendet.

- Speichern Sie die
Toolbox Sextante Linienendpunkte , Topologie, Linienendpunkte extrahieren -> endpt_0

- Aufheben der Teilung Ihrer Linien
Toolbox Sextante, Werkzeuge für

Linienebenen , Verbinden benachbarter Linien - Speichern Sie die Toolbox
Sextante Linienendpunkte , Topologie, Extrahieren der Endpunkte von line -> endpt_1

- die Punkte Ende durch Pseudo - Knoten „benachbarte Linien verbinden“ entfernt

Werkzeuge, Abfragen, Spatial Query -
Source - Schicht „endpt_0“
Relation „erbohrt“
Maskenebene "endpt_1"

aktivieren oder auf Ergebnis ergänzen klicken

Lösung zum Auffinden von Pseudoknoten (und anderen Topologieregeln) in QGIS 1.9.0 master (2.0) - Plugin "Topology Checker". Links: Topology Checker Plugin auf YouTube , Topology Checker Plugin auf GitHub

In PostGIS können Sie eine geänderte Version der Abfrage verwenden, um in diesem Thema behandelte Dangles zu finden , da Pseudoknoten Knoten sind, die 2 Linienfolgen abfangen, und Dangles Knoten sind, die 1 Linienfolge abfangen.

WITH nodes AS 
(SELECT ST_StartPoint(geom) AS pt FROM
linestring_table UNION ALL 
SELECT ST_EndPoint(geom) AS pt FROM
linestring_table) 
SELECT pt FROM nodes
GROUP BY pt HAVING count(*) = 2;