[Insight-users] Fwd: How to use the class MRIBiasCorrection
Oliver Trebbe
otrebbe at uni-muenster.de
Mon Feb 16 10:21:20 EST 2009
Hi Luis,
can´t find the papers...
404 not found
(link broken?)
would like to have a look at them...
Regards
Oliver
Luis Ibanez schrieb:
>
>
> ---------- Forwarded message ----------
> From: *Luis Ibanez* <luis.ibanez at kitware.com
> <mailto:luis.ibanez at kitware.com>>
> Date: Sun, Feb 15, 2009 at 3:08 PM
> Subject: Re: How to use the class MRIBiasCorrection
> To: orientation f <fuzengliang435 at gmail.com
> <mailto:fuzengliang435 at gmail.com>>
>
>
>
> Hi Orientation f,
>
> Please read the REAME file in the
>
> InsightApplications/MRIBiasCorrection
>
> and the associated papers:
>
> http://www.ia.unc.edu/public/styner/docs/tmi00.pdf
> http://www.ia.unc.edu/public/styner/docs/StynerTR97.pdf
>
> If you have the BrainWeb datasets,
> (that you can download from
>
> http://public.kitware.com/pub/itk/Data/BrainWeb/
>
> you can run the bias field estimator with the following command:
>
> ./BiasFieldEstimator --input ~/data/BrainWeb/brainweb165a10f17.mha
> --class-mean 80 130 --class-sigma 10 10 --use-log yes --input-mask
> brainMask.mha --degree 1 --grow 1.01 --shrink 0.95 --max-iteration
> 20000 --init-step-size 0.1 --schedule 2 2 2
>
> The brainMask.mha image can be generated by thresholding the input image
>
> brainweb165a10f17.mha with a threshold of 80, in order to get a binary
> mask
> that roughly covers the brain.
>
>
> This will produce a list of coefficients as output,
> for example, I get the following output:
>
>
> Loading images...
> Images loaded.
> Setting multires schedule :1
> Estimating the bias field...
> --input /home/ibanez/data/BrainWeb/brainweb165a10f17.mha --input-mask
> brainMask.mha --class-mean 80 130 --class-sigma 10 10 --use-log yes
> --degree 1 --grow 1.01 --shrink 0.95 --volume-max-iteration 20
> --inter-slice-max-iteration 4000 --init-step-size 0.1 --coefficients
> 0.00383953 0.000627722 0.00206439 0.00339799
>
>
> You can then take the list of coefficients and use them as input
> to the ImageCorrector tool, with a command like the following:
>
> ./ImageCorrector --input
> /home/ibanez/data/BrainWeb/brainweb165a10f17.mha --output
> correctedImage.mha --degree 1 --use-log yes --coefficients 0.00383953
> 0.000627722 0.00206439 0.00339799
>
>
> which will give you a corrected image.
>
>
> Please not that here, for simplicity I'm using only a polynomial of
> degree 1
> (e.g. a linear function.), In practice you may want to use degree = 3.
>
>
> Regards,
>
>
> Luis
>
>
> --------------------------------------------------------------------------------------------------------
>
>
> On Thu, Dec 25, 2008 at 3:31 AM, orientation f
> <fuzengliang435 at gmail.com <mailto:fuzengliang435 at gmail.com>> wrote:
>
> Hi,
> I want to use the class MRIBiasCorrection for knee MR image
> inhomogeneity correction. But the class includes so many methods
> which confuse me about how to use it for correction. So, is there
> any instruction( as the ITK Software Guide for some class) for
> this class usage?
> Thanks a lot
>
>
>
> ------------------------------------------------------------------------
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ITK FAQ at: http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>
--
Oliver Trebbe
Department of Neurology
University of Muenster
48129 Muenster, Germany
More information about the Insight-users
mailing list