How to access the points in a polygon?

[maxfreu]

I just don't get it (and can't find the relevant passage in the docs), please help:

simplepolygon = ArchGDAL.createpolygon([(0.,0.), (0.,1.), (1.,1.)])
AG.getpoint(simplepolygon,1)  # CE Failure Code 6, incompatible geometry for operation
simplepolygon[1]  # no getindex
first(simplepolygon)  # no iterate

In the end I want to convert the polygons to their julia counterpart from Meshes.jl.

[visr]

It looks like you can use getgeom for this:

ArchGDAL.jl/src/geointerface.jl

Lines 69 to 72 in e401e00

	elseif gtype in multilinetypes || gtype in polygontypes
	Vector{Vector{Float64}}[
	GeoInterface.coordinates(getgeom(g, i - 1)) for i in 1:ngeom(g)
	]

In the end I want to convert the polygons to their julia counterpart from Meshes.jl.

A completely renewed https://github.com/JuliaGeo/GeoInterface.jl/ was just released, with a pending PR in #290. If Meshes would also be open to add this interface (it has 0 dependencies), conversion could be a simple convert.

[maxfreu]

Omg, just now I read that getgeom indexing starts at 0... And I wondered why I always get a NULL Geometry. Shouldn't it simply raise an out of bounds error? Anyway, thanks for the hint & I'm looking forward to the new implementation!

[visr]

That would probably be a good improvement, yes. I assume it would be a cheap check.

[maxfreu]

The initial issue has been solved, so I think this can be closed :)