Blog dataconsulting.pl

Czasem funkcje wbudowane (st_sharedpatsh, st_intersect) są niewystarczające. Część danych geograficznych zapisanych w systemach GIS zawiera informacje o przecięciach tych obiektów w formie zdefiniowania punktu wspólnego dla obu obiektów w miejscu przecięcia (tak np. działa format OSM). PostGIS nie potrafi prosto rozróżnić np. linii które się przecinają i nie są zdefiniowane w oparciu o wspólny węzeł oraz tych które wspólny węzeł posiadają. Funkcje st_shared_paths, st_intersect nie są w tym przypadku wystarczające, a funkcji st_sharedpoints nie ma (za to napiszemy swoją!). Przykładowo:


SELECT st_asText('LINESTRING(-10 0, 10 0)'::geometry);

SELECT st_asText('LINESTRING(0 -10, 0 10)'::geometry);


Te dwie linie przecinają się, ale nie posiadają wspólnego węzła. Z punktu widzenia PostGIS’a wynik sprawdzenia ich przecięcia zwróci wartość TRUE:


SELECT st_Intersects('LINESTRING(-10 0, 10 0)'::geometry, 'LINESTRING(0 -10, 0 10)'::geometry);


Czy da się te sytuacje w PostGIS rozróżnić? Oczywiście!

Continue reading →