eecolidar/postgis

AHN2_Kadaster

使用Kadaster分类丰富AHN2数据

工作流

项目工作流示意图。该示意图的PDF版本也位于仓库中。

PostGIS服务器

PostGIS服务器的Docker镜像可在https://hub.docker.com/u/eecolidar/中找到。

可通过以下命令拉取镜像：

docker pull docker.xuanyuan.run/eecolidar/postgis

安装Docker

https://docs.docker.com/install/安装Docker。

更改Docker卷存储路径

sudo systemctl stop docker
sudo mkdir /data/local/docker-data

在**/etc/docker/daemon.json**中添加以下行：

{
	"data-root": "/data/local/docker-data"
}

启动服务器

sudo docker run -d --name eecolidar-postgis -e POSTGRES_PASSWORD=mysecretpassword -p 5432:5432 eecolidar/postgis

配置防火墙

如果需要配置防火墙，以下命令在Ubuntu上为PostGIS服务器开放端口：

sudo ufw allow 5432/tcp

连接PostGIS服务器

使用pgAdmin（推荐）

安装https://www.pgadmin.org并按照https://www.pgadmin.org/docs/pgadmin4/3.x/connect_to_server.html操作。

使用psql

psql -h 服务器IP -p 5432 -U postgres

从Docker镜像使用psql

也可使用PostgreSQL Docker镜像连接服务器：

docker pull docker.xuanyuan.run/postgres
docker run --rm -ti docker.xuanyuan.run/postgres /usr/bin/psql --help

PostgreSQL命令示例

列出数据库：

\l 或 \list

连接数据库：

\c database_name 或 \connect database_name

显示数据库中的表：

\dt+

查询表中数据：

SELECT * FROM table_name;

或

SELECT column, column2….
FROM docker.xuanyuan.run/table;

统计行数：

SELECT COUNT (*)
FROM docker.xuanyuan.run/table_name;

检查数据库（PostGIS）扩展

要使用PostGIS功能，需安装相应扩展。可通过以下命令检查数据库的现有扩展：

\connect database_name
\x
\dx postgis*

GML数据

显示GML文件信息

显示基本信息：

ogrinfo -fields=YES  TOP10NL_01O.gml

显示所有图层和要素：

ogrinfo -fields=YES -al TOP10NL_01O.gml

将GML转换为CSV

ogr2ogr -f "CSV" Top10NL_000001.csv Top10NL_000001.gml

将GML文件导入服务器

ogr2ogr -f "PostgreSQL" PG:"host=localhost user=postgres dbname=eecolidar password=mysecretpassword" Top10NL_000001.gml

导入到新表：

ogr2ogr -f "PostgreSQL" PG:"host=localhost user=postgres dbname=template_postgis password=mysecretpassword" Top10NL_000001.gml -nln newtablename

边界框

导入PostGIS数据库

边界框信息存储在数据库的bboxes表中。要导入Tile02O的边界框文件：

ogr2ogr -f "PostgreSQL" PG:"host=localhost user=postgres dbname=Tile02O password=mysecretpassword"  Tile02O_bounding_boxes_geom.csv -nln bboxes

psql -h localhost -U postgres --dbname='Tile02O' -c "
  ALTER TABLE public.bboxes ADD COLUMN wkb_geometry geometry;
  UPDATE public.bboxes SET wkb_geometry  = ST_GeomFromText(bboxes.geom, 28992);
"

TOP10NL中的类名

1: Gebouw（建筑物） 2: Inrichtingselement（设施元素） 3: Terrein (Polygon)（区域（多边形）） 4: Waterdeel（水域） 5: GeografischGebied (Point)（地理区域（点）） 6: FunctioneelGebied（功能区域） 7: Plaats（地点） 8: RegistratiefGebied（注册区域） 9: Hoogte（高度） 10: Relief (Line String)（地形（线串）） 11: Wegdeel（道路）

数据库操作

显示列详情

\d+ tabl_name;

显示前10行

SELECT * FROM table_name LIMIT 10;

查询交集

统计Tile02O的交集数量：

sql
SELECT COUNT(*)
FROM docker.xuanyuan.run/bboxes, waterdeel
WHERE ST_Intersects(bboxes.wkb_geometry, waterdeel.wkb_geometry);

查找Tile02O特定子瓦片的交集（以ogc_fid作为子瓦片标识符）：

sql
SELECT bboxes.ogc_fid, bboxes.name, waterdeel.ogc_fid AS waterdeel_ogc_fid, waterdeel.lokaalid AS waterdeel_lokaalid, ST_INTERSECTION(bboxes.wkb_geometry, waterdeel.wkb_geometry), ST_GeometryType(ST_INTERSECTION(bboxes.wkb_geometry, waterdeel.wkb_geometry))
FROM docker.xuanyuan.run/bboxes, waterdeel
WHERE (ST_GeometryType(waterdeel.wkb_geometry) = 'ST_Polygon' AND ST_Intersects(bboxes.wkb_geometry, waterdeel.wkb_geometry) AND bboxes.ogc_fid = 1)
GROUP BY bboxes.ogc_fid, waterdeel.ogc_fid

数据

下载top10NL数据

信息页面

直接下载链接。

依赖项/工具

首先创建虚拟环境

python -m venv venv
source ./venv/bin/activate
pip install --upgrade pip
pip install -r requirements.txt

资源

PostGIS

PostGIS文档可在此获取： https://postgis.net/documentation/

第4节提供数据管理和查询文档，第3节提供常见问题的即时解答。

PDAL

PDAL文档可在此获取： [***]

PDAL模块和过滤器可通过JSON声明格式组合成管道。PDAL已用于AHN2数据的合并。