SQL 计算两点之间的距离（纬度、经度）

在本文中，我们将介绍如何使用SQL计算两个地理位置之间的距离。计算地理距离是一项常见的任务，特别是在需要根据地理位置筛选和排序数据时。我们将使用标准的大圆距离公式（Haversine formula）来计算两个经纬度坐标之间的距离。

阅读更多：SQL 教程

大圆距离公式

大圆距离公式是一种计算两个地球上任意两点之间距离的数学模型。该公式基于大圆上两点之间的弧长。在使用该公式之前，我们需要确保经纬度坐标已转换为弧度。

以下是大圆距离公式的表达式：

d = 2 * R * asin(sqrt(sin((lat2 – lat1) / 2) * sin((lat2 – lat1) / 2) + cos(lat1) * cos(lat2) * sin((lon2 – lon1) / 2) * sin((lon2 – lon1) / 2)))

其中：
– d是两个点之间的距离（单位为米）
– R是地球的半径（单位为米）
– lat1, lon1是第一个点的纬度和经度（以弧度为单位）
– lat2, lon2是第二个点的纬度和经度（以弧度为单位）

SQL 实现

在SQL中，我们可以使用函数和运算符来实现上述的大圆距离公式。以下是一个示例SQL查询，用于计算两个坐标之间的距离：

SELECT 
    id, 
    (6371000 * acos(cos(radians(lat1)) * cos(radians(lat2)) * 
    cos(radians(lon2) - radians(lon1)) + sin(radians(lat1)) * 
    sin(radians(lat2)))) AS distance
FROM 
    locations

在上面的查询中，我们假设有一个名为locations的表，其中包含id、lat1、lon1、lat2和lon2等列。使用上述公式，我们计算了每对坐标之间的距离，并将结果作为distance列返回。

请注意，查询中的纬度和经度应以弧度表示，因此我们使用radians函数将角度转换为弧度。

示例

假设我们有一个包含多个地理位置的表，并且我们想要计算每个位置与给定坐标（纬度：40.7128，经度：-74.0060）之间的距离。该坐标对应于纽约市的经纬度。以下是一个示例表和查询：

-- 创建示例表
CREATE TABLE locations (
    id INT,
    name VARCHAR(50),
    lat FLOAT,
    lon FLOAT
);

-- 插入示例数据
INSERT INTO locations VALUES
    (1, 'Location A', 41.8781, -87.6298),
    (2, 'Location B', 34.0522, -118.2437),
    (3, 'Location C', 29.7604, -95.3698),
    (4, 'Location D', 32.7767, -96.7970),
    (5, 'Location E', 40.7128, -74.0060);

-- 计算距离
SELECT 
    id,
    name,
    (6371000 * acos(cos(radians(lat)) * cos(radians(40.7128)) * 
    cos(radians(-74.0060) - radians(lon)) + sin(radians(lat)) * 
    sin(radians(40.7128))))) AS distance
FROM 
    locations;