最近在做一个项目，有一个功能想要实现类似于查询附近的人的功能。由于项目的原因数据库只能使用 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)')
数据就不用转来转去了。
只是项目已有数据，不能这样做了。