之前用的是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;
搜索结果
以上为个人经验,希望能给大家一个参考,也希望大家多多支持。如有错误或未考虑完全的地方,望不吝赐教。