From 138b7777c43d5606720c33dd6e461e67f6ef115d Mon Sep 17 00:00:00 2001 From: rafaqz Date: Wed, 27 Dec 2023 01:45:47 +0100 Subject: [PATCH] da --- ext/DimensionalDataInterpolations.jl | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) 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