zoukankan      html  css  js  c++  java
  • 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便于编辑样式。

    更改网格大小为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)

  • 相关阅读:
    Sqli-labs less 12
    Sqli-labs less 13
    Sqli-labs less 14
    Python3之collections模块
    Python3之 contextlib
    Python3之sqlalchemy
    Python3之RabbitMQ
    Python3之redis使用
    Python3之Memcache使用
    python自动化开发学习 进程, 线程, 协程
  • 原文地址:https://www.cnblogs.com/shitao/p/2109620.html
Copyright © 2011-2022 走看看