用postgresql(postgis) 查找指定范围内的兴趣点
在地图中一个任意经,纬度下点击.想要找到以这个点为中心指定半径下包含多少个兴趣点.
查看postgis手册找到了这么一个函数(ST_Distance(geometry g1, geometry g2);)
官方文档上是这么说的:
ST_Distance — For geometry type Returns the 2-dimensional cartesian minimum distance (based on spatial ref) between two geometries in projected units. For geography type defaults to return spheroidal minimum distance between two
geographies in meters.
即给定两个空间点计算两点间的距离.计算结果的单位与你的空间数据的参考系有关.
如果你使用的是4326(wgs84)这个坐标系的话他是以度为单位的.要想转成米为单位的话还得做一下转换.下面会提到:
GEOCS代表的是地理坐标系，也就是以经纬度表示的坐标系统,例如4326
PROJCS代表的投影坐标系，它是通过一种算法把球面坐标系转成平面坐标系，以便计算，一般是以米为单位表示, 例如26986
因此，在求两点之间的距离时，由于存的数据都是经纬度，因此它参考的是GEOCS，要想得到以米为单位的结果，首先要把它转成PROJCS，可以通过ST_Transform来实现
查看postgis手册 geometry ST_Transform(geometry g1, integer srid);
第一个参数是原来的几何对像. 第二个参数为要把他转换到这个投影所代表的做标系下.
这时我们只要找一个单位是米的投影做标系把他转换过去就好了
例:
SELECT ST_Distance(
ST_Transform(ST_GeomFromText('POINT(-87.734087560562 43.770129071141)',4326),26986),
ST_Transform(ST_GeomFromText('POINT(-87.747382933006 43.759234252055)', 4326),26986)
);
这个查出来的结果即是以米为单位的两点间的距离了
现在在说如何查找一定范围内的点
这里用到了postgis里的这第一个函数:boolean ST_DWithin(geometry g1, geometry g2, double precision distance_of_srid);
第一个参数为参考对像.第二个参数为目标对像.第三个参数为距离(同样如果是地理坐标系单位是度.投影做标系单位是米)
即以g1为中心,半径为distance_of_srid,这个范围内包不包含g2,如果包含反回true,否则即为假
如了给一个完整的例子.查找以(-87.71 43.741)为中心半径1516米范围内的兴趣点,之后按与这个中心点由近到远的顺序排列结果
SELECT t.feat_id,astext(t.geometry) FROM gis_site t
WHERE ST_DWithin(
ST_Transform(GeomFromText('POINT(-87.71 43.741)',4326),26986),
ST_Transform(t.geometry,26986), 1516)
ORDER BY ST_Distance(GeomFromText('POINT(-87.71 43.741)',4326), t.geometry);
结果:
24;"POINT(-87.718330082111 43.753078987035)"
17;"POINT(-87.726085716036 43.736952192682)"
18;"POINT(-87.726085716036 43.736952192682)"
找到了三个点