MySQL ST_FrechetDistance 怎样比较两条轨迹的形状差异

如果业务要比较两条 GPS 轨迹、路线折线或设备运动路径，单看起点和终点的距离往往不够：两条路线可能从同一点出发、在中途绕行后又回到相近终点。MySQL 的 ST_FrechetDistance() 更适合回答“沿着各自轨迹前进时，两条折线最难跟随的差异有多大”。

官方地址：https://dev.mysql.com/doc/refman/8.4/en/spatial-relation-functions-object-shapes.html

它返回的是离散 Fréchet 距离：只使用几何中明确写出的点，不把线段之间的每个连续点自动补进计算。本文用两条 LineString 演示调用方式，并把最容易踩到的采样、SRS、单位和异常边界放在同一条实践路径里。

先准备同一 SRS 的两条 LineString

ST_FrechetDistance(g1, g2[, unit]) 的两个几何参数要能在同一空间参考系下比较。为了先看懂算法含义，下面使用 SRID 0 的平面坐标；这时数值可以理解为当前坐标单位中的距离，不要直接把它当成米。

-- 用少量显式顶点构造两条待比较轨迹
SET @route_a = ST_GeomFromText('LINESTRING(0 0, 0 5, 5 5)', 0);
SET @route_b = ST_GeomFromText('LINESTRING(0 1, 0 6, 3 3, 5 6)', 0);

-- 先确认两个几何的类型与 SRID，避免把不匹配的输入直接送入距离函数
SELECT ST_GeometryType(@route_a) AS type_a,
       ST_GeometryType(@route_b) AS type_b,
       ST_SRID(@route_a) AS srid_a,
       ST_SRID(@route_b) AS srid_b;

这里故意让第二条轨迹多一个拐点，便于观察“形状跟随差异”。两条值都应是 LineString，并且 SRID 一致；如果一条是 LineString、另一条是 Polygon，或者两条几何的 SRS 不同，就不属于这个函数的正常输入范围。

调用 ST_FrechetDistance 读取轨迹差异

输入准备好后，直接把两条几何传给函数。返回值是一个双精度数，数值越小通常表示在当前采样方式和坐标单位下越接近；它不是百分比，也不是“相似度 0 到 100 分”。

-- 计算两条轨迹的离散 Fréchet 距离
SELECT ST_FrechetDistance(@route_a, @route_b) AS frechet_distance;

-- 用业务阈值筛选形状差异较小的轨迹对
SELECT CASE
         WHEN ST_FrechetDistance(@route_a, @route_b) 

这个判断的阈值必须由坐标单位、设备误差、采样频率和业务容忍度共同决定。示例中的 3 只是 SQL 写法演示，不能直接搬到经纬度或真实道路数据中。

与只比较首尾点不同，离散 Fréchet 距离会关注两条折线沿路径前进时的整体匹配关系。因此，一条路线中间明显偏离、但终点仍相近时，距离值仍可能把这个偏离反映出来。

理解“离散”带来的采样影响

MySQL 手册明确说明，针对 LineString，函数只考虑几何中显式列出的点，线段之间的点不参与计算。这一点决定了轨迹入库方式：同一条真实路线，如果一次每 5 秒采样、一次每 60 秒采样，得到的折线顶点集合不同，结果也可能不同。

-- 两条折线在端点相近，但第二条显式记录了一个偏离点
SET @coarse = ST_GeomFromText('LINESTRING(0 0, 10 0)', 0);
SET @detour = ST_GeomFromText('LINESTRING(0 0, 5 4, 10 0)', 0);

-- 只比较当前几何中已经记录的顶点，不会自动采样线段内部
SELECT ST_FrechetDistance(@coarse, @detour) AS sampled_shape_gap;

因此，在写入数据库前应先约定采样策略：固定时间间隔、固定距离间隔，或者按轨迹转向点保留。比较结果时同时保存采样规则，后续才知道两个数值是否可比。不要通过随意增加顶点来“调到”想要的结果；增加的点必须来自真实采样或明确的几何预处理。

