SQL

Объясните мне как работает st_snap в postgis??

Прочитал документацию , пробую сделать примеры, но не получается.
 SELECT ST_AsText( 

    ST_Snap(poly,line, ST_Distance(poly,line)*10) 

  ) AS polysnapped 

FROM (SELECT 

  ST_GeomFromText('POLYGON((30.59280395507813 60.100456776182845,30.45547485351563 59.992118742381244,30.60104370117188 59.93300042374631,30.82901000976563 60.06209914960289,30.69442749023438 60.13466696253271,30.59280395507813 60.100456776182845))') As poly, 

       ST_GeomFromText('LINESTRING(30.778198242187504 60.098403034443926,30.85372924804688 59.89995826181929)') As line 

  ) As foo
Почему он мне возвращает тот же полигон без привязки линии?
Что за параметры передаются в скобках? я так понимаю, первый - основной полигон, второй - линия, полигон, который я хочу привязать, третий - это радиус который учитывается при привязке?
Что значит as polysnapped, as foo?
Как мне привязать например, вместо линии, другой полигон?
Олег
Олег
482
Столько вопросов...

Что значит as polysnapped, as foo?

Это псевдонимы, используемые для доступа. В частности, polysnapped — это псевдоним (новое имя) для поля, определяемого выражением
 ST_AsText(ST_Snap(...))
foo — псевдоним (новое имя) для подзапроса
 SELECT  

  ST_GeomFromText('...') As poly,  

  ST_GeomFromText('...') As line
То есть фактически твой запрос возвращает поле polysnapped из таблицы (виртуальной) foo. Само имя foo никак не используется, но его указание требуется синтаксисом.


Как мне привязать например, вместо линии, другой полигон?

В запросе для foo указать вместо линии другой полигон (я переименовал поля для наглядности):
 SELECT ST_AsText( ST_Snap(poly1, poly2, ST_Distance(poly1, poly2) * 10) ) AS polysnapped  

FROM (SELECT

  ST_GeomFromText('POLYGON(...)') As poly1,  

  ST_GeomFromText('POLYGON(...)') As poly2

) AS foo


Почему он мне возвращает тот же полигон без привязки линии?

Ну так в документации же сказано: «If no snapping occurs then the input geometry is returned unchanged». А вот почему no snapping occurs, можно только предполагать. Скорее всего, как опять же сказано в документации, «Too much snapping can result in invalid topology being created... This can result in some potential snaps being omitted». А вот это too much snapping вполне могло произойти от *10 в твоей привязке, в то время как в примерах *1.01 и * 1.25. Как говорится, почувствуйте разницу.


Что за параметры передаются в скобках? я так понимаю, первый - основной полигон, второй - линия, полигон, который я хочу привязать, третий - это радиус который учитывается при привязке?

Да. Сперва из форматированного текста функция ST_GeomFromText создаёт объекты, которые затем предаются в функцию ST_Snap.
Дима Пас
Дима Пас
12 091
Лучший ответ
Олег Искренне благодарю за такой развернутый ответ!!
Еще вопрос, может, вы поможете.
Я использую эту функцию чтобы получить одну общую сплошную границу, между полигонами, но получается абсолютно не то.
Вот пример, 2 полигона, хочу чтобы их границы соединились, но чтобы не получился один сплошной полигон:Вот, что у мня получается:Что я делаю не так?
Хочу добиться вот такого эффекта, именно так будут выглядеть границы моих полигонов:
было:стало:Может быть, я использую не те функции?
Дима Пас Сперва уточню ответ.

я так понимаю, первый - основной полигон, второй - линия, полигон, который я хочу привязать, третий - это радиус который учитывается при привязке?

Немного не так.

Первый параметр — то, что привязывается, то есть обрабатываемая геометрия. В результате привязки она может (и должна, вообще-то) измениться.

Второй параметр — то, к чему привязывается, вершины этой геометрии служат в качестве сетки. В результате привязки она не изменяется.

Третий параметр — расстояние, на котором для элементов обрабатываемой геометрии (первый параметр) видны потенциальные узлы привязки или, иначе, на какое расстояние могут перемещаться элементы обрабатываемой геометрии.
Дима Пас
Почему он мне возвращает тот же полигон без привязки линии?

Дело в том, что данном конкретном случае линия и многоугольник пересекаются, то есть расстояние между ними равно нулю. Если ноль умножить на что угодно, всё равно получится ноль, то есть третий параметр (расстояние привязки) равен нулю, а это как раз и приводит к неизменности обрабатываемой геометрии.

В таком случае вместо 
 ST_Distance(poly1, poly2)
лучше использовать
 ST_Distance(ST_Points(poly1), ST_Points(poly2))
то есть вычислять расстояние строго по вершинам, без учёта рёбер.
Дима Пас
Что я делаю не так?

Судя по результату, ты задаёшь слишком большое расстояние привязки.
Дима Пас Я набросал небольшой запрос. Поработай с ним, увидишь, как влияет расстояния привязки на результат.
 SELECT mult, 

       ST_AsText( ST_Snap(poly1, poly2, ST_Distance(ST_Points(poly1), ST_Points(poly2)) * mult) ) AS poly1_snapped_to_poly2,   

       ST_AsText( ST_Snap(poly2, poly1, ST_Distance(ST_Points(poly1), ST_Points(poly2)) * mult) ) AS poly2_snapped_to_poly1 

FROM (SELECT  

  ST_GeomFromText('POLYGON((...))') As poly1,  

  ST_GeomFromText('POLYGON((...))') As poly2 

) AS geoms, 

(VALUES (1.0), (1.1), (1.25), (1.5), (1.75), (2.0), (2.5), (3.0), (4.0), (5.0)) AS mults (mult);

В POLYGON((...)) вместо многоточий надо поставить свои координаты.

Множители подбирал на глазок по твоей карте.

Похожие вопросы