PostGIS: разбить полигон в заданном соотношении

11

Есть ли удобный способ разбить полигон на nчасти, более или менее равные по размеру в PostGIS?

Адам Матан
источник
Равная площадь или равный размер (то есть одинаковая ширина и высота)?
Энтони-GISCOE-
Площадь вне зависимости от пропорций.
Адам Матан
1
похожие посты gis.stackexchange.com/questions/5300/…
Кирк Куйкендалл
Я удивлен, что нет решения для этого в пространственном R
Брэд Несом

Ответы:

8

Это старая проблема без простого решения. Единственный подход, с которым я столкнулся, - это создать функцию, в которой вы дадите заголовок, а количество деталей и компьютер будут проверены, пока он не получит равные области. В AutoCAD есть функция LISP. В Postgis он работает так же, вот выдержка из PostGIS в действии от Мэннинга, этот код разбивает многоугольник на две равные части:

WITH RECURSIVE
ref(the_geom, env) AS (
SELECT the_geom,
ST_Envelope(the_geom) As env,
ST_Area(The_geom)/2 As targ_area,
1000 As nit
FROM us.states
WHERE state = 'Idaho'
),

T(n,overlap) AS (
VALUES (CAST(0 As Float),CAST(0 As Float))
UNION ALL
SELECT n + nit, ST_Area(ST_Intersection(the_geom, ST_Translate(env, n+nit, 0)))
FROM T CROSS JOIN ref
WHERE ST_Area(ST_Intersection(the_geom, ST_Translate(env, n+nit, 0)))> ref.targ_area
) ,  

bi(n) AS
(SELECT n
FROM T
ORDER BY n DESC LIMIT 1)  

SELECT bi.n,
ST_Difference(the_geom, ST_Translate(ref.env, n,0)) As geom_part1,
ST_Intersection(the_geom, ST_Translate(ref.env, n,0)) As geom_part2
FROM bi CROSS JOIN ref;
Pablo
источник
2

Одним из подходов может быть разделение многоугольника на треугольники, каждый с определенной площадью. Тогда нужно было бы попытаться сгруппировать эти (смежные) треугольники обратно в многоугольники (более или менее) размера области / n. Это будет своего рода настраиваемая версия проблемы «подмножество сумм» или «ранца» (и я не знаю, как начать с этого с PostGIS).

Энтони-ГИСКО-
источник
Можете ли вы дать больше информации об алгоритмах?
Маджид Ходжати