JEMRIS
2.8.1
opensource MRI simulations

In the following, a 2D spatiallyselective excitation experiment [Pauly 1989] will be demontrated in detail in three different ways. The first will reflect the flexibility of the symbolic mathematics evaluation capabilities of JEMRIS. In the second part, the use and implementation of external gradient and RF pulses will be demonstrated, and in a third part, the framework will be extended to boost the performance of the simulation process.
File
menu select Open
sequence. MyEPI.xml
. PaulyEtAl.xml
The paper describes the spatiallyselective small tip angle excitation with the application of two gradients accompanied by an RF pulse.
To start off, let us delete the hard rf pulse and replace it with a 2D selective excitation block consisting of an analytic RF pulse, and two analytic gradient pulses in read and phase direction.
HRF
. Exc
. ERF
, for excitation RF pulse. Exc
by clicking on the corresponding modul symbol and afterwards on Exc
. EGX
for analytic gradient on the X axis. EGY
for analytic gradient on the Y axis. Axis
of the latter to GY
.Let us now have a look at the analytical solution for the rf pulse and the kspace trajectory presented by the author, which we want to implement:
The paper investigates the case where beta is 2, n is 8, and A is 1 while the excitation is played out for 2 ms. The following steps will associate the excitation block to the according formulas.
ERF
. Duration
to 2. Observe
to ERF
,Duration.What the last step essentially means, is that some attribute of RF pulse needs to know about its own duration. You can see, the magnitude of the RF pulse is time dependent and varies over the 2 ms pulse duration. The concept of observing oneselfs attribute may appear confusing on a first glance, but is due to the flexibility of the attribute implmentation.
Let us now type in the constants involved in order of apearance. They can be accessed by depending attributes as c1
through c5
.
1
,.01,1,2,8 for gamma, alpha, A, beta, and n.Here, some explanation is needed as to why gamma is set to 1. Simply put, JEMRIS does not simulate the nuclear physics deeper than to the point, where the relaxation times for proton imaging are given. In that case the field strength dependant part is intrinsically built into the sample. So, in order to simulate the effects associated with faster T2 decay at high fields, one changes the value for T2 accordingly.
Let us now type in the formula of the time dependant rf magnitude.
Shape
to c1*c2*c3/a1*exp
( c4^2 * (1T/a1)^2) * ((2*Pi*c5*(1T/a1))^2 + 1)^0.5. TPOIs
to 200
.You should now be looking at the GUI in the following state:
We are not finished yet. The Gradients are next.
EGX
.We need to introduce the necessary parameters, all stemming from ERF
. As there are the duration, A, and n.
Observe
to ERF
,Duration/ERF,Constant3/ERF,Constant5. Duration
to a1
. Shape
to a2*
(1T/a1)*cos(2*Pi*a3*T/a1). Diff
to 1. EGY
. Duration
to a1
. Shape
to a2*
(1T/a1)*sin(2*Pi*a3*T/a1). Diff
to 1.The not obvious step in the last instruction block is the setting of the value for Diff
. In this step, the flexibility of the symbolic math library is demonstrated. The formula, that was typed in for the Shape
describes the kspace trajectory's time evolution. However, we are interested in the shape of the corresponding gradient. Which is just the first derivative with respect to time or T
. Alternatively you may also reset Diff to 0 and set Shape
to the formula for the gradients as described in the paper:
It is easily comprehendable that the differentiation will cost simulation time.
In Part 2, we will expand the framework by two very easily coded classes to demonstrate the extensibility of the framework. Additionally, as the use of symbolic mathematics is very appealing and should be the first choice for implementing new ideas, there are some performance issues to be considered. If, however, the computation is shifted to the level of the framework in C++, major speed ups may be gained depending on the complexity of the mathematics involved.
You may now simulate the sequence as described in the first steps tutorial. But! be very careful here. The simulation involves some computation power and, depending on the sample size, some RAM.
Step on to: In depth tutorial  Part 2
[Pauly 1989] Pauly JM, Nishimura DG, Macovski A. A kSpace Analysis of SmallTipAngle Excitation. Journal of Magnetic Resonance, 81, 4356 (1989)
 last change 17.06.2016 
Tony Stoecker 
Imprint 