Mysql|mysql 5.7.28 地理位置计算详解


mysql 5.7.28 空间地理位置计算

  • 1. 创建测试表
  • 2. 插入数据
  • 3. 计算两个坐标点之间的距离
  • 4. 计算指定的坐标点是否位于多边形(面)中
      • 4.1 第一种方式
      • 4.2 第二种方式
  • 5. 计算两个多边形(面)是否相交 (st_intersects )
  • 6. 计算两个多边形(面)是否包含(st_contains )
  • 7. 计算一个多边形(面)是否在另一个多边形(面)内(st_within)
  • 8. 遇到的异常
          • 1. MySQL: Invalid GIS data provided to function st_geometryfromtext

1. 创建测试表 创建polygon表,该表用于存储多边形(面)数据:
CREATE TABLE `polygon` ( `id` int(10) unsigned NOT NULL AUTO_INCREMENT, `polygongeo` polygon NOT NULL, `name` varchar(255) DEFAULT NULL, PRIMARY KEY (`id`) ) ENGINE=MyISAM AUTO_INCREMENT=2 DEFAULT CHARSET=utf8;

创建spots表,该表用于存储坐标点数据:
CREATE TABLE `spots` ( `id` int(11) NOT NULL, `name` varchar(255) DEFAULT NULL COMMENT '地点名称', `spot` point DEFAULT NULL COMMENT '经纬度点', PRIMARY KEY (`id`) ) ENGINE=InnoDB DEFAULT CHARSET=utf8;

2. 插入数据 polygon表:
INSERT INTO polygon ( polygongeo, NAME ) VALUES( POLYGONFROMTEXT ( 'POLYGON((36.6882437874 117.1774291992,36.6915474355 117.2364807129,36.6607078682 117.1499633789,36.6882437874 117.1774291992))' ), '济南市历下区第一区域' ); INSERT INTO polygon ( polygongeo, NAME ) VALUES( POLYGONFROMTEXT ( 'POLYGON((36.8588150600 117.0421600300, 36.8565487700 117.2181129500, 36.7488396900 117.1599519300, 36.7642261500 117.0372676800, 36.8588150600 117.0421600300))' ), '济南市历城区第一区域' );

spots表:
INSERT INTO spots ( `name`, `spot` ) VALUES('山东省济南市济南高新技术产业开发区遥墙街道马官寨村',ST_GeomFromText ( 'POINT(36.8370421400 117.1580314600)' )); INSERT INTO spots ( `name`, `spot` ) VALUES('山东省济南市济南高新技术产业开发区舜华路街道齐鲁软件园',ST_GeomFromText ( 'POINT(36.6718611000 117.1339130400)' ));

3. 计算两个坐标点之间的距离
select a.name,b.name,truncate(st_distance(a.spot,b.spot) * 111195 / 1000,2) as distance from (select * from spots where name='山东省济南市济南高新技术产业开发区舜华路街道齐鲁软件园') a, (select * from spots where name='山东省济南市济南高新技术产业开发区遥墙街道马官寨村') b

以上sql以spots表中的两个坐标点为起始点,使用st_distance函数进行两点之间的距离计算。结果如下:
地点一 地点二 距离(km)
山东省济南市济南高新技术产业开发区舜华路街道齐鲁软件园 山东省济南市济南高新技术产业开发区遥墙街道马官寨村 18.56
备注:mysql内置的 st_distance 函数计算的结果单位是度,需要乘111195(地球半径6371000*PI/180)是将值转化为米,除以1000换算为千米
4. 计算指定的坐标点是否位于多边形(面)中 4.1 第一种方式
利用mysql内置函数 MBRWithin(点,面) 计算,具体如下:
select AsText(p.polygongeo) from polygon p,(select spot from spots where name='山东省济南市济南高新技术产业开发区遥墙街道马官寨村') s whereMBRWithin(s.spot,p.polygongeo)

