最近在做一个项目,有一个功能想要实现类似于查询附近的人的功能。由于项目的原因数据库只能使用 postgresql
,空间查询就使用了 postgis
来实现。
具体业务像这样:业务需要返回附近距自己 1000
米的人的一个列表,返回列表中要带上距离,就是说某人距你多少米?
使用定位方式是什么?
我们使用 gps
定位发现定位精度,时间,误差和抖动都需要发大量的时间去处理。
于是采用了百度的定位 SDK
,百度的上面几处都处理不错。
有一个小问题,从百度拿到的坐标是 经度和纬度,这个坐标系应该是地理坐标,由于也是使用 gps
,其基础坐标应该是 wgs84
。所以拿到的值应该是:wgs84
+ 国家偏移 + 百度偏移 = 百度坐标
后面两项是没有太好的办法解决,我们拿百度坐标与 google
的大致对比了一下,误差 50
米以内,业务可以接受。
而且这个后面两项是系统偏移,不是随机的,这样的话,计算距离这种相对计算,应该可以消除部分系统偏移的影响。
再说说postgis
计算距离
postgis
是一个地理空间查询引擎,计算地理距离这很方便。
- 开始我们在表中建一列用于存地理类型的
geometry
。
SELECT AddGeometryColumn ('public','chest','position',4326,'POINT',2, false);
- 然后保存的时候,将 经纬度和 经纬度变成
geometry
进行保存,也就是三个字段,其中保存geometry
是用下面的sql
:
ST_GeomFromText('POINT(${lng} ${lat})', 4326)
sql
中的 4326
是指的wsg84
的系统。
- 然后对保存后的值进行距离查询,使用
ST_distance
进行,后面有二个坑等着我们:
第一个:geometry
的问题
postgis
中 geomery
使用 ST_distance
计算出来的单位竟然是弧度
,不是 米
,找了一段时间问题,发现是没有使用投影坐标引起的。
于是在计算的时候,先将 点
变成 投影坐标
,再来计算。
ST_Transform(geo, 2346)
这里的 2346
是中国西安 80 高斯克里格投影
的编号,是分带的,使用的中央经线在南京。
加上这个,我们能正确计算出来距离,也是 米
。
然后,一段时间后,把应用发到别的城市,也就是苏州,出现了问题
第二个: 投影带号的问题
原因是苏州不在 2346
中的带号里面,计算就出错了!
这样就有问题了,全国这么多的带号,不可能写在程序里,动态计算啊?再说万一出国了,没有中国的投影怎么办?
回头再去找方案,这个计算距离应该要支持最低就是全国,最好是全世界。
还真找到了,就是 geography
。
解释一下: geography
和 geometry
都是 postgis
中的数据类型,翻绎过来:地理图形和几何图形,用 postgis
的说法就是, geography
就是使用 wgs84
坐标系 的图形。
postgis
可以直接使用 geography
进行计算,支持全球。于是这个问题,终于有解决方案了,就是把我们的 geometry
转成 geography
,然后进行距离计算。
下面是我们使用 sql 语句
select
(ST_distance(position::geography, ST_GeomFromText('POINT(121 32)', 4326)::geography) as distance
from chest
where ST_dwithin(position::geography, ST_GeomFromText('POINT(121 32})', 4326)::geography, 1000)
orderBy distance desc
其中:position
是我存在数据库中的字段,使用的是 geometry
的类型, (121,32 )
是我输入的中心点。
::geography
就是postgis 中的 转换类型语法,把 geometry
转成 geography
,然后计算,这样全球都能使用。
总结
这些问题的解决,算是把以前的知识在复习了一下。 在项目中,我们也可以直接把坐标存成 geography
,使用如下的 sql
:
ST_GeographyFromText('SRID=4326;POINT(-110 30)')
数据就不用转来转去了。
只是项目已有数据,不能这样做了。