胰腺癌放疗图像引导锥束CT运动伪影修正技术研究.docx

胰腺癌放疗图像引导锥束CT运动伪影修正技术研究
罗 辰，任 青，刘建强，牛田野
Development of Motion Artifact Correction Solutions for the Cone-beam CT Images during Pancreatic Cancer Image-guided Radiotherapy
LUO Chen, REN Qing, LIU Jianqiang, and NIU Tianye
在线阅读 View online: https://doi.org/10.15953/j.ctta.2022.066

您可能感兴趣的其他文章
Articles you may be interested in

一种锥束CT动态Bowtie的蒙特卡罗模拟研究
Monte Carlo Simulation for a Dynamic Bowtie Filter with Cone-beam CT CT理论与应用研究. 2018, 27(3): 315-327
工业X射线CT中基于深度学习的射束硬化伪影抑制方法
Deep Learning Based Beam Hardening Artifact Reduction in Industrial X-ray CT CT理论与应用研究. 2018, 27(2): 227-240
可同时抑制多种图像伪影的最优骨校正
An Optimal Bone Correction Capable of Simultaneously Suppressing Various Types of Image Artifacts CT理论与应用研究. 2018, 27(3): 301-314
口腔锥形束CT剂量评估模式的研究现状
Status of Study in Dental Cone Beam Computed Tomography Dose Estimation Modes CT理论与应用研究. 2019, 28(1): 1-11
PET/CT影像辅助CT引导下肺穿刺活检的应用
Application of PET/CT Image Assisted CT Guided Lung Biopsy CT理论与应用研究. 2020, 29(4): 503-510
能谱CT单能量成像结合MAR技术降低金属植入物伪影的体模研究
Utility of Spectral CT Monochromatic Imaging with Metal Artifacts Reduction (MAR) for the Reduction of Metal Artifacts of Embolization Coil Implants
CT理论与应用研究. 2019, 28(5): 529-539

第 31 卷 第 6 期 CT 理论与应用研究 Vol.31, No.6
2022 年 12 月（761-771） CT Theory and Applications Dec., 2022
罗辰, 任青, 刘建强, 等. 胰腺癌放疗图像引导锥束 CT 运动伪影修正技术研究[J]. CT 理论与应用研究, 2022, 31(6): 761-771. DOI:10.15953/j.ctta.2022.066.
LUO C, REN Q, LIU J Q, et al. Development of motion artifact correction solutions for the cone-beam CT images during pancreatic cancer image-guided radiotherapy[J]. CT Theory and Applications, 2022, 31(6): 761-771. DOI:10.15953/j.ctta.2022.066. (in Chinese).

胰腺癌放疗图像引导锥束 CT 运动伪影修正技术研究
罗辰1,2✉ , 任青3, 刘建强4, 牛田野2,5✉
1.浙江机电职业技术学院自动化学院, 杭州 310053 2.深圳湾实验室, 广东 深圳 518107
3.中能医用加速器系统（广东）有限公司, 广东 东莞 523808 4.江苏康众数字医疗科技股份有限公司, 江苏 苏州 215123
5.浙江大学医学院, 杭州 310058
摘要：基于二维平板探测器技术的锥束 CT 系统广泛应用于放射治疗前的摆位验证，但在以胰 腺癌为代表的腹腔肿瘤放疗应用中会遇到呼吸运动和肠道蠕动导致的严重阴影与条状伪影问 题，难以从 CT 影像中分辨病灶区域。由于胰腺等柔性脏器在呼吸运动作用下会存在非刚性形 变，体表运动监测结果和实际脏器运动间存在难以量化的偏差，也难以监测以肠道蠕动为代 表的不规则运动，故腹腔锥束 CT 运动伪影修正问题缺乏有效解决方案。本文基于生物动力学 理论和人体生理学常识，提出全新的无需运动监测或体内标记物植入的放疗图像引导锥束 CT 运动伪影修正方法，设计基于伪影图像特征、融合各类 CT 图像域处理算法的伪影修正策略。 将该策略应用于临床腹腔锥束 CT 影像处理后所得结果表明锥束 CT 图像质量有显著改善，典 型软组织区域的平均 CT 数误差从 90HU 降到 30HU，肠道空腔边界和周边软组织信息得到部分 恢复。本文开发的伪影修正策略无需呼吸门控或增加投影数，避免标记物植入手术，可集成 到现有工作流，为胰腺癌放疗图像引导提供肿瘤定位信息。
关键词：锥束CT；运动伪影修正；放疗图像引导；数字图像处理
DOI:10.15953/j.ctta.2022.066 中图分类号：TP 391.41；R 814 文献标识码：A

