diff --git a/nitransforms/linear.py b/nitransforms/linear.py
index fcee7432..dc231085 100644
--- a/nitransforms/linear.py
+++ b/nitransforms/linear.py
@@ -492,7 +492,7 @@ def apply(
                 (
                     dataobj
                     if dataobj is not None
-                    else np.asanyarray(spatialimage.dataobj[..., t], dtype=input_dtype)
+                    else spatialimage.dataobj[..., t].astype(input_dtype, copy=False),
                 ),
                 yvoxels.T,
                 output=output_dtype,