Introduction to PostGIS 之 Spatial clustering(二)

上一篇中主要介绍了ST_SnapToGrid,现在我们用前几篇文中的features表来做下应用:

简要介绍下用到的方法:

geometry ST_Collect(geometry[] g1_array);把一系列geometries组合为单个ST_Geometry(旧版本称为GEOMETRYCOLLECTION).

geometry ST_Multi(geometry g1);把g1转化为一个MULTI* geometry。

和前几篇方法一样,连接数据库,打开查询可视化框,添加如下SQL:

SELECT ST_Multi(ST_Collect(the_geom)) As the_geom , ST_astext(ST_Multi(ST_Collect(the_geom)))
FROM features
GROUP BY ST_SnapToGrid(ST_Transform(the_geom,2327),2000,2000)

      把features中的点集聚到2000*2000 平方米的格子中,ST_Transform(the_geom,2327)是把坐标系统从4326转换为2327,因为4326的测量单位是度,2327(这里我随便找了个xi‘an 80的坐标系统)应该是米。ST_astext(ST_Multi(ST_Collect(the_geom))没多大用,只是在QGIS便于编辑样式。

Introduction to PostGIS 之 Spatial clustering(二)_第1张图片

更改网格大小为5000*5000:

SELECT ST_Multi(ST_Collect(the_geom)) As the_geom , ST_astext(ST_Multi(ST_Collect(the_geom)))
FROM features
GROUP BY ST_SnapToGrid(ST_Transform(the_geom,2327),5000,5000)

Introduction to PostGIS 之 Spatial clustering(二)_第2张图片

你可能感兴趣的:(cluster)