Add geointerface_geomtypes as a method (#189)

* Add geointerface_geomtypes

Should get GeoInterface.convert(GeometryBasics, ...) working.  Needs testing.

* test converts to GeometryBasics

Author: asinghvi17

src/geointerface.jl

@@ -20,6 +20,15 @@ GeoInterface.geomtrait(::MultiPolygon) = MultiPolygonTrait()
 GeoInterface.geomtrait(::Ngon) = PolygonTrait()
 GeoInterface.geomtrait(::AbstractMesh) = PolyhedralSurfaceTrait()
 
+# GeoInterface calls this method in `GeoInterface.convert(GeometryBasics, ...)`
+geointerface_geomtype(::GeoInterface.PointTrait) = Point
+geointerface_geomtype(::GeoInterface.MultiPointTrait) = MultiPoint
+geointerface_geomtype(::GeoInterface.LineStringTrait) = LineString
+geointerface_geomtype(::GeoInterface.MultiLineStringTrait) = MultiLineString
+geointerface_geomtype(::GeoInterface.PolygonTrait) = Polygon
+geointerface_geomtype(::GeoInterface.MultiPolygonTrait) = MultiPolygon
+geointerface_geomtype(::GeoInterface.PolyhedralSurfaceTrait) = Mesh
+
 GeoInterface.geomtrait(::Simplex{Dim,T,1}) where {Dim,T} = PointTrait()
 GeoInterface.geomtrait(::Simplex{Dim,T,2}) where {Dim,T} = LineStringTrait()
 GeoInterface.geomtrait(::Simplex{Dim,T,3}) where {Dim,T} = PolygonTrait()

test/geointerface.jl

@@ -70,15 +70,15 @@ end
     multipolygon_json = GeoJSON.read(multipolygon_str)
     multipolygon_hole_json = GeoJSON.read(multipolygon_hole_str)
 
-    point_gb = GeoInterface.convert(Point, point_json)
-    point_3d_gb = GeoInterface.convert(Point, point_3d_json)
-    linestring_gb = GeoInterface.convert(LineString, linestring_json)
-    polygon_gb = GeoInterface.convert(Polygon, polygon_json)
-    polygon_hole_gb = GeoInterface.convert(Polygon, polygon_hole_json)
-    multipoint_gb = GeoInterface.convert(MultiPoint, multipoint_json)
-    multilinestring_gb = GeoInterface.convert(MultiLineString, multilinestring_json)
-    multipolygon_gb = GeoInterface.convert(MultiPolygon, multipolygon_json)
-    multipolygon_hole_gb = GeoInterface.convert(MultiPolygon, multipolygon_hole_json)
+    point_gb = GeoInterface.convert(GeometryBasics, point_json)
+    point_3d_gb = GeoInterface.convert(GeometryBasics, point_3d_json)
+    linestring_gb = GeoInterface.convert(GeometryBasics, linestring_json)
+    polygon_gb = GeoInterface.convert(GeometryBasics, polygon_json)
+    polygon_hole_gb = GeoInterface.convert(GeometryBasics, polygon_hole_json)
+    multipoint_gb = GeoInterface.convert(GeometryBasics, multipoint_json)
+    multilinestring_gb = GeoInterface.convert(GeometryBasics, multilinestring_json)
+    multipolygon_gb = GeoInterface.convert(GeometryBasics, multipolygon_json)
+    multipolygon_hole_gb = GeoInterface.convert(GeometryBasics, multipolygon_hole_json)
 
     @test point_gb === Point{2, Float64}(30.1, 10.1)
     @test point_3d_gb === Point{3, Float64}(30.1, 10.1, 5.1)