¿Encontrar pseudo nodos en software SIG gratuito?

El software gvSIG OA Digital Edition 2010 tiene herramientas de topología para encontrar pseudo nodos en geometría lineal. Configuré la tolerancia de clúster 0.00002 y el número máximo de errores -10000 para la geometría lineal de conteo de enlaces 20000. Pero resultado infructuoso.

¿Hay alguna solución que encuentre pseudo nodos en el software SIG gratuito?

Necesito superponer pseudo nodos (una solución a este problema: usar la topología de herramientas de ArcInfo, pero la prioridad para mí es usar software libre). La geometría lineal creó varios usuarios en QGIS 1.8.0 en la base de datos PostGIS (v. 2.0.1).

Agregue una nueva imagen: 12 características lineales con tres pseudo nodos en A (línea 4/5), B (línea 6/7), C (línea 9/10). Los pseudo nodos deberían ser puntos en su lugar: dos entidades lineales con intersección en un punto (nodo) deberían ser una entidad lineal (línea 4/5 - línea 4, ...).

¿Es posible hacer una solicitud en PostGIS, lo que dará como resultado una capa de pseudo nodos?

Agregue una nueva imagen de pseudo nodos de ejemplos: si recibo pseudo nodos de capa de punto de capa lineal (rectificaciones azules), corregí los siguientes errores en la capa lineal: A - agregue geometría faltante, B - línea ajustada en la intersección, C - elimine el pseudo nodo.

Aquí una solución genérica, que puede mejorar con PostGIS o cualquier otro software compatible con OGC.

NOTA: como dije antes , un concepto clave en FOSS y GIS es la estandarización : las mejores soluciones adoptan estándares, como OGC .

Su problema es "encontrar pseudo nodos" ... Pero creo que es un poco más, "encontrar no pseudo nodos y unir líneas de pseudo nodos". Mi solución puede usarse para ambos.

Oferta de estándares OGC:

• ST_Boundary (geom) : para detectar los nodos de las líneas

• ST_Dump (geom) : para poner cada nodo individual en un registro de tabla SQL.

• ST_DWithin, ST_Equals, ST_SnapToGrid, ST_Snap se pueden usar para la tolerancia al cambio. Estoy usando ST_DWithin.

Podemos suponer que su problema principal se puede especificar con estos objetos y propiedades,

• solo hay segmentos de línea (de una tabla en línea) ), representados por una geometría LINESTRING ... No probé con MULTILNE, si tiene geometrytype = MULTIPOINT, puede dividir y emitir MULTILINEs con ST_Dump y ST_LineMerge;

• cada segmento de línea tiene un gid (ID de geometría) y un idline (ID de color) .

Entonces, el primer paso es obtener los nodos que provienen de unir líneas,

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;

NOTA: usar cachés porque son más rápidos que las vistas. Use "EXPLAIN SELECT ..." para verificar el tiempo de CPU, puede llevar mucho tiempo.

Aquí se detectan ciclos y líneas continuas (del mismo color) como ncolors=1puntos, y los pseudo nodos por ncolors=2puntos, por lo que tiene una capa con esos puntos.

Su tabla de "nodos buenos" está con los "puntos delimitadores" originales y sin "pseudonodos".

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 ha creado una herramienta Line Cleaner que parece hacer lo que usted desea.

Line Cleaner limpia redes simplificando geometrías complejas, cíclicas, muy cortas y de longitud cero, y eliminando pseudo-nodos y vértices insignificantes. Lo más significativo es que, en la fase de limpieza, puede garantizar que las coincidencias de características se puedan considerar automáticamente

El código fuente se puede encontrar en GitHub.

Solución no libre: transformador FME + MRF + SmartCleaner

Solución gratuita GRASS v.clean (El último QGIS 1.8.0 con herramientas GRASS es la forma más fácil de usarlo) y otras herramientas de limpieza de topología

Aquí hay pasos para encontrar sus pseudo nodos usando OpenJump y un SIG gratuito.
QGIS y gvSIG tienen el complemento Sextante, por lo que estos mismos pasos también deberían funcionar,
la unión espacial puede ser ligeramente diferente.
Usé la versión 1.2 para probar.

- guarde la
caja de herramientas Sextante de puntos finales de línea , topología, extraiga puntos finales de líneas -> endpt_0

- separe las líneas de la
caja de herramientas Sextante, herramientas para capas de línea, unir líneas adyacentes

- guarde la caja de
herramientas Sextante de puntos finales de línea sin dividir , topología, extraiga los puntos finales de línea -> endpt_1

- los puntos finales eliminados por "Unir líneas adyacentes" son pseudo nodos

Herramientas, Consultas, Consulta espacial,
Capa de origen "endpt_0"
Relación "Interseca"
Habilite la capa "endpt_1"

o haga clic en Complemento Resultado

Solución para encontrar pseudo nodos (y otras reglas de topología) disponibles en QGIS 1.9.0 master (2.0) - plugin "Topology Checker". Enlaces: complemento de Topology Checker en YouTube , complemento de Topology Checker en GitHub

Con PostGIS, puede usar una versión modificada de la consulta para encontrar los colgantes discutidos en este tema , ya que los pseudonodos son nodos que interceptan 2 cadenas lineales y los colgantes son nodos que interceptan 1 cadena lineal.

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;