使用MySQL空间扩展来select圆内的点

我有一个名为flags的表,其中包含一个称为coordinates的列,充满了MySQL的“点”。 我需要执行一个查询,根据100米半径的经度和纬度位置获得一个圆圈内的所有标志。

从使用的angular度来看,这是基于用户的位置。 例如,手机会给用户的经纬度位置,然后将其传递给API的这一部分。 然后由API在用户周围创build一个半径为100米的不可见圆圈,然后返回此圆圈中的标志。

这是API的这一部分我不知道如何创build,因为我不确定如何使用SQL来创build这个不可见的圆,并只在这个半径内select点。

这可能吗? 有没有一个MySQL的空间function ,可以帮助我做到这一点?

我相信Buffer()函数可以做到这一点,但我找不到任何文档如何使用它(例如SQL)。 理想情况下,我需要一个答案,告诉我如何使用这个function或最接近它。 在将这些坐标存储为地理空间点的地方,我应该使用地理空间函数来做我要求的最大效率。

标志表:

  • ID
  • 坐标
  • 名称

示例行:

1 | [几何 – 25B] | Tenacy AB

旗标表有经纬度位置和东,北(UTM)

用户的位置只是标准的纬度/经度,但我有一个库可以将此位置转换为UTM

MySQL中没有支持纬度/经度距离计算的地理空间扩展函数。

你要求在地球表面接近圆。 在你的问题中,你提到你的flags表中的每行都有纬度/长度的值,还有在几个不同的UTM区域之一中的通用横向墨卡托 (UTM)投影值。 如果我正确记住我的英国军械测量地图,UTM对于在这些地图上查找物品很有用。

计算UTM中相同区域中两点之间的距离是一件简单的事情:笛卡尔距离就是一个窍门。 但是,当点位于不同区域时,该计算不起作用。

因此,对于您的问题中描述的应用程序,有必要使用大圆距离 ,这是使用海峡或其他合适的公式计算的。

MySQL增加了地理空间扩展,支持将各种平面形状(点,多义线,多边形等)表示为几何图元的方法。 MySQL 5.6实现了一个未st_distance(p1, p2)距离函数st_distance(p1, p2) 。 但是,这个函数返回笛卡尔距离。 所以这是完全不适合基于经度和纬度的计算。 在温带地区,由于纬度线在两极之间越来越近,所以一度纬度的纬度(南北)几乎是东经(东西)的两倍。

所以,一个圆形的近似公式需要使用真正的经度和纬度。

在你的应用程序中,你可以在给定的点的十个法定英里内find所有的flagslatpoint,longpoint

  SELECT id, coordinates, name, r, units * DEGREES( ACOS( COS(RADIANS(latpoint)) * COS(RADIANS(X(coordinates))) * COS(RADIANS(longpoint) - RADIANS(Y(coordinates))) + SIN(RADIANS(latpoint)) * SIN(RADIANS(X(coordinates))))) AS distance FROM flags JOIN ( SELECT 42.81 AS latpoint, -70.81 AS longpoint, 10.0 AS r, 69.0 AS units ) AS p ON (1=1) WHERE MbrContains(GeomFromText ( CONCAT('LINESTRING(', latpoint-(r/units),' ', longpoint-(r /(units* COS(RADIANS(latpoint)))), ',', latpoint+(r/units) ,' ', longpoint+(r /(units * COS(RADIANS(latpoint)))), ')')), coordinates) 

如果要search20公里范围内的点,请更改查询的这一行

  20.0 AS r, 69.0 AS units 

对此,例如

  20.0 AS r, 111.045 AS units 

r是要search的半径。 units是地球表面纬度的距离单位(英里,公里,弗隆,无论你想要的)。

这个查询使用边界纬度/长度与MbrContains一起排除绝对离开始点太远的点,然后使用大圆距离公式来生成剩余点的距离。 所有这些的解释可以在这里find 。 如果您的表使用MyISAM访问方法并具有空间索引,则MbrContains将利用该索引来快速search。

最后,上面的查询select矩形内的所有点。 要缩小到只有圆圈中的点,并按照近似值sorting,请按如下方式包装查询:

  SELECT id, coordinates, name FROM ( /* the query above, paste it in here */ ) AS d WHERE d.distance <= dr ORDER BY d.distance ASC 

