3

我将大约 1000 条路径从某处导入我的 postgis 数据库到线串字段。

（编辑）我的桌子是这样的

+---------+---------------+------------------+
| id(int) | name(varchar) | path(LINESTRING) |
+---------+---------------+------------------+
| 123     | foo           | 000002...        |
| 124     | bar           | 000002...        |

我遇到的问题是每条路径都被分成块，并且在某些情况下这些块被混淆了。

假设线串在点号 50 和 70 处被分割：

  1. 块 A：点 1-50
  2. 块 B：点 51-70
  3. 块 C：点 71-100

当我将它迁移到我的数据库中时，它们混合在一起，所以结果线串可能会像这样结束：

  1. 块 A：点 1-50
  2. 块 C：点 71-100
  3. 块 B：点 51-70

所以这会产生从 50 到 71 的跳跃和另一个从 100 到 51 的跳跃

（编辑）当我导入这些分成块的路径时，我认为它们是有序的，但事实是有些是混合的，这使得我的一些线串与它们的点一样排列，就像第二个例子一样。

我希望能够重新排序这些块（点），所以我想构建一个 SQL 查询来检测哪些路径有混合点，然后我可以手动（使用使用 openlayers 制作的工具）重新排列它们。

最好有一个 SQL 更新查询来解决这个问题，但我认为检测更容易（我认为有大约 5% 或更少的路径有错误）

EDIT3：我认为检测脚本可以检查路径是否包含一对太远的连续点。也许从包含最长段的路径排序路径的 SQL 会很好。

如何创建一个函数来获取线串中最大段的长度？

这里我展示一个例子：这就是它在数据库中的样子 不好的方式

这就是我希望它被修复的方式 一个好的（固定的）方式

EDIT4：就像我在 EDIT3 上计划的那样，可以编写一个函数来查找线串中两个连续点之间的最长距离，使用ST_NPoints()和ST_PointN()遍历线串的点，然后可以进行查询以订购路径与那最长的距离。这个距离太长的线串很可能会出现所描述的问题。这样我就可以检测到它们，并手动修复它们。

检测 SQL 的结果将是这样的：

                                             |ordered by this|
+---------+---------------+------------------+---------------+
| id(int) | name(varchar) | path(LINESTRING) |  msbtcp(int)  |
+---------+---------------+------------------+---------------+
| 123     | foo           | 000002...        | 1000          |
| 124     | bar           | 000002...        | 800           |

*msbtcp 将是函数的结果：max_separation_between_two_consecutive_points(path)

4

2 回答 2

3

这听起来有点令人费解，但是如果您刚刚在线串中两个连续点之间的最大距离之后：

CREATE OR REPLACE FUNCTION max_distance_in_linestring(line geometry) RETURNS float as $BODY$
DECLARE
    i integer;
    n integer;
    d float;
    m float;
BEGIN
    d := 0;
    n := ST_NPoints(line);
    i := 2;
    LOOP
        EXIT WHEN i >= n;
        m := ST_Distance(ST_PointN(line,i-1),ST_PointN(line,i));
        -- use for lon,lats:
        -- m := ST_Distance(ST_PointN(line,i-1)::geography,ST_PointN(line,i)::geography);
        IF m > d THEN
            d := m;
        END IF;
        i := i + 1;
    END LOOP;
    RETURN d;
END;
$BODY$
LANGUAGE plpgsql;

SELECT max_distance_in_linestring('LINESTRING(0 0, 1 1, 2 2)'::geometry);
SELECT max_distance_in_linestring('LINESTRING(0 0, 4 3, 2 2)'::geometry);

您可能希望将 ST_PointN 调用重新投射到 ::geography 以获取以米为单位的距离。

SQL 将是这样的：

SELECT
  name, path 
FROM
  paths
ORDER BY
  max_distance_in_linestring(path) DESC
于 2012-04-25T11:13:20.367 回答
0

那么你想知道的是段是如何链接的，所以你需要将每个段的第一个和最后一个点与所有其他段进行比较：

SELECT foo.gid as segment_a, bar.gid as segment_b
    FROM 
        segments AS foo, 
        (SELECT the_geom, gid FROM segments) AS bar
    WHERE 
        bar.gid != foo.gid AND ( -- avoid same segments  
        ST_DWithin( 
                ST_GeometryN(foo.the_geom,  ST_NumGeometries(foo.the_geom)) , -- last from foo
                ST_GeometryN(bar.the_geom,  ST_NumGeometries(bar.the_geom)) , -- last from bar
                0.00005 ) OR -- precision, depends of your SRID
        ST_DWithin( 
                ST_GeometryN(foo.the_geom,  1) , -- first from foo (start = 1 index)
                ST_GeometryN(bar.the_geom,  ST_NumGeometries(bar.the_geom)) , -- last from bar
                0.00005 ) OR -- precision, depends of your SRID
        ST_DWithin( 
                ST_GeometryN(foo.the_geom,  ST_NumGeometries(foo.the_geom)) , -- last from foo
                ST_GeometryN(bar.the_geom,  1) , -- first from bar
                0.00005 ) OR -- precision, depends of your SRID
        ST_DWithin( 
                ST_GeometryN(foo.the_geom,  1) , -- first from foo 
                ST_GeometryN(bar.the_geom,  1) , -- first from bar
                0.00005 ) ) -- precision, depends of your SRID

瞧……！您可以将此链接保存在新表中，并使用 ST_Union 连接链接的段。

于 2012-04-18T15:56:57.073 回答