使用SQL Spatial在线与可变距离缓冲区交点创建点要素

完全可以通过SQL Spatial实现你描述的全流程空间运算，不需要借助第三方GIS工具即可完成。

整套操作可以拆分为4步，所有步骤都有标准SQL空间函数支撑：

• 提取线要素的起、终点：主流空间数据库均提供原生函数支持，PostGIS、MySQL可使用ST_StartPoint(geom)提取线起点、ST_EndPoint(geom)提取线终点，SQL Server对应函数为STStartPoint()、STEndPoint()，Oracle Spatial可通过SDO_GEOM.SDO_START_POINT、SDO_GEOM.SDO_END_POINT实现。
• 生成可变距离缓冲区：调用缓冲区生成函数时，直接将每条线对应的自定义缓冲距离作为参数传入即可，不需要设置全局固定缓冲半径，比如PostGIS的ST_Buffer(point_geom, dynamic_distance)、SQL Server的point_geom.STBuffer(dynamic_distance)都支持按要素属性动态传参。
• 计算缓冲区与原线的几何交集：使用空间相交计算函数获取两个几何的重叠部分，以线的起/终点为圆心生成的缓冲区，和原线的交集一定是原线上距离起/终点长度等于缓冲半径的一段短线。
• 提取目标点要素：对交集返回的短线段，提取远离起/终点的那个端点，就是你需要的交点，直接写入点要素表即可。

-- 假设源线表名为road_line，字段包含：fid(要素唯一ID)、shape(线几何字段)、buf_dis(每条线对应的可变缓冲距离)
SELECT
  fid,
  -- 生成起点侧的目标交点
  ST_StartPoint(ST_Intersection(ST_Buffer(ST_StartPoint(shape), buf_dis), shape)) AS start_side_point,
  -- 生成终点侧的目标交点
  ST_EndPoint(ST_Intersection(ST_Buffer(ST_EndPoint(shape), buf_dis), shape)) AS end_side_point
FROM road_line
-- 提前过滤缓冲距离大于线总长、无法生成有效交点的异常数据
WHERE ST_Length(shape) > buf_dis;

注意：不同空间数据库的SQL Spatial函数命名存在细微差异，但上述运算逻辑在PostGIS、SQL Server Spatial、MySQL Spatial、Oracle Spatial中均能实现，不存在能力层面的限制。如果缓冲距离大于对应线要素的总长度，缓冲区会完全覆盖整条线，此时无法生成符合要求的交点，运算前建议增加长度判断过滤这类数据。

内容的提问来源于stack exchange，提问作者Jorge