Postgis使用
一、安装与配置
先查看目前yum源中有哪些版本可用，并安装合适的版本
yum search postgis
安装Postgresql 10对应的版本PostGIS
yum install -y postgis24_10
二、创建GIS数据库
PostGIS是Postgresql的一个扩展，创建数据库并加载PostGIS扩展：
create database geo;
create extension postgis;
三、几何对象
几何对象输入
PostGIS支持很多几何类型：点、线、多边形、复合几何体等，并提供大量实用的相关函数。对于PostGIS而言，所有几何对象都有一个公共父类Geometry，所以在实际使用中，通常PostGIS的空间数据类型使用统一的Geometry类型，无论是点、折线还是多边形，都可以放入Geometry类型字段中。
PostGIS支持多种空间对象创建方式，大体上可以分为几类：
文本格式（WKT，Well-Known-Text）
二进制格式（WKB，Well-Known-Binary）
GeoJSON等编码
返回几何类型的函数
注：所有PostGIS相关函数都是以st开头
例如创建几何点（1,2），可以通过以下四种方式，得到的结果都一样
SELECT
'Point(1 2)'::GEOMETRY AS wkt,
'0101000000000000000000F03F0000000000000040':: GEOMETRY AS wkb,
st_geomfromgeojson('{"type":"Point","coordinates":[1,2]}') AS geo_json,
ST_Point(1,2) AS func;
几何对象存储
PostGIS的几何类型对Postgresql内建的几何类型使用了不同的存储方式。以点为例，使用内置的Point与ST_Point创建一个点，返回结果如下：
SELECT Point(1,2),ST_Point(1,2)
Postgresql中的Point只是一个包含两个Double的结构体（16字节），但PostGIS的点类型ST_Point则采用了不同的存储方式（21字节），除了两个坐标分量，他还包括了一些额外元数据：例如几何对象的实际类型、参考系的ID等。
如果需要人类可读的格式，则可以用ST_AsText输出WKT
SELECT st_astext(ST_Point(1,2))
几何对象运算
（1）两点间距离（几何距离）
SELECT ST_Point(1,1) <-> ST_Point(2,2)
注：上述方法只能用于几何距离计算，不能用于坐标距离计算。
（2）地理坐标距离计算
地理坐标距离计算需要引入地理坐标系（4326号坐标系，指代WGS84国际标准GPS坐标系），只有当两个点用同一坐标系时才能进行计算。
SELECT st_geomfromtext('POINT(116.321367 39.966956)',4326) :: geography <-> st_geomfromtext('POINT(116.315346 39.997398)',4326) :: geography
计算出A（116.321367，39.966956），B（116.315346，39.997398）两点间的距离为3.4km
ST_GeomFromText(character-string[, srid]) 根据字符串表示构造几何
参数
SELECT st_geomfromtext('POINT(116.321367 39.966956)',4326)
st_length可以计算某条路（折线）的总长度
st_area可以计算出某个闭合区域的面积
四、应用场景
场景一
给出一个中心点，找出该点周围一定距离 范围内所有符合条件的对象。例如，找出以用户为中心，周围1公里内所有的公交站，并按距离远近排序。
假设用户正在A地铁站:(116.321367,39.966956)，使用PostGIS计算距离
CREATE TABLE stations(name TEXT,position geography,
);
SELECT name,ST_Point(116.321367,39.966956)::geography <-> position::geography AS distance FROM stations WHERE st_point(116.321367,39.966956)::geography <-> position::geography < 500
ORDER BY st_point(116.321367,39.966956)::geography <-> position::geography;
场景二
电子围栏：用点和距离画圆圈人是一种常见场景，另外一种场景是，判断一个点落在了哪些地理围栏中
例如有车辆和用户的位置坐标，现在希望从坐标得到用户所处的城市(或者区域、商圈等)，又比如共享单车的禁停区检测，无人机的禁飞区识别，都是这种场景。
--兴趣区域（AOI，Area of Interest）
CREATE TABLE aoi(name TEXT,bound GEOMETRY
)
--检测A地铁站中心点所属的商圈
SELECT name FROM aoi WHERE ST_Contains(bound,ST_Point(116.458855,39.909863));
--或
SELECT name FROM aoi WHERE ST_Within(ST_Point(116.458855,39.909863),bound);
ST_Contains与ST_Within介绍
添加链接描述
boolean ST_Contains(geometry A, geometry B);
如果geometry B完全在geometry A里面则返回TRUE
添加链接描述
boolean ST_Within(geometry A, geometry B);
如果geometry A完全在geometry B里面则返回TRUE
ST_Contains 和ST_Within是相反的,即 ST_Contains(A,B) = ST_Within(B,A).
SELECT st_contains(st_setsrid(st_geomfromtext('POLYGON ((113.49628347923993 34.82665078796983, 113.6050170152778 34.88282716965251, 113.6427533999493 34.815514563935054, 113.5744551359546 34.77978979669181, 113.49008246199834 34.79650279857806, 113.49628347923993 34.82665078796983))'),4326),st_setsrid(st_point(113.619469,34.823071), 4326))
SELECT st_within(st_setsrid(st_point(113.619469,34.823071), 4326),st_setsrid(st_geomfromtext('POLYGON ((113.49628347923993 34.82665078796983, 113.6050170152778 34.88282716965251, 113.6427533999493 34.815514563935054, 113.5744551359546 34.77978979669181, 113.49008246199834 34.79650279857806, 113.49628347923993 34.82665078796983))'),4326))
--locaation为地理坐标表，digital_railing为电子围栏区域表
select * from locaation c,digital_railing d
where st_within(st_setsrid(st_point(c.lng,c.lat), 4326), d.geometry) and d.name = 'XXX'
参考书籍《Postgresql实战》
