Difference between revisions of "2013 Summer Project Week:Epilepsy Surgery"

From NAMIC Wiki
Jump to: navigation, search
(Created page with '__NOTOC__ <gallery> Image:PW-MIT2013.png|Projects List Image:genuFAp.jpg|Scatter plot of the original FA data through the genu of the corpus…')
 
 
(42 intermediate revisions by 2 users not shown)
Line 2: Line 2:
 
<gallery>
 
<gallery>
 
Image:PW-MIT2013.png|[[2013_Summer_Project_Week#Projects|Projects List]]
 
Image:PW-MIT2013.png|[[2013_Summer_Project_Week#Projects|Projects List]]
Image:genuFAp.jpg|Scatter plot of the original FA data through the genu of the corpus callosum of a normal brain.
 
Image:genuFA.jpg|Regression of FA data; solid line represents the mean and dotted lines the standard deviation.
 
 
</gallery>
 
</gallery>
 
==Instructions for Use of this Template==
 
#Please create a new wiki page with an appropriate title for your project using the convention 2013_Winter_Project_Week:<Project Name>
 
#Copy the entire text of this page into the page created above
 
#Link the created page into the list of projects for the project event
 
#Delete this section from the created page
 
#Send an email to tkapur at bwh.harvard.edu if you are stuck
 
  
 
==Key Investigators==
 
==Key Investigators==
* UNC: Isabelle Corouge, Casey Goodlett, Guido Gerig
+
* [http://www5.usp.br/en/ USP] - Luiz Murta
* Utah: Tom Fletcher, Ross Whitaker
 
  
 
<div style="margin: 20px;">
 
<div style="margin: 20px;">
<div style="width: 27%; float: left; padding-right: 3%;">
+
<div style="width: 47%; float: left; padding-right: 3%;">
  
 
<h3>Objective</h3>
 
<h3>Objective</h3>
We are developing methods for analyzing diffusion tensor data along fiber tracts. The goal is to be able to make statistical group comparisons with fiber tracts as a common reference frame for comparison.
+
This project will investigate the presence and location of the epileptogenic focus in temporal lobe by analyzing patterns of texture in magnetic resonance imaging (MRI) after segmentation using anisotropic diffusion filters anomalous and geodesic active contour.
  
 +
<li>The Temporal Lobe Epilepsy (TLE) with hippocampal sclerosis affects other mesial temporal lobe structures, showing the change in signal intensity of white matter (WM) and gray matter (GM) in Magnetic Resonance Imaging (MRI) in patients with TLE.
 +
<li>Such abnormalities of signal intensity changes in the temporal lobes, or loss of marking between WM and GM are referred as blurring in literature.
 +
<li>TLE is the most common type of refractory epilepsy in adults, and thus strong candidate to surgery.
 +
<br><br>
 +
<b>Purpose</b>
 +
<li>To provide a better quality of life for those with TLE, a surgical procedure is usually proposed for extracting the brain region that is the focus of epilepsy.
 +
<li>This project aims to implment a semi-automated method capable of locating lesions, that is difficult to be detected by visual inspection, is built to help experts identify tissues presenting blurring in temporal lobes.
 +
<li>The proposed software tool was developed in C + +, using three toolkits to support: VTK, ITK and QT.
 +
<li>To provide segmentation of the temporal lobes, we used the Geodesic Active Contour method combined with the anomalous anisotropic diffusion filter.
  
 +
</div>
  
 +
<div style="width: 50%; float: left;">
  
 +
<h3>Progress</h3>
 +
Examples:<br>
 +
[[File:ex_a.png]]<br>
 +
Normal MRI  at mesial temporal lobe<br>
  
 +
[[File:ex_b.png]]<br>
 +
MRI containing blurring phenomena on right side as indicated by the yellow arrow<br>
  
</div>
+
</div> <br>
  
<div style="width: 27%; float: left; padding-right: 3%;">
+
<div style="margin: 20px;">
 +
<div style="width: 60%; float: left; padding-right: 3%;">
  
<h3>Approach, Plan</h3>
+
<h3>Methods</h3>
 +
[[File:param.png]]
 +
[[File:imageUI.png]] <br>
 +
Default adjustable parameters, and original image,<br>
 +
<br>
 +
<b>Segmented TL córtex</b>
 +
[[File:segmB.png]]
 +
[[File:segmR.png]]
  
Our approach for analyzing diffusion tensors is summarized in the IPMI 2007 reference below.  The main challenge to this approach is <foo>.
 
  
Our plan for the project week is to first try out <bar>,...
 
  
 
</div>
 
</div>
  
<div style="width: 40%; float: left;">
+
<div style="width: 37%; float: left;">
  
<h3>Progress</h3>
+
<h3>Results</h3>
Software for the fiber tracking and statistical analysis along the tracts has been implemented. The statistical methods for diffusion tensors are implemented as ITK code as part of the [[NA-MIC/Projects/Diffusion_Image_Analysis/DTI_Software_and_Algorithm_Infrastructure|DTI Software Infrastructure]] project. The methods have been validated on a repeated scan of a healthy individual. This work has been published as a conference paper (MICCAI 2005) and a journal version (MEDIA 2006). Our recent IPMI 2007 paper includes a nonparametric regression method for analyzing data along a fiber tract.
+
<b>Classifier</b>
 +
<li>Once segmented, the images are analyzed using texture descriptors from the intensity histogram and local pixels dependence measures proposed by Haralick were extracted.
 +
Subsequently, the set of feature vectors from the segmented images is evaluated using three classifiers in [http://www.cs.waikato.ac.nz/ml/weka/ WEKA] environment.
 +
<li>Images of 32 patients with TLE were used, 15 of those exhibiting blurring phenomenon in TL and 17 do not.
 +
<li>In total, a set of 98 planar images was segmented, and 51 are classified in "with blurring" category and 47 "without blurring" category.
 +
<br>
 +
<li>Of all three classifiers:
 +
<br>-- artificial neural network;
 +
<br>-- nearest neighbour;
 +
<br>-- and decision tree
 +
<li> tested, the one based on "J48" decision tree had the best, and most interesting results.
 +
<li>32 feature have been used to constructo descriptors vectors of segmented images:
 +
<br>-- 24 were extracted from the co-occurrence matrix using Haralick texture descriptors and
 +
<br>-- 8 were intensity statistics obtained from histogram.  
 +
<li>When evaluating the performance of the three classifiers selected for the learning process, as a basis of 70 training samples and 28 test.
 +
<li>It was found that, for decision tree, precision and area under the ROC curve were 71.42 % and 69.30 %, respectively.
 +
<br>
 +
[[File:tlNoblur.png]]
 +
[[File:histNoblur.png]] no blurring
 +
[[File:tlBlurred.png]]
 +
[[File:histBlurred.png]] blurring
 +
<br>
 +
<br>
 +
<b>Decision Tree</b>
 +
<br><br>
 +
D_Hist_MeanIntensity <= 0.453<br>
 +
|  D_Hist_stdDev <= 0.138: with_blurring (5.0)<br>
 +
|  D_Hist_stdDev > 0.138<br>
 +
|  |  E_Hist_Kurtosis <= 0.497<br>
 +
|  |  |  E_COOmeanHomogeneity_dist3 <= 0.425<br>
 +
|  |  |  |  D_COOmeanHomogeneity_dist2 <= 0.318: without_blurring (9.0/2.0)<br>
 +
|  |  |  |  D_COOmeanHomogeneity_dist2 > 0.318: with_blurring (5.0)<br>
 +
|  |  |  E_COOmeanHomogeneity_dist3 > 0.425<br>
 +
|  |  |  |  E_COOmeanHomogeneity_dist3 <= 0.895: without_blurring (26.0/1.0)<br>
 +
|  |  |  |  E_COOmeanHomogeneity_dist3 > 0.895: with_blurring (3.0/1.0)<br>
 +
|   |  E_Hist_Kurtosis > 0.497: with_blurring (5.0)<br>
 +
D_Hist_MeanIntensity > 0.453<br>
 +
|  E_Hist_Kurtosis <= 0.48: with_blurring (12.0)<br>
 +
|  E_Hist_Kurtosis > 0.48<br>
 +
|  |  E_COOmeanHomogeneity_dist1 <= 0.432: with_blurring (3.0)<br>
 +
|  |  E_COOmeanHomogeneity_dist1 > 0.432: without_blurring (2.0)<br>
 +
 Number of Leaves  : 9<br>
  
 +
</div> <br>
  
</div>
+
<h3>Conclusions</h3>
</div>
 
  
==Delivery Mechanism==
+
<li>As a result, from preselected 32 texture feature descriptors, 6 were shown to be relevant for classification of TL tissues.
 +
<li>Thus, we observe that using approach of texture analysis of segmented images, can helps the specialist to identify signal change MRIs in the temporal lobes of patients with TLE.
 +
<li>Once trained, classifier is easily included in this system, and can help surgery decision.
 +
<li>Although these preliminary results are promising, the method still needs to be validated considering TL separatelly and in a larger image database.
  
This work will be delivered to the NA-MIC Kit as a (please select the appropriate options by noting YES against them below)
+
<h3>To do list:</h3>
  
#ITK Module
+
<li>As a nest step, I´ll integrate the code to 3DSlicer.
#Slicer Module
+
<li>Try new schems of texture analisys to this application.
##Built-in
+
<li>Implement few machine learning algorithms, in addition to image processing in C++.
##Extension -- commandline
+
<li>Speed up some time consumming algorithms with GPU programming.
##Extension -- loadable
+
<li>Extend these and other ideas to cotical dysplasia: 2014 !
#Other (Please specify)
 
  
 
==References==
 
==References==
*Fletcher P, Tao R, Jeong W, Whitaker R. [http://www.na-mic.org/publications/item/view/634 A volumetric approach to quantifying region-to-region white matter connectivity in diffusion tensor MRI.] Inf Process Med Imaging. 2007;20:346-358. PMID: 17633712.
+
*Shaker, M. & Soltanian-Zadeh, H., 2008. Voxel-Based Morphometric Study of Brain Regions from Magnetic Resonance Images in Temporal Lobe Epilepsy. Image Analysis and Interpretation, 2008. SSIAI 2008. IEEE Southwest Symposium on, 209-212.
* Corouge I, Fletcher P, Joshi S, Gouttard S, Gerig G. [http://www.na-mic.org/publications/item/view/292 Fiber tract-oriented statistics for quantitative diffusion tensor MRI analysis.] Med Image Anal. 2006 Oct;10(5):786-98. PMID: 16926104.
 
* Corouge I, Fletcher P, Joshi S, Gilmore J, Gerig G. [http://www.na-mic.org/publications/item/view/1122 Fiber tract-oriented statistics for quantitative diffusion tensor MRI analysis.] Int Conf Med Image Comput Comput Assist Interv. 2005;8(Pt 1):131-9. PMID: 16685838.
 
* Goodlett C, Corouge I, Jomier M, Gerig G, A Quantitative DTI Fiber Tract Analysis Suite, The Insight Journal, vol. ISC/NAMIC/ MICCAI Workshop on Open-Source Software, 2005, Online publication: http://hdl.handle.net/1926/39 .
 

Latest revision as of 15:01, 21 June 2013

Home < 2013 Summer Project Week:Epilepsy Surgery

Key Investigators

  • USP - Luiz Murta

Objective

This project will investigate the presence and location of the epileptogenic focus in temporal lobe by analyzing patterns of texture in magnetic resonance imaging (MRI) after segmentation using anisotropic diffusion filters anomalous and geodesic active contour.

  • The Temporal Lobe Epilepsy (TLE) with hippocampal sclerosis affects other mesial temporal lobe structures, showing the change in signal intensity of white matter (WM) and gray matter (GM) in Magnetic Resonance Imaging (MRI) in patients with TLE.
  • Such abnormalities of signal intensity changes in the temporal lobes, or loss of marking between WM and GM are referred as blurring in literature.
  • TLE is the most common type of refractory epilepsy in adults, and thus strong candidate to surgery.

    Purpose
  • To provide a better quality of life for those with TLE, a surgical procedure is usually proposed for extracting the brain region that is the focus of epilepsy.
  • This project aims to implment a semi-automated method capable of locating lesions, that is difficult to be detected by visual inspection, is built to help experts identify tissues presenting blurring in temporal lobes.
  • The proposed software tool was developed in C + +, using three toolkits to support: VTK, ITK and QT.
  • To provide segmentation of the temporal lobes, we used the Geodesic Active Contour method combined with the anomalous anisotropic diffusion filter.
  • Progress

    Examples:
    Ex a.png
    Normal MRI at mesial temporal lobe

    Ex b.png
    MRI containing blurring phenomena on right side as indicated by the yellow arrow


    Methods

    Param.png ImageUI.png
    Default adjustable parameters, and original image,

    Segmented TL córtex SegmB.png SegmR.png


    Results

    Classifier

  • Once segmented, the images are analyzed using texture descriptors from the intensity histogram and local pixels dependence measures proposed by Haralick were extracted. Subsequently, the set of feature vectors from the segmented images is evaluated using three classifiers in WEKA environment.
  • Images of 32 patients with TLE were used, 15 of those exhibiting blurring phenomenon in TL and 17 do not.
  • In total, a set of 98 planar images was segmented, and 51 are classified in "with blurring" category and 47 "without blurring" category.
  • Of all three classifiers:
    -- artificial neural network;
    -- nearest neighbour;
    -- and decision tree
  • tested, the one based on "J48" decision tree had the best, and most interesting results.
  • 32 feature have been used to constructo descriptors vectors of segmented images:
    -- 24 were extracted from the co-occurrence matrix using Haralick texture descriptors and
    -- 8 were intensity statistics obtained from histogram.
  • When evaluating the performance of the three classifiers selected for the learning process, as a basis of 70 training samples and 28 test.
  • It was found that, for decision tree, precision and area under the ROC curve were 71.42 % and 69.30 %, respectively.
    TlNoblur.png HistNoblur.png no blurring TlBlurred.png HistBlurred.png blurring

    Decision Tree

    D_Hist_MeanIntensity <= 0.453
    | D_Hist_stdDev <= 0.138: with_blurring (5.0)
    | D_Hist_stdDev > 0.138
    | | E_Hist_Kurtosis <= 0.497
    | | | E_COOmeanHomogeneity_dist3 <= 0.425
    | | | | D_COOmeanHomogeneity_dist2 <= 0.318: without_blurring (9.0/2.0)
    | | | | D_COOmeanHomogeneity_dist2 > 0.318: with_blurring (5.0)
    | | | E_COOmeanHomogeneity_dist3 > 0.425
    | | | | E_COOmeanHomogeneity_dist3 <= 0.895: without_blurring (26.0/1.0)
    | | | | E_COOmeanHomogeneity_dist3 > 0.895: with_blurring (3.0/1.0)
    | | E_Hist_Kurtosis > 0.497: with_blurring (5.0)
    D_Hist_MeanIntensity > 0.453
    | E_Hist_Kurtosis <= 0.48: with_blurring (12.0)
    | E_Hist_Kurtosis > 0.48
    | | E_COOmeanHomogeneity_dist1 <= 0.432: with_blurring (3.0)
    | | E_COOmeanHomogeneity_dist1 > 0.432: without_blurring (2.0)
     Number of Leaves : 9

  • Conclusions

  • As a result, from preselected 32 texture feature descriptors, 6 were shown to be relevant for classification of TL tissues.
  • Thus, we observe that using approach of texture analysis of segmented images, can helps the specialist to identify signal change MRIs in the temporal lobes of patients with TLE.
  • Once trained, classifier is easily included in this system, and can help surgery decision.
  • Although these preliminary results are promising, the method still needs to be validated considering TL separatelly and in a larger image database.

    To do list:

  • As a nest step, I´ll integrate the code to 3DSlicer.
  • Try new schems of texture analisys to this application.
  • Implement few machine learning algorithms, in addition to image processing in C++.
  • Speed up some time consumming algorithms with GPU programming.
  • Extend these and other ideas to cotical dysplasia: 2014 !

    References

    • Shaker, M. & Soltanian-Zadeh, H., 2008. Voxel-Based Morphometric Study of Brain Regions from Magnetic Resonance Images in Temporal Lobe Epilepsy. Image Analysis and Interpretation, 2008. SSIAI 2008. IEEE Southwest Symposium on, 209-212.