Difference between revisions of "Projects:ConformalFlatteningRegistration"
From NAMIC Wiki
| Line 5: | Line 5: | ||
'''Objective''' | '''Objective''' | ||
| − | The goal of this project is for better visualizing and computation of neural activity from fMRI brain imagery. Also, with this technique, shapes can be mapped to | + | The goal of this project is for better visualizing and computation of neural activity from fMRI brain imagery. Also, with this technique, shapes can be mapped to spheres for shape analysis, registration or other purposes. Our technique is based on conformal mappings which map genus-zero surface: in fMRI case cortical or other surfaces, onto a sphere in an angle preserving manner. |
The explicit transform is obtained by solving a partial differential equation. Such transform will map the original surface to a plane(flattening) and then one can use classic stereographic transformation to map the plane to a sphere. | The explicit transform is obtained by solving a partial differential equation. Such transform will map the original surface to a plane(flattening) and then one can use classic stereographic transformation to map the plane to a sphere. | ||
| − | The process of the algorithm is | + | The process of the algorithm is briefly given below: |
# The conformal mapping <span class="texhtml">''f''</span> is defined on the originla surface <span class="texhtml">Σ</span> as <math>\triangle f = (\frac{\partial}{\partial u} - i\frac{\partial}{\partial v})\delta_p</math>. In that <span class="texhtml">''u''</span> and <span class="texhtml">''v''</span> are the conformal coordinates defined on the surface and the <span class="texhtml">δ<sub>''p''</sub></span> is a Dirac function whose value is non-zero only at point <span class="texhtml">''p''</span>. By solving this partial differential equation the mapping <span class="texhtml">''f''</span> can be obtained. | # The conformal mapping <span class="texhtml">''f''</span> is defined on the originla surface <span class="texhtml">Σ</span> as <math>\triangle f = (\frac{\partial}{\partial u} - i\frac{\partial}{\partial v})\delta_p</math>. In that <span class="texhtml">''u''</span> and <span class="texhtml">''v''</span> are the conformal coordinates defined on the surface and the <span class="texhtml">δ<sub>''p''</sub></span> is a Dirac function whose value is non-zero only at point <span class="texhtml">''p''</span>. By solving this partial differential equation the mapping <span class="texhtml">''f''</span> can be obtained. | ||
| Line 25: | Line 25: | ||
[[Image:Brain.PNG|[[Image:Brain.PNG|Image:Brain.PNG]]]] ===> [[Image:Brain-flat.PNG|[[Image:Brain-flat.PNG|Image:Brain-flat.PNG]]]] | [[Image:Brain.PNG|[[Image:Brain.PNG|Image:Brain.PNG]]]] ===> [[Image:Brain-flat.PNG|[[Image:Brain-flat.PNG|Image:Brain-flat.PNG]]]] | ||
| − | |||
''References'' | ''References'' | ||
| − | + | * S. Angenent, S. Haker, A. Tannenbaum, and R. Kikinis. On the Laplace-Beltrami operator and brain surface flattening. IEEE Trans. on Medical Imaging, Vol. 18, pp. 700-711, 1999 <br /> [2] Y. Gao, J. Melonakos, and A. Tannenbaum. Conformal Falttening ITK Filter. to appear in In 2006 MICCAI Open-source Workshop <br /> | |
'''Links''' | '''Links''' | ||
| − | |||
| − | |||
| − | |||
Revision as of 13:21, 4 September 2007
Home < Projects:ConformalFlatteningRegistrationConformal Flattening Registration
Back to NA-MIC_Collaborations
Objective
The goal of this project is for better visualizing and computation of neural activity from fMRI brain imagery. Also, with this technique, shapes can be mapped to spheres for shape analysis, registration or other purposes. Our technique is based on conformal mappings which map genus-zero surface: in fMRI case cortical or other surfaces, onto a sphere in an angle preserving manner.
The explicit transform is obtained by solving a partial differential equation. Such transform will map the original surface to a plane(flattening) and then one can use classic stereographic transformation to map the plane to a sphere.
The process of the algorithm is briefly given below:
- The conformal mapping f is defined on the originla surface Σ as [math]\triangle f = (\frac{\partial}{\partial u} - i\frac{\partial}{\partial v})\delta_p[/math]. In that u and v are the conformal coordinates defined on the surface and the δp is a Dirac function whose value is non-zero only at point p. By solving this partial differential equation the mapping f can be obtained.
- To solve that equation on the discrete mesh representation of the surface, finite element method(FEM) is used. The problem is turned to solving a linear system D'x = b. Since b is complex vector, the real and imaginary parts of the mapping f can be calculated separately by two linear system.
- Having the mapping f, the original surface can be mapped to a plane.
- Further, the plane can be mapped to a sphere by the stereographic projection.
Also, in the work of Multiscale Shape Segmentation, conformal flattening is used as the first step for remeshing the surface.
Progress
This algorithm is now written into SandBox as itkConformalFlatteningFilter. More test is being made and then hopefully can be integrated into the ITK CVS repository.
===> 
===> 
References
- S. Angenent, S. Haker, A. Tannenbaum, and R. Kikinis. On the Laplace-Beltrami operator and brain surface flattening. IEEE Trans. on Medical Imaging, Vol. 18, pp. 700-711, 1999
[2] Y. Gao, J. Melonakos, and A. Tannenbaum. Conformal Falttening ITK Filter. to appear in In 2006 MICCAI Open-source Workshop
Links
