We investigated combining UTE sequences with different TR without requiring knowledge of T1, to enable increasing the number of short TE scans for T2* quantification. Many short T2* tissues have multiple compartments with ultra-short and somewhat longer T2 values. To quantify both a substantial number of images with ultra-short TE and images with a substantial maximal TE are required. The large maximal TE requires relatively large TR and hence long scan times, while the ultra-short TE scans have to be acquired separately. Hence, being able to combine images with different TR would be beneficial for such studies.
UTE imaging is used to study tissues that have a compartment with very short T∗2(~1ms) and longer T∗2(~20-30ms). To quantify both, several images with short (<3ms) TE and a few with longer (~30ms) TE have to be acquired. This requires multiple (multi-echo) UTE acquisitions, conventionally acquired with identical TR>max TE. However, to acquire sufficient images with very short TE this scheme acquires unnecessarily many images with a long TE. Hence, we propose reducing the TR of some of the multi-echo scans, while preserving T1 weighting by appropriately choosing the flip angle (FA).
The UTE-SPGR pseudo steady state signal is given by [1] S(TR,FA,TE)=S(TR,FA)f(TE)(Eq.1) with S(TR,FA)=(1−E1)sin(FA)1−E1cos(FA), E1=exp(−TRT1).
We want to use images with different echo time (TE), flip angle(FA) and repetition time (TR) to characterize f(TE), without having to quantify T1. Hence we want a relation between TR1, FA1, TR2, and FA2 of two acquisitions such that the ratio of signals is independent of T1. Unfortunately, no closed form analytic relation obtaining that can be derived. Hence, a local relation with T1 is solved for small FA. A series expansion of that solution in FA1, TR2, and TR2, assuming small FA and TR<<T1, is given by FA2=FA1√TR2TR1(Eq.2) with an intensity scaling C=TR1tan(FA2/2)TR2tan(FA1/2)≈S(TR2,FA2)S(TR1,FA1)(Eq.3) between the images. Note that Eq.[2] is independent of T1 and hence images satisfying this relation satisfy it for tissues with any T1.
Results and discussion
Figure 1 shows equation 3 and the numerical solution. Up to FA_1=39°, the approximation holds to good accuracy; <0.5° difference in FA. Using the numerical optimized values, the errors as quantified with \frac{S(TR_1,FA_1 ) C-S(TR_2,FA_2)}{S(TR_1,FA_1 ) C+S(TR_2,FA_2) } are less than 0.06% for all FA and T_1.
Figure 2 shows the ratio of the phantom images. It shows that this ratio is independent of the T_1 and approximately 4% below the expected C. This small mismatch could be due to imperfect spoiling as the phantom does not have short T_2^*. However, the difference is constant throughout the image and hence it can be compensated by an additional calibration.
Figure 3 shows an in-vivo validation of the relation. The scaled axial views of a single UTE image of the two acquisitions differing in FA and TR are essentially equal; the difference image shows mainly zero mean noise inside the tendon, which is indicated by a contour. Note the factor 12.5 difference in intensity scale. Some residual motion might cause the main edge artifacts seen in the difference image.