postgresql 计算距离的实例(单位直接生成米)

数据库 2024/12/27 佚名

3 2 1

之前用的是ST_Distance 函数，但是貌似需要进行一次单位的转换，而且网上有说那种转换不是特别准确，现在暂时将该算法记录在此：

select st_distance(ST_GeomFromText('POINT(120.451737 36.520975)',900913),ST_GeomFromText('POINT(120.455636 36.520885)',900913))*60*1.852;

这里的计算方式倒是可以换坐标系，但是，测试了两个坐标系都没有起作用。而且该种方式转换过单位后跟arcgis计算出的结果相差甚远，最终决定使用下面的方式；

今天发现了另外一种方式来计算距离，这种方式可以直接生成单位为米的结果：

select ST_Length(Geography(ST_GeomFromText('LINESTRING(120.451737 36.520975,120.455636 36.520885)')));

这种方式的不便在于：

1.要把点转换成线或者其他的图形而不是点；

2.geography函数现在只支持4326坐标系，不能换成其他的。

追加：

上面的方式是计算点到点的距离，但是如果要想知道某一个点到某条线的距离是不是在某个范围内，又该如何计算呢；如下：

select ST_Contains(St_Astext(ST_Buffer(geography(geomfromtext('MULTILINESTRING((线的坐标点))')),25.00{以米为单位的距离})),st_astext(geography(geomfromtext('POINT(121.37805 37.54142)')))) as result

使用类似上面的方式，就可以输入以米为单位的距离判断某个点是否在某个距离范围内；

补充：postgresql 搜索指定距离内的记录按近到远排序并返回距离

实例如下：

CREATE TABLE mylocation ( 
 id SERIAL PRIMARY KEY,
 geom GEOMETRY(Point, 4326),
 name VARCHAR(128),
 x double precision,
 y double precision
); 
 
INSERT INTO mylocation (geom,name,x,y) VALUES (
 ST_GeomFromText('POINT(0.0001 0)', 4326),'zhangsan',0.0001,0
);
INSERT INTO mylocation (geom,name,x,y) VALUES (
 ST_GeomFromText('POINT(0.001 0)', 4326),'zhangsan',0.001,0
);
INSERT INTO mylocation (geom,name,x,y) VALUES (
 ST_GeomFromText('POINT(0.001 0)', 4326),'zhangsan',0.001,0
);
INSERT INTO mylocation (geom,name,x,y) VALUES (
 ST_GeomFromText('POINT(0.1 0)', 4326),'zhangsan',0.1,0
); 
 
SELECT id, name,geom,x,y,  ST_DistanceSphere(
           geom,
           ST_GeometryFromText('POINT(0 0)')) distance
FROM mylocation
WHERE ST_DWithin(
 geom, 
 ST_GeomFromText('POINT(0 0)', 4326),
 0.001
)ORDER BY distance asc;;

查询语句下面距离单位为m

SELECT id, name,geom,x,y,  ST_DistanceSphere(
           geom,
           ST_GeometryFromText('POINT(0 0)')) distance
FROM mylocation
WHERE ST_DWithin(
 geom::geography, 
 ST_GeomFromText('POINT(0 0)', 4326)::geography,
 1000
) ORDER BY distance asc;

搜索结果