在PostGIS中进行点数据的等值线提取
说明
介绍在PostGIS中从点数据提取等值线。
关键字： raster、point、PostGIS、等值线
环境准备
- Postgresql版本：
PostgreSQL 14.0, 64-bit - PostGIS版本：
POSTGIS="3.3.2" - QGIS版本：
3.28.3-Firenze（验证用）
基本步骤
一、数据准备
略.
二、二维点数据转换为三维点
略.
三、将三维点数据转换为栅格
略.
以上三步参见《在PostGIS中进行点数据的栅格化》
四、从栅格数据中提取等值线
关键SQL如下：
Create table _ResultTableName asWITH c AS (SELECT (ST_Contour(rast, 1, fixed_levels => 100,polygonize=>false)).*FROM _TheRastTableName )SELECT geom, id, value FROM c;
其中ST_Contour详细参数说明参见:
https://postgis.net/docs/RT_ST_Contour.html
按照说明参数polygonize=true时，可以生成等值面，但在实际使用时依旧会生成线，所以用了另外一种方法提取等值面，参加其他文档。