B1 simulator : simulate B1 field for MRI RF coils |
|||
These Matlab routines include demonstration and individual Matlab functions to simulate B1 field generated by the RF coil (including array). Given the 3D spatial disbribution of current on the array coil, Biot-Savart's law is used to calculate the magnetic field associated with the current distribution on the RF coil. A newer version using quasi-static Maxwell equations to model the magnetic field patterns is coming soon. |
|||
| Codes | |||
| Inside Martinos Center: | |||
| If you are inside Martinos Center in MGH, you can add the following path to your starup.m file to access all the routines. | |||
addpath('~fhlin/matlab/toolbox/fhlin_toolbox'); |
|||
| If you are outside Martinos Center in MGH, please download the codes to your local computer. | |||
| b1sim_dc_core.m | calculate the quasistatic B1 field given the coil geometry and FOV using Biot-Savart's law | ||
| Example | |||
| 8-channel head array coil: | |||
| Date: | Jan-10-2005 | ||
Purpose: |
Simulate and display the B1 field for an 8-channel array coil. |
||
| Description:: | The spatial distribution of the current on an 8-element RF array coil is simulated. Each element of the array is a circular surface coil. Using Biot-Savart's law, the program calculates the 3D spatial distribution of the magnetic field. | ||
| Code: | Matlab code: b1sim_dc_array_circle8.m Matlab code: show_b1sim_dc_array_circle8.m Etc: array_cicle8.txt (definition of the coil geometry) Etc: outer_skull642.tri (data file of head outer skull 3D mesh) Etc: inverse_read_tri.m (the file i/o routine to read triangulation 3D mesh) |
||
| Example: | b1_sim_dc_array_circle8; %calculate the B1 field show_b1sim_dc_array_circle8; %display the array coil geometry |
||
|
| The coil geometry of the 8-channel head array. All eight element of the array coil wrap around the head. The center of the RF array coil shows the outer skull 3D model. One calculated B1 field (the combinatio of x- and y- directional components) from the combination of all 8 channels is shown at right. | |||
 |
| The combination of x- and y-components of the B1 field calculated for each coil of the array on an axial slice. | |||