shapes = AG.registerdrivers() do
AG.read("Nonmotorized_Transportation_Plan.shp") do dataset
[
AG.getgeom(feature) for feature in AG.getlayer(dataset, 0)
if GDAL.getgeometryref(feature.ptr) != C_NULL
]
end
end
from there, to project them, the following should work:
AG.createcoordtrans(AG.importEPSG(2927), AG.importEPSG(4326)) do coordtrans
AG.transform!.(shapes, coordtrans)
end
using https://pcjericks.github.io/py-gdalogr-cookbook/projection.html#reproject-a-layer as an example. e.g. pull out geometries from a shapefile in ArchGDAL using http://gisdata-piercecowa.opendata.arcgis.com/datasets/nonmotorized-transportation-plan
from there, to project them, the following should work:
alternatively: