如何在大量2002类型线数据中查找首尾重合的线对？

核心思路：直接匹配端点坐标 #

针对2002类型线数据，最精准的方式是提取每条线的起点和终点，然后查找与其他线端点完全重合的线对。避免用邻近分析（如SDO_NN），因为邻近分析会包含距离接近但未完全重合的点，而你需要的是完全重合的端点。

步骤1：提取线的端点 #

使用SDO_GEOM.SDO_POINT_AT_POSITION函数提取线的起点（位置参数为1）和终点（位置参数为-1）：

• 起点：SDO_GEOM.SDO_POINT_AT_POSITION(geom, 1)
• 终点：SDO_GEOM.SDO_POINT_AT_POSITION(geom, -1)

步骤2：编写SQL查询匹配线对 #

假设你的线表名为lines，几何字段为geom，主键为id，以下SQL可直接找出所有首尾重合的线对，并标注匹配类型：

SELECT
  a.id AS line_id_1,
  b.id AS line_id_2,
  CASE
    WHEN SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, 1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, 1), 'mask=EQUAL') = 'TRUE' THEN '起点-起点重合'
    WHEN SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, 1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, -1), 'mask=EQUAL') = 'TRUE' THEN '起点-终点重合'
    WHEN SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, -1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, 1), 'mask=EQUAL') = 'TRUE' THEN '终点-起点重合'
    WHEN SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, -1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, -1), 'mask=EQUAL') = 'TRUE' THEN '终点-终点重合'
  END AS match_type
FROM lines a
JOIN lines b ON a.id != b.id
WHERE
  SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, 1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, 1), 'mask=EQUAL') = 'TRUE'
  OR SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, 1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, -1), 'mask=EQUAL') = 'TRUE'
  OR SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, -1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, 1), 'mask=EQUAL') = 'TRUE'
  OR SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, -1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, -1), 'mask=EQUAL') = 'TRUE'
AND a.id < b.id; -- 避免重复输出(A,B)和(B,A)

优化建议 #

• 确保geom字段已创建空间索引，否则大数量数据查询会极慢：

  CREATE INDEX lines_sidx ON lines(geom) INDEXTYPE IS MDSYS.SPATIAL_INDEX;
  

• 如果需要排除自身匹配，a.id != b.id已经处理，加上a.id < b.id可避免重复配对。

关于SDO_UTIL.CONCAT_LINES的说明 #

你提到的SDO_UTIL.CONCAT_LINES不是用来查找线对的工具，它的作用是合并两条首尾相连的线为单条线。找到匹配的线对后，可使用它完成合并操作，示例：

SELECT SDO_UTIL.CONCAT_LINES(a.geom, b.geom) AS merged_geometry
FROM lines a
JOIN lines b ON a.id = 123 AND b.id = 456 -- 替换为实际匹配的线ID
WHERE SDO_RELATE(SDO_GEOM.SDO_POINT_AT_POSITION(a.geom, -1), SDO_GEOM.SDO_POINT_AT_POSITION(b.geom, 1), 'mask=EQUAL') = 'TRUE';

内容的提问来源于stack exchange，提问作者Pierre de la Verre