MRI Contrast Calculator: TR, TE, Flip Angle and Scan Time

Compare spin-echo and spoiled gradient-echo MRI signals, tissue contrast, Ernst angles and scan time. Explore TR and TE curves with a precession and relaxation animation.

#Ernst angle #MRI #scan time #spin echo #spoiled gradient echo #T1 relaxation #T2 relaxation

This MRI contrast calculator helps students, technologists and engineers connect pulse-sequence settings to the signal that reaches the echo. Change repetition time (TR), echo time (TE) and, for spoiled gradient echo, the excitation flip angle. Compare two tissues without needing a research simulation package.

Start with the default spin-echo example: white matter and gray matter at 1.5 T, TR 500 ms and TE 15 ms. The signal cards report relative amplitudes on one common scale. A positive A − B means tissue A is brighter; a negative value means tissue B is brighter. The shaded circles are signal swatches, not reconstructed anatomy. Relative proton density defaults to one for both tissues so that relaxation differences can be studied separately.

Choose either 1.5 T or 3 T to use the corresponding measurements from Stanisz and colleagues. These presets describe excised specimens at 37 °C, including animal tissues; they are teaching examples rather than universal values for living patients. Select custom T1 and T2 to explore a phantom, another publication or a hypothetical tissue. Changing field strength selects relaxation data only: the calculator does not predict the field-dependent change in equilibrium magnetization, coil sensitivity or signal-to-noise ratio.

For T1 weighted spin echo, explore shorter TR with short TE. To emphasize T2 differences, lengthen TR and then increase TE. The TR and TE charts hold the other settings fixed and show both signals. Their sweep ranges keep the spin-echo approximation within TE ≤ TR/10. Contrast depends on the tissue pair; a weighting label alone cannot predict which tissue is brighter.

Switch to spoiled gradient echo to investigate flip angle and T2 decay. Enter T2 independently: the illustrative starting values are not Stanisz measurements and do not update with field strength. The Ernst angle maximizes one tissue's steady-state signal at fixed TR. It generally differs from the angle giving the largest difference between two tissues. The flip-angle chart deliberately uses TE = 0 to isolate excitation and T1 saturation.

Scan time assumes conventional Cartesian sampling with one acquired phase-encoding line per TR. Set phase lines and averages, then either sequential 2D slice packages or acquired 3D partitions. Interleaved slices within a single package do not each multiply the time. Preparation pulses, calibration, triggering, partial Fourier, acceleration and turbo echo trains are outside this estimate.

Press ▶ to watch tissue A precess, lose transverse magnetization and recover longitudinal magnetization after one isolated pulse. Precession is slowed and its displayed phase is illustrative; relaxation follows the chosen time constants. This educational tool explains contrast mechanisms and timing trade-offs. It does not select clinical protocols, assess patient safety, calculate contrast-agent doses or replace a scanner's timing and safety system.

How it is calculated

The conventional short-TE spin-echo approximation is

SSE=ρ(1−e−TR/T1)e−TE/T2.S_{SE}=\rho\left(1-e^{-TR/T_1}\right)e^{-TE/T_2}.

It assumes ideal 90° excitation and 180° refocusing with TE small compared with TR; the change in longitudinal recovery caused by the refocusing pulse is neglected. See Sung, UCLA M219 (2024), slide 25.

For a completely spoiled gradient echo in steady state,

SGRE=ρsin⁡α 1−E11−E1cos⁡α e−TE/T2∗,E1=e−TR/T1.S_{GRE}=\rho\sin\alpha\,\frac{1-E_1}{1-E_1\cos\alpha}\,e^{-TE/T_2^*},\qquad E_1=e^{-TR/T_1}.

