03 PostGIS如何有效比较2D/3D两种维度下的几何是否相同

最新推荐文章于 2024-01-09 22:14:05 发布
最新推荐文章于 2024-01-09 22:14:05 发布
阅读量686 收藏 1
点赞数 1
分类专栏： # PostGIS 文章标签： 3d mybatis java postgresql
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44188072/article/details/130722040
版权

一 背景

近期遇到一些需求，使用PostGIS存储几何数据，做的过程中需求产生过一些变化，大概如下：

  1. 服务端采用PostGIS存储几何数据，几何类型为Point、LineString、Polygon、MultiPoint、MultiLineString、MultiPolygon等常见的几种几何类型，维度为二维，不包含Z坐标；
  2. 创建几何时采用几何信息的WKT格式进行创建，例如“POINT (x y)”；
  3. 创建时如果几何存在则返回指定记录id（涉及到几何对比）；
  4. 新需求，几何允许Z坐标的存储；

在要求存储Z坐标之前，解决的主要思路为：

  1. 服务端接受几何WKT数据构建几何对象，实体中几何属性采用org.postgis.geometry或者org.locationtech.jts.geom；
  2. 自定义配置MyBatisPlus对postgis的几何（postgis-Jdbc）与自定义类型几何（org.postgis.geometry或者org.locationtech.jts.geom)的转换；
  3. 采用MyBatis/MyBatisPlus进行几何数据记录数据的CRUD，插入几何记录前调用在Mapper.xml层自定义的是否有指定几何存在的方法（PostGIS的ST_Equals(geo1,geo2)函数）来判断是否有存在的几何记录；

整个思路在未添加Z维度坐标之前是正确的，但是在添加Z坐标维度后，有些地方需要进行调整，主要为：

  1. 实体中几何属性不建议采用GeoTools（org.locationtech.jts.geom）中组织的几何，因为其不具有Z坐标的支持，可使用postgis的几何（postgis-Jdbc）替换（此部分不是重点，仅提醒一下）；

image.png
以点类型对比，虽然GeoTools中的Point类中Z信息也存在于Coordinate中，但Point类并没有很好的描述，所以如果存在Z甚至M等多维度数据，建议使用postgis-jdbc中对几何的组织方式。

      	<!-- postgis数据库驱动 -->
        <dependency>
            <groupId>net.postgis</groupId>
            <artifactId>postgis-jdbc</artifactId>
        </dependency>
  1. 数据库中建立几何表，采用3D几何类型创建
    1. Point -> PointZ
    2. Polygon -> PolygonZ
    3. LineString -> LineStringZ
    4. MultiPoint ->MultiPointZ
    5. MultiLineString -> MultiLineStringZ
    6. MultiPolygon -> MultiPolygonZ
-- ----------------------------
-- 添加PostGIS扩展
-- ----------------------------
-- CREATE EXTENSION postgis;


DROP TABLE IF EXISTS "public"."data_point";
CREATE TABLE data_point (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(Point, 4326)
);
CREATE INDEX data_point_geom_index ON data_point USING GIST (geometry);

DROP TABLE IF EXISTS "public"."data_pointz";
CREATE TABLE data_pointz (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(PointZ, 4326)
);
CREATE INDEX data_pointz_geom_index ON data_pointz USING GIST (geometry);



DROP TABLE IF EXISTS "public"."data_line";
CREATE TABLE data_line (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(LineString, 4326)
);
CREATE INDEX data_line_geom_index ON data_line USING GIST (geometry);

DROP TABLE IF EXISTS "public"."data_linez";
CREATE TABLE data_linez (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(LineStringZ, 4326)
);
CREATE INDEX data_linez_geom_index ON data_linez USING GIST (geometry);



DROP TABLE IF EXISTS "public"."data_polygon";
CREATE TABLE data_polygon (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(Polygon, 4326)
);
CREATE INDEX data_polygon_geom_index ON data_polygon USING GIST (geometry);

DROP TABLE IF EXISTS "public"."data_polygonz";
CREATE TABLE data_polygonz (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(PolygonZ, 4326)
);
CREATE INDEX data_polygonz_geom_index ON data_polygonz USING GIST (geometry);



DROP TABLE IF EXISTS "public"."data_point_multi";
CREATE TABLE data_point_multi (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(MultiPoint, 4326)
);
CREATE INDEX data_point_multi_geom_index ON data_point_multi USING GIST (geometry);

