Разбить полигон на ячейки

[Svility]

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

IMG_2301.PNG (158.01 КБ) 10299 просмотров

[trir]

ST_SquareGrid, ST_Intersection

[Svility]

ST_ SquareGrid

А если версия postgis 2.5 ?

[drrrrr]

А если версия postgis 2.5 ?

Лет 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 - из квадратов';

[Svility]

Код: Выделить всё

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;

А шаг сетки <1км как сделать?

[freeExec]

Использовать метрическую проекцию, я думаю