普通计算两点之间距离：
st_distance(point(x1, y1), point(x2, y2))
如果需要转为米，需要 *111195，而且由于地球是球形，实际结果会有偏差

计算GEO坐标之间的距离（单位：米）：
st_distance_sphere(point(lon1, lat1), point(lon2, lat2))

注意，GEO 坐标顺序为 lon, lat

补充说明：
关于 point(x, y) 的参数，有的认为应该是 lat, lon，也有的人认为应该是 lon, lat。
如果是矩形的坐标系，两个参数的顺序其实没什么区别，但用在球面上就可能会产生一些误差了。
实际编程中，大多使用 LatLon 的顺序，但是我翻阅了 MySQL 官方文档，找到了这么一句话：

Point X and Y coordinates are interpreted as longitude and latitude

所以，正确的顺序应该是 lon（经度）,lat（纬度）

PostgreSQL 也支持这个功能，但是，需要安装 postgis 扩展，此扩展比 postgresql 本身还大，共70个依赖（postgresql 才 19个），直接安装 postgis 时会自动安装 postgresql。

安装完postgis后(Alpine命令：apk add postgis)，需要使用 postgres 账号执行操作，其他账号无权限，例如：

docker exec -ti -u postgres 容器ID psql
\c youdatabase
create extension postgis;

添加扩展后，会往当前库里面创建一个 spatial_ref_sys 表并插入 8500 行 GIS 配置数据。

调用也稍微不一样，方法是 st_distancesphere：

select st_distancesphere(point(0, 0)::geometry, point(1, 1)::geometry)

需要声明类型，否则会提示参数类型不符。

如果不想安装 postgis，也可以直接用 postgresql 内置的两点距离计算方法 point_distance(point(0, 0), point(1, 0)) * 104567（单位：米），由于两点角度不一样，104567 是个估算参数，例如，在横向测量可以使用 97955，纵向可以使用 111195，45度角可以使用101500。使用折中参数的话，距离计算误差基本在 10% 以内（5公里以内的商家列出来是 4.5 ~ 5.5 公里之间的高椭圆形）。

实际上，根据纬度的不同，横向测量参数也有变化，计算方法实际应该是：111195*cos(纬度*3.14159/180)，例如长沙地区是111195 * cos(28.2244 / 57.295827) ≈ 97974。

可以使用 gis 函数对比出该缩放比例：

select *, pt2/pt1 as scalerate from (
    select point_distance(point(116.420146, 39.897679), point(116.415122, 39.847263)) as pt1,
    st_distancesphere(point(116.420146, 39.897679)::geometry, point(116.415122, 39.847263)::geometry) as pt2
) a

也可以这样，少了平方和平方根的计算：

create or replace function distance_sphere
(in lon1 numeric, in lat1 numeric, in lon2 numeric, in lat2 numeric)
returns numeric return acos(
    sin(radians(lat1)) * sin(radians(lat2)) +
    cos(radians(lat1)) * cos(radians(lat2)) * cos(radians(lon1 - lon2))
) * 6371004