PostGIS教程十八:线性参考

目录

一、创建线性参考

二、函数列表


    附:关于线性参考的概念,可以参考这篇文章:ArcGIS中的线性参考/动态分段技术(一):Linear Referencing背景

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

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

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

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

一、创建线性参考

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

 
  1. -- Simple example of locating a point half-way along a line

  2. SELECT ST_LineLocatePoint('LINESTRING(0 0, 2 2)', 'POINT(1 1)');

  3. -- Answer 0.5

  4.  
  5. -- What if the point is not on the line? It projects to closest point

  6. -- 即做(0, 2)点到线串(0 0, 2 2)的垂线,使用对应的垂足点来求比例

  7. SELECT ST_LineLocatePoint('LINESTRING(0 0, 2 2)', 'POINT(0 2)');

  8. -- Answer 0.5

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

 
  1. -- 下面所有的SQL都是用来创建新的事件表的

  2. CREATE TABLE nyc_subway_station_events AS

  3. -- 我们首先需要找到一组可能最接近的候选者

  4. -- streets, 按id和distance排列...

  5. WITH ordered_nearest AS (

  6. SELECT

  7. ST_GeometryN(streets.geom,1) AS streets_geom,

  8. streets.gid AS streets_gid,

  9. subways.geom AS subways_geom,

  10. subways.gid AS subways_gid,

  11. ST_Distance(streets.geom, subways.geom) AS distance

  12. FROM nyc_streets streets

  13. JOIN nyc_subway_stations subways

  14. ON ST_DWithin(streets.geom, subways.geom, 200)

  15. ORDER BY subways_gid, distance ASC

  16. )

  17. -- 我们使用'distinct on'PostgreSQL功能为每各唯一的街道gid获取第一条街道(最近的街道)。

  18. -- 然后,我们可以将这条街道信息置入ST_LinLocatePoint函数,使其沿着它的候选地铁站来计算

  19. SELECT

  20. DISTINCT ON (subways_gid)

  21. subways_gid,

  22. streets_gid,

  23. ST_LineLocatePoint(streets_geom, subways_geom) AS measure,

  24. distance

  25. FROM ordered_nearest;

  26.  
  27. -- 主码对于可视化软件很有用

  28. ALTER TABLE nyc_subway_station_events ADD PRIMARY KEY (subways_gid);

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

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

 
  1. -- Simple example of locating a point half-way along a line

  2. SELECT ST_AsText(ST_LineInterpolatePoint('LINESTRING(0 0, 2 2)', 0.5));

  3.  
  4. -- Answer POINT(1 1)

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

 
  1. -- New view that turns events back into spatial objects

  2. CREATE OR REPLACE VIEW nyc_subway_stations_lrs AS

  3. SELECT

  4. events.subways_gid,

  5. ST_LineInterpolatePoint(ST_GeometryN(streets.geom, 1), events.measure)AS geom,

  6. events.streets_gid

  7. FROM nyc_subway_station_events events

  8. JOIN nyc_streets streets

  9. ON (streets.gid = events.streets_gid);

    查看街道上的原始地铁站点(红星)和事件点(蓝色圆圈),你可以看到事件是如何被直接捕捉到最近的街道线的(事件点全部位于街道线上)。

_images/lrs1.jpg

    注意线性参考函数的一个令人惊讶的用法与线性参考模型无关。如上所示,可以使用这些函数将点捕捉到线性要素(即可以使用线性参考来实现捕捉功能)。对于像GPS轨迹或其他预期参考线性网络的输入这样的用例,捕捉是一个方便的功能。

二、函数列表

猜你喜欢

转载自blog.csdn.net/xujingzhong0077/article/details/90168277