Introduction to PostGIS 之线性参考、地理编码

今天学习了postGIS中的两个很有用的方法：

• ST_Line_Interpolate_Point

        geometry ST_Line_Interpolate_Point(geometry a_linestring, float a_fraction)：返回线上内插的一个点

 

返回线上20%处得一个点

• ST_Line_Locate_Point

        float ST_Line_Locate_Point(geometry a_linestring, geometry a_point)：返回a_point离a_linestring在上最近的点位置（这个位置表示为一个0到1浮点数）。点具体的位置我们可以使用前面介绍的ST_Line_Interpolate_Point取出。

 

 

返回给定点到线的最近点 

       这两个方法对于地理编码（Geocoding）是很有用的，下面我就用上一篇导入的数据来做下试验：

       首先我这里用到了QGIS这个软件使查询结果可视化（另外OpenJUMP也可以查询可视化）

      1.首先添加图层

      2.添加完图层之后，就开始执行我们的查询了，这里用到了QGIS中的一个插件工具RT Sql Layer

 这个插件工具需要另外下载安装：1）插件—>Fetch Python Plugins..

                                       

                                          2)点击Fetch Python Plugins后找到RT Sql Layer  安装

      3.安装完插件后  点击RT Sql Layer 连接到数据库island，然后点击Query Builder，在弹出的窗体中输入如下查询语句：

SELECT ST_Line_Interpolate_Point(
        ln.the_geom,
        ST_Line_Locate_Point(
            ln.the_geom,
            pt.the_geom
        )
    )
FROM road as ln, features as pt;

 

       4.得到如下结果：

从图上就会看到 路上多了很多点，但是这个结果并不是我所要的，因为这样的结果是把路两边的很远的features都拉过来了。所以要做下限制，只返回离道路一定距离的点，执行以下SQL：

SELECT DISTINCT ON (pt.gid)

pt.name AS pt_name,
pt.gid AS pt_id,
ST_Line_Interpolate_Point(
ln.the_geom,
ST_Line_Locate_Point(ln.the_geom, pt.the_geom)
) As snapped_point
FROM
features AS pt INNER JOIN
road AS ln
ON
ST_DWithin(pt.the_geom, ln.the_geom, 0.004)
ORDER BY
pt.gid;

这样的结果就比较理想了，这里用到了ST_DWithin方法，这个方法也是很常用的，执行一些地物查找就会用到它。

注：如果你所用的POSTGIS是1.5+的话，ST_ClosestPoint将是一个更好的方法：

geometry ST_ClosestPoint(geometry g1, geometry g2)，这个就简洁了许多，g1不局限于线，g2不局限于点。

你可以用ST_ClosestPoint来替代ST_Line_Interpolate(ST_Line_Locate_Point ..  效果是一样的。