postgis 中的距离计算

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

数据就不用转来转去了。

只是项目已有数据，不能这样做了。

参考文档

PostGIS-04.md

Geography