Difference between revisions of "Projects:ShapeAnalysisWithOvercompleteWavelets"

From NAMIC Wiki
Jump to: navigation, search
 
(66 intermediate revisions by 4 users not shown)
Line 1: Line 1:
Back to [[NA-MIC_Collaborations|NA-MIC_Collaborations]], [[Algorithm:MIT|MIT Algorithms]]
+
Back to [[Algorithm:MIT|MIT Algorithms]]
 +
__NOTOC__
 +
= Shape Analysis with Overcomplete Wavelets =
  
In this work, we extend the Euclidean wavelets on the sphere. The resulting over-complete spherical wavelets are invariant to rotation of the parameterization of the original spherical image. We apply the over-complete spherical wavelet to cortical folding development and show significantly consistent results as well as improved sensitivity compared with the previous method of using bi-orthogonal spherical wavelet. In particular, we are able to detect developmental asymmetry in the left and right hemispheres.
+
In this work, we extend the Euclidean wavelets to the sphere. The resulting over-complete spherical wavelets are invariant to the rotation of the spherical image parameterization. We apply the over-complete spherical wavelet to cortical folding development and show significantly consistent results as well as improved sensitivity compared with the previously used bi-orthogonal spherical wavelet. In particular, we are able to detect developmental asymmetry in the left and right hemispheres.
  
 
= Description =
 
= Description =
Line 7: Line 9:
 
Bi-orthogonal spherical wavelets have been shown to be
 
Bi-orthogonal spherical wavelets have been shown to be
 
powerful tools in the segmentation and shape analysis of 2D
 
powerful tools in the segmentation and shape analysis of 2D
closed surfaces [1,2], but unfortunately they suffer from aliasing
+
closed surfaces, but unfortunately they suffer from aliasing
 
problems and are therefore not invariant under rotations of
 
problems and are therefore not invariant under rotations of
 
the underlying surface parameterization. See the toy example in the figure below.
 
the underlying surface parameterization. See the toy example in the figure below.
Line 14: Line 16:
 
Instead, we propose to use the over-complete spherical wavelets. These over-complete spherical wavelets are based on filter bank theory, directly extending the ideas of Euclidean steerable pyramid to the sphere. We demonstrate the theoretical advantage of over-complete
 
Instead, we propose to use the over-complete spherical wavelets. These over-complete spherical wavelets are based on filter bank theory, directly extending the ideas of Euclidean steerable pyramid to the sphere. We demonstrate the theoretical advantage of over-complete
 
wavelets over bi-orthogonal wavelets. We also show that over-complete spherical wavelets allow us to
 
wavelets over bi-orthogonal wavelets. We also show that over-complete spherical wavelets allow us to
build more stable cortical folding development models, and detect a wider array of regions of folding development in a newborn dataset.
+
build more stable cortical folding development models, and detect a wider array of regions of folding development in a newborn dataset. The use of spherical wavelet transform in cortical
 +
shape analysis allows us to study cortical folds of different
 +
spatial scales, which are difficult to analyze by cortical
 +
folding analysis methods based on local features such as
 +
curvature and sulcal depth measurements.
  
 +
== Experimental Results ==
 +
 +
'''Comparison with Bi-orthogonal Wavelets'''
 +
 +
The two images on the left of the figure below show the cortical folding speed detected by the bi-orthogonal wavelets. Notice, how the detection changes wildly with different parameterizations. On the other hand, the results of over-complete wavelets are stable across different parameterizations. Furthermore, the over-complete wavelets are much more sensitive than the bi-orthogonal wavelets in detecting changes.
 +
 +
<center>
 
<table>
 
<table>
<tr height="200">
+
<tr>
 
<td>
 
<td>
[[Image:RobustnessOfWaveletAnalysis1.png|thumb|center|height="80"|width="80"|Bi-orthogonal wavelets: Original Parameterization]]
+
[[Image:RobustnessOfWaveletAnalysis1.png|thumb|center|150px|Bi-orthogonal wavelets: Original Parameterization]]
 
</td>
 
</td>
 
<td>
 
