Разбить полигон на ячейки
Добавлено: 09 авг 2022, 20:13
Добрый день. Подскажите, с помощью какого набора команд можно разбить полигон на ячейки нужного размера? см. вложение
Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
Лет 10 назад понадобилась планшетная сетка, пришлось написать функцию по этому поводу. Подойдет функция и для вашего случая.
Код: Выделить всё
create or replace function
make_grid(xmin int,
ymin int,
xmax int,
ymax int,
step int,
point boolean default false)
returns table (geom geometry)
as $$
begin
if point then
return query
select st_point(x,y) from
(select * from generate_series(xmin, xmax, step) x,
generate_series(ymin, ymax, step) y) gridcoord;
else
return query
select ST_MakeEnvelope(x, y, x+step, y+step, -1) from
(select * from generate_series(xmin, xmax, step) x,
generate_series(ymin, ymax, step) y) gridcoord;
end if;
end;
$$ language plpgSQL;
comment on function make_grid(int, int, int, int, int, boolean) is
'Функция возвращает таблицу геометрий, образующих сетку из точек
или полигонов (квадратов).
аргументы (Xmin, Ymin, Xmax, Ymax, шаг сетки, флаг default false )
флаг - true - сетка строится из точек, false - из квадратов';
Код: Выделить всё
with mgrid as (
select st_transform(geom, 4326) geom
from my_poligon where id_=1
)
select
make_grid(
st_xmin(geom)::int,
st_ymin(geom)::int,
st_xmax(geom)::int,
st_ymax(geom)::int,
1,
false
) geom
from mgrid;