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 where  MBRWithin(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 where  MBRWithin(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

空间地理位置参数不完整,造成该问题的原因可能有:
1、面状WKT字符串本身不完整、缺失内容。
2、面状WKT字符串的面要素中首尾点没有重合。
3、面状WKT字符串存在岛洞现象,或者自相交现象。

你可能感兴趣的:(Mysql,mysql,数据库)