结果如下:
搜索结果:
POLYGON((36.85881506 117.04216003,36.85654877 117.21811295,36.74883969 117.15995193,36.76422615 117.03726768,36.85881506 117.04216003))
以上结果表示该坐标点位于多边形(面)内
执行如下sql:
select AsText(p.polygongeo) from polygon p,(select spot from spots where name='山东省济南市济南高新技术产业开发区舜华路街道齐鲁软件园') s whereMBRWithin(s.spot,p.polygongeo)

结果如下:
搜索结果:
POLYGON:null
以上结果表示该坐标点位于多边形(面)外
4.2 第二种方式
利用mysql内置函数 MBRContains(面,点) 计算,具体如下:
SET @g1 = geomFromText('POLYGON((121.474243 31.234504,121.471775 31.233348,121.470724 31.23155,121.471603 31.230229,121.472655 31.230357,121.475777 31.232045,121.474243 31.234504))'); SET @g2 = geomFromText('POINT(121.471775 31.233348)'); select MBRContains(@g1,@g2);

结果如下:
结果: 1
结论:该坐标点位于该多边形内
备注:如果返回值为 0 ,说明该坐标点不在多边形内,返回值为 1 ,说明坐标点位于多边形内。
5. 计算两个多边形(面)是否相交 (st_intersects ) 利用mysql内置函数 st_intersects 判断两个空间地理位置多边形(面)是否相交,具体如下:
SET @g1 = geomFromText('POLYGON((121.474243 31.234504,121.471775 31.233348,121.470724 31.23155,121.471603 31.230229,121.472655 31.230357,121.475777 31.232045,121.474243 31.234504))'); SET @g2 = geomFromText('POLYGON((121.474243 31.234804,121.471775 31.233948,121.471724 31.23155,121.471903 31.230229,121.472655 31.230157,121.475777 31.231045,121.474243 31.234804))'); SELECT st_intersects(@g1,@g2);

结果如下:
结果: 1
结论:两个多边形相交
备注:如果返回值为 0 ,说明两个多边形未相交,返回值为 1 ,说明两个多边形相交。
6. 计算两个多边形(面)是否包含(st_contains ) 利用mysql内置函数 st_contains 进行空间地理位置多边形(面)包含计算,具体如下:
SET @g1 = geomFromText('POLYGON((121.474243 31.234504,121.471775 31.233348,121.470724 31.23155,121.471603 31.230229,121.472655 31.230357,121.475777 31.232045,121.474243 31.234504))'); SET @g2 = geomFromText('POLYGON((121.474243 31.234504,121.471775 31.233348,121.470724 31.23155,121.471603 31.230229,121.474243 31.234504))'); SELECT st_contains (@g1,@g2);

结果如下:
结果:1
结论:多边形@g1 完全包含 @g2
备注:该函数的参数区分先后顺序,即结果为0:参数一 不包含 参数二,结果为1:参数一完全包含参数二
7. 计算一个多边形(面)是否在另一个多边形(面)内(st_within) 利用mysql内置函数 st_within进行空间地理位置多边形(面)包含计算,具体如下:
SET @g1 = geomFromText('POLYGON((121.474243 31.234504,121.471775 31.233348,121.470724 31.23155,121.471603 31.230229,121.472655 31.230357,121.475777 31.232045,121.474243 31.234504))'); SET @g2 = geomFromText('POLYGON((121.474243 31.234504,121.471775 31.233348,121.470724 31.23155,121.471603 31.230229,121.474243 31.234504))'); SELECT st_within(@g1,@g2);

结果如下:
结果:1
结论:多边形@g1 在 @g2 内
备注:该函数的参数区分先后顺序,即结果为0:参数一 不在 参数二内,结果为1:参数一在参数二内
8. 遇到的异常 1. MySQL: Invalid GIS data provided to function st_geometryfromtext
【Mysql|mysql 5.7.28 地理位置计算详解】空间地理位置参数不完整,造成该问题的原因可能有:
1、面状WKT字符串本身不完整、缺失内容。
2、面状WKT字符串的面要素中首尾点没有重合。
3、面状WKT字符串存在岛洞现象,或者自相交现象。

    推荐阅读