<td>
[[Image:RobustnessOfWaveletAnalysis2.png|thumb|center|height="80"|width="80"|Bi-orthogonal wavelets: Different Parameterization]]
+
[[Image:RobustnessOfWaveletAnalysis2.png|thumb|center|150px|Bi-orthogonal wavelets: Different Parameterization]]
 
</td>
 
</td>
 
<td>
 
<td>
[[Image:RobustnessOfWaveletAnalysis3.png|thumb|center|height="80"|width="80"|Over-complete wavelets: Original Parameterization]]
+
[[Image:RobustnessOfWaveletAnalysis3.png|thumb|center|154px|Over-complete wavelets: Original Parameterization]]
 
</td>
 
</td>
 
<td>
 
<td>
[[Image:RobustnessOfWaveletAnalysis4.png|thumb|center|height="80"|width="80"|Over-complete wavelets: Different Parameterization]]
+
[[Image:RobustnessOfWaveletAnalysis4.png|thumb|center|154px|Over-complete wavelets: Different Parameterization]]
 
</td>
 
</td>
 
</tr>
 
</tr>
 
</table>
 
</table>
 +
</center>
  
As seen in the figures above, we see that the bi-orthogonal wavelets analysis of cortical folding is
+
'''Overcomplete Wavelets Global Shape Analysis'''
  
 +
As seen in the figure below, we find that larger folds (lower wavelet scales) develop earlier but slower. This is consistent with the previous postmortem study[3]. Furthermore, we find that the fastest folding development occurs at a younger age on the left hemisphere than on the right.
  
[1] P. Yu, P. E. Grant, Y. Qi, X. Han, et al., "Cortical surface shape
+
<center>
analysis based on spherical wavelets," IEEE Transaction on Medical
+
<table>
Imaging, vol. 26, pp. 582-97, 2007.
+
<tr>
 +
<td>
 +
[[Image:speed_bar_plot.png|thumb|center|222px|Maximum growth rate (1/week) for each wavelet scale]]
 +
</td>
 +
<td>
 +
[[Image:age_bar_plot.png|thumb|center|222px|Age (weeks) of fastest growth for each wavelet scale]]
 +
</td>
 +
</tr>
 +
</table>
 +
</center>
  
[2] D. Nain, S. Haker, A. Bobick, and A. Tannenbaum, "Multiscale 3-d
+
'''Overcomplete Wavelets Local Shape Analysis'''
shape representation and segmentation using spherical wavelets," IEEE
 
Transaction on Medical Imaging, vol. 26, pp. 598-618, 2007.
 
  
= Key Investigator =
+
The figure below shows the results of the regional analysis. Consistent with the global development result discussed above, regions that develop earlier (darker blue) also grow more slowly (more red). For example, the lateral side of the parietal lobe on the left hemisphere develops earlier than the right hemisphere, but at a slower speed.
  
MIT: B.T. Thomas Yeo, Peng Yu, Wanmei Ou, Ellent Grant, Bruce Fischl, Polina Golland
+
Also consistent with the global analysis, we find that larger folds develop earlier but slower. On the lateral side, the pre- and post-central gyri develop the fastest during 30-31 weeks on both hemispheres while smaller structures such as supramarginal and angular gyri develop the fastest at a much later time, as indicated at frequency level 3.
  
= Publication =
+
<center>
 
+
<table>
B.T.T. Yeo, W. Ou, P. Golland. On the Construction of Invertible Filter Banks on the 2-Sphere. Yeo, Ou and Golland. Accepted to the IEEE Transactions on Image Processing [In Press]
+
<tr>
 +
<td>
 +
Maximum Growth Rate (Left)
 +
</td>
 +
<td>
 +
Maximum Growth Rate (Right)
 +
</td>
 +
<td>
 +
Age of Maximum Growth (Left)
 +
</td>
 +
<td>
 +
Age of Maximum Growth (Right)
 +
</td>
 +
</tr>
 +
<tr>
 +
<td>
 +
