diff --git a/ext/DimensionalDataInterpolations.jl b/ext/DimensionalDataInterpolations.jl index cdf4ce451..9285b5d03 100644 --- a/ext/DimensionalDataInterpolations.jl +++ b/ext/DimensionalDataInterpolations.jl @@ -12,7 +12,8 @@ function interp(A::AbstractDimArray; isempty(otherdims(to, dims(A))) || throw(DimensionMismatch("Cannot interpolate over dimensions not in the source array")) # TODO handle permutations shared_dims = commondims(A, to) - dest_dims = dims(Rasters._extent2dims(commondims(to, shared_dims); size, res, crs=crs(A)), shared_dims) + # TODO move `_extent2dims` to DimensionalData` + dest_dims = dims(Rasters._extent2dims(commondims(to, shared_dims); size, res), shared_dims) other_dims = otherdims(A, to) degrees = if isnothing(degree) map(_ -> Linear(), dims(A)) @@ -41,7 +42,7 @@ function interp(A::AbstractDimArray; error("interpolate is only implemented for regular grids") end - return Raster(data, dest_dims) + return DimArray(data, dest_dims) end end