#
###############################################################################
+source `dirname $0`/midp_common.sh
+
################# BLUTDIR #####################
registration_blutdir()
blutdir_params="--levels $nb_levels --metric $metric --optimizer $optimizer --samples $nb_samples --spacing $spacing,$spacing,$spacing --bins $hist_bins --maxIt $nb_iter --interp $interpolator --verbose"
cmd="clitkBLUTDIR -r $reference -t $target -m $mask_ref --targetMask $mask_targ --vf $vf -o $result $blutdir_params"
$cmd > $log
+
+ abort_on_error registration_blutdir $? clean_up_registration
}
################# ELASTIX #####################
echo "Computing ELASTIX $reference -> $target ..."
exec_dir=`which elastix`
exec_dir=`dirname $exec_dir`
+ suffix=${nb_samples}_${nb_iter}_${nb_levels}
cat $exec_dir/params_BSpline.txt | sed -e "s+<NbIterations>+$nb_iter+" \
-e "s+<LabelsFile>++" \
-e "s+<HistBins>+$hist_bins+" \
-e "s+<Levels>+$nb_levels+" \
-e "s+<NbSamples>+$nb_samples+" \
-e "s+<SamplerType>+$sampling_algo+" \
- -e "s+<Spacing>+$spacing+" > params_BSpline.txt
+ -e "s+<Spacing>+$spacing+" > params_BSpline_${suffix}.txt
vf_dir=`dirname $vf`
vf_base=`basename $vf .mhd`
result_base=`basename $result .mhd`
# image registration
- cmd="elastix -f $reference -m $target -fMask $mask_ref -mMask $mask_targ -out $result_dir -p params_BSpline.txt"
+ cmd="elastix -f $reference -m $target -fMask $mask_ref -mMask $mask_targ -out $result_dir -p params_BSpline_${suffix}.txt"
$cmd > /dev/null
+ abort_on_error registration_elastix $? clean_up_registration
# generate vector field
cmd="transformix -tp $result_dir/TransformParameters.0.txt -out $vf_dir -def all"
$cmd > /dev/null
+ abort_on_error registration_elastix $? clean_up_registration
# post-processing
mv $vf_dir/deformationField.mhd $vf