在肿瘤放射治疗中，治疗计划制定一般以患者在多排螺旋 CT 设备上的离线扫描结果即 计划 CT 影像为基准，因此在每个疗程剂量投放前需要确认患者及肿瘤靶区相对于加速器头 及多叶光栅的位置，为实现剂量投放前患者的现场 CT 扫描，当代放射治疗系统中普遍集成 了锥束 CT 设备，此时无需移动患者即可得到CT 影像；该影像用于指导患者摆位，并可提 供肿瘤定位信息。
锥束 CT 利用 X 射线球管发出的立体锥形束和高分辨大尺度平板探测器实现三维立体 CT 成像[1] ，保留了常规多排诊断 CT 解剖结构无畸变的成像特性，可真实反映肿瘤形态和位 置。锥束 CT 相对于诊断 CT 具有 3 个显著优势：结构简单，体积小，可集成进放疗加速器 治疗系统；探测器成像单元线度小，且成像系统空间分辨率高；总辐射剂量较低，可降低
收稿日期 ：2022-04-16。
基金项目 ：国家重点研发计划政府间国际科技创新合作重点专项（基于多自由度四维双能锥束 CT 的胰腺 癌放射治疗运动控制（2018YFE0114800））；浙江省教育厅 2021 年度高校国内访问工程师“校 企合作项目”（无需运动监测的CT 运动伪影修正方法研究（FG2021043））；恶性肿瘤发病机制 及转化研究教育部/北京市重点实验室开放基金（2022 开放-2）。

762 CT 理论与应用研究 31 卷

患者辐射损伤风险。因此，获得准确、高质量的锥束 CT 影像有助于提升放疗疗效，基于锥 束 CT 的图像引导放射治疗（imaged-guided radiation therapy，IGRT）技术是放疗研究 热点[2] ，锥束 CT 也成为 IGRT 最常用的成像手段。
患者每次接受放疗前现场锥束 CT 扫描获得的 CT 影像可提供解剖结构变化信息，包括患 者体位、脏器位置及肿瘤靶区位置等；但由于锥束 CT 扫描速度较诊断 CT 慢数十倍，扫描中 脏器运动不可忽略，则在锥束 CT 影像中会出现运动导致的伪影，影像质量受到显著影响。在 胰腺癌放疗过程中胰腺等柔性脏器的运动不仅会造成术前锥束 CT 影像中的严重伪影问题，进 而导致肿瘤区域定位精度变差，还会直接影响放疗剂量投放的精度。因此，锥束 CT 扫描过 程中胰腺等腹腔脏器运动造成的各类伪影问题已严重限制了该技术在胰腺癌 IGRT 中的应用。
扫描过程中的主要运动包括周期性的肺部呼吸运动和不规则的肠道蠕动。在这些来源 和特征不同的脏器运动以及柔性脏器非刚性形变的共同作用下，腹腔锥束 CT 影像中的伪影 主要呈现为腹腔各脏器附近的条状伪影和阴影伪影。此时从锥束 CT 影像中难以提取消化道 肿瘤的定位信息，且基于骨结构分辨和配准的放疗前患者几何摆位精度也受限。
图1对同一位患者的放疗术前锥束 CT（瓦里安公司 Trilogy 放疗机 OBI 系统[3] ）和诊 断 CT（西门子公司 64 排 Somatom Definition AS CT 设备）的单层横断面图像进行比较， 直观展示锥束扫描过程中腹腔脏器运动引起的图像伪影和成像质量恶化。腹腔锥束 CT 扫描 时长受限于所采用平板探测器的读出和数据传输速率，一般需要几十秒，其中Trilogy 系 统的OBI 设备完成患者腹腔扫描需 60 s 以上。伪影形态则主要呈现为肠道空腔附近的条状 和阴影伪影（图1（a）中红色箭头所示），此时难以区分肌肉、脂肪及其他软组织，且无 法获得脏器的准确位置和形状。相对地，诊断 CT 的腹腔扫描时间一般在 1 s 内，此时腹腔 脏器可近似为准静止状态，不会出现运动引发的伪影现象（图1（b））。

