问题

    PostGIS是否有方法能将一个Polygon面切割成若干份小的Polygon面，且每一份的面积差不多大？其实并没有现成的方法，但是通过灵活运用postgis函数可以快速实现这样的功能，总共只要简单的5步就可以了，下文具体说明。

二 步骤分解

2.1 准备测试数据

本文选择南京市行政区边界作为测试数据。

create table nanjing as

select name,geom from city where name=‘南京市’;

2.2 面内生成任意数量点

本文在南京区域内生成2000测试点。

CREATE TABLE nanjing_points AS

SELECT (ST_Dump(ST_GeneratePoints(geom, 2000))).geom

AS geom FROM nanjing;

2.3 点聚合成簇

使用ST_ClusterKMeans均值中心聚类分析，把散点聚簇分类。

CREATE TABLE nanjing_pts_clustered AS

SELECT geom, ST_ClusterKMeans(geom, 10) over () AS cluster

FROM nanjing_points;

2.4 提取簇的中心

CREATE TABLE nanjing_centers AS

SELECT cluster, ST_Centroid(ST_collect(geom)) AS geom

FROM nanjing_pts_clustered GROUP BY cluster;

2.5 使用voronoi算法生成面

CREATE TABLE nanjing_voronoi AS

SELECT (ST_Dump(ST_VoronoiPolygons(ST_collect(geom)))).geom AS geom

FROM nanjing_centers;

2.6 使用ST_Intersection方法切割

用voronoi算法构成的面，去切割原始的南京行政区边界：

CREATE TABLE nanjing_divided AS
SELECT ST_Intersection(a.geom, b.geom) AS geom
FROM nanjing a
CROSS JOIN nanjing_voronoi b;

三 总结

    Postgis中造数据