Pure Julia WKT2 to PROJJSON conversion

src/GeoIO.jl

@@ -53,7 +53,7 @@ import GeoInterface as GI
 import GeoFormatTypes as GFT
 import ArchGDAL.GDAL
 
-# ProjJSON CRS
+# PROJJSON CRS
 import JSON3
 
 # VTK extensions
@@ -112,10 +112,13 @@ function formats(io=stdout; sortby=:extension)
   pretty_table(io, sorted, alignment=:c, crop=:none, show_subheader=false)
 end
 
-# utilities
+# basic utilities
 include("utils.jl")
 
-# conversions
+# utilities for CRS strings
+include("crsstrings.jl")
+
+# utilities for geometry conversion
 include("conversion.jl")
 
 # extra code for backends

src/conversion.jl

@@ -2,6 +2,9 @@
 # Licensed under the MIT License. See LICENSE in the project root.
 # ------------------------------------------------------------------
 
+raw(coords::CRS) = coords.x, coords.y
+raw(coords::LatLon) = coords.lon, coords.lat
+
 # --------------------------------------
 # Minimum GeoInterface.jl to perform IO
 # --------------------------------------
@@ -23,9 +26,6 @@ GI.ncoord(::GI.PointTrait, p::Point) = CoordRefSystems.ncoords(crs(p))
 GI.getcoord(::GI.PointTrait, p::Point) = ustrip.(raw(coords(p)))
 GI.getcoord(trait::GI.PointTrait, p::Point, i) = GI.getcoord(trait, p)[i]
 
-raw(coords::CRS) = coords.x, coords.y
-raw(coords::LatLon) = coords.lon, coords.lat
-
 GI.ncoord(::GI.LineTrait, s::Segment) = CoordRefSystems.ncoords(crs(s))
 GI.ngeom(::GI.LineTrait, s::Segment) = nvertices(s)
 GI.getgeom(::GI.LineTrait, s::Segment, i) = vertex(s, i)