A nonparametric, entropy-minimizing MRI tissue classification algorithm implementation using ITK

Please use this identifier to cite or link to this publication: http://hdl.handle.net/1926/29
This paper focuses on the role of open-source software in the development of a novel magnetic resonance image (MRI) tissue classification algorithm. Specifically, we describe the of use existing classes in the Insight Segmentation and Registration Toolkit (ITK) and several new classes that were implemented to perform non-parametric density estimation and entropy minimization. These new classes also provide a general framework for nonparametric density estimation and related applications.
minus 1 File (114Kb)
There is no code review at this time.

minus Using ITK can clean up software project design, but Science Thrown Open still has its scientific needs. by Greg Harris on 09-09-2005 for revision #1
starstarstarstarstar expertise: 3 sensitivity: 5
The authors propose a (uniform?) map-distortion approach to getting crisp anatomy diagrams for MRI images of an organ whose positional variability with respect to tissue class information is aparently unknown to them. The basal ganglia and perhaps the cerebellum could well be solved this way, but not the cerebral cortex.

The authors propose the (novel?) hypothesis that tissue segmentation should be contingent on the deformability of a tissue-type probability vector atlas onto the acquisition structural image set; that bias field correction should be iterative and anatomy-based; that a stochastic approach to what sampling of features to use will lead to a convergent labeling of the image set with stable posterior class likelihood images; and that the MRI signal from adjacent voxels constitutes the right pool from which to draw feature samples. In this, the authors' interest in texture-driven "surface reconstruction" reported elsewhere appears to be what they have in mind (conjecture).

They claim to be evaluating tissue classifications with a "synthetic data ground truth" phantom image generation process, but didn't validate the appropriateness of their phantom, either. We didn't even find out if the phantom was based on an MRI scan that classified well.

Open Science:
There were no Open Science software submissions with this work-in-progress report.

The only way to reproduce the work corresponding to this report is to steal their idea and crank out code of our own. However, I already know that prior probability field deformation and bias field estimation are inadequate to solve the problem of classifying the human cerebral cortex volume as to gyri and sulci and as to gray matter wrapped around positionally variable white matter spines.

Use of Open Source Software:
The authors dutifully cited their use of ITK, but seemed to miss the point of what someone with experience in brain imaging and already using ITK would want in the way of significant details: knowing that pixel features are stored in an Image class instantiated with a PixelType of [type] Vector is hardly news; instead (in addition?) there should be a discussion of how the kernel size varies from early iterations to late ones: does it go from gross to fine on a schedule? What kind of signal processing filter does this in effect establish, and how does it compare to, say, an edge-preserving filter like itkGradientAnisotropicDiffusionImageFilter?

Open Source Contributions:
They provided a flow chart, and said such a scheme is open to user experimentation, but their code has not been released into ITK yet.

Code Quality:
We are left to expect they are planning to meet the ITK consortium standards for things added to ITK.

Applicability to other problems:
The Parzen method, sorting out worthwhile signal from characterizable error, is general, and suggests considerations of interest to anyone coping with numerically founded map or diagram fitting.

Suggestions for future work:
Statistical texture is a part of choosing the training classes in the stereotaxic cluster Mahalanobis metric classifier I contributed[1] to brains2, but we found that it was important to have the same prior probabilities for classes everywhere throughout the image. We currently have split bias field correction out as a pre-processing stage (that uses ITK) for our T1-T2 tissue classification workup.

Requests for additional information from authors:
The authors propose an implausible kind of measurement feedback based on stability of its outcome. Someone needs to make the case somewhere that MRI tissue classification is a discrete-outcome maximum likelihood kind of probability estimate, rather than a voxel composition fractions kind of probability estimate as histology would suggest.

Additional Comments:
Brain tissue classification is an important problem, but one that must be distinguished from the more anatomically positional problem of parcellation into the functioning units of the cerebral cortex and its supporting nuclei. An atlas has no business decreeing whether you may see gray matter at each of the anatomical coordinates you have imaged.

[1] Harris, G., Andreasen, N.C.,Cizadlo, T., Bailey, J.M., Bockholt, H.J., Magnotta, V.A., and Arndt, S. 1999. Improving Tissue Classification in MRI: A Three-Dimensional Multi-Spectral Discriminant Analysis Method With Automated Training Class Selection. Journal of Computer-Assisted Tomography, 23(1), 144-154.
minus Probably good software, but no evidence provided by Ingmar Bitter on 09-08-2005 for revision #1
starstarstarstarstar expertise: 3 sensitivity: 5

The paper presents a method to classify voxels in MRI brain images into brain tissue classes. It uses probabilistic brain atlases and non-parametric density estimation and entropy minimization to determine convergence of the classification.

That the method can classify brain tissues better than other state of the art classification methods.

There is no evidence for any of the claims. The reader is referred to a currently unpublished paper for the details on mathematics and experiments and comparisons.

Open Science:
The paper states that the code is based on and extents ITK. All new components are (or will be) contributed to ITK. However, the implementation is described as work in progress.
No input data or result images are available.

Reproducing the work is impossible without information on what to reproduce.

Use of Open Source Software:
The described software is open source, based on ITK, and the important classes are well described. However, there are no pointers to the complete application discussed in the paper. It seems that I am supposed to build such an application myself from the existing and new ITK classes.

Open Source Contributions:
I did not evaluate their contribution, nor double checked that the new classes are in the ITK CVS repository.

Code Quality:
I did not evaluate the code.

Applicability to other problems:
The methods appears to be quite general and can easily go beyond brain tissue classification. With few changes it does not even have to be restricted to MRI data.

Suggestions for future work:
Add evidence for the capabilities of your method.

Requests for additional information from authors:

Additional Comments:
The references have many capitalization errors in the paper titles.

The described additions to ITK are most likely very valuable, but the paper does not provide any evidence for its performance claims and it therefore not an example of open science. It does address a need in the community (if it really works) and there is some background material and probably other applications that can benefit from this work.

I'll be happy to change my review if the paper adds the missing evidence and makes it easy to reproduce the results.

Comment by Mahira rai Rai: Mahira Rai yellow
Add a new review
Quick Comments

Download All

Statistics more
Global rating: starstarstarstarstar
Review rating: starstarstarstarstar [review]
Code rating:
Paper Quality: plus minus

Information more
Keywords: MRI tissue classification, Parzen window, ITK, MRI segmentation
Export citation:


Linked Publications more
GeoInterp: Contour Interpolation with Geodesic Snakes
Geodesic Snakes GeoInterp: Contour Interpolation with Geodesic Snakes Geodesic Snakes
by Saboo R., Julian R., Stephen P.
Applying Image Mining to Support Gestational Age Determination Applying Image Mining to Support Gestational Age Determination
by Araujo A., Bellon O., Silva L., Vieira E., Cat M.

View license
Loading license...

Send a message to the author
ISSN 2327-770X
Powered by Midas