PostGIS实现圆形检索和矩形检索

找到圆心 POINT(1000 1000) 100米范围内的查询以下:

SELECT * FROM geotable
WHERE ST_DWithin(geocolumn, ‘POINT(1000 1000)’, 100.0);

又或者你须要的是矩形检索：ST_MakeBox2D(geometry pointLowLeft, geometry pointUpRight);

SELECT feature_id, feature_name, the_geom
FROM features
WHERE the_geom && ST_SetSRID(ST_MakeBox2D(ST_Point(-989502.1875, 528439.5625),
ST_Point(-987121.375 ,529933.1875)),2163)

 如下是详细的解释：

在地图中一个任意经,纬度下点击.想要找到以这个点为中心指定半径下包含多少个兴趣点.

查看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)"

找到了三个点