DROP TABLE IF EXISTS "public"."data_pointz_multi";
CREATE TABLE data_pointz_multi (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(MultiPointZ, 4326)
);
CREATE INDEX data_pointz_multi_geom_index ON data_pointz_multi USING GIST (geometry);



DROP TABLE IF EXISTS "public"."data_line_multi";
CREATE TABLE data_line_multi (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(MultiLineString, 4326)
);
CREATE INDEX data_line_multi_geom_index ON data_line_multi USING GIST (geometry);

DROP TABLE IF EXISTS "public"."data_linez_multi";
CREATE TABLE data_linez_multi (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(MultiLineStringZ, 4326)
);
CREATE INDEX data_linez_multi_geom_index ON data_linez_multi USING GIST (geometry);



DROP TABLE IF EXISTS "public"."data_polygon_multi";
CREATE TABLE data_polygon_multi (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(MultiPolygon, 4326)
);
CREATE INDEX data_polygon_multi_geom_index ON data_polygon_multi USING GIST (geometry);

DROP TABLE IF EXISTS "public"."data_polygonz_multi";
CREATE TABLE data_polygonz_multi (
  id SERIAL PRIMARY KEY,
  geometry GEOMETRY(MultiPolygonZ, 4326)
);
CREATE INDEX data_polygonz_multi_geom_index ON data_polygonz_multi USING GIST (geometry);
  1. 原本采用PostGIS的空间分析函数ST_Equals(geo1,geo2)来比较几何相等，但具有Z坐标后无法正确比较，即采用ST_Equals比较"POINT (1 1 1)"和“POINT (1 1 0)”是同一个几何

image.png

二 如何有效比较2D/3D两种维度下的几何是否相同

我将我尝试过的方案划分为了两种：一种是使用PostGIS空间分析函数进行的比较；另外一种是不使用空间分析函数，也就是直接将几何转为WKT字符串，进行字符串的比较，这种方案首先说可以，但是效率低、开销大，不建议使用。

1 利用空间分析函数的比较方式

涉及到的空间分析函数参考下表：

函数名称描述
ST_OrderingEquals(geo1, geo2)Exactly Equal通过逐点比较两个几何体来确定精确的相等性，以确保它们在位置上相同。
ST_Equals(geo1, geo2)Spatially Equal精确相等（Exactly Equal）并没有考虑几何图形的空间性质。有一个函数，恰当地命名为ST_Equals，可用于测试几何图形的空间相等性或等效性。
geo1 ~= geo2Bounds Equal在最坏的情况下，精确相等需要对几何体中的每个顶点进行比较，以确定相等。这可能很慢，并且可能不适合比较大量的几何形状。为了更快地进行比较，提供了相等边界运算符~=。这只对边界框（矩形）进行操作，确保几何图形占据相同的二维范围，但不一定占据相同的空间。

这部分官方提供了名为“Equality”可比较性的专题说明，建议浏览参考。
下列测试以Point、PointZ、MultiPoint、MultiPointZ等进行，测试数据如下：

-- Point
-- id 1
INSERT INTO data_point (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 1), 4326));
-- id 2
INSERT INTO data_point (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 1), 4326));
-- id 3
INSERT INTO data_point (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 2), 4326));

-- PointZ
-- id 1
INSERT INTO data_pointz (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 1, 1), 4326));
-- id 2
INSERT INTO data_pointz (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 1, 0), 4326));
-- id 3
INSERT INTO data_pointz (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 2, 0), 4326));

-- MultiPoint
-- id 1
INSERT INTO data_point_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((1 1),(2 2))',4326));
-- id 2
INSERT INTO data_point_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((2 2),(1 1))',4326));
-- id 3
INSERT INTO data_point_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((2 2),(1 2))',4326));

-- MultiPointZ
-- id 1
INSERT INTO data_pointz_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((1 1 1),(2 2 2))',4326));
-- id 2
INSERT INTO data_pointz_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((2 2 2),(1 1 1))',4326));
-- id 3
INSERT INTO data_pointz_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((2 2 2),(1 1 0))',4326));
-- id 4
INSERT INTO data_pointz_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((1 1 1),(3 3 3))',4326));

（1）ST_Equals(geo1, geo2)

考虑几何图形的空间性质，可用于测试几何图形的空间相等性或等效性，更符合我们对几何相等的直观理解（即包含相同的区域）。需要注意，无论是绘制多边形的方向、定义多边形的起点，还是使用的点数，在这里都不重要，重要的是几何包含相同的空间。
下列的示例将对一些特别的地方进行说明：

a.忽略Z维度（特别注意！特别注意！特别注意！）
SELECT id,ST_AsText(geometry) AS wkt
FROM data_pointz
WHERE 
ST_Equals(
	geometry, 
	ST_GeomFromText(
		'POINT (1 1 1)',
		 4326
	)
)

image.png
注意！虽然data_pointz.1和data_pointz.2的Z不相同，但仍然会判定相等，因为ST_Equals函数仅仅对二维比较

b.忽略几何构件顺序
SELECT id,ST_AsText(geometry) AS wkt
FROM data_point_multi
WHERE 
ST_Equals(
	geometry, 
	ST_GeomFromText('MULTIPOINT ((1 1),(2 2))',4326)
);

image.png
注意，data_point_multi.1和data_point_multi.2两个MultiPoint都是由Point(1 1)和Point(2 2)组成的，只不过顺序不同，但依然被判定相等，同理在LineString、Polygon中也是一样的。

（2）ST_OrderingEquals(geo1, geo2)

ST_OrderingEquals函数则比ST_Equals函数更加“苛刻”，无论是绘制多边形的方向、定义多边形的起点，还是使用的点数，比较时都会被考虑。

a.Z维度的有效比较
-- id 1
INSERT INTO data_pointz (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 1, 1), 4326));
-- id 2
INSERT INTO data_pointz (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 1, 0), 4326));
-- id 3
INSERT INTO data_pointz (geometry)
VALUES (ST_SetSRID(ST_MakePoint(1, 2, 0), 4326));

SELECT id,ST_AsText(geometry) AS wkt
FROM data_pointz
WHERE 
ST_OrderingEquals(
	geometry, 
	ST_GeomFromText(
		'POINT (1 1 1)',
		 4326
	)
)

image.png

b.几何构建顺序的有效比较
-- id 1
INSERT INTO data_pointz_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((1 1 1),(2 2 2))',4326));
-- id 2
INSERT INTO data_pointz_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((2 2 2),(1 1 1))',4326));
-- id 3
INSERT INTO data_pointz_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((2 2 2),(1 1 0))',4326));
-- id 4
INSERT INTO data_pointz_multi (geometry)
VALUES (ST_GeomFromText('MULTIPOINT ((1 1 1),(3 3 3))',4326));

SELECT id,ST_AsText(geometry) AS wkt
FROM data_pointz_multi
WHERE 
ST_OrderingEquals(
	geometry, 
	ST_GeomFromText('MULTIPOINT ((2 2 2),(1 1 1))',4326)
);

image.png

2 不太常规的比较方式

（1）转WKT进行纯字符串层面的比较


-- 比较方式4：字符串比较
SELECT id,ST_AsText(geometry) AS wkt
FROM data_point
WHERE 
ST_AsText(geometry) = ST_AsText(ST_GeomFromText('POINT(1 1)'));

image.png


-- 比较方式4：字符串比较
SELECT id,ST_AsText(geometry) AS wkt
FROM data_pointz
WHERE 
ST_AsText(geometry) = ST_AsText(ST_GeomFromText('POINT(1 1 1)'));

image.png

（2）利用PostGIS的距离求算函数判断是否有距离为0的几何存在，如果有则代表已经存在；

SELECT id,ST_AsText(geometry) AS wkt
FROM data_pointz
WHERE 
(ST_Distance(geometry,ST_GeomFromText('POINT (1 1 1)',4326))) = 0;

image.png

SELECT id,ST_AsText(geometry) AS wkt
FROM data_pointz
WHERE 
(ST_3DDistance(geometry,ST_GeomFromText('POINT (1 1 1)',4326))) = 0;

image.png

三 总结

类似插入记录时，如果几何存在则返回已存在记录id的需求还是很常见的，但是需要结合具体的业务场景以及数据情况（2D/3D/…)来判断采用什么类型的空间分析函数。
ST_OrderingEquals足够精确，无论是绘制多边形的方向、定义多边形的起点，还是使用的点数，比较时都会被考虑。我用GeoJson举一个例子：

