在 JTS 上处理几何端 WKBwriter 时出现的问题

Question

我刚刚开始学习几何，但遇到了一些非常奇怪的错误。我正在尝试在我的 java 代码中使用特定的 SRID 创建 WKB。

我已经完成了：

GeometryFactory gm = new GeometryFactory(new PrecisionModel,4326)
WKBWriter w = new WKBWriter(2,ByteOrderValues.LITTLE_ENDIAN)

Geometry geom: Geometry = gm.createPoint(new Coordinate(4,5))
Array[Byte] the_geom = w.write(geom)

println(geom)
println(geom.getSRID)
println(WKBWriter.toHex(the_geom.get))

并获得了

POINT (4 5)
4326
010100000000000000000010400000000000001440

但是在我的数据库中使用postgis：

GEOMETRYFROMTEXT('Point(4 5)',4326)) 

导致

0101000020E610000000000000000010400000000000001440

我在这里出错了？

Answer 1

当存在 SRID 时，PostGIS 不（仅）使用 WKB 标准，它使用 EWKB 支持 3d 和 SRID 的格式，其中 WKB 规范 没有。将 WKB 视为一种轻型二进制二维数组，它是一种非常简单的格式。它不包括二进制表示形式的高程或 SRID，在 OGC 标准空间参考系统中是元数据。以下是一些交互类型，可能会帮助您了解如何使用 WKB 和 EWKB：

moveable=> insert into my_table (some_other_geom) values (GeomFromText('POINT(4 5)'));
moveable=> insert into my_table (some_other_geom) values (GeomFromText('POINT(4 5)', 4326));

这为您提供：

010100000000000000000010400000000000001440
0101000020E610000000000000000010400000000000001440

正如您所看到的。如果您从十六进制值开始并想要添加 SRID 投影并存储为 EWKB，您可以执行以下操作的一些变体：

==>insert into my_table (some_other_geom) values (ST_GeomFromWKB(decode('010100000000000000000010400000000000001440', 'hex'), 4326));
INSERT 0 1

或：

moveable=> insert into my_table (some_other_geom) values (ST_GeomFromWKB(decode('010100000000000000000010400000000000001440', 'hex')));
INSERT 0 1

正如您所期望的，分别产生这些结果：

0101000020E610000000000000000010400000000000001440
010100000000000000000010400000000000001440

其中 null -1 默认 SRID 是与非 SRID 2d WKB 格式相同。总的来说，我认为使用 WKT、GeoJSON 和 KML 更安全，因为它们也是 OGC(?)¹ 开放标准，并且 SQL 版本有一些需要的更新待办的。正如 PostGIS 文档所建议的：“PostGIS 扩展格式目前是 OGC 格式的超集（每个有效的 WKB/WKT 都是有效的 EWKB/EWKT），但这在未来可能会有所不同，特别是如果 OGC 推出与我们的扩展冲突的新格式因此，您不应该依赖此功能！”

(1) 实际上不知道其中一半标准的状态如何，但 json、text 和 xml 似乎可能存在于任何给定应用程序的生命周期中。