Pre-allocated interpolation

I am using AffineTransform to transform a large number of volumes and the interpolation process allocated a huge amount of memory. I know there is a `interpolate!` function. However `interpolate!` would destroy the input matrix. Also it looks like it can only be used with `BSpline(Quadratic(InPlace()))`. Is there a way to output the interpolation to a pre-allocated matrix but also keep the input matrix? Thanks!