WellKnownGeometry
Reading and writing Well Known Text (WKT) and Well Known Binary (WKB) based on GeoInterface.jl. See this Wikipedia page for an explanation of Well-known text and binary geometry.
Given a GeoInterface compatible geometry, this package can generate the WKT and WKB representation of it. It also does the reverse, as it implements GeoInterface for WKT or WKB strings.
WellKnownText
Given a WKT string, we can retrieve the type and underlying coordinates, and thus convert it to other geometries using GeoInterface. Note that WKT strings are wrapped by GeoFormatTypes so we can distinguish them from any other strings.
using ArchGDAL
using GeoFormatTypes
using WellKnownGeometry
using GeoInterface
wkts = "POINT (30 10)"
wkt = GeoFormatTypes.WellKnownText(GeoFormatTypes.Geom(), wkts)
GeoInterface.geomtrait(wkt) # PointTrait()
GeoInterface.ncoord(wkt) # 2
GeoInterface.coordinates(wkt) # 2-element Vector{Float64}: 30.0 10.0
p = convert(ArchGDAL.IGeometry{ArchGDAL.wkbPoint}, wkt) # Geometry: POINT (30 10)
As ArchGDAL geometries implement GeoInterface, we can generate the WKT for it.
wkt = WellKnownGeometry.getwkt(p) # WellKnownText{GeoFormatTypes.Geom}(GeoFormatTypes.Geom(), "POINT (30.0 10.0)")
GeoFormatTypes.val(wkt) # "POINT (30.0 10.0)"
WellKnownGeometry.getwkt
— Functiongetwkt(geom)
Retrieve the Well Known Text (WKT) as GeoFormatTypes.WellKnownText
for a geom
that implements the GeoInterface. Use GeoFormatTypes.val
to get the String representation.
WellKnownBinary
Given a WKB byte string, we can retrieve the type and underlying coordinates, and thus convert it to other geometries using GeoInterface. Note that WKB byte strings are wrapped by GeoFormatTypes so we can distinguish them from any other byte strings.
using ArchGDAL
using GeoFormatTypes
using WellKnownGeometry
using GeoInterface
wkbs = [0x01, 0x01, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x3e, 0x40, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x24, 0x40]
wkb = GeoFormatTypes.WellKnownBinary(GeoFormatTypes.Geom(), wkbs)
GeoInterface.geomtrait(wkb) # PointTrait()
GeoInterface.ncoord(wkb) # 2
GeoInterface.coordinates(wkb) # 2-element Vector{Float64}: 30.0 10.0
p = convert(ArchGDAL.IGeometry{ArchGDAL.wkbPoint}, wkb) # Geometry: POINT (30 10)
As ArchGDAL geometries implement GeoInterface, we can generate the WKB for it.
wkb = WellKnownGeometry.getwkb(p) # WellKnownText{GeoFormatTypes.Geom}(GeoFormatTypes.Geom(), "POINT (30.0 10.0)")
GeoFormatTypes.val(wkb) # "POINT (30.0 10.0)"
WellKnownGeometry.getwkb
— Functiongetwkb(geom)
Retrieve the Well Known Binary (WKB) as GeoFormatTypes.WellKnownBinary
for a geom
that implements the GeoInterface. Use GeoFormatTypes.val
to get the Vector{UInt8} representation.
We thank Julia Computing for supporting contributions to this package.