Форумы GIS-Lab.info

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

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;