在地图中一个任意经,纬度下点击.想要找到以这个点为中心指定半径下包含多少个兴趣点.
查看postgis手册找到了这么一个函数(ST_distance(geometryg1,geometryg2);)
官方文档上是这么说的:
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手册geometryST_transform(geometryg1,integersrID);
第一个参数是原来的几何对像. 第二个参数为要把他转换到这个投影所代表的做标系下.
这时我们只要找一个单位是米的投影做标系把他转换过去就好了
例:
SELECT ST_distance(
ST_transform(ST_GeomFromText('POINT(-87.734087560562 43.770129071141)',4326),26986),
ST_transform(ST_GeomFromText('POINT(-87.747382933006 43.759234252055)',26986)
);
这个查出来的结果即是以米为单位的两点间的距离了
现在在说如何查找一定范围内的点
这里用到了postgis里的这第一个函数:booleanST_DWithin(geometryg1,geometryg2,double precisiondistance_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)',
ST_transform(t.geometry,1516)
ORDER BY ST_distance(GeomFromText('POINT(-87.71 43.741)',t.geometry);
结果:
24;"POINT(-87.718330082111 43.753078987035)"
17;"POINT(-87.726085716036 43.736952192682)"
18;"POINT(-87.726085716036 43.736952192682)"
找到了三个点
总结以上是内存溢出为你收集整理的用postgresql(postgis) 查找指定范围内的兴趣点全部内容,希望文章能够帮你解决用postgresql(postgis) 查找指定范围内的兴趣点所遇到的程序开发问题。
如果觉得内存溢出网站内容还不错,欢迎将内存溢出网站推荐给程序员好友。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)