[[Image:Lh.smoothwm.csw.lsq.rate.level0_lat.png|center|200px|Bi-orthogonal wavelets: Original Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Rh.smoothwm.csw.lsq.rate.level0_lat.png|center|200px|Bi-orthogonal wavelets: Different Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Lh.smoothwm.csw.lsq.age.level0_lat.png|center|200px|Over-complete wavelets: Original Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Rh.smoothwm.csw.lsq.rate.level0_lat.png|center|200px|Over-complete wavelets: Different Parameterization]]
 +
</td>
 +
</tr>
 +
<tr>
 +
<td>
 +
[[Image:Lh.smoothwm.csw.lsq.rate.level1_lat.png|center|200px|Bi-orthogonal wavelets: Original Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Rh.smoothwm.csw.lsq.rate.level1_lat.png|center|200px|Bi-orthogonal wavelets: Different Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Lh.smoothwm.csw.lsq.age.level1_lat.png|center|200px|Over-complete wavelets: Original Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Rh.smoothwm.csw.lsq.rate.level1_lat.png|center|200px|Over-complete wavelets: Different Parameterization]]
 +
</td>
 +
</tr>
 +
<tr>
 +
<td>
 +
[[Image:Lh.smoothwm.csw.lsq.rate.level2_lat.png|center|200px|Bi-orthogonal wavelets: Original Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Rh.smoothwm.csw.lsq.rate.level2_lat.png|center|200px|Bi-orthogonal wavelets: Different Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Lh.smoothwm.csw.lsq.age.level2_lat.png|center|200px|Over-complete wavelets: Original Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Rh.smoothwm.csw.lsq.rate.level2_lat.png|center|200px|Over-complete wavelets: Different Parameterization]]
 +
</td>
 +
</tr>
 +
<tr>
 +
<td>
 +
[[Image:Lh.smoothwm.csw.lsq.rate.level3_lat.png|center|200px|Bi-orthogonal wavelets: Original Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Rh.smoothwm.csw.lsq.rate.level3_lat.png|center|200px|Bi-orthogonal wavelets: Different Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Lh.smoothwm.csw.lsq.age.level3_lat.png|center|200px|Over-complete wavelets: Original Parameterization]]
 +
</td>
 +
<td>
 +
[[Image:Rh.smoothwm.csw.lsq.rate.level3_lat.png|center|200px|Over-complete wavelets: Different Parameterization]]
 +
</td>
 +
</tr>
 +
</table>
 +
</center>
  
P. Yu, B.T.T. Yeo, P.E. Grant, B. Fischl, P. Golland. Cortical Folding Development Study based on Over-Complete Spherical Wavelets.  In Proceedings of MMBIA: IEEE Computer Society Workshop on Mathematical Methods in Biomedical Image Analysis, 2007.
+
= Key Investigators =
  
 +
