High Speed Electric and Magnetic Field calculation Algorithm for fMRI-MEG Integrative Analysis of Neural Activity
Project/Area Number |
17560371
|
Research Category |
Grant-in-Aid for Scientific Research (C)
|
Allocation Type | Single-year Grants |
Section | 一般 |
Research Field |
Measurement engineering
|
Research Institution | Kyoto University |
Principal Investigator |
HAMADA Shoji Kyoto University, Graduate school of Engineering, Associate Professor, 工学研究科, 助教授 (20246656)
|
Co-Investigator(Kenkyū-buntansha) |
KOBAYASHI Tetsuo Kyoto University, Graduate school of Engineering, Professor, 工学研究科, 教授 (40175336)
|
Project Period (FY) |
2005 – 2006
|
Project Status |
Completed (Fiscal Year 2006)
|
Budget Amount *help |
¥3,500,000 (Direct Cost: ¥3,500,000)
Fiscal Year 2006: ¥600,000 (Direct Cost: ¥600,000)
Fiscal Year 2005: ¥2,900,000 (Direct Cost: ¥2,900,000)
|
Keywords | Electric and Magnetic Field Analysis / fMRI / MEG / Surface Charge Simulation Method / Fast Multipole Method / Voxel model / Integrative Analysis |
Research Abstract |
1.We have developed an 'equivalent multipole moment method with polar translations (EMMM-PT)'. This method enables us to numerically calculate exact solution described by Laplace's equation in multi spheres model. We also have developed a 'fast-multipole surface-charge-simulation method for voxel models (FMM-SCM-V)'. This method adopts the FMM that is well known as an efficient 0(N) solver, and it enables us to treat high-resolution human model composed of several to hundred million voxels based on MR images. 2.Numerical calculation procedures for magnitude of MEG coil signal are classified. One way traces 'Equivalent current dipole→Potential on boundary→Magnetic flux by Geselowitz formula', the other traces 'AC coil current→Induced electric field-Magnetic flux by Eaton's reciprocity theorem'. Analyses by EMMM-PT numerically proved the identification of these two ways in multi-media tissue models. 3.The accuracy of FMM-SCM-V is estimated by comparing the numerical results by FMM-SCM-V and EMMM-PT in the grey matter region of a seven spheres head model. The standard deviation of error vector of the electric field is about 0.81% (normalized by the maximum field in the grey matter region) in the case of 0.34mm voxel size. In this case, the number of voxels, the number of surface elements, and calculation time are 106,575,801 voxels, 4,936,874 elements, 59 min. 39 sec., respectively. 4.fMRI-MEG Integrative method is performed by utilizing 'an anomalously applied least squares' and 'linearly-constrained adaptive beamformer' based on the lead field matrix calculated by the Sarvas equation. It has been found that a non-adaptive version of the latter is equivalent to the former, and the both methods have advantages and disadvantages in the spatial and temporal resolution properties.
|
Report
(3 results)
Research Products
(25 results)