用 SRS 和单位让结果可解释

如果轨迹使用地理 SRS，例如经纬度坐标，距离会按照该 SRS 的语义处理；如果指定了单位参数，则结果按指定的长度单位返回。工程上最重要的不是记住某个数字，而是保证两条轨迹使用相同 SRS，并在结果字段旁记录单位。

-- 使用同一地理 SRS 创建两条轨迹，示例坐标为经度和纬度
SET @gps_a = ST_GeomFromText(
  'LINESTRING(116.3800 39.9000, 116.3900 39.9050)', 4326
);
SET @gps_b = ST_GeomFromText(
  'LINESTRING(116.3800 39.9010, 116.3910 39.9060)', 4326
);

-- 明确要求按米返回，让阈值含义更直观
SELECT ST_FrechetDistance(@gps_a, @gps_b, 'metre') AS distance_metre;

-- 同一比较也可以按英尺返回，但不要混用两个单位的阈值
SELECT ST_FrechetDistance(@gps_a, @gps_b, 'foot') AS distance_foot;

官方文档还说明，支持的单位可以从 INFORMATION_SCHEMA.ST_UNITS_OF_MEASURE 查询；未知单位会报错。对于 SRID 0，指定线性单位可能触发未知长度单位错误，所以不能只把 'metre' 机械地加到所有平面坐标上。

在入库查询中处理空值和异常边界

正式表查询时，几何列可能为空、类型不对、SRID 不一致，或者包含越界的地理坐标。不要把函数调用直接拼进大批量更新；先筛出可比较的记录，再按业务决定是否记录异常。

-- 只取两条轨迹都存在且类型为 LineString 的候选记录
SELECT a.track_id AS track_a,
       b.track_id AS track_b,
       ST_FrechetDistance(a.geom, b.geom, 'metre') AS distance_metre
FROM track a
JOIN track b ON a.track_id 

这段筛选只解决了明显的空值、类型和 SRID 条件，不能替代对几何有效性的业务检查。根据手册的通用规则，空参数或空几何可能返回 NULL，格式不正确会报 ER_GIS_INVALID_DATA，不同 SRS 会报 ER_GIS_DIFFERENT_SRIDS；如果几何本身无效，结果可能未定义或报错。

把距离值放进可复现的判定规则

一个稳妥的轨迹相似筛选，至少应把距离、单位、SRS 和采样策略一起落库。这样当阈值变化或重新采样时，可以解释为什么同一对设备的结果发生变化。

-- 保存比较口径，避免只留下一个无法解释的浮点数
INSERT INTO track_similarity (track_a, track_b, distance_value, distance_unit, srid, sample_policy)
SELECT a.track_id,
       b.track_id,
       ST_FrechetDistance(a.geom, b.geom, 'metre'),
       'metre',
       ST_SRID(a.geom),
       '固定时间间隔采样'
FROM track a
JOIN track b ON a.track_id 

最后再根据业务阈值做分组或排序，而不是把距离直接翻译成“相似百分比”。如果需要观察路线形状差异，保留原始 LineString 和本次采样元数据，比只保留一个分类标签更有用。

常见问题

ST_FrechetDistance 支持 Polygon 吗？
在相同 Cartesian 或 geographic SRS 下，官方说明该函数只支持 LineString；不要把它当成任意几何类型的通用距离函数。

为什么轨迹看起来接近，距离却变大？
先检查显式顶点和采样密度。离散算法只看几何中明确记录的点，采样规则不同会改变比较对象。

SRID 0 能不能直接使用 metre？
不要直接假设可以。SRID 0 没有已知长度单位，指定线性单位可能触发未知长度单位错误；应使用有明确单位语义的 SRS，或先完成坐标转换。

这个返回值越小就一定代表道路完全相同吗？
不一定。它只描述给定两组显式顶点在当前 SRS 和采样策略下的离散形状差异，还要结合时间、方向、道路拓扑等业务条件判断。