（a） （b）
图 1 同一位患者的腹腔锥束 CT（a）和诊断 CT（b）扫描横断面图像（显示窗：[-500, 500] HU）
Fig.1 Comparison between the transverse images from abdominal cone-beam CT and diagnostic CT scans of one patient (Display window: [-500, 500] HU)
CT 扫描过程中腹腔脏器运动主要来源于胸腔呼吸运动和腹腔肠道蠕动，其中呼吸运动 为间接影响，此时需要考虑呼吸作用下腹腔柔软脏器的非刚性形变；而肠道蠕动为非刚性、 不可预测的不规则运动，此时难以找到运动规律。因此，腹腔锥束 CT 伪影修正的最大难点 是腹腔脏器的不规则蠕动和非刚性形变情况下无法获得准确的运动参量信息。由于缺乏准 确定位肿瘤靶区的方法， 目前胰腺癌等放疗中会增加靶区范围[4] ，此时会有更大比例的正 常组织受到辐照，其中可能包含小肠等危机器官（organ at risk，OAR），进而导致剂量利 用率降低和术后不良反应风险增加。另外在强子治疗如质子放疗中，带电粒子的绝大部分

6 期 罗辰等：胰腺癌放疗图像引导锥束 CT 运动伪影修正技术研究 763

能量沉积在布拉格峰，此时若没有准确的肿瘤靶区定位及放疗计划，绝大部分能量均会作 用于正常组织，引发严重后果。
目前 IGRT 中脏器运动引发的锥束 CT 伪影的相关研究集中于解决呼吸运动造成的伪影， 而针对非刚性、不规则的肠道蠕动的研究十分有限。为抑制呼吸运动，早期曾采用训练患 者在诊断扫描和治疗过程中屏息的方法[5] ，但其临床实用性差，实际应用中一般仅能抑制 最大运动幅度，即确保不会出现剧烈运动。随着光学成像技术的发展，可利用双目成像等 具备三维坐标测量功能的设备对患者体表贴附的标识点进行呼吸跟踪，并以此为基础进行 锥束 CT 投影图的呼吸时相分组，从而得到分呼吸时相的锥束 CT 影像[6-8] 。这种相位相关方 法也被称为四维锥束 CT（four-dimensional cone-beam CT，4D CBCT），已逐步应用于各 类肺癌放射治疗中的胸腔图像引导[9-10] 。需要注意的是 4D CBCT 中投影数是常规锥束 CT 的 数倍，相应辐射剂量和扫描时间有显著提升；同时分时相后相位划分不准确引发的投影图 不连续性以及稀疏角度重建均会造成图像质量的显著降低[11] 。另外 2018 年 Steiner 等[12]发 现现有4D CBCT 技术应用于肺部肿瘤定位时会低估肿瘤的实际运动幅度，该误差部分来源 于体表标识点和肺部肿瘤定位的差异，但和重建算法的选取无关。
近年来，4D CBCT 技术也进一步推广应用于腹腔的肝癌[13] 以及胃肠道肿瘤如胃癌[14]、 食管癌[15]及胰腺癌[16]放射治疗中的肿瘤辅助定位。考虑到腹腔脏器的不规则运动和在呼吸 等运动作用下柔性脏器的形变，一般采用在肿瘤区域附近植入金属标识点的方式来辅助定 位肿瘤。针对 CT 影像伪影和肿瘤定位问题，Marchant 等[17]于2011年提出了基于仿射变换 的投影图修正方案，该方案中首先手动提取投影图中各个标识点的定位，再设计模型计算 得出无运动情况下的标识点定位，从而获得投影图修正所需的变换参数，最终锥束 CT 影像 中伪影得到抑制，肿瘤区域定位精度提升；2012 年该课题组进一步开发集成了标识点自动 提取算法，提升了计算效率[18] ，但该方案中的有限个标识点无法准确描述肿瘤区域的复杂 形状变化，同时金属标识点自身会引发 CT 影像中的金属伪影，另外临床实践中植入操作增 加了手术风险和复杂度，且植入物可能在患者体内迁移。
为避免金属植入物问题，无需植入标识点的运动提取方法是近年的研究热点，主要包 括基于配准和迭代重建的两类解决方案。Wang 等[19]于2014 年提出了基于 4D CBCT 和二维在 线 X 射线透视图配准的肿瘤跟踪定位方法，利用 4D CBCT 生成二维数字重建 X 射线影像 （digital reconstructed radiograph，DRR），再将 DRR 和 X 射线透视图做二维配准获得 肿瘤定位；Zachiu 等[20]于2018 年设计了用于计划 CT 和锥束 CT 影像配准的 Evolution 算法， 该算法中仅需两个输入参数，大大降低了参数优化难度，其配准结果用于在线放疗图像引 导；2015 年 Liu 等[21]建立了五维呼吸模型，并应用于 4D CBCT 的迭代重建，无需呼吸跟踪 数据即可得到 4D CBCT 重建结果；2017 年 Zhang 等[22]则将全变分方法应用于4D CBCT 的迭 代重建；2021 年 Peterlik 等[23]在美国瓦里安公司最新的Halcyon 放疗机上进行基于迭代优 化算法的运动伪影修正算法研究，在 iCBCT 迭代重建中加入梯度权重，将平均 CT 数误差 （考虑 5 s 内的平移运动）从 136.0HU 降低至 56.0HU。 目前基于配准的解决方案和植入物 方案相比加入了柔性脏器非刚性形变因素的影响，有效提升了定位精度，但需要额外进行 X 射线透视或 CT 扫描，同时缺乏普适的配准参数优化方法，可以造成配准后的影像失真。 而基于迭代重建的解决方案则受限于简化模型的假设是否成立以及迭代重建计算效率问题。
腹腔锥束 CT 影像中肠道蠕动引发的伪影修正研究可追溯到 2007 年 Voogt 等[24]在放射

764 CT 理论与应用研究 31 卷

治疗中的计算机应用国际会议（International Conference on the Use of Computers in Radiation Therapy, ICCR）上报道的投影筛选方法，即 Flatex 算法[25] ，该算法主要解决 锥束 CT 影像中肠道内气泡附近的条状伪影问题，其中首先筛选对应较小幅度运动的部分投 影，再进行基于少量投影的解析重建；该算法的应用受限于人工筛选投影的效率，且稀疏 投影重建容易引入新的条状伪影。近年来没有出现后续或者相关研究报道，无法建立准确 运动模型的肠道蠕动已成为腹腔锥束 CT 伪影修正中的最大难题。
综上所述， 目前腹腔脏器运动引发的锥束 CT 伪影修正方法主要包括患者屏息训练、体 表呼吸跟踪、体内植入标识点辅助定位、基于图像配准的运动提取、迭代重建中运动建模、 投影图筛选等方式。这些方法中脏器运动参数的获取依赖于非直接测量或数学模型假设， 且未同时覆盖呼吸运动和肠道蠕动两大伪影来源，伪影消除效果受限；另外其运动测量估 计或配准提取过程对诊疗方案的要求更为苛刻，且更多的投影数也对应了单次扫描中总辐 射剂量的提升。
针对上述问题，本文一方面从锥束 CT 解析重建理论出发，挖掘重建图像中的实、频域 信息，设计包含一系列信号及图像处理算法的伪影修正工作流，对包含伪影的腹腔锥束 CT 影像进行处理，避免运动监测的精度和复杂度问题以及运动模型估计的误差问题。另一方 面，本文按照伪影的形态特征将腹腔锥束 CT 影像中的伪影划分为阴影和条状伪影并分别进 行处理，实现来源和类型不同的呼吸运动和肠道蠕动的同时处理，避免常规伪影修正方法 中无法直接获取肠道运动模型的难题。
1 运动伪影修正图像处理流程
针对胰腺癌放疗图像引导锥束 CT 运动伪影修正问题，本文设计的图像处理流程图如 图 2所示。首先进行数字模体仿真实验和临床诊断影像的特征分析，结合腹腔器官动力学 和人体生理学经验知识，建立运动伪影的图像特征和运动种类及定量参数间的关系，并以 此为基础选取降噪、分割、正弦图插值、空间滤波等操作的初始参数。接着对腹腔锥束 CT 影像进行预处理：先进行 CT 图像域降噪，其中注意保持边界；再进行腹腔区域分割，即 在 CT 影像中去除病床、衣物等干扰信息。然后在经过粗分割的腹腔区域中根据各器官的一 般线衰或 CT 值进行第 2 次图像分割，分离出高对比度的肠道空腔区域（最暗）和类骨结构 区域（最亮）；在此基础上针对肠道空腔区域邻域进行基于正弦图插值的阴影伪影修正，并 针对肠道空腔和类骨结构区域的更大邻域进行基于空间滤波的条状伪影修正。另外在运动 伪影已经严重影响感兴趣区域分割结果的情况下，需要从修正后的CT 影像中重新提取感兴 趣区域并进行上述操作。
参数设置
参数优化

预处理
阴影伪影修正
条状伪影修正
（降噪、分割）
（正弦图插值）
（空间滤波）

腹腔锥束 CT 影像 感兴趣区域 阴影伪影修正 CT 影像 运动伪影修正 CT 影像
图 2 胰腺癌放疗图像引导锥束 CT 图像域运动伪影修正策略
Fig.2 Proposed image-domain motion artifact correction strategy
for CBCT in IGRT of pancreatic cancer

6 期 罗辰等：胰腺癌放疗图像引导锥束 CT 运动伪影修正技术研究 765

预处理的目标是为后续感兴趣区域分割提供更高质量的 CT 影像，并排除病床、衣物影 像等干扰因素。预处理流程主要包括降噪和分割两部分。为在降噪的同时保持感兴趣区域 的边界，采用基于块匹配的滤波方法[26]（block-matching 3D filtration，BM3）进行降噪。 由于腹腔区域的外轮廓（皮肤）一般不会受到运动伪影影响，边界锐度较好，降噪后的 CT 影像利用主动轮廓法经过少量迭代运算即可快速收敛到腹腔区域轮廓（一般近似椭圆）。在 预处理中如果发现锥束 CT 影像中散射污染严重，则需要进行基于模板图的快速图像域散射 修正，这里的模板图就是上述主动轮廓法分割的结果。
阴影伪影修正算法流程如图 3所示，核心操作是归一化后的正弦图插值。首先利用阈 值或分水岭等常规方法进行肠道空腔区域分割，并通过图像形态学膨胀操作获得拓展的感 兴趣区域，再通过二值化操作获得区域标记，并将该区域重新赋值为常数，该常数选取为 该区域所有非零像素值的中位数，得到模板图。接着通过正向投影算法获得肠道空腔区域、 感兴趣区域、以及模板图对应的正弦图，在处理单层 CT 图像时可简化为扇束位形，直接应 用 Radon 变换。其中考虑到正弦图上直接插值会因为数值不连续导致新的伪影出现，在插 值前需要进行正弦图的逐点归一化操作，其归一化参数因子来源于模板图对应的正弦图。 采用二维线性插值方法处理归一化后的正弦图，并在插值结果中补回归一化参数因子即可 利用Radon 逆变换或解析重建得到阴影伪影修正后的 CT 影像。
投影变换 正弦图插值（归一化方法） 阴影伪影补偿
分割

腹腔锥束 CT 影像 肠道空腔 空腔对应正弦图 阴影伪影图像组分 阴影伪影修正 CT 影像 膨胀 反投影变换

感兴趣区域拓展 感兴趣区域正弦图 插值后的正弦图
区域标记

模板图 模板正弦图（归一化）
图 3 胰腺癌放疗图像引导锥束 CT 图像域阴影伪影修正流程
Fig.3 Proposed shading artifact correction flowchart for CBCT in IGRT of pancreatic cancer
条状伪影修正算法流程如图 4所示，除了获得虚拟投影数据的投影变换外，主要包括 空间带通滤波以及正弦图的布尔量操作。首先需要对条状运动伪影的形态特征进行统计分 析，得出其代表性的特征是在空间频率域呈现中频分布特征；接着设计空间窄带滤波器提 取条状伪影对应形态图，并利用正弦图布尔量操作进一步提取条状伪影对应区域；最后在 修正前图像上直接扣除条状伪影对应图像即可得到条状伪影修正结果。
2 CT 图像处理结果与分析
2.1 CT 影像数据来源说明
本文处理的临床 CT 影像数据来源于 2015 至 2018 年间浙江大学医学院附属邵逸夫医院

766 CT 理论与应用研究 31 卷

放疗科胰腺癌患者在放疗前为辅助摆位而进行的日常锥束 CT 扫描，同一患者在锥束 CT 扫 描前 2 天内进行的计划 CT 扫描获得的高质量 CT 影像用于参考。其中锥束 CT 扫描采用美国 瓦里安公司 Trilogy 放疗机上配备的 OBI锥束 CT 系统的胸腹腔半扇大视野模式，计划 CT 扫描采用德国西门子公司的 Somatom 系列 CT 机的腹腔模式。本文处理的影像数据来源包 括 2 例患者的腹腔锥束 CT 和计划 CT 扫描数据，其中锥束 CT 的主要扫描参数如表 1所示。
投影变换 正弦图布尔量操作 条状伪影补偿
滤波

中频组分 中频对应正弦图 条状伪影图像组分 条状伪影修正 CT 影像
反投影变换
分割

感兴趣区域拓展 高对比度亮度 亮区对应正弦图 高对比度区域
中频正弦图

肠道空腔 空腔对应正弦图
图 4 胰腺癌放疗图像引导锥束 CT 图像域条状伪影修正流程
Fig.4 Proposed streaking artifact correction flowchart for CBCT
in IGRT of pancreatic cancer
2.2 图像处理结果
经过本文设计的运动伪影修正工作流处理，阴影和条状等形态的运动伪影得到显著抑 制，且噪声水平有明显改善。
图 5 是锥束 CT 单层横断面图像经过流程处理的中间结果以及和计划CT 参考图的对比。 预处理降噪显著抑制噪声水平，对于图 5（f）给出的感兴趣区域进行分析，发现降噪后峰 值信噪比（peak signal-to-noise ratio，PSNR）提升约 15 dB。进一步经过阴影伪影和条 状伪影修正，相较于经过配准的计划 CT 参考图，结构相似性 （structural similarity， SSIM）从 0.844 提升至 0.939，感兴趣 区域 CT 数 的平均绝对值误差（ Mean absolute error，MAE）从约 92HU 降低至约 34HU，而均方根误差（root mean square error，RMSE）

表 1 Trilogy 放疗机OBI锥束 CT 系统的扫描参数 Table 1 Scan parameters of the OBI CBCT sys- tem on Trilogy radiotherapy device

扫描参数名
参数值
源到探测器距离（SID）/mm
1501
源到转轴距离（SAD）/mm
1000
探测器水平偏移量/mm
148
扫描电压/kVp
125
扫描电流/mA
80
扫描模式
腹腔、半扇
投影角数量
663


表 2 腹腔锥束 CT 运动伪影修正后图像
质量指标变化
Table 2 Changes in image quality
indicators before and after abdominal CBCT motion artifact correction

图像质量指标
修正前
修正后
整体图像 PSNR/dB
41.562
61.110
感兴趣区域PSNR/dB
47.593
62.263
整体图像 SSIM
0.844
0.939
感兴趣区域 CT 数MAE/HU
92.348
33.762
感兴趣区域 CT 数 RMSE/HU
98.986
42.873

6 期 罗辰等：胰腺癌放疗图像引导锥束 CT 运动伪影修正技术研究 767

从约 99HU 降低至约 43HU。经过处理，处理前的完全淹没在噪声和伪影中的肠道空腔边界 和肠道周边软组织结构得到部分恢复。表 2 列出了锥束 CT 运动伪影修正前后的关键图像指 标的对比结果。

（a）未处理图像 （b）预处理降噪后 （c）阴影伪影修正后

（d）条状伪影修正后 （e）计划 CT 参考图 （f）感兴趣区域选取
图 5 腹腔锥束 CT 运动伪影修正效果（横断面视图）
Fig.5 Motion artifact correction results for one abdominal CBCT scan (transverse view)
图 6 是同一位患者腹腔锥束 CT 矢状位视图中的运动伪影修正效果，红色虚线框指示了 伪影问题较为严重的肠道空腔及周边区域。经过修正，肠道空腔边界清晰可见，且周边软 组织中的高频条状伪影得到显著抑制。

（a）未处理图像 （b）预处理降噪后 （c）阴影伪影修正后 （d）条状伪影修正后
图 6 腹腔锥束 CT 运动伪影修正效果（矢状位视图）
Fig.6 Motion artifact correction results for one abdominal CBCT scan (sagittal view)
图 7 是另一位患者腹腔锥束 CT 横断面和矢状位视图中的运动伪影修正效果，红色虚线 框指示了伪影问题较为严重的肠道空腔及周边区域。降噪处理后，皮下脂肪和肌肉之间的 组织边界更清晰。进一步修正后，肠道空腔边界得到部分恢复，且条状伪影得到显著抑制。
3 讨论
经过数字模体仿真实验发现腹腔锥束 CT 影像中运动伪影的重要来源是肠道空腔的蠕 动，包括放缩以及平移两种主要运动形式。通过模体图像和临床 CT 图像分析，运动伪影形

768 CT 理论与应用研究 31 卷

态也包含阴影和条状两种，且均分布在肠道空腔周边。因此，本文设计的运动伪影修正方 法的关键步骤是肠道空腔区域的准确分割，故分割操作是在预处理降噪后进行。但在运动 显著、运动伪影严重的情况下，降噪后通过阈值等方法进行分割时仍可能出现较大偏差， 此时需要在第 1 轮伪影修正操作结束后再从处理后的结果得到更为准确的肠道空腔区域分 割模板，并在预处理降噪图像上重新应用新分割结果进行伪影修正。此外若不同肠道空腔 间隔很近时，不同空腔运动引入的伪影将出现强耦合，也将显著影响空腔区域的分割精度， 此时可先结合计划CT 高质量参考图像非刚性配置到锥束 CT 的方式辅助判断空腔区域的实 际位置、形状和数量，对空腔区域分割结果手动修正后再进行伪影修正。

（a）未处理图像 （b）预处理降噪后 （c）阴影伪影修正后 （d）条状伪影修正后
图 7 腹腔锥束 CT 运动伪影修正效果（横断面和和矢状位视图）
Fig.7 Motion artifact correction results for another abdominal CBCT
scan (transverse and sagittal views)
本文提出的处理流程中针对阴影和条状伪影分别设计了修正操作，这两项操作的顺序 不能颠倒。阴影伪影修正中主要涉及到正弦图插值操作，而条状伪影修正则是空间滤波操 作。若先进行滤波操作会导致图像的高频信息缺失，影响插值操作对象的数据完整性。
正弦图插值操作中归一化是关键步骤，在正弦图插值操作中归一化方案的选取分别进 行了中位数、平均值以及最大值和最小值的测试，其中利用基于中位数的模板图进行归一 化后的插值修正效果最优。中位数方案中归一化正弦图中的数值相对于 1 的分布更为对 称，有助于获得更准确的插值结果。当肠道空腔区域靠近皮肤（人体边界）时，阴影伪影 修正中的关键正弦图插值修复操作中需要进行外插，则插值精度和正弦图修复效果将受到 影响。此时需要对于插值操作所选用的区域进行非对称设计，即不包含皮肤以及皮肤以外 的区域，而在空腔另一侧加入更多的非空腔区域。因此，通过形态学操作中的参数优化调 整可以避免外插问题，保证当肠道空腔区域靠近皮肤边界时的伪影修正效果。
由于运动估计难，肠道蠕动为自发不规则运动等问题， 目前仍缺乏针对腹腔锥束 CT 扫 描中的运动伪影问题的有效修正方法，经过以腹腔锥束 CT 运动伪影修正为主题的文献检索， 2021 年 Peterlik 等[23]应用梯度权重迭代优化方法，对于不同幅度的运动分别将平均 CT 数 误差从 72.1HU 和 136HU 降低为 36.1HU 和 56HU，但难以实现肠道空腔边界的还原以及肠 道周边软组织结构的恢复。本文提出的方法对于肠道周边的伪影修正效果更优，且根据针

