PostGIS教程十八：线性参考

最新推荐文章于 2025-03-15 18:30:16 发布

不睡觉的怪叔叔

最新推荐文章于 2025-03-15 18:30:16 发布

阅读量4.8k

点赞数 8

分类专栏：数据库

数据库专栏收录该内容

35 篇文章

订阅专栏

一、创建线性参考

二、函数列表

附：关于线性参考的概念，可以参考这篇文章：ArcGIS中的线性参考/动态分段技术（一）：Linear Referencing背景

线性参考是一种表示要素的方法，这些要素可以通过引用一个基本的线性要素来描述。使用线性参照建模的常见示例包括：

高速公路资源，使用沿着公路网中的英里作为参照。
公路维护作业，参照的是沿着公路网的一对英里测量。
水产库存，其中鱼的存在位置被记录为距离上游的一段位置之间。
河流的水文特征，以河流的某一个点到另一个点作为参考。

线性参考模型的优点是，从属空间观测信息不需要与基准空间观测信息分开记录，并且对基础观测信息进行更新时，从属观测信息将自动更新从而追踪新几何图形。

注意：ESRI的线性参照约定是有一个线性空间要素的基表和一个非空间的事件表，其中非空间的事件表包括对空间要素的外键引用和沿参照要素的测量值。我们将使用术语"事件表(event table)"来表示我们构建的非空间表。

一、创建线性参考

如果要参照线性网络的现有点表，请使用ST_LineLocatePoint函数，该函数接受线串和点，并返回该点沿线串的线性参考比例。

-- Simple example of locating a point half-way along a line
SELECT ST_LineLocatePoint('LINESTRING(0 0, 2 2)', 'POINT(1 1)');
-- Answer 0.5

-- What if the point is not on the line? It projects to closest point
-- 即做(0, 2)点到线串(0 0, 2 2)的垂线，使用对应的垂足点来求比例
SELECT ST_LineLocatePoint('LINESTRING(0 0, 2 2)', 'POINT(0 2)');
-- Answer 0.5

我们可以使用ST_LineLocatePoint函数根据nyc_subway_stations创建相对于街道的"事件表"。

-- 下面所有的SQL都是用来创建新的事件表的
CREATE TABLE nyc_subway_station_events AS
-- 我们首先需要找到一组可能最接近的候选者
-- streets, 按id和distance排列...
WITH ordered_nearest AS (
SELECT
  ST_GeometryN(streets.geom,1) AS streets_geom,
  streets.gid AS streets_gid,
  subways.geom AS subways_geom,
  subways.gid AS subways_gid,
  ST_Distance(streets.geom, subways.geom) AS distance
FROM nyc_streets streets
  JOIN nyc_subway_stations subways
  ON ST_DWithin(streets.geom, subways.geom, 200)
ORDER BY subways_gid, distance ASC
)
-- 我们使用'distinct on'PostgreSQL功能为每各唯一的街道gid获取第一条街道（最近的街道）。
-- 然后，我们可以将这条街道信息置入ST_LinLocatePoint函数，使其沿着它的候选地铁站来计算
SELECT
  DISTINCT ON (subways_gid)
  subways_gid,
  streets_gid,
  ST_LineLocatePoint(streets_geom, subways_geom) AS measure,
  distance
FROM ordered_nearest;

-- 主码对于可视化软件很有用
ALTER TABLE nyc_subway_station_events ADD PRIMARY KEY (subways_gid);

一旦我们有了一个事件表，将其转换回一个空间视图是很有趣的，这样我们就可以将事件相对于派生出它们的原始点进行可视化。

要从线性参考比例值得到对应点，我们使用ST_LineInterpolatePoint函数，下面是关于我们前面的简单例子的逆过程：

-- Simple example of locating a point half-way along a line
SELECT ST_AsText(ST_LineInterpolatePoint('LINESTRING(0 0, 2 2)', 0.5));

-- Answer POINT(1 1)

而且，我们可以将nyc_subway_station_events表连接回nyc_streets表，并使用measure属性生成空间事件点（这个示例中是地铁站点），而无需引用原始nyc_subway_stations表。

-- New view that turns events back into spatial objects
CREATE OR REPLACE VIEW nyc_subway_stations_lrs AS
SELECT
  events.subways_gid,
  ST_LineInterpolatePoint(ST_GeometryN(streets.geom, 1), events.measure)AS geom,
  events.streets_gid
FROM nyc_subway_station_events events
JOIN nyc_streets streets
ON (streets.gid = events.streets_gid);