MySQL 存储地理信息通常使用 GEOMETRY 数据类型或其子类型（如 POINT, LINESTRING, POLYGON 等）。为了支持这些数据类型，MySQL 提供了 SPATIAL 索引，这允许我们执行高效的地理空间查询。
1. 创建支持地理信息的表
首先，我们需要一个包含 GEOMETRY 或其子类型列的表。以下是一个示例，展示如何创建一个包含 POINT 类型的表：
CREATE TABLE locations (
id INT AUTO_INCREMENT PRIMARY KEY,
name VARCHAR(255) NOT NULL,
position POINT NOT NULL,
SPATIAL INDEX(position) -- 为位置列创建空间索引
) ENGINE=InnoDB;
2. 插入地理信息数据
我们可以使用 GeomFromText() 或 PointFromText() 函数插入地理数据。以下是如何插入一个点的示例：
INSERT INTO locations (name, position)
VALUES ('Location A', GeomFromText('POINT(10 20)'));
-- 或者使用 PointFromText
INSERT INTO locations (name, position)
VALUES ('Location B', PointFromText('POINT(30 40)'));
3. 查询地理信息数据
我们可以使用 MBRContains(), Distance_Sphere(), ST_Distance_Sphere() 等函数来查询地理数据。以下是一些示例：
3.1查找指定矩形区域内的位置
-- 查找位置在 (0, 0) 到 (20, 20) 矩形区域内的所有位置
SELECT * FROM locations
WHERE MBRContains(
GeomFromText('POLYGON((0 0, 20 0, 20 20, 0 20, 0 0))'),
position
);
3.2查找距离特定点一定距离内的位置
注意：这里使用了 Distance_Sphere() 函数，它基于地球是完美球体的假设。对于更精确的计算，我们可以使用 ST_Distance_Sphere() 并指定地球半径。
-- 查找距离 (15, 15) 点 10 公里内的所有位置
-- 假设地球半径为 6371 公里（平均半径）
SELECT *, (6371 * acos(cos(radians(15))
* cos(radians(X(position)))
* cos(radians(Y(position)) - radians(15))
+ sin(radians(15))
* sin(radians(X(position))))) AS distance_km
FROM locations
HAVING distance_km < 10;
3.3使用 ST_Distance_Sphere() 查找距离
这是一个更精确的距离计算示例，它使用 ST_Distance_Sphere() 函数并指定地球的平均半径。
-- 查找距离 (15, 15) 点 10 公里内的所有位置
SELECT *, ST_Distance_Sphere(point(15, 15), position, 6371) AS distance_km
FROM locations
HAVING distance_km < 10;
注意：上述查询中的距离计算是基于 Haversine 公式的简化版本，它假设地球是一个完美的球体。在实际应用中，我们可能需要使用更复杂的算法来考虑地球的不规则形状。
此外，我们还可以使用 MySQL 的其他地理空间函数和操作符来执行更复杂的地理空间查询和操作。
4.查询地理信息进阶示例
我们可以探讨一个更复杂的示例，该示例涉及POLYGON地理数据类型，并使用ST_Contains函数来检查一个点是否位于多边形内部。同时，我们也会使用ST_Distance_Sphere函数来计算点与多边形中心点的距离。
4.1创建表并插入数据
首先，我们创建一个包含POLYGON列的表，并插入一些多边形数据。
CREATE TABLE polygons (
id INT AUTO_INCREMENT PRIMARY KEY,
name VARCHAR(255) NOT NULL,
shape POLYGON NOT NULL,
SPATIAL INDEX(shape)
) ENGINE=InnoDB;
INSERT INTO polygons (name, shape)
VALUES ('Polygon A', GeomFromText('POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))'));
INSERT INTO polygons (name, shape)
VALUES ('Polygon B', GeomFromText('POLYGON((20 20, 30 20, 30 30, 20 30, 20 20))'));
-- 创建一个包含点的表
CREATE TABLE points (
id INT AUTO_INCREMENT PRIMARY KEY,
name VARCHAR(255) NOT NULL,
position POINT NOT NULL,
SPATIAL INDEX(position)
) ENGINE=InnoDB;
INSERT INTO points (name, position)
VALUES ('Point 1', GeomFromText('POINT(5 5)'));
INSERT INTO points (name, position)
VALUES ('Point 2', GeomFromText('POINT(25 25)'));
4.2查询点是否在多边形内部，并计算距离
现在，我们可以编写一个查询来检查点是否位于多边形内部，并计算这些点与多边形中心点的距离。
-- 假设我们想要检查'Point 1'和'Point 2'是否分别位于'Polygon A'和'Polygon B'内部
-- 并计算它们与各自多边形中心点的距离
-- 首先，我们需要计算每个多边形的中心点
SET @polygonA_center = ST_Centroid(GeomFromText('POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))'));
SET @polygonB_center = ST_Centroid(GeomFromText('POLYGON((20 20, 30 20, 30 30, 20 30, 20 20))'));
-- 然后，我们可以使用这些中心点与点表中的点进行比较和距离计算
SELECT
p.name AS point_name,
p.position,
CASE
WHEN ST_Contains(pg.shape, p.position) THEN 'Inside'
ELSE 'Outside'
END AS location_status,
ST_Distance_Sphere(p.position, CASE pg.name WHEN 'Polygon A' THEN @polygonA_center ELSE @polygonB_center END, 6371) AS distance_km
FROM
points p
JOIN
polygons pg ON (
(p.name = 'Point 1' AND pg.name = 'Polygon A') OR
(p.name = 'Point 2' AND pg.name = 'Polygon B')
);
这个查询首先计算了两个多边形的中心点，并使用JOIN语句将点表与多边形表连接起来。它使用ST_Contains函数来检查点是否位于多边形内部，并使用ST_Distance_Sphere函数来计算点与对应多边形中心点的距离（以公里为单位）。注意，我们使用了CASE语句来根据点的名称选择正确的多边形中心点进行计算。
这个查询将返回每个点的名称、位置、是否在多边形内部的状态以及与对应多边形中心点的距离。