{
    "type": "FeatureCollection",
    "name": "水位监测",
    "crs": {},
    "features": [
        {
            "type": "Feature",
            "properties": {
                "time": "time 1",
                "waterLevel": [
                    1.1,
                    2.2,
                    3.3
                ]
            },
            "geometry": {
                "type": "MultiPoint",
                "coordinates": [
                    [1,1],
                    [2,2],
                    [3,3]
                ]
            }
        },
        {
            "type": "Feature",
            "properties": {
                "time": "time 2",
                "waterLevel": [
                    1.11,
                    2.22,
                    3.33
                ]
            },
            "geometry": {
                "type": "MultiPoint",
                "coordinates": [
                    [1,1],
                    [2,2],
                    [3,3]
                ]
            }
        }
    ]
}

数据包含两个要素，要素为MultiPoint类型，包含三个坐标点，每个坐标点为水位监测仪器的的空间位置，属性“waterLevel”则对应三个点的具体水位值。类似这样的场景，我们就不能仅仅考虑几何是否包含（占据）相同的空间。
假设某些场景中并不需要考虑Z坐标以及几何的构建顺序，则ST_Equals就是最佳选择，因为它的效率优于前者，

billy_gisboy
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
3D-PostGIS3DExplorer.zip
09-18
3D-PostGIS3DExplorer.zip,学习可视化使用postgis三维功能创建三维模型！,3D建模使用专门的软件来创建物理对象的数字模型。它是3D计算机图形的一个方面，用于视频游戏，3D打印和VR，以及其他应用程序。
判断postgis 图层的几何类型
weixin_35749796的博客
01-09 229
可以使用 ST_GeometryType 函数来判断 PostGIS 图层的几何类型。例如，如果你想要检查表中名为 "geom_column" 的列的几何类型，可以使用以下 SQL 语句： SELECT ST_GeometryType(geom_column) FROM table_name; ...
engine中空间运算接口ITopologicalOperator与空间关系接口IRelationalOperator
rrrrssss00的专栏
08-12 9803
<br />示例:<br /> <br />ITopologicalOperator sourcePolygon= (ITopologicalOperator)outerPolygon; IPolygon resultPolygon = (IPolygon)sourcePolygon.Difference(cutPolygons);<br /> <br />IRelationalOperator re = (IRelationalOperator)outerPolygon; if (!re.Contai
PostGIS空间数据库三维数据存储和查询
weixin_43138044的博客
06-12 1713
在插入数据时，使用 ST_GeomFromText 函数将文本格式的几何对象转换为 PostGIS 的几何类型。在插入数据时，使用 ST_GeomFromText 函数将文本格式的几何对象转换为 PostGIS 的几何类型。在插入数据时，使用 ST_GeomFromText 函数将文本格式的几何对象转换为 PostGIS 的几何类型。如果需要创建 LineStringZM 类型的三维线段，可以使用类似的方式进行操作。如果需要创建 PointZM 类型的三维点，可以使用类似的方式进行操作。
判断二维数组中元素是否互不相同
qq_61706112的博客
11-17 2517
#include<stdio.h> #define MaxSize 20 //打印二维数组 void printArr(int arr[MaxSize][MaxSize],int m,int n) { int i,j; for(i=0;i<m;++i) { for(j=0;j<n;++j) printf("%d ",arr[i][j]); printf("\n"); } } //判断二维数组中是否有相同元素存在 bool judgeSameElem(int.
postgis/postgresql实验报告
03-08
PostGIS是一个基于开源数据库管理系统PostgreSQL的空间数据库扩展，它允许将空间数据与其他非空间数据一起存储在同一个数据库中，并提供了一组用于处理空间数据的函数和查询。本实验主要介绍了如何在PostgreSQL中...
PostGIS空间数据库扩展到PostgreSQL [mirror] - PostGIS / PostGIS
01-27
Debbie: Winnie: Dronie: Travis: ...Official chat room is the #postgis channel on the network. To participate, point your preferred to: irc://irc.freenode.net/#postgis Or try a web IRC
12.（地图工具篇）PostgreSQL+PostGIS实现两坐标点之间最短路径查询算法函数.zip
09-26
下载如有问题，可私信博主。下载前建议先查看博客内容，其地址为：https://blog.csdn.net/QQ98281642/article/details/120483897
pg2b3dm:从PostGIS转换为b3dm切片的工具
04-28
pg2b3dm 用于将3D几何图形从PostGIS转换为 / b3dm tile的工具。 该软件最初是py3dtiles（ ）的端口，用于生成b3dm tile。 可以在Cesium JS或MapBox GL JS中可视化生成的3D Tiles。 与py3dtiles的区别： 性能改善； ...
ArcGIS小技术 ——2D数据转3D的数据（要素类添加Z、M）
alun550的博客
04-07 5261
2D数据转3D的数据（要素类添加Z、M）原始数据方法一 —— 复制要素方法二 —— 依据属性实现要素转3D方法三 —— 插值shape         原始数据         方法一 —— 复制要素 在复制要素的时候，设置环境 enable z和enable m，arcpy同理。 仅仅开启z   ...
PostGIS教程二十一：3-D
不睡觉的怪叔叔的博客
03-05 4314
目录 一、3-D几何图形 二、3-D函数 三、N-D索引 注意：本节介绍许多PostGIS2.0及更高版本才支持的功能。 一、3-D几何图形 到目前为止，我们一直在处理2-D几何图形（二维几何图形），只有X和Y坐标。但是PostGIS支持所有几何图形类型额外的维度，对于每个坐标，另外还能支持用于表示高度信息的"Z"维度以及用于添加额外附加信息的"M"维度（通常为时间、道...
A full Web 3D GIS, using PostGIS and X3D output
weixin_30594001的博客
08-10 178
A full Web 3D GIS, using PostGIS and X3D output Presentation | Presented Olivier Courtin, Camptocamp France SAS David Jonglez, Camptocamp France SAS This presentation shows how P...
PostGIS学习教程二十：3-D
最新发布
数字人生
01-09 1232
对于3-D和4-D几何图形，额外的维度将作为几何图形中每个顶点的额外坐标添加，并且几何图形类型将得到增强，以指示如何解译额外的维度信息。构建索引后，可以在查询中使用&&&索引操作符，&&&和&&是相同的语义——边界框相交——区别在于，&&&使用几何图形的所有维度来应用这个语义。一般来说，N-D索引只比2-D索引执行速度稍慢一些，所以只使用N-D索引，因为N-D查询将提高查询的多维度选择性。要使用N-D索引搜索nyc_streets表，只需将&&2-D索引运算符替换为&&&3-D索引运算符。
基于PostGIS的高级应用(3)--线性参考
weixin_34315189的博客
08-30 1010
一 线性参考干啥用的   如果直接写个“高大上”的定义结果往往是一脸懵逼的，也不知道为什么要定义这么一个概念。其实线性参考技术在我们生活中是非常常见的，比如打开高德，百度地图的App，查看实时路况，道路被不同路况的颜色动态分段显示了；高速中发生交通事故，电视广播中常常对地点描述为“距离xx高速入口xx公里处”...
PostGIS:计算一个二维面在三维线上的投影
专注GIS软件开发。
04-10 219
上述方法只适用于地上三维对象的计算。如果是地下三维对象，应当将ST_Extrude对象的全部Z值减去高差的1/2。主要思路：将二维面对象转换为POLYHEDRALSURFACE对象。
PostgreSQL 3D City 应用
weixin_33862514的博客
01-27 700
标签 PostgreSQL , CityJSON , CityGML , 3D City , 3D City Database , Geometrical Model , Semantic Model , COLLADA , VRML , X3D , IFC , gbXML , KML , OGC (Open Geospatial Consortium)...
PostgreSQL + PostGIS + SFCGAL 优雅的处理3D数据
weixin_34124651的博客
10-28 1670
标签 PostgreSQL , SFCGAL , 3D 背景 PostGIS 整合 SFCGAL，优雅的处理3D空间数据。 例子 PDF: 3D and exact geometries for PostGIS , FOSDEM PGDay http://www.sfcgal.org/ https://www.tuicool.com/arti...
PostGIS教程七：几何图形（Geometry）
南方有岛
05-11 4047
目录 一、介绍 二、元数据表 三、表示真实世界的对象 3.1、点（Points） 3.2、线串（Linestring） 3.3、多边形（Polygon） 3.4、集合（Collection） 四、几何图形输入和输出 五、从文本转换 一、介绍 在前面的章节中，我们已经往数据库中加载了数据，现在让我们来先看一些简单的例子。 在pgAdmin中，再次选择nyc数据...
mybatis几何geometry类型转换
qq_40224400的博客
01-06 6318
利用mybatis从postgresql中读取geometry类型字段到实体类中，可以多种接受方式： 1. 直接利用String去接受wkt 2. 利用mybatis在写入实体类时将wkt转换为geometry类型 第一种方式：非常简单粗暴，不过在代码中需要用到geometry类型时，都会写上如下代码： try { ...... JtsUtils.createLineStringByWKT(w...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值