* MIT: [[http://people.csail.mit.edu/ythomas/ | B.T. Thomas Yeo]], Peng Yu, Wanmei Ou, Polina Golland.
 +
* Harvard: Ellent Grant, Bruce Fischl.
  
 +
= Publications =
  
[http://www.na-mic.org/Special:Publications?text=Projects%3ARegistrationRegularization&submit=Search&words=all&title=checked&keywords=checked&authors=checked&abstract=checked&sponsors=checked&searchbytag=checked| NA-MIC Publications Database]
+
[http://www.na-mic.org/publications/pages/display?search=Projects%3AShapeAnalysisWithOvercompleteWavelets&submit=Search&words=all&title=checked&keywords=checked&authors=checked&abstract=checked&sponsors=checked&searchbytag=checked| NA-MIC Publications Database on Shape Analysis With Overcomplete Wavelets]
  
= Links =
+
[[Category: Shape Analysis]]

Latest revision as of 20:18, 11 May 2010

Home < Projects:ShapeAnalysisWithOvercompleteWavelets
Back to MIT Algorithms

Shape Analysis with Overcomplete Wavelets

In this work, we extend the Euclidean wavelets to the sphere. The resulting over-complete spherical wavelets are invariant to the rotation of the spherical image parameterization. We apply the over-complete spherical wavelet to cortical folding development and show significantly consistent results as well as improved sensitivity compared with the previously used bi-orthogonal spherical wavelet. In particular, we are able to detect developmental asymmetry in the left and right hemispheres.

Description

Bi-orthogonal spherical wavelets have been shown to be powerful tools in the segmentation and shape analysis of 2D closed surfaces, but unfortunately they suffer from aliasing problems and are therefore not invariant under rotations of the underlying surface parameterization. See the toy example in the figure below.

Bump on the sphere (left side). When the north pole is right under the bump, both the bi-orthogonal and overcomplete wavelets detect the bump (first column on the right). When the north pole is rotated away from the bump, only the overcomplete wavelet detects the bump (second column on the right).

Instead, we propose to use the over-complete spherical wavelets. These over-complete spherical wavelets are based on filter bank theory, directly extending the ideas of Euclidean steerable pyramid to the sphere. We demonstrate the theoretical advantage of over-complete wavelets over bi-orthogonal wavelets. We also show that over-complete spherical wavelets allow us to build more stable cortical folding development models, and detect a wider array of regions of folding development in a newborn dataset. The use of spherical wavelet transform in cortical shape analysis allows us to study cortical folds of different spatial scales, which are difficult to analyze by cortical folding analysis methods based on local features such as curvature and sulcal depth measurements.

Experimental Results

Comparison with Bi-orthogonal Wavelets

The two images on the left of the figure below show the cortical folding speed detected by the bi-orthogonal wavelets. Notice, how the detection changes wildly with different parameterizations. On the other hand, the results of over-complete wavelets are stable across different parameterizations. Furthermore, the over-complete wavelets are much more sensitive than the bi-orthogonal wavelets in detecting changes.

Bi-orthogonal wavelets: Original Parameterization
Bi-orthogonal wavelets: Different Parameterization
Over-complete wavelets: Original Parameterization
Over-complete wavelets: Different Parameterization

Overcomplete Wavelets Global Shape Analysis

As seen in the figure below, we find that larger folds (lower wavelet scales) develop earlier but slower. This is consistent with the previous postmortem study[3]. Furthermore, we find that the fastest folding development occurs at a younger age on the left hemisphere than on the right.

Maximum growth rate (1/week) for each wavelet scale
Age (weeks) of fastest growth for each wavelet scale

Overcomplete Wavelets Local Shape Analysis

The figure below shows the results of the regional analysis. Consistent with the global development result discussed above, regions that develop earlier (darker blue) also grow more slowly (more red). For example, the lateral side of the parietal lobe on the left hemisphere develops earlier than the right hemisphere, but at a slower speed.

Also consistent with the global analysis, we find that larger folds develop earlier but slower. On the lateral side, the pre- and post-central gyri develop the fastest during 30-31 weeks on both hemispheres while smaller structures such as supramarginal and angular gyri develop the fastest at a much later time, as indicated at frequency level 3.

Maximum Growth Rate (Left)

Maximum Growth Rate (Right)

Age of Maximum Growth (Left)

Age of Maximum Growth (Right)

Bi-orthogonal wavelets: Original Parameterization
Bi-orthogonal wavelets: Different Parameterization
Over-complete wavelets: Original Parameterization
Over-complete wavelets: Different Parameterization
Bi-orthogonal wavelets: Original Parameterization
Bi-orthogonal wavelets: Different Parameterization
Over-complete wavelets: Original Parameterization
Over-complete wavelets: Different Parameterization
Bi-orthogonal wavelets: Original Parameterization
Bi-orthogonal wavelets: Different Parameterization
Over-complete wavelets: Original Parameterization
Over-complete wavelets: Different Parameterization
Bi-orthogonal wavelets: Original Parameterization
Bi-orthogonal wavelets: Different Parameterization
Over-complete wavelets: Original Parameterization
Over-complete wavelets: Different Parameterization

Key Investigators

  • MIT: [| B.T. Thomas Yeo], Peng Yu, Wanmei Ou, Polina Golland.
  • Harvard: Ellent Grant, Bruce Fischl.

Publications

NA-MIC Publications Database on Shape Analysis With Overcomplete Wavelets