Pure Julia WKT2 to PROJJSON conversion (#156)

* add testing against json schema

* half-complete implementation for processing GEOGCRS

* use DeepDiffs.jl for semantic comparison between nested Dicts

* add processing for GEOGCRS(Variant 1: ENSEMBLE)

* add processing for GEOGCRS(Variant 2: DATUM)

* handle UNIT optionally occurring in coordinate_system or each axis

* testing utils

* testing

* add processing for GEODCRS

* add processing for DYNAMIC entries is GEOG and GEOD.

These end up effecting the "datum" json items

* uncommitted changes

* add processing for PROJCRS

* replace the actual projjsonstring function

* remove ArchGDAL dependency

* Revert "remove ArchGDAL dependency"

This reverts commit 33fac70.
turns out ArchGDAL is still needed for other functionality

* solid cleanup round

* refactor testing utils

* side step EPSGs unsupported by the WKT dataset

* Fix 8% of cases where axis.direction should not start with uppercase

* add support for meridian and prime_meridian optional nodes

* support custom units for measurments (Clarke's foot??)

* full unit support

* ellipsoid cleanup

* git revert hot fix

* finalize testsuite

* testset fix

* reorganize wkt2json.jl

* rename wktdict2jsondict to wkt2json

* simple test to make sure our main callsite is tested

* cleanup test/jsonutils.jl

* remove test dependency on DeepDiffs as it still had few false-negatives

DeepDiffs and PrettyPrinting are still vastly useful when developing

* remove import DeepDiffs

* reintroduce changes after merging with master

* turns out spatialref is still needed for extra/gdal.jl

* change file locations

* rename vars and helper functions to be more inline with style

* adjust and rename itembykey (to finditem)

* various improvements to comments

* delete commented code

* run JuliaFormatter

* various comments, error msgs, and code cleanups

* one round of cleanups for test/jsonutils.jl

* cleanup debug_json (now check_projjson) test/jsonutils.jl

* test/jsonutils.jl renames

* style renames

* refactor wkt2json_datum and change it's encapsulation

* mass consistency renames

* error msg improvements

* remove 80% of remaining underscores to fit project style

* Move test/jsonutils.jl to outside testfiles

* Rename test/jsonutils.jl to test/testutils.jl

* Move all test utilities to test/testutils.jl

* Miscellaneous cleanup

* More refactoring

* Fix typo

* Rename src/projjson.jl to src/crsstrings.jl

* Refactor src/conversion.jl

* Refactor tests

* Refactor tests

* Minor adjustments

* follow style, remove function return type

* minor adjustments as per PR review

* organize and elaborate on testing edge cases

* remove helper debugging functions

* Final cleanup

* Final code style fixes

* Move functions around

* Additional cleanup

* jsonroundtrip is not needed anymore. Do not leave any SubStrings in the jsondict

I found the reason `jsonroundtrip` was needed for `JSONSchema.validate`
to work properly. Turns out there is a bug in JSONSchema that makes it
faultily deal with the underlying String behind a SubStrings

---------

Co-authored-by: Júlio Hoffimann <[email protected]>

Author: Omar-Elrefaei

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)