通用操作

1)客户端连接实例

2)初次装载PostGIS扩展模块

创建扩展：

create extension postgis;

查看版本：

select postgis_version();

select postgis_full_version();

3)空间数据写入数据库表

首先创建带Geometry字段的表，SQL参考：

create table testg ( id int, geom geometry )

distributed by (id);

该SQL表示插入的空间数据不区分几何类型，几何类型包括Point / MultiPoint / Linestring / MultiLinestring / Polygon

/ MultiPolygon等。如果在创建表时已知Geometry类型和SRID，SQL参考：

create table test ( id int, geom geometry(point, 4326) )

distributed by (id);

Geometry类型指定Point，SRID为4326。SRID不指定默认为0。

有关SRID可参考 SRID

数据插入，SQL参考：

-- without srid

insert into testg values (1, ST_GeomFromText('point(116 39)'));

-- with srid

insert into test values (1, ST_GeomFromText('point(116 39)', 4326));

JDBC Java程序参考：

import java.sql.Connection;

import java.sql.DriverManager;

import java.sql.Statement;

public class PGJDBC {

public static void main(String args[]) {

Connection conn = null;

Statement stmt = null;

try{

Class.forName("org.postgresql.Driver");

//conn = DriverManager.getConnection("jdbc:postgresql://:3432/","", "");

conn.setAutoCommit(false);

stmt = conn.createStatement();

String sql = "INSERT INTO test VALUES (1001, ST_GeomFromText('point(116 39)', 4326) )";

stmt.executeUpdate(sql);

stmt.close();

conn.commit();

conn.close();

} catch (Exception e) {

System.err.println(e.getClass().getName() + " : " + e.getMessage());

System.exit(0);

}

System.out.println("insert successfully");

}

}

4)空间索引管理

创建空间索引

create index idx_test_geom on test using gist(geom);

idx_test_geom为自定义索引名，test为表名，geom为Geometry列名。

查看表有哪些索引

select * from pg_stat_user_indexes

where relname='test';

查看索引大小

select pg_indexes_size('idx_test_geom');

索引重构

reindex index idx_test_geom;

删除索引

drop index idx_test_geom;

5)典型空间查询SQL

矩形范围查询

-- without srid

select st_astext(geom) from testg

where ST_Contains(ST_MakeBox2D(ST_Point(116, 39),ST_Point(117, 40)), geom);

-- with srid

select st_astext(geom) from test

where ST_Contains(ST_SetSRID(ST_MakeBox2D(ST_Point(116, 39),ST_Point(117, 40)), 4326), geom);

ST_MakeBox2D算子生成一个Envelope。

几何缓冲范围查询

-- without srid

select st_astext(geom) from testg

where ST_DWithin(ST_GeomFromText('POINT(116 39)'), geom, 0.01);

-- with srid

select st_astext(geom) from test

where ST_DWithin(ST_GeomFromText('POINT(116 39)', 4326), geom, 0.01);

ST_DWithin用法参考：ST_DWithin

多边形相交判定(在内部或在边界上)

-- without srid

select st_astext(geom) from testg

where ST_Intersects(ST_GeomFromText('POLYGON((116 39, 116.1 39, 116.1 39.1, 116 39.1, 116 39))'), geom);

-- with srid

select st_astext(geom) from test

where ST_Intersects(ST_GeomFromText('POLYGON((116 39, 116.1 39, 116.1 39.1, 116 39.1, 116 39))', 4326), geom);

ST_*算子对大小写不敏感，更多用法可参考 PostGIS官方资料

说明

AnalyticDB PG 6.0不完全兼容PostGIS功能集，例如不支持 create extension postgis_topology，不推荐用Geography类型创建表(非要用，SRID默认为0或4326)。