Feasibility of Longitudinal Relaxation Rate Mapping with Non-Cartesian Sampling and Compressed Sensing on a 1.5T MR-Linac

Quantitative mapping of the longitudinal relaxation rate (R1=1/T1) is a major building block for several multiparametric MRI protocols intended for adaptive radiation therapy planning. The implementation of these protocols is challenging in anatomical sites that experience large physiological motion. To implement and validate a motion-resolved quantitative T1 mapping method on a 1.5T MR-Linac that combines non-Cartesian k-space sampling trajectories with compressed sensing (CS) reconstruction techniques. Four 3D non-Cartesian k-space trajectories were evaluated: radial and stack-of-stars sampling using half- and full-spoke coverage. A variable flip angle acquisition was performed using the spoiled gradient-echo sequence, and T1 mapping was validated using two standard phantoms. Gradient delay timing was optimized empirically to minimize trajectory-induced artifacts. Eight compressed sensing reconstruction strategies were tested using spatial and spatiotemporal regularization operators. Reconstructions were evaluated across multiple implementation parameters and ranked based on spatial resolution, bias, and variability. In vivo studies included one healthy volunteer and one patient undergoing radiotherapy to a target in the kidney. Motion-resolved imaging was performed using respiratory self-gating and phase-sorted reconstruction. All non-Cartesian trajectories demonstrated high repeatability and low longitudinal bias in phantom studies, with coefficients of variation below 3.3%. Radial half-spoke sampling achieved the shortest scan times and highest agreement with Cartesian benchmarks. Reconstruction methods incorporating spatiotemporal regularization maintained spatial resolution and quantitative accuracy across undersampling factors up to 20-fold. In human subjects, non-Cartesian T1 mapping provided improved accuracy and reduced variability in mobile abdominal tissues compared to Cartesian acquisitions, particularly in the kidney cortex and medulla, where motion artifacts led to overestimation and higher variance in the reference method. T1 mapping using non-Cartesian trajectories and compressed sensing reconstruction is feasible on a 1.5T MR-Linac. The proposed approach enables accurate, motion-resolved quantitative imaging within clinically practical acquisition times. These results support integration of quantitative T1 mapping into adaptive MR-guided radiotherapy workflows and establish a foundation for future development of multiparametric imaging and response-adaptive treatment strategies.