This program implements the motion compensation algorithm described in Wollny G, Kellman P, Santos A, Ledesma-Carbayo M-J, "Automatic Motion Compensation of Free Breathing acquired Myocardial Perfusion Data by using Independent Component Analysis" Medical Image Analysis, 2012, DOI:10.1016/j.media.2012.02.004.
input perfusion data set
output perfusion data set
file name base for registered fiels
save cropped set to this file
save the features images resulting from the ICA and some intermediate images used for the RV-LV segmentation with the given file name base to PNG files. Also save the coefficients of the initial best and the final IC mixing matrix.
save synthetic reference images
save intermediate registered images
verbosity of output, print messages of given level and higher priorities. Supported priorities starting at lowest level are:
|info:||Low level messages|
|trace:||Function call trace|
|fail:||Report test failures|
|fatal:||Report only fatal errors|
print copyright information
print this help
print a short help
print the version number and exit
ICA components 0 = automatic estimation
don't strip the mean from the mixing curves
segment and scale the crop box around the LV (0=no segmentation)
skip images at the beginning of the series e.g. because as they are of other modalities
maximum number of iterations in ICA
|delta-peak:||difference of the peak enhancement images|
|delta-feature:||difference of the feature images|
minimal mean frequency a mixing curve can have to be considered to stem from brething. A healthy rest breating rate is 12 per minute. A negative value disables the test.
Maxiumum number of threads to use for processing,This number should be lower or equal to the number of logical processor cores in the machine. (-1: automatic estimation).
Optimizer used for minimization. For supported plug-ins see Plugin type: minimizer/singlecost
optimizer used for refinement after the main optimizer was called. For supported plug-ins see Plugin type: minimizer/singlecost
start coefficinet rate in spines, gets divided by --c-rate-divider with every pass
cofficient rate divider for each pass
start divcurl weight, gets divided by --divcurl-divider with every pass
divcurl weight scaling with each new pass
image cost. For supported plug-ins see Plugin type: 2dimage/fullcost
Register the perfusion series given in 'segment.set' by using automatic ICA estimation. Skip two images at the beginning and otherwiese use the default parameters. Store the result in 'registered.set'.
mia-2dmyoica-nonrigid -i segment.set -o registered.set -k 2