Rigid body transformation for big datasets

Inputs are a volume that can be uint8 and the transformation parameters (rotations, translations, scaling). Output is the linearly interpolated volume using the transformations. The transformation matrix is also returned. Not very fast, but does not require the grid as in intern, and works with integer. It is meant to transform a large volume of integer values.

% Big_Transform: rigid body transformation of a big dataset
%
% [Vi,M] = Big_Transform(V,p)
%
% V: Volume to transform, can be integer uint8.
% p: [ti tj tk ai aj ak si sj sk];
% (translation, rotation, scaling)
% Vi: new interpolated volume using trilinear interpolation
% M: transformation matrix in homogeneous coordinates
% if only one output variable: returns the transformation matrix only