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 I have a geometry in a specific CRS I can get the crs information via ArchGDAL.getspatialref(geom) but then I can't use this crs for the reprojection.
It would be good to be able to use the crs that I get from ArchGDAL in the relevant ArchGDAL functions.
I can get the projstring via ArchGDAL.toProj4 but then to use it in the reproject call I have to wrap it in the ProjString type.
julia> deu = GADM.get("DEU", depth=0)
(geom = ArchGDAL.IGeometry{ArchGDAL.wkbMultiPolygon}[Geometry: MULTIPOLYGON (((7.9864206310001647.5555305480001, ...02)))], GID_0 = ["DEU"], COUNTRY = ["Germany"])
julia> crs = ArchGDAL.getspatialref(deu.geom[1])
Spatial Reference System:+proj=longlat +datum=WGS84 +no_defs
# This works: but could be made easier
julia> projdeu = ArchGDAL.reproject(only(deu.geom), ProjString(ArchGDAL.toPROJ4(crs)), ProjString(newproj))
Geometry: MULTIPOLYGON (((17.199717191748-0.16543096960377, ...47)))
julia> projdeu = ArchGDAL.reproject(only(deu.geom), crs, ProjString(newproj))
ERROR: MethodError: no method matching reproject(::ArchGDAL.IGeometry{ArchGDAL.wkbMultiPolygon}, ::ArchGDAL.ISpatialRef, ::ProjString)
Closest candidates are:reproject(::ArchGDAL.AbstractGeometry, ::GeoFormat, ::GeoFormat; kwargs...)
@ ArchGDAL ~/Documents/NFDI4Earth/deliverables/D2.5_1-5_Overview_DC_tech/dev/ArchGDAL/src/spatialref.jl:155reproject(::T, ::GeoFormat, ::Nothing; kwargs...) where T
@ ArchGDAL ~/Documents/NFDI4Earth/deliverables/D2.5_1-5_Overview_DC_tech/dev/ArchGDAL/src/spatialref.jl:100reproject(::Union{Tuple{var"#s244", var"#s244"} where var"#s244"<:Number, Tuple{var"#s243", var"#s243", var"#s243"} where var"#s243"<:Number, AbstractVector{<:Number}}, ::GeoFormat, ::GeoFormat; kwargs...)
@ ArchGDAL ~/Documents/NFDI4Earth/deliverables/D2.5_1-5_Overview_DC_tech/dev/ArchGDAL/src/spatialref.jl:114...
Stacktrace:
[1] top-level scope
@ REPL[160]:1
The text was updated successfully, but these errors were encountered:
When I have a geometry in a specific CRS I can get the crs information via ArchGDAL.getspatialref(geom) but then I can't use this crs for the reprojection.
It would be good to be able to use the crs that I get from ArchGDAL in the relevant ArchGDAL functions.
I can get the projstring via ArchGDAL.toProj4 but then to use it in the reproject call I have to wrap it in the ProjString type.
The text was updated successfully, but these errors were encountered: