bugfix extend

src/methods/crop_extend.jl

@@ -174,6 +174,9 @@ end
 function _extend_to(A::AbstractRaster, to::DimTuple;
     filename=nothing, suffix=nothing, touches=false
 )
+    others = otherdims(to, A)
+    # Allow not specifying all dimensions
+    to = (set(dims(A), map(=>, dims(A, to), to)...)..., others...)
     # Calculate the range of the old array in the extended array
     ranges = _without_mapped_crs(A) do A
         _without_mapped_crs(to) do to
@@ -192,6 +195,18 @@ function _extend_to(A::AbstractRaster, to::DimTuple;
         suffix, parent=parent(A), missingval=missingval(A),
         name=name(A), metadata=metadata(A)
     )
+    # Input checks
+    map(dims(A, to), dims(newA, to)) do d1, d2
+        if lookup(d1) isa Union{AbstractSampled,NoLookup}
+            b1, b2 = bounds(d1), bounds(d2)
+            b1[1] >= b2[1] || throw(ArgumentError("Lower bound of $(basetypeof(d1)) lookup of `$(b2[1])` are not larger than the original `$(b1[1])`"))
+            b1[2] <= b2[2] || throw(ArgumentError("Upper bound of $(basetypeof(d2)) lookup of `$(b2[2])` is not larger than the original `$(b1[2])`"))
+        elseif lookup(d1) isa Categorical
+            map(lookup(d1)) do x 
+                x in d2 || throw(ArgumentError("category $x not in new dimension"))
+            end
+        end
+    end
     # The missingval may have changed for disk-based arrays
     if !isequal(missingval(A), missingval(newA))
         A = replace_missing(A, missingval(newA))