-
Notifications
You must be signed in to change notification settings - Fork 27
Open
Description
When looking into a polygon, the constituent geometries print as linear rings but possess LineStringTrait
s. This is a problem "upstream" in GeometryOps when reconstructing geometries from ArchGDAL inputs.
Consider the following MWE:
import GeoInterface as GI, ArchGDAL as AG
polygon = GI.Polygon([GI.LinearRing([(cos(t), sin(t)) for t in LinRange(0, 2pi, 100)])])
GI.geomtrait(GI.getgeom(polygon, 1)) # LinearRingTrait
ag_polygon = GI.convert(AG, polygon)
GI.geomtrait(GI.getgeom(ag_polygon, 1)) # LineStringTrait
@which GI.geomtrait(GI.getgeom(ag_polygon, 1))
reveals:
ArchGDAL.jl/src/geointerface.jl
Lines 365 to 369 in 8d57512
function GeoInterface.geomtrait( | |
geom::Union{map(T -> AbstractGeometry{T}, linetypes)...}, | |
) | |
return GeoInterface.LineStringTrait() | |
end |
which seems incorrect...
yeesian
Metadata
Metadata
Assignees
Labels
No labels