6 期 罗辰等：胰腺癌放疗图像引导锥束 CT 运动伪影修正技术研究 769

对不同腹腔锥束 CT 影像数据集分别的定量指标评估，对于软组织区域的平均 CT 数误差修 正能力和 Igor Peterlik 的方法在同一数量级。
在上述讨论基础上，本文提出的处理流程可以满足伪影更严重、肠道空腔区域位形更 复杂的情况下的伪影修正需求，有助于在胰腺癌图像引导放疗任务中提升锥束 CT 影像质量。
4 结论
胰腺癌放疗图像引导锥束 CT 成像通常需要 30 s 以上的扫描时间，呼吸运动和肠道蠕 动叠加到柔软脏器后造成显著的位移和非刚性形变，在锥束 CT 影像中对应以阴影和条状形 态为主的运动伪影问题，且由于肠道蠕动无法进行准确的实时运动监测导致难以实现基于 实时轨迹的伪影修正。相应地，本文建立的锥束 CT 运动伪影修正工作流包括： ① 三维 CT 图像块匹配降噪； ② 建立正向投影模型，获得虚拟投影图； ③ 基于肠道分割和正弦图插 值的阴影伪影修正； ④ 基于频域空间滤波和插值的条状伪影修正； ⑤ 基于图像处理并行 加速的算法效率优化； ⑥ 修正参数优化和效果评估。
本文针对胰腺癌放疗图像引导锥束 CT 运动伪影修正难题，提出一种基于数字图像处理 技术的高效图像域处理方案。经过本文参数优化后的预处理降噪、阴影伪影修正、条状伪 影修正等一系列操作处理，伪影和噪声问题得到全面改善，肠道空腔的边缘和周边软组织 解剖细节得到部分恢复，将有助于提升放疗前锥束 CT 扫描的图像质量，方便胰腺肿瘤的定 位。由于本文提出的伪影修正方案完全在图像域进行图像处理操作，不依赖额外先验知识， 且无需修改成像系统硬件、调整成像协议、或增加剂量，可直接嵌入现有临床扫描工作流 中作为高效后处理方案，适应性和实用性强，将驱动锥束 CT 在放疗科现有摆位验证常规应 用外的拓展。
参考文献
[1] SIEWERDSEN J H, JAFFRAY D A. Cone-beam computed tomography with a flat-panel imager: Magnitude and effects of X-ray scatter[J]. Medical Physics, 2001, 28(2): 220-231.
[2] JAFFRAY D A, SIEWERDSEN J H, WONG J W, et al. Flat-panel cone-beam computed tomography for image-guided radiation therapy[J]. International Journal of Radiation Oncology, Biology, Physics, 2002, 53(5): 1337-1349.
[3] SABEL M. First results with the varian on-board imager (TM)[J]. Strahlentherapie und Onkologie, 2005: 181.
[4] ABBAS H, CHANG B, CHEN Z J. Motion management in gastrointestinal cancers[J]. Journal of Gastrointestinal Oncology, 2014, 5(3): 223-235.
[5] MURPHY M J, MARTIN D, WHYTE R, et al. The effectiveness of breath-holding to stabilize lung and pancreas tumors during radiosurgery[J]. International Journal of Radiation Oncology, Biology, Physics, 2002, 5