You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
When looking into a polygon, the constituent geometries print as linear rings but possess LineStringTraits. 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 inLinRange(0, 2pi, 100)])])
GI.geomtrait(GI.getgeom(polygon, 1)) # LinearRingTrait
ag_polygon = GI.convert(AG, polygon)
GI.geomtrait(GI.getgeom(ag_polygon, 1)) # LineStringTrait
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:
@which GI.geomtrait(GI.getgeom(ag_polygon, 1))
reveals:ArchGDAL.jl/src/geointerface.jl
Lines 365 to 369 in 8d57512
which seems incorrect...
The text was updated successfully, but these errors were encountered: