Jak uprościć routowalną sieć?

Mam wykres sieci, który muszę uprościć w celu zmniejszenia liczby krawędzi . Pomysł polegałby na połączeniu węzłów znajdujących się blisko siebie i usunięciu łączących krótkich krawędzi.

Jak można to osiągnąć w PostGIS lub GRASS? Czy są jakieś lepsze sposoby automatycznego uproszczenia takiej sieci?

Próbowałem już funkcji ST_SnapToGrid, ale nie jestem zadowolony z wyników (szary = oryginalny, czarny = przyciągnięty):

Najbliższe, do jakiego do tej pory doszedłem, to:

v.clean input=roads output=snap5rmline tool=snap,rmline thresh=5 

Zatrzaskuje drogi z tolerancją 5 metrów i usuwa wszystkie linie zerowej długości. Nie jest to optymalne rozwiązanie, ponieważ wydaje się przyciągać raczej losowo do niektórych wierzchołków.

Czy próbowałeś już GRASS v.generalize ?

v.generalizepozwala wybrać algorytm generalizacji z atrybutem metody . Istnieje kilka: douglas, douglas_reduction, lang, redukcja, reumann, boyle, przesuwne_averaging, ważenie odległości, chaiken, pustelnik, węże, sieć, przemieszczenie .

I dodatkowe parametry jak threshold, degree_thresh, angle_thresh(w zależności od wybranego algorytmu) może pomóc uzyskać dokładny wynik.

Oto samouczek .

Nie zrobiłem tego, ale myślę, że mogę zasugerować kierunek.

1. Utwórz topologię za pomocą PostGIS dla swojego wykresu.
2. Znajdź wszystkie węzły z tylko dwiema krawędziami.
3. Wylecz krawędzie.

ST_ModEdgeHeal połączy jedną krawędź w drugą. ST_NewEdgeHeal zastąpi oba nowymi krawędziami.

Podręcznik topologii PostGIS

@ w mroku, widzę, że napisałeś narzędzie do zagęszczania linii w Sextante. Dlatego sugeruję następujący algorytm, aby uniknąć „losowego” przyciągania jednego z twoich punktów.

Wybierz segmenty linii, których chcesz się pozbyć na podstawie ich długości.

Dla każdego z tych segmentów utwórz punkt w punkcie środkowym

Usuń mały segment

Teraz możesz używać ST_Snap w PostGIS (patrz przykład tutaj )

EDYCJA: pamiętaj, że w twoim przypadku możesz najpierw użyć v.net , aby usunąć pseudo-węzły (węzeł łączący tylko dwie linie)

Przekazywanie, w jaki sposób Michaël Michaud przeanalizował to na liście programistów OpenJUMP:

> Interesting question. There has been so much research works about
> generalization that it is surely not a simple task. I have tested the
> following approach with mitigated results :
> - make the layer planar with the noder plugin
> - detect small cycles with the graph plugin
> - merge adjacent cycles
> - create a point inside each cycles surface
> - remove network segments along these cycles (query + DE-9IM)
> - detect roads touching the cycles
> - project the center points of cycles to the road endpoints ==> replace small roundabout by simple nodes
> 
> Possible improvements (probably need a dedicated plugin)
> - make it work on non planar graph (or just remove bridges/tunnels from the graph)
> - replace small edges between two nodes with degree 3+ by a single node,    not only small cycles
> - process iteratively starting with smallest edges/cycles