190 likes | 208 Views
Total Variation Based Magnetic Field Mapping. Zhuo Zheng Joint work with C. Anand and T. Terlaky Advanced Optimization Lab. Outline. Motivation. Phase shift in the ideal case. Field mapping methods. Phase shift in the presence of noise. Total variation. Conclusions. Motivation.
E N D
Total Variation Based Magnetic Field Mapping Zhuo Zheng Joint work with C. Anand and T. Terlaky Advanced Optimization Lab
Outline • Motivation Phase shift in the ideal case • Field mapping methods Phase shift in the presence of noise Total variation • Conclusions
Motivation • Magnetic fields determine the resonant frequency • The field inhomogeneity is inevitable • It causes geometric distortions and artifacts in images • SSFP imaging is quite sensitive to the resonant frequency
Correction Methods Steps Involved By acquiring two images at different TEs Use of the field map to compensate in the reconstruction Measure spatial variation Field Compensation Reasonably good Images
Previous Model Revisited • Steady-state model established
Search And Position (SAP) • Perform pixel-wise evaluation to locate a certain values • No noise in the ideal case • It results in a slightly different steady-state matrix S • We introduce a new P in the presence of the field variation
Search And Position (SAP) Modeling in the rotating frame of reference
Search And Position (SAP) Suppose we have a ,We need to know where it locates in the image
Noise Corruption Case • Measurement becomes • A least square program has been formulated Subject to some constraints. Numerical results show that it works !
Total Variation Approach • First proposed by Rudin, Osher and Fatemi • Total variation based noise removal algorithms • Denoising can be solved as a constrained optimization program Objective: TV norm of the denoised image Constraint: deviation from the observation
Total Variation Approach Model Equation: Optimization Problem:
Total Variation Approach However, we are concerned about the signal instead of the image Therefore, our model should be:
Total Variation Approach However, we are concerned about the signal instead of the image Therefore, our model should be: We formulate the nonconstrained problem and derive the optimality condition.
Total Variation Approach Therefore, our model should be: We formulate the nonconstrained problem and derive the optimality condition.
Total Variation Approach Euler-Lagrange Equation: Substitute into the expression of field map term and tissue density, we try to solve a hard optimization program.
Euler-Lagrange Equation: Substitute into the expression of field map term and tissue density, we try to solve a hard optimization program.
Conclusions • We explore several methods to generate a field map • Total variation based method is promising in our case • Better formulation has to be called for