mia-2dmyoicapgt(1)
โจ NAME โจ
okay so bestie, meet mia-2dmyoicapgt(1) โ it runs a registration of a series of 2D images. No cap, it's literally the moment. ๐
๐ธ SYNOPSIS ๐ธ
mia-2dmyoicapgt -i <in-file> -o <out-file> [options]
๐ฆ DESCRIPTION ๐ฆ
So mia-2dmyoicapgt implements a two passs motion compensation algorithm. (Yes, "passs" with three s's, we don't question it.) First a linear registration is run based on a variation of
Gupta et al. "Fully automatic registration and segmentation of first-pass myocardial perfusion MR image sequences", Academic Radiology 17, 1375-1385
as 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,
followed by a non-linear registration based on
Chao Li and Ying Sun, "Nonrigid Registration of Myocardial Perfusion MRI Using Pseudo Ground Truth", In Proc. Medical Image Computing and Computer-Assisted Intervention MICCAI 2009, 165-172, 2009.
Note that for this nonlinear motion correction a preceding linear registration step is usually required. This version of the program may run all registrations in parallel. ๐
๐ OPTIONS ๐
โจ Pseudo Ground Thruth estimation โจ
-A,--alpha=0.1โ spacial neighborhood penalty weight ๐ (the little weight that keeps your neighbors in check, bestie)-B,--beta=4โ temporal second derivative penalty weight โณ (it's giving "second derivative but make it penalty")-T,--rho-thresh=0.85โ correlation threshold for neighborhood analysis ๐ (the correlation cutoff for when the neighborhood gets analyzed)
๐ File-IO ๐
-i,--in-file=(input, required); stringโ input perfusion data set ๐พ (required, bestie, you gotta feed it something)-o,--out-file=(output, required); stringโ output perfusion data set ๐พ (required, where the magic lands)-r,--registered=โ File name base for the registered images. Image type and numbering scheme are taken from the input images as given in the input data set. ๐ผ๏ธ (the base name for your registered pics, and it steals the type and numbering from the input data set)--save-cropped=(output); stringโ save cropped set to this file, the image files will use the stem of the name as file name base โ๏ธ (crop it, save it, and the stem becomes the file name base)--save-feature=(output); stringโ save segmentation feature images and initial ICA mixing matrix ๐ง (feature images plus the initial ICA mixing matrix, slay)--save-refs=(output); stringโ for each registration pass save the reference images to files with the given name base ๐ (every registration pass gets its reference images saved)--save-regs=(output); stringโ for each registration pass save intermediate registered images ๐ (intermediate registered images for each pass)
๐ Help & Info ๐
-V,--verbose=warningโ verbosity of output, print messages of given level and higher priorities. Supported priorities starting at lowest level are:traceโ Function call tracedebugโ Debug outputinfoโ Low level messagesmessageโ Normal messageswarningโ Warningsfailโ Report test failureserrorโ Report errorsfatalโ Report only fatal errors
--copyrightโ print copyright information ยฉ๏ธ-h,--helpโ print this help ๐-?,--usageโ print a short help โ--versionโ print the version number and exit ๐ข
๐งฌ ICA ๐งฌ
--fastica=internalโ FastICA implementationto be used ๐งช For supported plugins see PLUGINS:fastica/implementation-C,--components=0โ ICA components 0 = automatic estimation ๐ข (zero means automatic, iykyk)--normalizeโ normalized ICs โจ--no-meanstripโ don't strip the mean from the mixing curves ๐ซ (mean stays, periodt)-s,--segscale=0โ segment and scale the crop box around the LV (0=no segmentation) ๐ซ (zero means no segmentation, bestie)-k,--skip=0โ skip images at the beginning of the series e.g. because as they are of other modalities โญ๏ธ (skip the first few, maybe they're other modalities)-m,--max-ica-iter=400โ maximum number of iterations in ICA ๐-E,--segmethod=featuresโ Segmentation method ๐งฉdelta-featureโ difference of the feature imagesdelta-peakโ difference of the peak enhancement imagesfeaturesโ feature images
-b,--min-breathing-frequency=-1โ 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. A value 0.0 forces the series to be identified as acquired with initial breath hold. ๐ซ (breathing vibes only; negative disables, 0.0 forces initial breath hold)
โ๏ธ Processing โ๏ธ
--threads=-1โ 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). ๐งต (max threads, and-1means automatic estimation)
๐ Registration ๐
-L,--linear-optimizer=gsl:opt=simplex,step=1.0โ Optimizer used for minimization of the linear registration The string value will be used to construct a plug-in. ๐งฎ For supported plugins see PLUGINS:minimizer/singlecost--linear-transform=affineโ linear transform to be used The string value will be used to construct a plug-in. ๐งต For supported plugins see PLUGINS:2dimage/transform-O,--non-linear-optimizer=gsl:opt=gd,step=0.1โ Optimizer used for minimization in the non-linear registration. The string value will be used to construct a plug-in. ๐งฎ For supported plugins see PLUGINS:minimizer/singlecost-a,--start-c-rate=16โ start coefficinet rate in spines, gets divided by --c-rate-divider with every pass. ๐--c-rate-divider=2โ Cofficient rate divider for each pass. โ-d,--start-divcurl=10000โ Start divcurl weight, gets divided by --divcurl-divider with every pass. โ๏ธ--divcurl-divider=2โ Divcurl weight scaling with each new pass. ๐-R,--reference=-1โ Global reference all image should be aligned to. If set to a non-negative value, the images will be aligned to this references, and the cropped output image date will be injected into the original images. Leave at -1 if you don't care. In this case all images with be registered to a mean position of the movement ๐ฏ (set-1if you don't care, and everything aligns to a mean position; non-negative means align to this reference and inject the cropped output image date into the originals)-w,--imagecost=image:weight=1,cost=ssdโ image cost, do not specify the src and ref parameters, these will be set by the program. The string value will be used to construct a plug-in. ๐ธ For supported plugins see PLUGINS:2dimage/fullcost-l,--mg-levels=3โ multi-resolution levels ๐ช-p,--linear-passes=3โ linear registration passes (0 to disable) ๐-P,--nonlinear-passes=3โ non-linear registration passes (0 to disable) ๐
๐ PLUGINS: 1d/splinebc ๐
mirrorโ Spline interpolation boundary conditions that mirror on the boundary ๐ช- (no parameters)
repeatโ Spline interpolation boundary conditions that repeats the value at the boundary ๐- (no parameters)
zeroโ Spline interpolation boundary conditions that assumes zero for values outside 0๏ธโฃ- (no parameters)
๐ PLUGINS: 1d/splinekernel ๐
bsplineโ B-spline kernel creation , supported parameters are: ๐งตd = 3; int in [0, 5]โ Spline degree.
omomsโ OMoms-spline kernel creation, supported parameters are: ๐งถd = 3; int in [3, 3]โ Spline degree.
๐ PLUGINS: 2dimage/cost ๐
lnccโ local normalized cross correlation with masking support., supported parameters are: ๐w = 5; uint in [1, 256]โ half width of the window used for evaluating the localized cross correlation.
lsdโ Least-Squares Distance measure ๐- (no parameters)
miโ Spline parzen based mutual information., supported parameters are: ๐cut = 0; float in [0, 40]โ Percentage of pixels to cut at high and low intensities to remove outliers.mbins = 64; uint in [1, 256]โ Number of histogram bins used for the moving image.mkernel = [bspline:d=3]; factoryโ Spline kernel for moving image parzen hinstogram. For supported plug-ins see PLUGINS:1d/splinekernelrbins = 64; uint in [1, 256]โ Number of histogram bins used for the reference image.rkernel = [bspline:d=0]; factoryโ Spline kernel for reference image parzen hinstogram. For supported plug-ins see PLUGINS:1d/splinekernel
nccโ normalized cross correlation. ๐- (no parameters)
ngfโ This function evaluates the image similarity based on normalized gradient fields. Various evaluation kernels are available., supported parameters are: ๐งeval = ds; dictโ plugin subtype. Supported values are:sqโ square of differencedsโ square of scaled differencedotโ scalar product kernelcrossโ cross product kernel
ssdโ 2D imaga cost: sum of squared differences, supported parameters are: ๐งฎautothresh = 0; float in [0, 1000]โ Use automatic masking of the moving image by only takeing intensity values into accound that are larger than the given threshold.norm = 0; boolโ Set whether the metric should be normalized by the number of image pixels.
ssd-automaskโ 2D image cost: sum of squared differences, with automasking based on given thresholds, supported parameters are: ๐ญrthresh = 0; doubleโ Threshold intensity value for reference image.sthresh = 0; doubleโ Threshold intensity value for source image.
โจ PLUGINS: 2dimage/fullcost โจ
image ๐
Generalized image similarity cost function that also does multi-resolution processing, bestie โ the actual similarity measure gets handed in as an extra parameter. Here's what you can pass:
cost = ssd; factoryโ the cost function kernel. For supported plug-ins seePLUGINS:2dimage/cost๐debug = 0; boolโ save the intermediate results so you can snoop on them for debugging ๐ต๏ธโ๏ธref =(input, io)โ the reference image. For supported file types seePLUGINS:2dimage/iosrc =(input, io)โ the study image. For supported file types seePLUGINS:2dimage/ioweight = 1; floatโ the weight of the cost function โ๏ธ
labelimage ๐ท๏ธ
Similarity cost function that maps the labels of two images, and it handles label-preserving multi-resolution processing too. Parameters, go:
debug = 0; int in [0, 1]โ write the distance transforms out to a 3D image ๐maxlabel = 256; int in [2, 32000]โ the maximum number of labels to considerref =(input, io)โ the reference image. For supported file types seePLUGINS:2dimage/iosrc =(input, io)โ the study image. For supported file types seePLUGINS:2dimage/ioweight = 1; floatโ the weight of the cost function โ๏ธ
maskedimage ๐ญ
Generalized masked image similarity cost function, and yes, it also handles multi-resolution processing. Heads up though: the masks you feed it should be densely filled regions when you're doing multi-resolution processing, because otherwise the mask info can literally get lost when the image gets downscaled ๐ญ. The reference mask and the transformed mask of the study image get combined with a binary AND. The actual similarity measure comes in as an extra parameter, and here's the menu:
cost = ssd; factoryโ the cost function kernel. For supported plug-ins seePLUGINS:2dimage/maskedcostref =(input, io)โ the reference image. For supported file types seePLUGINS:2dimage/ioref-mask =(input, io)โ the reference image mask (binary). For supported file types seePLUGINS:2dimage/iosrc =(input, io)โ the study image. For supported file types seePLUGINS:2dimage/iosrc-mask =(input, io)โ the study image mask (binary). For supported file types seePLUGINS:2dimage/ioweight = 1; floatโ the weight of the cost function โ๏ธ
๐ PLUGINS: 2dimage/io ๐
bmp ๐ผ๏ธ
BMP 2D-image input/output support, periodt. This plug-in reads and writes binary images and 8-bit gray scale images, and it gives you read-only support for 4-bit gray scale images. The color table gets straight-up ignored and the pixel values are taken as literal gray scale values ๐
- Recognized file extensions:
.BMP,.bmp - Supported element types:
- binary data, unsigned 8 bit
datapool ๐๏ธ
Virtual IO to and from the internal data pool โจ
- Recognized file extensions:
.@
dicom ๐ฅ
2D image io for DICOM
- Recognized file extensions:
.DCM,.dcm - Supported element types:
- signed 16 bit, unsigned 16 bit
exr ๐๏ธ
a 2dimage io plugin for OpenEXR images
- Recognized file extensions:
.EXR,.exr - Supported element types:
- unsigned 32 bit, floating point 32 bit
jpg ๐ธ
a 2dimage io plugin for jpeg gray scale images
- Recognized file extensions:
.JPEG,.JPG,.jpeg,.jpg - Supported element types:
- unsigned 8 bit
png ๐ฆ
a 2dimage io plugin for png images
- Recognized file extensions:
.PNG,.png - Supported element types:
- binary data, unsigned 8 bit, unsigned 16 bit
raw ๐งฑ
RAW 2D-image output support
- Recognized file extensions:
.RAW,.raw - Supported element types:
- binary data, signed 8 bit, unsigned 8 bit, signed 16 bit, unsigned 16 bit, signed 32 bit, unsigned 32 bit, floating point 32 bit, floating point 64 bit
tif ๐จ๏ธ
TIFF 2D-image input/output support
- Recognized file extensions:
.TIF,.TIFF,.tif,.tiff - Supported element types:
- binary data, unsigned 8 bit, unsigned 16 bit, unsigned 32 bit
vista ๐
a 2dimage io plugin for vista images
- Recognized file extensions:
.-,.V,.VISTA,.v,.vista - Supported element types:
- binary data, signed 8 bit, unsigned 8 bit, signed 16 bit, unsigned 16 bit, signed 32 bit, unsigned 32 bit, floating point 32 bit, floating point 64 bit
๐ PLUGINS: 2dimage/maskedcost ๐
lncc ๐
local normalized cross correlation with masking support, and it's giving parameters:
w = 5; uint in [1, 256]โ half width of the window used for evaluating the localized cross correlation
mi ๐ง
Spline parzen based mutual information with masking. The parameters:
cut = 0; float in [0, 40]โ percentage of pixels to cut at high and low intensities, so the outliers get yeeted โ๏ธmbins = 64; uint in [1, 256]โ number of histogram bins used for the moving imagemkernel = [bspline:d=3]; factoryโ spline kernel for the moving image parzen hinstogram. For supported plug-ins seePLUGINS:1d/splinekernelrbins = 64; uint in [1, 256]โ number of histogram bins used for the reference imagerkernel = [bspline:d=0]; factoryโ spline kernel for the reference image parzen hinstogram. For supported plug-ins seePLUGINS:1d/splinekernel
ncc ๐
normalized cross correlation with masking support.
- (no parameters) โ yeah, literally none, she's simple like that ๐
ssd โ
Sum of squared differences with masking.
- (no parameters) ๐ซถ
๐ PLUGINS: 2dimage/transform ๐
affine ๐
Affine transformation (six degrees of freedom), and the parameters are:
imgboundary = mirror; factoryโ image interpolation boundary conditions. For supported plug-ins seePLUGINS:1d/splinebcimgkernel = [bspline:d=3]; factoryโ image interpolator kernel. For supported plug-ins seePLUGINS:1d/splinekernel
rigid ๐
Rigid transformations (i.e. rotation and translation, three degrees of freedom), with parameters:
imgboundary = mirror; factoryโ image interpolation boundary conditions. For supported plug-ins seePLUGINS:1d/splinebcimgkernel = [bspline:d=3]; factoryโ image interpolator kernel. For supported plug-ins seePLUGINS:1d/splinekernelrot-center = [[0,0]]; 2dfvectorโ relative rotation center, i.e.<0.5,0.5>corresponds to the center of the support rectangle ๐ฏ
rotation ๐
Rotation transformations (i.e. rotation about a given center, one degree of freedom), with parameters:
imgboundary = mirror; factoryโ image interpolation boundary conditions. For supported plug-ins seePLUGINS:1d/splinebcimgkernel = [bspline:d=3]; factoryโ image interpolator kernel. For supported plug-ins seePLUGINS:1d/splinekernelrot-center = [[0,0]]; 2dfvectorโ relative rotation center, i.e.<0.5,0.5>corresponds to the center of the support rectangle ๐ฏ
spline ๐
Free-form transformation that can be described by a set of B-spline coefficients and an underlying B-spline kernel. Parameters:
anisorate = [[0,0]]; 2dfvectorโ anisotropic coefficient rate in pixels, nonpositive values will be overwritten by theratevalue..imgboundary = mirror; factoryโ image interpolation boundary conditions. For supported plug-ins seePLUGINS:1d/splinebcimgkernel = [bspline:d=3]; factoryโ image interpolator kernel. For supported plug-ins seePLUGINS:1d/splinekernelkernel = [bspline:d=3]; factoryโ transformation spline kernel.. For supported plug-ins seePLUGINS:1d/splinekernelpenalty = ; factoryโ transformation penalty term. For supported plug-ins seePLUGINS:2dtransform/splinepenaltyrate = 10; float in [1, inf)โ isotropic coefficient rate in pixels
translate โก๏ธ
Translation only (two degrees of freedom), and the parameters:
imgboundary = mirror; factoryโ image interpolation boundary conditions. For supported plug-ins seePLUGINS:1d/splinebcimgkernel = [bspline:d=3]; factoryโ image interpolator kernel. For supported plug-ins seePLUGINS:1d/splinekernel
vf ๐งต
This plug-in implements a transformation that defines a translation for each point of the grid defining the domain of the transformation, and the parameters are:
imgboundary = mirror; factoryโ image interpolation boundary conditions. For supported plug-ins seePLUGINS:1d/splinebcimgkernel = [bspline:d=3]; factoryโ image interpolator kernel. For supported plug-ins seePLUGINS:1d/splinekernel
๐ PLUGINS: 2dtransform/splinepenalty ๐
okay bestie so this one is all about the divcurl penalty on the transformation ๐ซ and here are the parameters it supports:
curlโcurl = 1; float in [0, inf), which is the penalty weight on curl. yes, curl like the vibes, but numbers. ๐divโdiv = 1; float in [0, inf), the penalty weight on divergence. twin to the one above, fr fr.normโnorm = 0; bool. set it to1if the penalty should be normalized with respect to the image size ๐ (read the note carefully, babes, it's aboolnot a float).weightโweight = 1; float in (0, inf), the weight of penalty energy. note the parens: it's a strictly positive range, no zero allowed, iykyk. ๐ฅ
โจ PLUGINS: fastica/implementation โจ
the FastICA algorithm said "pick your fighter" and honestly? both slay ๐. here's the tea:
internalโ this is the MIA implementation of the FastICA algorithm. homegrown, iconic, no cap. ๐
(no parameters)
itppโ this is the IT++ implementation of the FastICA algorithm. guest star energy, still ate. ๐ซถ
(no parameters)
๐ PLUGINS: minimizer/singlecost ๐
Okay bestie, this is the other minimizer plugin family โ the single-cost one. Let's get into it ๐ซถ
๐ธ gdas ๐ธ
Gradient descent with automatic step size correction. supported parameters are:
ftolr=0;double in [0, inf)โ stop if the relative change of the criterion is below. ๐max-step=2;double in (0, inf)โ the maximal absolute step size.maxiter=200;uint in [1, inf)โ stopping criterion: the maximum number of iterations. ๐min-step=0.1;double in (0, inf)โ the minimal absolute step size.xtola=0.01;double in [0, inf)โ stop if the inf-norm of the change applied to x is below this value. ๐
๐ธ gdsq ๐ธ
Gradient descent with quadratic step estimation, supported parameters are:
ftolr=0;double in [0, inf)โ stop if the relative change of the criterion is below. ๐gtola=0;double in [0, inf)โ stop if the inf-norm of the gradient is below this value.maxiter=100;uint in [1, inf)โ stopping criterion: the maximum number of iterations. ๐scale=2;double in (1, inf)โ fallback fixed step size scaling.step=0.1;double in (0, inf)โ initial step size.xtola=0;double in [0, inf)โ stop if the inf-norm of x-update is below this value.
๐ธ gsl ๐ธ
optimizer plugin based on the multimin optimizers of the GNU Scientific Library (GSL) https://www.gnu.org/software/gsl/, supported parameters are:
eps=0.01;double in (0, inf)โ gradient based optimizers: stop when|grad| < eps, simplex: stop when simplex size< eps. ๐iter=100;uint in [1, inf)โ maximum number of iterations.opt=gd;dictโ specific optimizer to be used. Supported values are:simplexโ Simplex algorithm of Nelder and Meadcg-frโ Flecher-Reeves conjugate gradient algorithmcg-prโ Polak-Ribiere conjugate gradient algorithmbfgsโ Broyden-Fletcher-Goldfarb-Shannbfgs2โ Broyden-Fletcher-Goldfarb-Shann (most efficient version)gdโ Gradient descent.
step=0.001;double in (0, inf)โ initial step size.tol=0.1;double in (0, inf)โ some tolerance parameter.
๐ธ nlopt ๐ธ
Minimizer algorithms using the NLOPT library, for a description of the optimizers please see http://ab-initio.mit.edu/wiki/index.php/NLopt_Algorithms, supported parameters are:
ftola=0;double in [0, inf)โ stopping criterion: the absolute change of the objective value is below this value.ftolr=0;double in [0, inf)โ stopping criterion: the relative change of the objective value is below this value.higher=inf;doubleโ higher boundary (equal for all parameters). โฌ๏ธlocal-opt=none;dictโ local minimization algorithm that may be required for the main minimization algorithm. Supported values are:gn-directโ Dividing Rectanglesgn-direct-lโ Dividing Rectangles (locally biased)gn-direct-l-randโ Dividing Rectangles (locally biased, randomized)gn-direct-noscalโ Dividing Rectangles (unscaled)gn-direct-l-noscalโ Dividing Rectangles (unscaled, locally biased)gn-direct-l-rand-noscaleโ Dividing Rectangles (unscaled, locally biased, randomized)gn-orig-directโ Dividing Rectangles (original implementation)gn-orig-direct-lโ Dividing Rectangles (original implementation, locally biased)ld-lbfgs-nocedalโ Noneld-lbfgsโ Low-storage BFGSln-praxisโ Gradient-free Local Optimization via the Principal-Axis Methodld-var1โ Shifted Limited-Memory Variable-Metric, Rank 1ld-var2โ Shifted Limited-Memory Variable-Metric, Rank 2ld-tnewtonโ Truncated Newtonld-tnewton-restartโ Truncated Newton with steepest-descent restartingld-tnewton-precondโ Preconditioned Truncated Newtonld-tnewton-precond-restartโ Preconditioned Truncated Newton with steepest-descent restartinggn-crs2-lmโ Controlled Random Search with Local Mutationld-mmaโ Method of Moving Asymptotesln-cobylaโ Constrained Optimization BY Linear Approximationln-newuoaโ Derivative-free Unconstrained Optimization by Iteratively Constructed Quadratic Approximationln-newuoa-boundโ Derivative-free Bound-constrained Optimization by Iteratively Constructed Quadratic Approximationln-neldermeadโ Nelder-Mead simplex algorithmln-sbplxโ Subplex variant of Nelder-Meadln-bobyqaโ Derivative-free Bound-constrained Optimizationgn-isresโ Improved Stochastic Ranking Evolution Strategynoneโ don't specify algorithm
lower=-inf;doubleโ lower boundary (equal for all parameters). โฌ๏ธmaxiter=100;int in [1, inf)โ stopping criterion: the maximum number of iterations. ๐opt=ld-lbfgs;dictโ main minimization algorithm. Supported values are:gn-directโ Dividing Rectanglesgn-direct-lโ Dividing Rectangles (locally biased)gn-direct-l-randโ Dividing Rectangles (locally biased, randomized)gn-direct-noscalโ Dividing Rectangles (unscaled)gn-direct-l-noscalโ Dividing Rectangles (unscaled, locally biased)gn-direct-l-rand-noscaleโ Dividing Rectangles (unscaled, locally biased, randomized)gn-orig-directโ Dividing Rectangles (original implementation)gn-orig-direct-lโ Dividing Rectangles (original implementation, locally biased)ld-lbfgs-nocedalโ Noneld-lbfgsโ Low-storage BFGSln-praxisโ Gradient-free Local Optimization via the Principal-Axis Methodld-var1โ Shifted Limited-Memory Variable-Metric, Rank 1ld-var2โ Shifted Limited-Memory Variable-Metric, Rank 2ld-tnewtonโ Truncated Newtonld-tnewton-restartโ Truncated Newton with steepest-descent restartingld-tnewton-precondโ Preconditioned Truncated Newtonld-tnewton-precond-restartโ Preconditioned Truncated Newton with steepest-descent restartinggn-crs2-lmโ Controlled Random Search with Local Mutationld-mmaโ Method of Moving Asymptotesln-cobylaโ Constrained Optimization BY Linear Approximationln-newuoaโ Derivative-free Unconstrained Optimization by Iteratively Constructed Quadratic Approximationln-newuoa-boundโ Derivative-free Bound-constrained Optimization by Iteratively Constructed Quadratic Approximationln-neldermeadโ Nelder-Mead simplex algorithmln-sbplxโ Subplex variant of Nelder-Meadln-bobyqaโ Derivative-free Bound-constrained Optimizationgn-isresโ Improved Stochastic Ranking Evolution Strategyauglagโ Augmented Lagrangian algorithmauglag-eqโ Augmented Lagrangian algorithm with equality constraints onlyg-mlslโ Multi-Level Single-Linkage (require local optimization and bounds)g-mlsl-ldsโ Multi-Level Single-Linkage (low-discrepancy-sequence, require local gradient based optimization and bounds)ld-slsqpโ Sequential Least-Squares Quadratic Programming
step=0;double in [0, inf)โ initial step size for gradient free methods.stop=-inf;doubleโ stopping criterion: function value falls below this value.xtola=0;double in [0, inf)โ stopping criterion: the absolute change of all x-values is below this value.xtolr=0;double in [0, inf)โ stopping criterion: the relative change of all x-values is below this value.
โจ EXAMPLE โจ
Okay so picture this ๐ญ โ you want to register the perfusion series given in segment.set by first using automatic ICA estimation to run the linear registration and then the PGT registration. You skip two images at the beginning and otherwise use the default parameters. And you store the result in registered.set. Iconic, no notes. ๐
mia-2dmyoicapgt -i segment.set -o registered.set -k 2
๐ AUTHOR(s) ๐
Gert Wollny
๐ COPYRIGHT ๐
This software is Copyright (c) 1999โ2015 Leipzig, Germany and Madrid, Spain. It comes with ABSOLUTELY NO WARRANTY and you may redistribute it under the terms of the GNU GENERAL PUBLIC LICENSE Version 3 (or later). For more information run the program with the option --copyright.
fr fr, that warranty disclaimer is not a vibe, it's a warning โ no warranty, periodt. ๐
That's a wrap on the whole page bestie ๐ฆ this one was from USER COMMANDS, version v2.4.7, living in the mia-2dmyoicapgt(1) section โ go register something gorgeous and slay ๐
โจ๐ธ
๐ฆ where this page comes from
- ubuntu noble: mia-tools 2.4.7-13build4 ยท GPL-3.0-or-later
- ubuntu resolute: mia-tools 2.4.7-18 ยท GPL-3.0-or-later
- ubuntu stonking: mia-tools 2.4.7-19build1 ยท GPL-3.0-or-later
- ๐ girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms