0 前言

       MRI因良好的软组织对比度、较高空间分辨率、无创伤、无电离辐射等优势^{[1, 2]}，被认为是腰椎及间盘病变的首选影像学检测^{[3, 4, 5]}。在扫描检查期间，患者必须长时间保持静止以防止运动伪影^[6]。尤其对于有严重腰椎病变患者，难以配合完成检查。因此，在更短时间内获得满足临床诊断需求的图像是目前腰椎及间盘MRI检查所面临的重要问题。传统敏感度编码（sensitivity encoding, SENSE）利用多通道相控阵线圈并行采集减少相位编码的数目，从而缩短扫描时间。但SENSE采集过程中，会出现幻像伪影、混叠伪影、卷积伪影，影响图像质量^[7]。压缩感知（compressed sensing, CS）利用随机化稀疏采样，采样频率远小于奈奎斯特采样定律^{[8, 9]}，应用于腰椎MRI可以显著缩短扫描时间，减少图像伪影^{[10, 11]}。以往研究将CS应用于腰椎三维（three dimensional, 3D）T2WI序列以及二维（two dimensional, 2D）矢状位梯度回波（gradient recalled echo, GRE）和快速自旋回波（turbo spin echo, TSE）T2WI序列^{[12, 13]}可行性良好。然而，临床上2D矢状位T1WI、T2WI和轴位T2WI应用更加广泛。本研究首次将CS应用于腰椎2D序列中，旨在探讨不同加速因子（acceleration factor, AF）对2D腰椎MRI图像质量与扫描时间的影响，并优化AF。

1 材料与方法

1.1 一般资料

       于2021年6月至2022年3月在大连医科大学附属第一医院前瞻性招募32例（男12例，女20例）受检者，行腰椎MRI检查，年龄（45.28±14.11）岁。入组标准：（1）无椎体或盆腔恶性肿瘤、转移瘤及放化疗病史；（2）半年内无腰背部手术史。排除标准：（1）怀疑新发腰椎压缩骨折者；（2）严重脊柱后凸，难以卧床者；（3）腰椎钛合金内固定患者；（4）有MRI检查禁忌证者，如体内有起搏器、支架、金属植入物及患有幽闭恐惧症等。本研究遵照《赫尔辛基宣言》，经大连医科大学附属第一医院伦理委员会批准（批件号：PJ-KS-KY-2021-10），所有被试同意并签署知情同意书。

1.2 检查方法

       使用Philips Ingenia CX 3.0 T MRI（Philips Healthcare, Best, the Netherlands）设备，12通道内置体线圈。所有被试均行常规和5种不同AF的腰椎矢状位TSE-T1WI和T2WI、轴位TSE-T2WI序列，具体为SENSE AF=2以及CS AF=2、3、4、5（图1）。（1）矢状位T1WI扫描参数：TR 324 ms，TE 9 ms，体素0.89 mm×2.01 mm×4.00 mm，层间距1 mm，信号平均采集次数（number of signal average, NSA）1，无加速、SENSE 2、CS 2、CS 3、CS 4、CS 5扫描时间分别为162.0 s、84.0 s、84.0 s、58.3 s、45.4 s、37.6 s；（2）矢状位T2WI扫描参数：TR 1600 ms，TE 90 ms，体素0.80 mm×1.83 mm×4.00 mm，层间距1 mm，NSA 1，无加速、SENSE 2、CS 2、CS 3、CS 4、CS 5扫描时间分别为103.0 s、51.3 s、51.3 s、35.2 s、25.6 s、22.4 s；（3）轴位T2WI扫描参数：TR 3200 ms，TE 129 ms，体素0.60 mm×0.95 mm×4.00 mm，层间距0.4 mm，NSA 2，无加速、SENSE 2、CS 2、CS 3、CS 4、CS 5扫描时间分别为216.0 s、115.0 s、115.0 s、76.0 s、64.0 s、50.9 s。

点击查看大图

图1  男，24岁，腰椎第5间盘退化。图1A～1F分别为加速因子（AF）为无加速、敏感度编码（SENSE）AF=2和压缩感知（CS）AF=2、3、4、5的矢状位T1WI-快速自旋回波（TSE）序列；1G～1L分别为AF为无加速、SENSE AF=2和CS AF=2、3、4、5的矢状位T2WI-TSE序列；1M～1R分别为AF为无加速、SENSE AF=2和CS AF=2、3、4、5的轴位T2WI-TSE序列。可见随着AF增大，噪声逐渐增大。
Fig. 1  Male, 24 years old, lumbar disc degeneration. 1A-1F show sagittal T1WI turbo spin echo (TSE) sequences with acceleration factor (AF) being non accelerated, sensitivity encoding (SENSE) AF=2, compressed sensing (CS) AF=2, 3, 4, and 5 respectively; 1G-1L show sagittal T2WI-TSE sequences with AF being non accelerated, SENSE AF=2, CS AF=2, 3, 4, and 5, respectively; 1M-1R show respectively axial T2WI-TSE sequences with AF being no acceleration, SENSE AF=2, CS AF=2, 3, 4, and 5。It can be seen that with the increase of AF, the noise gradually increases.

1.3 数据测量及分析

1.3.1 数据测量

       将扫描所得图像传至ISP（Philips, IntelliSpace Portall Version7）工作站，由两名观察者（观察者1为从事MRI工作5年的主管技师，观察者2为从事MRI工作4年的技师）在矢状位T1WI、T2WI椎体最大层面划定感兴趣区（region of interest, ROI），测量L1～L5椎体、间盘信号强度（signal intensity, SI）和噪声强度标准差（standard deviation, SD），ROI大小分别为140～150 mm²、40～50 mm²，并在L3～L5水平划定三个ROI测量脑脊液SI和SD，大小20～30 mm²，取平均值。在轴位T2WI划定ROI测量L1～L5椎体、间盘SI和SD值，ROI大小140～150 mm²，并在L3～L5水平划定ROI测量右侧腰大肌SI和SD，大小40～50 mm²（图2），取平均值。计算信噪比（signal to noise ratio, SNR）和对比噪声比（contrast to noise ratio, CNR）^[5]。矢状位T1WI、T2WI：SNR=SI椎体或间盘/SD脑脊液，CNR=（SI椎体或间盘-SI脑脊液）/SD脑脊液；轴位T2WI：SNR=SI椎体或间盘/SD腰大肌，CNR=（SI椎体或间盘-SI脑脊液）/SD腰大肌。

点击查看大图

图2  男，24岁，L5椎间盘囊肿。2A：在矢状位T2WI-快速自旋回波（TSE）序列椎体最大层面划感兴趣区测量椎体、间盘和脑脊液的信号强度（SI）和噪声强度（SD）；2B：在轴位T2WI-TSE序列间盘划感兴趣区测量间盘和右侧腰大肌的SI和SD。
Fig. 2  Male, 24 years old, L5 intervertebral disc cyst. 2A: Region of interest (ROI) is drawn on the largest section of sagittal T2WI-TSE sequence vertebral body to measure the signal intensity (SI) and standard deviation (SD) of vertebral body, intervertebral disc and cerebrospinal fluid; 2B: Diagram measure SI and SD of the intervertebral disc and right psoas major muscle by drawing ROI in axial T2WI-TSE sequence.

1.3.2 主观评分

       由两名观察者根据椎体间盘终丝马尾形态及边界清晰度、图像运动伪影大小以及诊断价值用五分法对图像进行评分（表1）。4分及4分以上满足诊断要求。

点击查看表格

表1  主观评分标准表
Tab. 1  Table of subjective scoring standards

1.4 统计学分析

       使用统计软件SPSS 21.0对数据进行分析。使用组内相关系数（intra-class correlation coefficient, ICC）、Kappa检验分析两观察者测量数据及主观评分的一致性，结果≥0.75为一致性良好。对各序列不同组间图像的SNR、CNR和主观评分采用单因素ANOVA分析，P＜0.05为差异具有统计学意义。

2 结果

2.1 主观评分一致性检验

       两观察者测量数据和主观评分一致性良好（ICC：0.878～0.997，Kappa：0.763～0.948；表2）。选择观察者1的测量数据和主观评分进行后续分析。

点击查看表格

表2  两观察者测量数据及主观评分一致性比较
Tab. 2  Comparison of the consistency of measurement data and subjective scores between the two observers

2.2 两组ADC值及主观评分的差异性检验

       ANOVA检验矢状位TSE-T1WI、T2WI、轴位TSE-T2WI不同AF间椎体和间盘的SNR、CNR及主观评分差异有统计学意义（P＜0.05；表3）。两两比较结果：当CS=4时，矢状位T1WI椎体和间盘的SNR、矢状位T2WI椎体和间盘的SNR、CNR和主观评分较常规序列差异有统计学意义（P＜0.05）；当CS=3时，轴位T2WI椎体SNR较常规序列差异有统计学意义（P＜0.05；图3）。

点击查看大图

图3  不同加速序列椎体信噪比（SNR）、椎体对比噪声比（CNR）、间盘SNR、间盘CNR箱式图。3A～3D：矢状位T1WI序列；3E～3H：矢状位T2WI序列；3I～3L：轴位T2WI序列。当压缩感知（CS）=4时，矢状位T1WI序列椎体和间盘的SNR、矢状位T2WI序列椎体和间盘的SNR、CNR与常规序列差异有统计学意义；当CS=3时，轴位T2WI序列椎体SNR与常规序列差异有统计学意义。SENSE为敏感度编码。
Fig. 3  Box diagram of signal to noise ratio (SNR) of vertebral body, contrast to noise ratio (CNR) of vertebral body, SNR of intervertebral disc and CNR of intervertebral disc in different acceleration sequence. 3A-3D: T1WI sagittal positions; 3E-3H: T2WI sagittal positions; 3I-3L: T2WI axial positions. When compressed sensing (CS)=4, SNR of sagittal T1WI vertebral body and intervertebral disc, SNR and CNR of sagittal T2WI vertebral body and intervertebral disc are significantly different from those of conventional sequences; When CS=3, the SNR of vertebral body on axial T2WI is statistically different from that of conventional sequence. SENSE is sensitivity encoding.

点击查看表格

表3  不同AF客观评价以及主观评分的差异性比较
Tab. 3  Comparison of differences in objective evaluation and subjective scores of different AFs

3 讨论

       本研究将CS应用于2D腰椎矢状位T1WI、T2WI和轴位T2WI三个序列，优化AF。结果表明当AF为3、3、2时，椎体SNR、CNR和主观评分与常规序列差异无统计学意义，扫描时间缩短56.65%，扫描效率提高一倍以上，相关研究国内尚未有报道。通过CS使腰椎MRI时间显著缩短，不仅可以提高检查成功率，还可提高设备的社会经济效益，具有较好的临床应用前景。

3.1 腰椎MRI检查常规序列的优势及局限性

       随着年龄的增长，腰椎退行性病变的发病率逐渐增高^{[14, 15]}，常出现慢性腰背疼、下肢麻木疼痛等症状，严重者引发压缩骨折^[16]。计算机断层成像（computed tomography, CT）可以较好显示腰椎形态特征，广泛用于脊柱外科手术计划和术后随访。但CT对椎间盘的显示较差，且电离辐射无法避免^{[17, 18]}。MRI因具有良好的软组织对比度、较高空间分辨率、无创伤、无电离辐射等优势，通过多序列、多方位扫描，用以评估椎体及间盘的形态学改变，被认为是腰椎及间盘病变的首选影像学检测^{[19, 20]}。腰椎MRI检查是临床最常见的检查部位之一，但因其检查时间相对较长^[21]，对于新发压缩骨折、严重脊椎后凸的患者，因其难以长时间保持仰卧位姿势，常会因不自主运动导致图像质量产生伪影。在保证图像质量的前提下缩短扫描时间，不仅可以显著提高检查成功率，还可以有效提升MR设备的社会经济效益。

3.2 CS原理及应用优势

       CS是由信号稀疏分解和逼近理论进一步发展形成的一种信号处理技术^{[22, 23]}。依据MR信号具有稀疏性和不相干性的特点，通过数字化随机稀疏采样的方式获得离散样本的信号^[24]，然后采用非线性重建算法迭代重建获得稀疏信号，实现了部分K空间数据采集成像^{[25, 26]}。由于采样频率远小于奈奎斯特采样定律，扫描时间明显缩短^{[27, 28]}。该技术最先由LUSTIG等^[29]应用于临床。目前CS的临床应用逐渐推广，各部位不同序列最佳AF优化成为了国内外研究热点^{[30, 31]}。浦仁旺等^[32]将CS应用于3D屏气胰胆管成像，对比传统序列，当AF为24时，图像质量仍能满足临床诊断需求，扫描时间仅为15 s。林青等^[33]通过比较心脏MRI CS序列与常规序列，结果表明采用CS序列的图像与常规序列图像无明显差异，心功能分析准确可靠，且检查时间明显缩短。王学东等^[34]利用结合CS的mDixon Quant序列进行肝脏脂肪定量，定量结果不受影响，且扫描时间明显缩短。

3.3 CS应用于腰椎2D成像的优势及可行性

       常规腰椎MRI检查包括矢状位T1WI、T2WI和轴位T2WI三个序列，一次检查时间约10 min，已成为临床最常见的检查之一。但对于部分受检者，因慢性腰背疼等原因难以忍受长时间检查，缩短扫描时间可以有效提高检查成功率。BRATKE等^[35]将CS与腰椎3D-T2WI序列结合，结果表明当应用AF=4.5时，扫描时间缩短39%，图像质量与应用SENSE的序列无明显差异。MORITA等^[12]将CS应用于腰椎3D-T2WI-VISTA序列，可以在保证图像质量的前提下明显缩短扫描时间。然而，由于3D序列扫描时间相对较长，2D序列在临床的应用更加广泛。所以本研究将CS与腰椎2D MRI序列相结合，探究CS对2D腰椎MRI图像的影响，以推广其临床应用价值。BRATKE等^[13]将CS应用于腰椎矢状位2D TSE-T2WI与GRE-T2WI序列扫描，结果表明TSE序列较GRE序列伪影减少，且可以使用更高的AF。对于TSE序列，AF为3时图像质量满足诊断需求，扫描时间缩短43%。本研究根据前者的研究，创新性地将CS应用于腰椎矢状位T1WI、T2WI和轴位T2WI三个序列。当AF为3时，矢状位T2WI序列椎体、间盘SNR、CNR和主观评分与常规序列无明显差异，扫描时间缩短65.83%（103.0 s vs. 35.2 s），与前者研究结果一致，但节省了更多的扫描时间。对于矢状位T1WI和轴位T2WI序列，AF分别为3和2时，扫描时间缩短64.01%（162.0 vs. 58.3 s）、46.76%（216.0 vs. 115.0 s），且对陈旧性椎体压缩骨折、许莫氏结节、腰间盘突出、椎体脂肪化等腰椎退行性病变显示良好，图像质量满足诊断需求。

3.4 本研究的局限性

       本研究存在以下局限性：（1）样本量较小，有待进一步扩大样本量；（2）椎体退行性病变因患者不同年龄、性别而异^[36]，因患者量较小，未将患者年龄、性别等因素纳入考虑范围内；（3）因扫描时间较长，将怀疑新发压缩骨折患者、严重脊椎后凸受检者排除，但后续临床扫描中显示，应用CS的2D序列对以上病变的显示良好，且大部分患者能够配合完成检查。

4 结论

       综上所述，对比传统序列，结合CS的腰椎2D序列可以在保证图像质量的前提下，有效缩短扫描时间，提高检查成功率。临床推荐AF分别为3、3、2行腰椎矢状位T1WI、T2WI和轴位T2WI序列扫描，对腰间盘突出、陈旧性椎体压缩骨折等退行性病变显示清晰，部分新发压缩骨折患者、严重脊椎后凸患者也可配合完成检查。且总扫描时间缩短56.65%，扫描效率提高一倍以上，可以显著提高MR设备的社会经济效益，减少患者预约等待时间，具有一定的推广应用价值。