ImageJ2 script to compose displacement fields

composeDfields.groovy

#@ Dataset da
#@ Dataset db

dap = Views.moveAxis( da, 2, 0 );
dbp = Views.moveAxis( da, 2, 0 );
spacing = [1, 1, 1] as double[];

dfA = new DisplacementFieldTransform( dap, spacing );
dfB = new DisplacementFieldTransform( dbp, spacing );

seq = new RealTransformSequence();
seq.add( dfA );
seq.add( dfB );

sz = dap.dimensionsAsLongArray();
itvl = new FinalInterval( sz[1..<sz.size()] as long[] );
compRaw = DisplacementFieldTransform.createDisplacementField( seq, itvl, spacing );
comp = Views.moveAxis( compRaw, 0, 2 );

ImageJFunctions.show( comp );


import net.imglib2.*;
import net.imglib2.view.*;
import net.imglib2.realtransform.*;
import net.imglib2.type.numeric.real.*;
import net.imglib2.img.display.imagej.ImageJFunctions;
import net.imglib2.util.*;