The fixed point before excitation follows from Mz−=M0(1−E1)+E1cos⁡α Mz−M_z^-=M_0(1-E_1)+E_1\cos\alpha\,M_z^-. Maximizing the signal with respect to angle gives αE=arccos⁡E1\alpha_E=\arccos E_1. See Hargreaves 2012, Eq. 5 (https://doi.org/10.1002/jmri.23742) and Karakuzu 2022, Eqs. 2.5–2.6 (attached s4, pages 56–58). The Ernst result applies to each tissue separately, not to maximum tissue contrast.

Signed contrast is ΔS=SA−SB\Delta S=S_A-S_B; absolute contrast is ∣ΔS∣|\Delta S|. All amplitudes use the same arbitrary reference. The flip-angle chart evaluates spoiled GRE at TE = 0 on a 0.25° grid; its contrast peak is only a sampled demonstration.

Single-echo Cartesian scan time is t=TR Nphase Navg Nextrat=TR\,N_{phase}\,N_{avg}\,N_{extra}, after converting TR to seconds. NextraN_{extra} is sequential 2D packages or acquired 3D partitions. Sung, slide 35, gives the independent 1000 ms × 256 × 1 = 4:16 example.

The isolated-pulse animation starts at equilibrium: Mz/M0=1−(1−cos⁡α)e−t/T1M_z/M_0=1-(1-\cos\alpha)e^{-t/T_1} and Mxy/M0=sin⁡α e−t/T2M_{xy}/M_0=\sin\alpha\,e^{-t/T_2}. It uses a uniform field, no refocusing and an arbitrary slowed precession rate. Its pulse is 90° for spin echo, or the selected flip angle for GRE; it is separate from the steady-state echo calculation.

Assumptions

  • One uniform, monoexponential compartment per tissue, identical receiver sensitivity and a common relative-signal scale. Proton density inputs are user assumptions; equal defaults isolate relaxation contrast.
  • T1 and T2 preset means are the “This study” columns of Stanisz et al. (2005), Table 1, printed p. 509 (attached s2, page 3). Measurements were in vitro at 37 °C: mouse liver/muscle/heart, rat kidney/spinal cord, bovine cartilage/brain/optic nerve, and human blood. These are not patient reference ranges.
  • The 1.5 T / 3 T choice changes the preset relaxation times only. It does not scale M0 with B0 or predict SNR.
  • T2 defaults of 40 and 50 ms are illustrative editable assumptions, not published tissue presets. GRE requires T2 ≤ T2. A physical Bloch compartment also requires T2 ≤ 2T1.
  • Ideal instantaneous RF pulses, complete spoiling and established GRE steady state. No inflow, diffusion, magnetization transfer, chemical shift, exchange, B0/B1 errors or receiver noise.
  • One acquired phase line per TR; each 3D partition is acquired, while 2D slices can be interleaved within one TR. Package count is supplied by the user. No timing overhead or acceleration is inferred.

Where it applies

Educational use for ideal spin echo and completely spoiled GRE with 0°–90° excitation. It does not model balanced SSFP, unspoiled GRE, inversion recovery, EPI or fast/turbo spin-echo contrast. For the approximate spin-echo formula use TE ≤ TR/10; longer ratios show a warning because the neglected refocusing effect on longitudinal recovery matters. TE = 0 and flip = 0 are mathematical limits, not practical protocols.

Displayed signal is relative amplitude, not scanner intensity, image quality, SNR or CNR. The swatches use a common grayscale and are not an image simulation. Clinical relaxation times depend on tissue condition, species, temperature, measurement technique and magnetic field; the preset means have experimental uncertainty and are not exact constants.

The scan-time estimate covers conventional single-echo Cartesian acquisitions only. It excludes preparation and dummy cycles, calibration, gating delays, phase oversampling, acceleration, partial Fourier and echo trains. Real package packing and hardware timing must come from the scanner. The animation uses a separate isolated pulse and intrinsic T2, rather than a free-induction T2* decay or a simulated repeated pulse train. Precession phase and speed are illustrative, not a Larmor-frequency prediction.

Use for learning, not diagnosis, clinical protocol prescription, implant/SAR assessment or contrast-agent dosing.

References

  • Stanisz GJ et al. T1, T2 relaxation and magnetization transfer in tissue at 3T. Magn Reson Med 54:507–512 (2005), Table 1. https://doi.org/10.1002/mrm.20605 . Public research-group copy: https://www.ultrasoundandmriforcancertherapy.ca/wp-content/uploads/2017/03/2005-3..pdf
  • Hargreaves B. Rapid gradient-echo imaging. J Magn Reson Imaging 36:1300–1313 (2012), Eq. 5 and Ernst-angle discussion. https://doi.org/10.1002/jmri.23742 ; open manuscript: https://pmc.ncbi.nlm.nih.gov/articles/PMC3502662/
  • Karakuzu A. Bringing Quantitative Magnetic Resonance Imaging Under One Umbrella. PhD thesis, Polytechnique Montréal (2022), printed pp. 26–28, Eqs. 2.5–2.6 and Fig. 2.19. https://publications.polymtl.ca/10224/1/2022_AgahKarakuzu.pdf
  • Sung K. Imaging Sequences III. UCLA M219, 28 February 2024, slides 25 and 35. https://labs.dgsom.ucla.edu/file/549184/m219_lecture16_imaging_sequences_iii.key.pdf
  • Ernst RR, Anderson WA. Application of Fourier transform spectroscopy to magnetic resonance. Rev Sci Instrum 37:93–102 (1966). https://doi.org/10.1063/1.1719961

Inputs

InputSymbolUnitDefault
Signal modelseqSpin echo · ideal 90°–180°
Repetition time, TR
Time between successive excitations of the same slice or slab.
TRms500
Echo time, TE
Time from excitation to the centre of the measured echo. TE = 0 is an ideal mathematical limit.
TEms15
Excitation flip angle
The spoiled steady-state model applies from 0° to 90°. Spin echo uses fixed 90° excitation and 180° refocusing.
alphadeg20
Enter custom T1 and T2
Presets are excised-tissue measurements at 37 °C from Stanisz et al. (2005), Table 1; many are animal tissues.
manualfalse
Preset field strength
Selects measured relaxation times only. It does not model the change in equilibrium signal or SNR with magnetic field.
field1.5 T
Tissue Atissue_AWhite matter (bovine)
Tissue Btissue_BGray matter (bovine)
T1 of tissue AT1_A_inms884
T2 of tissue AT2_A_inms72
T1 of tissue BT1_B_inms1124
T2 of tissue BT2_B_inms95
Relative proton density A
A user-assumed relative equilibrium signal, including proton density. Equal defaults isolate relaxation contrast; these are not measured tissue density data.
rho_A1
Relative proton density B
Use the same reference scale as A. No automatic magnetic-field or receiver-gain scaling is applied.
rho_B1
T2* of tissue A
Illustrative editable assumption, not a Stanisz measurement or a field-dependent preset. Enter an independently known T2* for your example; it must not exceed T2.
T2s_Ams40
T2* of tissue B
Illustrative editable assumption, not a tissue preset. T2* includes additional dephasing and must not exceed T2.
T2s_Bms50
Cartesian acquisition
One phase-encoded line per TR. No turbo echo train, acceleration, partial Fourier or overhead is included.
acq2D · interleaved slice package
Acquired phase-encoding lines
The number actually acquired per image or partition; one line is acquired each repetition.
N_phase256
Signal averages (NEX / NSA)N_avg1
Acquired 3D partitions
An additional phase-encoding dimension. This is not the reconstructed slice count.
N_part128
Sequential 2D slice packages
Slices interleaved within one TR share a package. Enter packages acquired in separate passes, not the total slice count. Use the scanner's package count.
N_pack1
Time after one RF pulse · ▶ animate
An isolated pulse starting at equilibrium, not the repeated sequence. Intrinsic T2 decay in a uniform field; precession is slowed for display.
t_animms40

Results

With the default inputs above, the calc gives:

ResultValue
Tissue A · relative signal0.3507
Tissue B · relative signal0.3066
Contrast A − B · relative signal0.0441
Ideal scan time2.13 min
Ernst angle A · spoiled GRE55.39 deg
Ernst angle B · spoiled GRE50.14 deg

Checked against sources

Each case below was calculated by this calc and agrees with the value from its source (within the tolerance where one is shown).

Stanisz Table 1 · 1.5 T white/gray matter means

Source: Independent measured means from Stanisz et al. (2005), DOI:10.1002/mrm.20605, Table 1, printed p. 509. Use This study columns, not literature values. These checks verify data selection; they do not validate modeled echo amplitudes.

Inputs: the default inputs above.

ResultSourceThis calc
Tissue A longitudinal relaxation884 ms884 ms
Tissue A intrinsic transverse relaxation72 ms72 ms
Tissue B longitudinal relaxation1,124 ms1,124 ms
Tissue B intrinsic transverse relaxation95 ms95 ms

Stanisz Table 1 · 3 T white/gray matter means

Source: Independent measured means from Stanisz et al. (2005), Table 1, printed p. 509, DOI:10.1002/mrm.20605. In vitro bovine brain at 37 °C. This case checks the field-strength lookup, not echo-signal accuracy.

Inputs: Preset field strength = 3 T; the others as above.

ResultSourceThis calc
Tissue A longitudinal relaxation1,084 ms1,084 ms
Tissue A intrinsic transverse relaxation69 ms69 ms
Tissue B longitudinal relaxation1,820 ms1,820 ms
Tissue B intrinsic transverse relaxation99 ms99 ms

Published Bloch example · TR 30 ms, T1 800 ms

Source: Karakuzu (2022), PhD thesis, Fig. 2.19, printed p. 28 (PDF page 58), and qMRLab introductory tutorial: a Bloch-simulated signal maximum around 16° for TR = 30 ms and T1 = 800 ms. The analytic value is about 15.59°; tolerance allows the published whole-degree rounding, without changing the equation.

Inputs: Signal model = Spoiled gradient echo · steady state; Enter custom T1 and T2 = yes; Repetition time, TR = 30 ms; Echo time, TE = 0 ms; T1 of tissue A = 800 ms; Excitation flip angle = 16 deg; the others as above.

ResultSourceThis calc
Ernst angle A · spoiled GRE16 ± 0.5 deg15.593 deg

UCLA worked timing example · 4 min 16 s

Source: Kyung Sung, UCLA M219 Imaging Sequences III, 28 February 2024, slide 35: TScan = 1000 ms × 256 × 1 = 4:16 [mm:ss], one echo per TR. Independently published worked scan-time result. The original PDF was readable via web search but returned HTTP 403 to the attachment downloader, so s5 is a source link.

Inputs: Repetition time, TR = 1,000 ms; the others as above.

ResultSourceThis calc
Ideal acquisition time256 s256 s

Related Resources

Responses

Your email address will not be published. Required fields are marked *

L ading...