用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)"
找到了三个点
用postgresql(postgis) 查找指定范围内的兴趣点相关推荐
- 根据GPS经纬度查找指定范围内的对象
项目上需要根据当前的经纬度查询指定范围内离这个坐标最近的对象. 主要是思路.请大家指正 public static void CheckGps() { double temp = ...
- linux查找某时间段文件夹,shell查找指定时间段内的文件
#!/bin/bash #20170905 输入参数格式 echo "显示"$1"的备份文件" date_0=$1 date_1=`expr $date_0 + ...
- PostGIS 查询某点周围指定范围内的兴趣点
问题描述 在我们生活中,想必很多人都使用过一个功能就是查找附近XX米内的美食.景点等信息.那么该功能是如何实现的呢?本文将带你一探究竟. 具体问题:查询给定点(如:113.678 34.796)周围1 ...
- linux shell find命令 查找指定时间范围内的文件
例如查找`2013-08-08`到`2013-09-01`号之间的文件,使用如下命令即可: find /log/ -name 'production.log-2013*' -newermt '2013 ...
- linux查找特定日期之后的文件,Linux Find命令查找指定时间范围内的文件的例子
例如查找`2013-08-08`到`2013-09-01`号之间的文件,使用如下命令即可: 复制代码 代码如下: find /log/ -name 'production.log-2013*' -ne ...
- linux 指定范围内查找文件,Linux Find命令查找指定时间范围内的文件的例子
例如查找`2013-08-08`到`2013-09-01`号之间的文件,使用如下命令即可: 复制代码 代码如下: find /log/ -name 'production.log-2013*' -ne ...
- linux中使用egrep取值ip地址,egrep命令_Linux egrep 命令用法详解:在文件内查找指定的字符串...
egrep命令用于在文件内查找指定的字符串.egrep执行效果与grep -E相似,使用的语法及参数可参照grep指令,与grep的不同点在于解读字符串的方法.egrep是用extended regu ...
- VBA查找单元格内指定文字并标红加粗
VBA查找单元格内指定文字并标红加粗 功能: 选中想要查找的单元格,在函数内点击运行,实现将指定的文字标红加粗的功能 bug: 指定文字会区分大小写,如何让它不区分大小写? 改进 将是否加粗,和标记的 ...
- 使用simple transformation查找xml file内某个节点的attribute是否存在指定value
Created by Jerry Wang on Jun 05, 2014 下列report实现通过simple transformation查找xml 文件内下列路径的节点ds其attribute ...
最新文章
- 人脸检测对齐--Joint Face Detection and Alignment using Multi-task Cascaded Convolutional Networks
- C#学习笔记——密封类与密封方法
- spring中的quartz调度问题
- PHP之preg_replace()与ereg_replace()正则匹配比较讲解
- Centos之压缩和解压缩命令
- python怎么输出小数部分_python 输出小数控制
- 会linux基本命令是脚本语言吗,如何理解Linux Shell和基本Shell脚本语言?
- T4生成多文件时,不生成自己
- fork、getpid、getppid函数
- leetcode 633. 平方数之和(双指针)
- Kaggle数据竞赛——伪标签的使用
- java中finaljava中this_Java中this,static,final,const用法详解
- 修改服务器的AJP监听地址,修改服务器的AJP监听地址
- 前后端交互总结2:使用PHP进行表单数据上传与更新
- 2017.3.18 PPT汇报--总结
- PS怎么把人物扣的更干净_PS抠图技巧
- 挑战微软 + GitHub!谷歌联手 Replit,升级 AI 编程“神器”:曾拒绝微软 10 亿美元的收购...
- 移动端vue仿朋友圈项目总结
- Apple苹果产品MFi设计及标准汇总
- 1到n的数按字典序排序