postgis和地理围栏

postgis

postgis是pg数据库的一个插件，除原数据类型外(int varchar)、新增了空间数据类型(geography和geometry)。比如我们新建一张道路表road(字段有名称varchar、建设时间timestamp、地理位置geometry)，可以将道路名字、建设时间存进去，同时把道路的位置坐标和使用的坐标系存到一个geometry类型字段。再建一张建筑表存建筑位置，这样就能利用postgis提供的函数分析马路附近有哪些建筑这类问题。

geometry是几何类型，基础是平面，平面上两点之间的最短路径是直线。几何函数（面积、距离、长度、交点等）是使用线性向量和笛卡尔平面计算的。这简化了实施并运行得更快，但对于地球球面以上的数据是不准确的。它可以存地理坐标系数据，只是计算大范围计算精度有局限，更适合选择合适的投影坐标系存平面坐标。

geography地理数据类型基于球面模型。球体上两点之间的最短路径是一个大圆弧。地理上的函数（面积、距离、长度、交点等）是使用球体上的弧计算的。通过考虑世界的球体形状，函数可提供更准确的结果。

人员设备demo表

用4326二维地理坐标系(只有经纬度)，制作一份地图，导入pg，表结构如下

CREATE TABLE public.demo_building (

id int8 NOT NULL,

geom public.geometry(multipolygonzm, 4326) NULL,

"name" varchar(10) NULL,

CONSTRAINT demo_building_pkey PRIMARY KEY (id)

);

CREATE TABLE public.demo_person (

id int8 NOT NULL,

geom public.geometry(pointzm, 4326) NULL,

"name" varchar(10) NULL,

CONSTRAINT demo_person_pkey PRIMARY KEY (id)

);

查询设备最近的人员

select

person ."name" ,

ST_Distance(building.geom, person.geom)

from

demo_building building,

demo_person person

where

building."name" = '设备1'

ORDER BY ST_Distance(building.geom, person.geom) asc

-- 基于索引的KNN查询

select

person ."name"

from

demo_person person

order by person.geom <-> (select building.geom from demo_building building where building."name" = '设备1')

qgis使用“要素之间最短线”,源图层为building，目标图层为person，可绘制最短线，打开属性表可看距离

判断人员是否在围栏内

--距离原点指定0.5范围内的人员，圆形围栏

select

person.geom ,

person ."name" ,

person .geom <-> 'SRID=4326;POINT Z(-0.764265422077925 -0.0342836850649348 0)'::geometry as dist

from

demo_person person

where person .geom <-> 'SRID=4326;POINT Z(-0.764265422077925 -0.0342836850649348 0)'::geometry <= 0.5

--在指定多边形范围内的人，在矿工作区内的人,多边形围栏

select

*

from

demo_person person

where

ST_Contains((select building .geom from demo_building building where building ."name" = '矿工作区')::geometry,person.geom);

--在指定多边形边上的人，在高坡工作区边上的人，围栏边上

select

*

from

demo_person person

where

ST_Intersects((select building .geom from demo_building building where building ."name" = '高坡工作区')::geometry,person.geom);

聚合区域人员数量

--查询指定范围内的人员，在矿工作区内的人员数量

select count(1) from demo_person person where ST_Within(person.geom,(select building .geom from demo_building building where building ."name" = '矿工作区')::geometry);