postgis将multi和geometrycollection转换成为多行数据，自定义函数

Posted 2023-02-14 JerFer

工作中用到st_intersection 函数，返回的内容五花八门，简单图形，多部件，甚至于也有GeometryCollection ，这样的数据参与下一次计算的时候不是很方便，所以做了一个拆分，然后过滤想要的纬度的图形。

CREATE OR REPLACE FUNCTION geometry_to_array(in_geom "public"."geometry",demon int4)
  RETURNS "public"."geometry" AS $BODY$  
declare 
    result "public"."geometry"[];  
		n int4;
	tmp "public"."geometry";
BEGIN  
		n= st_numgeometries(in_geom);
		FOR var IN 1..n LOOP
			tmp =st_geometryn(in_geom,var) ;
			if st_dimension(tmp) = demon then 
				result[var]=st_geometryn(in_geom,var);
			end if ;
		END LOOP;
   RETURN result;  
END;  
$BODY$
  LANGUAGE plpgsql

后面发现其实PostGIS本身是提供了类似功能的。

# st_buffer是为了避免交集出现自相交的情况
(ST_Dump(st_intersection(st_buffer(t.shape,0.00), st_buffer(s.shape,0.00)))).geom