这假定表中的坐标作为POINT()数据types存储在标记为“点”的列中。 函数X(point)和Y(point)分别从点值中提取纬度和经度值。

 SET @lat = the latitude of the point SET @lon = the longitude of the point SET @rad = radius in Kilometers to search from the point SET @table = name of your table SELECT X(point),Y(point),*, ( 6373 * acos ( cos ( radians( @lat ) ) * cos( radians( X(point) ) ) * cos( radians( Y(point) ) - radians( @lon ) ) + sin ( radians( @lat ) ) * sin( radians( X(point) ) ) ) ) AS distance FROM @table HAVING distance < @rad 

如果你想以英里来做,用3959replace常量6373

对于那些想要减less查询语法的用户,下面是一个用户定义的MySQL函数的常见实现,用于实现基于Haversine公式的距离函数。

 CREATE FUNCTION HAVERSINE ( coord1 POINT, coord2 POINT ) RETURNS DOUBLE DETERMINISTIC BEGIN DECLARE dist DOUBLE; SET rlat1 = radians( X( coord1 ) ); SET rlat2 = radians( X( coord2 ) ); SET rlon1 = radians( Y( coord1 ) ); SET rlon2 = radians( Y( coord2 ) ); SET dist = ACOS( COS( rlat1 ) * COS( rlon1 ) * COS( rlat2 ) * COS( rlon2 ) + COS( rlat1 ) * SIN( rlon1 ) * COS( rlat2 ) * SIN( rlon2 ) + SIN( rlat1 ) * SIN( rlat2 ) ) * 6372.8; RETURN dist; END 

UPDATE

使用ST_Distance_Sphere()使用lat / long来计算距离

http://dev.mysql.com/doc/refman/5.7/en/spatial-convenience-functions.html#function_st-distance-sphere

因为缓冲区是一个多边形,MySQL <5.6中的多边形操作被实现为“最小边界矩形”(MBR),所以在MySQL <5.6中缓冲区并不会有多大帮助。

自MySQL 5.6以来, 完整的非MBR st_*操作得以实现 。 但是对于你来说,最好的解决scheme就是使用无证函数st_distance

 select * from waypoints where st_distance(point(@center_lon, @center_lat), coordinates) <= radius; 

很难find,因为它没有logging:-)但是在这个博客上提到,他的作者也填写了所提到的bugreport 。 有一些警告(援引博客):

坏消息是:

1)所有function仍然只使用平面系统坐标。 不支持不同的SRID。

2)空间索引(RTREE)仅支持MyISAM表格。 可以使用InnoDB表的函数,但不会使用空间键。

点1)表示距离的单位将与坐标单位相同(在WGS84的情况下为度数)。 如果需要以米为单位的距离,则必须使用具有对应于米的单位的预计协调系统(例如UTM或类似的)。

所以,如果你不想要这些警告,或者如果MySQL <5.6,你将不得不编写自己的自定义距离函数。

希望我的版本帮助

 SELECT * FROM `locator` WHERE SQRT(POW(X(`center`) - 49.843317 , 2) + POW(Y(`center`) - 24.026642, 2)) * 100 < `radius` 

详情在这里http://dexxtr.com/post/83498801191/how-to-determine-point-inside-circle-using-mysql

来自: https : //gis.stackexchange.com/questions/31628/find-points-within-a-distance-using-mysql

 SELECT id, ( 6371 * acos ( cos ( radians(78.3232) ) * cos( radians( lat ) ) * cos( radians( lng ) - radians(65.3234) ) + sin ( radians(78.3232) ) * sin( radians( lat ) ) ) ) AS distance FROM markers HAVING distance < 30 ORDER BY distance LIMIT 0 , 20; 

(记住要replace所有的常量,这个例子是以公里为单位的)

您可以使用:

 SELECT name, lat, lng FROM vw_mytable WHERE ST_Contains(ST_Buffer( ST_GeomFromText('POINT(12.3456 34.5678)'), (0.00001*1000)) , mypoint) = 1 

上面里面的expression式:0.00001 * 1000给你一个直径1000米的圆,它应用在一个视图上,名称列只是一个标签,点mypoint是我的点列的名称,lat是在里面计算的用ST_X(mytable.mypoint)查看,用ST_Y(mytable.mypoint)查看lng,他们只是显示lat和lng的文字值。 它会给你所有属于圆的坐标。