Liu Lianli, Shen Liyue, Johansson Adam, Balter James M, Cao Yue, Vitzthum Lucas, Xing Lei
Department of Radiation Oncology, Stanford University, Palo Alto, California, USA.
Department of Electrical Engineering, Stanford University, Palo Alto, California, USA.
Med Phys. 2024 Apr;51(4):2526-2537. doi: 10.1002/mp.16845. Epub 2023 Nov 28.
Volumetric reconstruction of magnetic resonance imaging (MRI) from sparse samples is desirable for 3D motion tracking and promises to improve magnetic resonance (MR)-guided radiation treatment precision. Data-driven sparse MRI reconstruction, however, requires large-scale training datasets for prior learning, which is time-consuming and challenging to acquire in clinical settings.
To investigate volumetric reconstruction of MRI from sparse samples of two orthogonal slices aided by sparse priors of two static 3D MRI through implicit neural representation (NeRP) learning, in support of 3D motion tracking during MR-guided radiotherapy.
A multi-layer perceptron network was trained to parameterize the NeRP model of a patient-specific MRI dataset, where the network takes 4D data coordinates of voxel locations and motion states as inputs and outputs corresponding voxel intensities. By first training the network to learn the NeRP of two static 3D MRI with different breathing motion states, prior information of patient breathing motion was embedded into network weights through optimization. The prior information was then augmented from two motion states to 31 motion states by querying the optimized network at interpolated and extrapolated motion state coordinates. Starting from the prior-augmented NeRP model as an initialization point, we further trained the network to fit sparse samples of two orthogonal MRI slices and the final volumetric reconstruction was obtained by querying the trained network at 3D spatial locations. We evaluated the proposed method using 5-min volumetric MRI time series with 340 ms temporal resolution for seven abdominal patients with hepatocellular carcinoma, acquired using golden-angle radial MRI sequence and reconstructed through retrospective sorting. Two volumetric MRI with inhale and exhale states respectively were selected from the first 30 s of the time series for prior embedding and augmentation. The remaining 4.5-min time series was used for volumetric reconstruction evaluation, where we retrospectively subsampled each MRI to two orthogonal slices and compared model-reconstructed images to ground truth images in terms of image quality and the capability of supporting 3D target motion tracking.
Across the seven patients evaluated, the peak signal-to-noise-ratio between model-reconstructed and ground truth MR images was 38.02 ± 2.60 dB and the structure similarity index measure was 0.98 ± 0.01. Throughout the 4.5-min time period, gross tumor volume (GTV) motion estimated by deforming a reference state MRI to model-reconstructed and ground truth MRI showed good consistency. The 95-percentile Hausdorff distance between GTV contours was 2.41 ± 0.77 mm, which is less than the voxel dimension. The mean GTV centroid position difference between ground truth and model estimation was less than 1 mm in all three orthogonal directions.
A prior-augmented NeRP model has been developed to reconstruct volumetric MRI from sparse samples of orthogonal cine slices. Only one exhale and one inhale 3D MRI were needed to train the model to learn prior information of patient breathing motion for sparse image reconstruction. The proposed model has the potential of supporting 3D motion tracking during MR-guided radiotherapy for improved treatment precision and promises a major simplification of the workflow by eliminating the need for large-scale training datasets.
从稀疏样本进行磁共振成像(MRI)的容积重建对于三维运动跟踪是可取的,并有望提高磁共振(MR)引导的放射治疗精度。然而,数据驱动的稀疏MRI重建需要大规模训练数据集进行先验学习,这在临床环境中获取既耗时又具有挑战性。
通过隐式神经表示(NeRP)学习,利用两个静态3D MRI的稀疏先验辅助,研究从两个正交切片的稀疏样本进行MRI的容积重建,以支持MR引导放疗期间的三维运动跟踪。
训练一个多层感知器网络来参数化特定患者MRI数据集的NeRP模型,该网络将体素位置和运动状态的4D数据坐标作为输入,并输出相应的体素强度。通过首先训练网络来学习具有不同呼吸运动状态的两个静态3D MRI的NeRP,通过优化将患者呼吸运动的先验信息嵌入到网络权重中。然后通过在内插和外推运动状态坐标处查询优化后的网络,将先验信息从两个运动状态扩展到31个运动状态。从先验增强的NeRP模型作为初始化点开始,我们进一步训练网络以拟合两个正交MRI切片的稀疏样本,并通过在3D空间位置查询训练后的网络获得最终的容积重建。我们使用具有340 ms时间分辨率的5分钟容积MRI时间序列对7例肝细胞癌腹部患者进行了评估,该序列使用黄金角径向MRI序列采集并通过回顾性排序进行重建。从时间序列的前30秒分别选择两个具有吸气和呼气状态的容积MRI进行先验嵌入和增强。其余4.5分钟的时间序列用于容积重建评估,我们在其中将每个MRI回顾性下采样到两个正交切片,并在图像质量和支持三维目标运动跟踪的能力方面将模型重建图像与真实图像进行比较。
在评估的7例患者中,模型重建的MR图像与真实MR图像之间的峰值信噪比为38.02±2.60 dB,结构相似性指数测量值为0.98±0.01。在整个4.5分钟时间段内,通过将参考状态MRI变形为模型重建和真实MRI估计的大体肿瘤体积(GTV)运动显示出良好的一致性。GTV轮廓之间的95%分位数豪斯多夫距离为2.41±0.77 mm,小于体素尺寸。在所有三个正交方向上,真实情况与模型估计之间的平均GTV质心位置差异小于1 mm。
已开发出一种先验增强的NeRP模型,用于从正交电影切片的稀疏样本重建容积MRI。仅需一次呼气和一次吸气的3D MRI来训练模型,以学习患者呼吸运动的先验信息用于稀疏图像重建。所提出的模型具有在MR引导放疗期间支持三维运动跟踪以提高治疗精度的潜力,并有望通过消除对大规模训练数据集的需求而大大简化工作流程。