Postgis随机点在多边形内

Question

如果我在Postgis中有一个多边形,我怎么能找到 - 计算这个多边形内的随机点？

Answer 1

@Mike引用的链接没有代码,但来自Dr.JTS的好线索:"点密度"地图......"基本上这涉及创建一组N个随机放置的点,它们位于给定的多边形内".一个函数做到这一点:输入是多边形,输出是随机点.

这些链接具有相同的SQL/PostGIS功能RandomPoint(Geometry):sorokine 2011 和osgeo.org/postgis/wiki.第二个链接(wiki)更完整,解释和显示示例,以及RandomPointsInPolygon(geometry,integer)解决问题的函数.

将solion扩展为输入"每个区域的点密度"或点之间的平均距离:

 CREATE OR REPLACE FUNCTION RandomPointsInPolygon(
     geom geometry,                -- the polygon
     avg_dist float DEFAULT 20.0,  -- average of 20 units between points 
     min_pts integer DEFAULT 1,    -- min. of points
     max_pts integer DEFAULT 1000  -- max. of points
 ) RETURNS SETOF geometry AS 
 $$
   SELECT CASE WHEN npts=1 THEN ST_Centroid($1) 
               ELSE RandomPointsInPolygon($1,npts) 
          END
   FROM (
     SELECT CASE WHEN d<$3 THEN $3 WHEN d>$4 THEN $4 ELSE d END AS npts
     FROM (SELECT (st_area($1)/(pi()*($2/2.0)^2))::integer AS d) AS t
   ) AS t2;
 $$ LANGUAGE SQL;

Answer 2

Postgis 2.3.0 及更高版本有一个新函数，可以将点生成到多边形ST_GeneratePoints中。