基于系统建模的低剂量CT重建研究

【摘要】X射线计算机断层成像(CT)以其高的时间分辨率、空间分辨率及对比度分辨率成为现代影像学的杰出代表,在临床诊断和治疗中广泛使用。但过量照射X光可诱发癌症、白血病或其他遗传性疾病,因此CT的剂量安全问题成为业界关注的焦点。如何以最小的代价和最低的X射线剂量获得最佳的CT图像质量,已经成为临床最为迫切的需求。降低X射线剂量的方法有多种,如通过降低管电流和缩短曝光时间来降低X射线球管的毫安秒(mAs)或者减少扫描角向采样数目等直接性的措施来降低X谢线剂量。但与此同时,上述扫描模式下解析重建的图像质量将严重退化,图像中的噪声和伪影影响了临床诊断和应用。近年来,针对低剂量CT图像优质重建的研究方兴未艾,主要集中于对解析重建的低剂量CT图像进行后处理滤波、对低剂量CT投影数据进行恢复,以及包括代数迭代和统计迭代的迭代重建算法。其中,针对迭代重建算法的研究尤其热烈。相对于解析重建算法,迭代重建算法在CT固有的物理系统建模方面具有较大的优势,如可精确模拟系统成像几何,X线多能光谱、线束硬化、散射、噪声等。因此,迭代重建算法能保持或者提高图像空间分辨力的同时抑制伪影和噪声。进一步的,统计迭代重建算法能基于光子统计学建立多种更精确的噪声模型,在降低图像噪声方面有更佳的表现。常规的,在迭代重建过程中,准确建模投影数据的测量方程是得到优质重建的前提和基础,而先验知识的合理利用对于求解的稳定性及获得高精度重建图像具有非常重要的意义。目前,在疾病的临床诊断和治疗过程中,常用到反复的CT扫描,如在灌注CT成像、4D-CT成像、CT图像引导的活体组织穿刺检查以及图像引导放疗等过程中。以图像引导放疗为例,除计划阶段CT扫描外,在整个放射治疗阶段,在分次治疗前需对患者进行CBCT扫描来进行定位,在此情况下,患者所接受的累加辐射剂量较常规CT检查扫描将会很高,反复CT扫描增加了罹患肿瘤的风险。在上述反复CT扫描过程中,对于同一个病人,先后两次CT扫描的图像之间除了几何位置和某些不自主运动引起的不同外,大部分解剖结构是相同的。换句话说,两次扫描获取的CT图像之间存在大量的结构冗余信息。传统的针对低剂量CT的研究方法仅仅引入图像自身的局部邻域的约束作为先验信息,未能考虑同一病人先前扫描获取的图像所能提供的先验信息来导引当前CT的图像的重建。鉴于此,本论文针对基于先前图像的先验信息的合理引入、先验信息引入后算法的优化等问题进行了深入的研究,同时本文对基于平板探测器的CBCT(ConeBeamCT)成像系统的噪声特性进行了进一步实验分析和验证,为后续迭代重建建立了理论和系统模型基础。归纳起来,本文的主要工作有：1)为合理的引入先验图像信息,同时突破传统先验单纯依赖于目标图像局部邻域内的像素灰度信息的约束,本文利用非局部的思想通过在固定大小的搜索窗内检测基于图像块匹配的像素相似性而形成正则化项来实现基于先验图像的先验信息的有效引入。新方法在能够引入先验图像信息的同时弱化了对当前待重建目标图像和先验图像间配准精度的要求。数值仿真,物理体模以及临床数据实验表明,新方法能够提高当前低剂量图像的重建质量,同时不引入伪结构信息。2)针对脑灌注CT成像的两大特点：(1)灌注扫描前的平扫描标准剂量图像的分辨率高且噪声低,可以为低剂量灌注序列增强图像的重建提供极为丰富的形态结构冗余信息；(2)对同一感兴趣层面反复扫描所得的增强序列图像除增强信息外,各图像之间的相对形变量较小；本文提出一种标准剂量平扫描图像导引的低剂量CT脑灌注序列图像重建方法。新方法首先利用利用标准剂量平扫描图像丰富的冗余信息实现基于改进的非局部均值算法的低剂量增强扫描图像的优质恢复,然后对恢复后的图像做去卷积迭代重建。数值仿真体模和病人数据实验表明,该算法能有效抑制图像中噪声,提高图像信噪比,从而使得相应血流动力学参数的计算更加准确。3)为解决先验图像和当前待重建图像之间的结构不匹配问题,本文提出了一种利用当前测量所得投影数据与先验图像配准的方法获取与当前待重建图像相近或类似的先验图像用于基于先验图像的CBCT少角度重建。具体来说,本文算法是通过当前测量投影数据与先验图像在当前成像几何下的前向投影数据的匹配来实现图像域变形场的估计,然后将该变形场作用于先验图像来获取配准的先验图像,用于基于先验图像的重建策略中。特别的,本文采用的重建算法为PICCS(PriorImageCompressConstrainedSensing)算法。本文提出的算法能够有效避免图像空间配准时FDK算法重建的图像中噪声和条形伪影对配准精度的影响。XCAT仿真实验以及临床数据实验结果表明,本文提出的方案能够获取优质先验图像,重建结果优于传统PICCS算法。4)鉴于CBCT成像系统中平板探测器与常规扇形束CT探测器物理设计上的不同,本文实验性的针对瓦里安TrueBeam系统中的平板探测器采集信号进行噪声相关性研究。通过对固定角度下反复测量的的探测数据的分析得出投影数据在探测器单元之间存在噪声相关性。实验结果表明,二维探测器八邻域内的噪声相关性是明显可见的,且一阶邻域内的相关系数明显高于二阶邻域内相关系数。同时,该噪声相关性与扫描所用X线剂量无关。本文将该噪声相关性应用于基于投影数据噪声模型的惩罚加权最小二乘恢复算法中,重建图像结果表明,该噪声相关性的引入能够提高重建图像质量。该实验研究为CBCT图像重建提供了更精准的系统建模。
【作者】张华；
【导师】陈武凡；马建华；
【作者基本信息】南方医科大学，生物医学工程，2014，博士
【关键词】低剂量CT；先验图像；迭代重建；CBCT；噪声相关性；

【参考文献】
[1]张宁宁.降解半纤维素嗜热菌的筛选及耐热木聚糖酶的性质[D].福建农林大学,植物病理学,2010,硕士.
[2]王冬.计算机光盘驱动器高速旋转盘片的解析动力学研究[D].吉林大学,机械设计及理论,2014,硕士.
[3]潘连玉阳.基于CAN总线的矿用自卸车电控系统研究[D].长安大学,机械电子工程,2013,硕士.
[4]张洋洋.增值税“扩围”对我国金融业的影响及对策研究[D].首都经济贸易大学,财政学,2014,硕士.
[5]展秋红.突发公共卫生事件应急管理研究[D].内蒙古大学,公共管理,2014,硕士.
[6]梁俊华.南中国海PAHs的分布、生物泵输出及时间序列研究[D].厦门大学,海洋化学,2014,硕士.
[7]张佳宁.第三人利益合同中受益第三人请求权研究[D].兰州大学,民法,2012,硕士.
[8]苏春雅.扬州市城区典型水体溶解氧分布特性及充氧技术研究[D].扬州大学,市政工程,2011,硕士.
[9]郭培军.中国专利审查流程服务系统的设计与实现[D].北京交通大学,2014.
[10]刘杰.应用蛋白质组学方法分析集胞藻PCC6803对己烷的耐受机理[D].天津大学,制药工程,2013,硕士.
[11]翟婷婷.黄鹤楼卷烟包装的“即时性主题”设计应用研究[D].湖北工业大学,设计艺术学,2014,硕士.
[12]唐武峰,崔鹏高,武守防.安钢100t电弧炉物料平衡及炉料结构的优化[J].河南冶金.2004(05)
[13]杨眉.基于机器视觉的贝类自动识别与定位方法研究[D].大连海洋大学,农业机械化工程,2014,硕士.
[14]胡金丹.用关联理论解读王尔德喜剧中的言语幽默[D].华中师范大学,外国语言学及应用语言学,2014,硕士.
[15]王建军.基于正负刚度并联的立式水泵隔振技术的研究[D].中北大学,工程力学,2013,硕士.
[16]张钰龙.户籍改革的积分制研究[D].安徽大学,社会学,2013,硕士.
[17]王俊.“陪伴式”实务模式的改善与推广：“网格化陪伴”实务模式[D].华中师范大学,社会工作,2014,硕士.
[18]孟晓娜.客户关系管理在我国银行中的应用及研究[D].天津工业大学,2005.
[19]孙超.论民事公益诉讼原告资格多元化[D].吉林大学,法律,2013,硕士.
[20]张群,彭奇志,滕颖.高校馆藏实体资源和虚拟资源的整合与服务[J].情报杂志,2005,07:127-128+131.
[21]邓德传.虚拟机资源分配的非合作博弈标价模型研究[D].杭州电子科技大学,计算机应用与技术,2011,硕士.
[22]韩丽萍.滨州市辖区全民体育健身现状与对策研究[D].山东师范大学,体育教育训练学,2003,硕士.
[23]缪轶.基于事故分析的电气设备制造企业EHS体系优化研究[D].中国地质大学（北京）,安全工程,2013,硕士.
[24]张琳悦.《书记员工作实务》精品课程建设项目管理研究[D].吉林大学,工业工程,2013,硕士.
[25]李国文.大掺量粉煤灰高强混凝土的试验研究[D].西北工业大学,结构工程,2004,硕士.
[26]邵祁峰.保障性住区的空间生产与社会再造[D].南京大学,城市规划与设计,2013,硕士.
[27]白雪.矿用挖掘机提升机构减速器动力学特性研究[D].吉林大学,机械设计及理论,2013,硕士.
[28]刘雄.低温NH_3-SCR脱硝MnO_x–CoO_x复合氧化物催化剂的研究[D].重庆大学,化学工程与技术,2014,硕士.
[29]张学彩.山东省民营高新技术企业创新型文化构建研究[D].济南大学,技术经济及管理,2012,硕士.
[30]司伟.椪柑采后蒂腐病菌生物学特性及其防治研究[D].华中农业大学,园林植物与观赏园艺,2011,硕士.
[31]刘金凤.中美高中力学教材中科学方法教育的比较研究[D].曲阜师范大学,课程与教学论,2013,硕士.
[32]徐晓蕾.海门方言古词语词义研究[D].华中师范大学,汉语言文字学,2014,硕士.
[33]胡祥萍.云制造环境下基于语义的制造资源建模与管理研究[D].北京交通大学,机械设计及理论,2013,硕士.
[34]张伟锋.高校交通智能化管理调研与探索——以浙江工业大学为例[J].中国水运(下半月),2012,04:58-59+245.
[35]金周,温治,司俊龙.烧结热过程数值仿真系统的研究与分析[J].冶金自动化,2008,05:38-41.
[36]赵雁军.同煤大唐热电公司资产管理信息平台设计与实现[D].电子科技大学,软件工程（专业学位）,2012,硕士.
[37]韩磊.论物的损害赔偿方式及实现[D].南京师范大学,法律,2012,硕士.
[38]张金海.广东国际光电采购中心项目投资效益研究[D].湖南大学,工商管理,2013,硕士.
[39]王祎好.谨慎性的“非对称性”视角下财务指标改进研究—从投资者的角度出发[D].兰州商学院,会计学,2014,硕士.
[40]任福平.科通公司仓库管理系统设计方案[D].兰州大学,工商管理（专业学位）,2013,硕士.
[41]李婧.论未成年人犯罪公诉案件办理模式之构建[D].广西师范大学,法律,2013,硕士.
[42]刘海珠.道路交通事故严重程度影响因素分析及预测模型建立[D].吉林大学,交通信息工程及控制,2014,硕士.
[43]朱凌云,吴宝明.医学数据挖掘的技术、方法及应用[J].生物医学工程学杂志,2003,03:559-562.
[44]颜丽.青霉素G酰化酶产生菌发酵条件研究[D].甘肃农业大学,营养与食品卫生,2013,硕士.
[45]段银祥.可控/“活性”自由基聚合法制备嵌段共聚物聚苯乙烯—聚丙交酯（PS-b-PLA）[D].湖南科技大学,化学工艺,2013,硕士.
[46]张丽.新中国的出版政策研究[D].广东海洋大学,马克思主义中国化研究,2013,硕士.
[47]缪新磊.基于超声辅助加工的超声振动系统优化设计及试验[D].扬州大学,机械制造及其自动化,2012,硕士.
[48]苏达.大功率LED封装散热性能的若干问题研究[D].浙江大学,2008.
[49]吴宇.小型移动水质监测系统的研究[D].浙江大学,2013.
[50]杨建虎.聚离子液体/铂纳米粒子催化剂的研制[D].苏州大学,2011.