ISSN 1004-4140
CN 11-3017/P

基于ASD-POCS框架的高阶TpV图像重建算法

闫慧文, 乔志伟

闫慧文, 乔志伟. 基于ASD-POCS框架的高阶TpV图像重建算法[J]. CT理论与应用研究, 2021, 30(3): 279-289. DOI: 10.15953/j.1004-4140.2021.30.03.01
引用本文: 闫慧文, 乔志伟. 基于ASD-POCS框架的高阶TpV图像重建算法[J]. CT理论与应用研究, 2021, 30(3): 279-289. DOI: 10.15953/j.1004-4140.2021.30.03.01
YAN Huiwen, QIAO Zhiwei. High Order TpV Image Reconstruction Algorithms Based on ASD-POCS Framework[J]. CT Theory and Applications, 2021, 30(3): 279-289. DOI: 10.15953/j.1004-4140.2021.30.03.01
Citation: YAN Huiwen, QIAO Zhiwei. High Order TpV Image Reconstruction Algorithms Based on ASD-POCS Framework[J]. CT Theory and Applications, 2021, 30(3): 279-289. DOI: 10.15953/j.1004-4140.2021.30.03.01

基于ASD-POCS框架的高阶TpV图像重建算法

基金项目: 

国家自然科学基金面上项目(62071281);山西省重点研发计划(201803D421012);山西省留学人员科技活动项目(RSC1622);山西省回国留学人员科研资助项目(2020-008)。

详细信息
    作者简介:

    闫慧文,女,山西大学计算机科学与技术专业硕士研究生,主要从事医学图像重建、图像处理等方面的研究,E-mail:578705307@qq.com;乔志伟*,男,博士,山西大学计算机与信息技术学院教授、博士生导师,主要从事电子顺磁共振成像、图像重建算法、高性能计算等方面的研究,E-mail:zqiao@sxu.edu.cn。

  • 中图分类号: O242;TP391.41

High Order TpV Image Reconstruction Algorithms Based on ASD-POCS Framework

  • 摘要: 总变差(TV)最小化模型目前已广泛应用于图像重建领域,其通过最小化一阶图像梯度大小变换的L1范数实现,能在稀疏投影采集下得到精确的重构。然而,TV模型是基于分段平滑的图像的假设提出的,有时会产生阶梯效应。研究发现,高阶总变差(HOTV)模型可以有效压制阶梯效应,提高重建精度。此外,TpV模型使用Lp范数来逼近L0范数,有望进一步提高稀疏重建能力。鉴于此,本文将HOTV模型与TpV模型结合,提出一种新的高阶TpV (HOTpV)重建模型,采用自适应梯度下降-投影到凸集(ASD-POCS)算法进行求解,分别在理想和有噪声条件下对灰度渐变仿真模体以及真实CT图像仿真模体进行稀疏重建实验。实验结果显示,相比于TV、TpV以及HOTV三种重建模型,HOTpV能得到精度最高的图像。
    Abstract: Total variation(TV) minimization model has been widely used in the field of image reconstruction. It can achieve accurate reconstruction under sparse projection acquisition by minimizing L1 norm of first-order image gradient size transformation. However, the TV model is based on the assumption of segmented smooth image, which sometimes leads to staircase effect. Researches show that the high-order Total Variation(HOTV) model can suppress the staircase effect effectively and improve the reconstruction accuracy. In addition, total p-variation(TpV, 0 < p ≤ 1) model uses Lpnorm to approximate L0 norm, which is expected to further improve the sparse reconstruction ability. In view of this, this paper combines HOTV model with TpV model, a new high-order TpV(HOTpV) reconstruction model is proposed, which is solved by adaptive steepest descent-projection onto convex sets(ASD-POCS) algorithm, sparse reconstruction experiments are carried out on grayscale gradual simulation phantom and real CT image simulation phantom under ideal and noisy conditions. The experimental results show that compared with TV, TpV and HOTV, HOTpV can get the highest accuracy image.
  • 肺动脉CTA作为一种无创血管成像技术,在诊断肺动脉栓塞、观察血管解剖具有较高的敏感度和准确度[1-2],是目前公认的诊断肺动脉栓塞的金标准。其中团注追踪法是肺动脉CTA最有效、最常使用的检查方式[3-4],由于肺动脉不同于体循环,除了受对比剂注射速率、对比剂剂量及对比剂浓度的影响,其还受到个体的身高、体重、肺动脉血流速度、右心泵血功能、呼吸等[5-6]多种因素的影响,想要得到理想的肺动脉CTA图像,关键在于把握准确的扫描时间窗。因此,在团注追踪时如何选择适当触发阈值非常关键。

    本文前瞻性纳入上海市金山区亭林医院肺动脉CTA共112例,观察不同触发阈值下对肺动脉CTA图像质量的影响,旨在选取最佳的触发阈值。

    前瞻性纳入上海市金山区亭林医院2021年12月至2023年4月临床怀疑肺动脉栓塞的患者112例作为研究对象,其中男性49例,女性63例,年龄在37~93岁,平均年龄为64.28岁。根据不同触发阈值随机将患者分配到A、B、C 3组中,A组触发阈值为120 HU共38例、B组触发阈值为200 HU共37例、C组触发阈值为250 HU共37例。3组间性别、年龄、身高、体重差异无统计学意义(表1),可以比较。

    表  1  3组不同触发阈值的一般情况比较
    Table  1.  Comparison of patient demographics across three different trigger thresholds
    项目组别统计检验
    A组B组C组FP
    男/女15/2214/2320/180.890.43
    年龄/岁62.82±12.4865.95±12.8364.81±18.865.060.06
    身高/mm166.07±9.57 163.94±7.61 164.48±11.03 1.240.39
    体重/kg68.08±14.6364.09±13.7560.57±16.671.350.17
    肺栓塞(阳性/阴性)8/297/3016/222.680.07
    下载: 导出CSV 
    | 显示表格

    纳入标准:①临床怀疑肺栓塞患者;②能顺利完成检查者。排除标准:①碘过敏患者;②严重的肾功能衰竭;③严重心衰不能平卧者;④怀孕 8~15周的妇女;⑤各种原因不能配合完成检查者。

    上述患者均通过医院伦理委员会的批准然后告知患者及家属存在的风险因素,并签署知情同意书后方可进行检查。

    联影uCT760-64排螺旋CT及其联影配套工作站,美德瑞达Salient(DC009DW)双筒高压注射器。对比剂选用350 mgI/mL碘海醇注射液(扬子江药业)。患者取仰卧位,双手上举交叉置于前额部,右侧肘静脉埋置18号静脉套管针。

    扫描范围从肺尖扫描至肺底部。注射方案:20 mL生理盐水+45 mL对比剂+30 mL生理盐水,注射速率均为4.5 mL/s[7-8]。ROI置于肺动脉根部,注射造影剂6 s后开始检测扫描,检测时间间隔为1 s,触发阈值根据分组分别为120、200和250 HU,阈值触发后2 s开始启动扫描。扫描螺距均为1.25,层厚5.0 mm,层间距5.0 mm;增强后重建层厚为1.25 mm,重建层间距1.25 mm。在工作站上做出三维VR、MIP及MPR图像。

    客观评判:分别测量上腔静脉、肺动脉主干及左右肺动脉、右肺静脉的CT值作为客观依据进行分析。

    主观评分:采用双盲5分值+附加减分法[2,7],由两位高年资主治医师对3组患者所得到的CT横断面图像、VR、MIP以及MPR图像进行综合评分。5分,肺动脉清晰,与周围组织界限清晰可辨,可以显示肺动脉第5~6级血管属支;4分,肺动脉清晰,与周围组织界限清晰可辨,可以显示肺动脉第4~5级血管属支;3分,肺动脉清晰,与周围组织界限清晰可辨,可以显示肺动脉第3~4级血管属支;2分,肺动脉显示较差,与周围组织对比较差,仅以显示肺动脉第1~2级血管属支;1分,肺动脉显示模糊,与周围组织界限不清,仅显示肺动脉主干。

    附加减分项:上腔静脉大量对比剂潴留产生线束硬化伪影干扰,-1分;肺静脉显影密度高,接近或超过肺动脉,干扰到肺动脉阅读,-1分。

    使用SPSS 19.0软件对所得数据进行统计学分析。采用单因素ANOVA检验,比较3组间的上腔静脉、肺动脉主干及左右肺动脉、右肺静脉CT值;3组间的图像质量主观的综合评分。以P<0.05为具有统计学意义。

    采用Kappa检验评价两名主治医师间对所得图像质量主观评价的一致性,0.0~0.20为一致性极低、0.21~0.40为一致性一般、0.41~0.60为的一致性中等、0.61~0.80为高度一致、0.81~1为完全一致。

    3组上腔静脉、肺动脉主干、右肺动脉、左肺动脉、右肺静脉CT值对比(表2)。采用单因素ANOVA检验3组数据,3组间的肺动脉主干、右肺动脉、左肺动脉CT值比较接近,无统计学差异。

    表  2  3组不同触发阈值的图像质量比较
    Table  2.  comparison of image quality across three different trigger thresholds
    指标组别统计检验
    A组B组C组FP
    肺动脉主干CT值360.34±104.74356.63±107.18344.30±131.030.200.82
    左肺动脉CT值333.53±95.67 349.38±103.52334.86±120.090.230.80
    右肺动脉CT值334.00±90.65 348.84±101.81348.46±121.140.240.79
    上腔静脉CT值701.08±434.61594.09±433.70440.57±263.674.350.02
    右肺静脉CT值184.82±77.27 226.50±66.95 252.05±71.64 8.250.00
    图像主观评分3.82±0.804.47±0.764.05±0.975.170.01
    下载: 导出CSV 
    | 显示表格

    3组间的上腔静脉、右肺静脉CT值比价有统计学差异。上腔静脉随触发阈值的提升而降低,而右肺静脉CT值随触发阈值的升高而下降。

    3组间的图像质量主观评分比较。A组:3.82±0.80、B组:4.47±0.76、C组:4.05±0.97。3组间的图像质量主观评分有统计学差异,B组>C组>A组,以B组触发阈值为200 HU时,评分最高。

    本组112例肺动脉CTA共检出肺动脉栓塞31例,采用单因素ANOVA检验3组间肺栓塞阳性率无统计差异(表2);采用Kappa检验评价两名医师间对所得112例图像质量主观评价的一致性结果为高度一致,Kappa=0.78。

    肺动脉CTA是被公认的诊断肺动脉栓塞的金标准[9]。高质量的肺动脉CTA图像不仅能提高肺栓塞的检出率,同时也能增强医师诊断信心,提升其诊断效能[10-11]。目前肺动脉CTA的检查方式有3种,第1种是经验法,即设置固定时间内启动扫描,其成像质量不稳定,目前已经逐步被淘汰。第2种是小剂量团注测试法,虽然成像质量较高,但其检查操作复杂,同时对比剂使用量多,检查时产生的辐射剂量也高于其他方式[12]。第3种是团注追踪法,其操作方法简单,对比剂使用量小,辐射剂量少,图像质量可靠,是目前临床上最常用的检查方式[1,6,13]。但是,想得到非常理想的的肺动脉CTA图像也并非易事,主要原因是肺循环时间短,而且受到个体因素影响较大,不同患者的肺循环达峰时间均有差异。如果扫描启动时间过早易造成对比剂在肺动脉干显影较淡,或不能充分灌注到5~6级分支血管,同时上腔静脉高浓度对比剂潴留,产生的线束硬化伪影影响右肺动脉干附壁血栓的观察[6,9,14];启动扫描过迟则造成肺静脉显影浓度高,干扰肺动脉的读取;因此选择适当的触发阈值来抓取合理的检查时间窗非常重要[6]

    在选择触发具体阈值时国内外学者均存在较大争议,多数学者认为触发阈值在80~150 HU间最合理[5]。奚佳汇[4]认为使用高于基线20 HU的低触发阈值,即触发阈值在70~80 HU间,可获取高质量的肺动脉CTA图像,其文中也承认所选各组病例并非都能得到理想的图像;曾玉平等[5]认为的触发阈值设置在100 HU时图像质量最佳,但是对于患者出现心影明显增大、右心功能不全、肺动脉高压等情况,应适当提高设定阈值才能获得优质图像;Palm等[3]认为,触发阈值应高出主肺动脉基线100 HU之后开始扫描,即触发阈值在150~160 HU左右,文中也提及检查时会出现上腔静脉造影剂潴留以及肺静脉显影现象。

    综上所述,之前的国内外多数学者的研究都存在一个倾向,就是认为选择较低的触发阈值可以取得较高质量的肺动脉CTA图像[3-5,15]。实际上,上述研究都基于这样一个现实,就是由于受到CT机型条件限制,后延迟时间较长,达到触发阈值还要延迟一段时间才能进行扫描[15]。一般各个机型不同,所设延迟时间也各不相同,国内外报道大多为4~7.5 s不等[3-5,15-16]。奚佳汇[4]认为触发阈值在70~80 HU间,后延迟时间在3 s;曾玉平等[5]认为的触发阈值设置在100 HU,后延迟描时间为9~11 s;Palm等[3]的研究触发阈值在150~160 HU左右,后延迟时间为5 s。

    随着CT技术的不断发展,后延迟时间也不断在缩短,达到目标阈值后可以在较短时间完成全肺动脉容积扫描,这样我们就可以实现以目标阈值做为触发阈值[13],以较短的后延迟时间(2 s)就可以抓取到肺目标阈值,并在对比剂达峰的时间段内完成扫描[14],以取得高质量的肺动脉CTA图像。目前,可以肯定的是国内外许多研究均认为,肺动脉主干及左右肺动脉的CT值需达到200~250 HU就可满足诊断[7,17-18]

    在正常情况下,肺动脉CTA触发阈值达到80~150 HU后,4 s启动扫描,由于肺循环时间在延迟3~4 s后肺动脉CT值是能达到预期目标。但是对比剂经肘静脉注入后到达肺动脉的峰值时间,除受个体的身高、臂长、体重的差异的影响外,还受右心功能以及肺动脉压力等的影响[10,14,19-20],特别在肺动脉主干栓塞、右心功能不全、肺动脉瓣膜病变、肺动脉高压等情况下,血流速度慢,需要时间会更长,对比剂不能快速充盈肺动脉主干,并灌注到5~6级分支血管[5-6,16],同时上腔静脉高浓度对比剂剂潴留产生的线束硬化伪影,影响右肺动脉干附壁血栓的观察,导致漏诊误诊(图1)。

    图  1  触发阈值为120 HU的VR、MIP、MPR及横断位图
    (a)~(d)患者体重65 kg,使用对比剂为45 mL,注射流速为4.5 mL/s,触发阈值为120 HU,肺动脉显示清晰,与周围组织界限清晰可辨,可以显示肺动脉第5-6级血管属支,5分;上腔静脉大量潴留并见少许放射状伪影(红箭头所指),干扰右肺动脉,-1分。综合评价:4分。
    Figure  1.  Trigger VR, MIP, MPR, and cross-sectional bitmaps with a threshold of 120 HU

    而适当调高触发阈值达到200 HU时,2 s后启动扫描,既能保证肺动脉主干CT值的浓度,满足临床诊断,又可以保证对比剂充分灌注到肺动脉5~6级分支,上腔静脉潴留少,虽然肺静脉显影,但是仍然较淡(图2),并不影响肺动脉观察,总体上是满足了临床需求,提高了肺动脉CTA的图像质量。本研究中B组37例病例中,肺动脉主干CT值为(356.63±107.18)HU、右肺动脉CT值为(349.38±103.52)HU、左肺动脉CT值为(348.84±101.81)HU,肺动脉5~6分支显影理想,而上腔静脉平均CT值为594.09,尚不形成线束硬化伪影,右肺静脉CT值较淡显影,为226.50 HU,主观综合评分最高。

    图  2  触发阈值为200 HU的VR、MIP、MPR及横断位图
    (a)~(d)患者体重65 kg,使用对比剂为45 mL,注射流速为4.5 mL/s,触发阈值为200 HU,肺动脉显示清晰,与周围组织界限清晰可辨,可以显示肺动脉第5~6级血管属支,5分。上腔静脉潴留少,未见少许放射状伪影,肺静脉淡,不干扰读片。综合评分:5分。
    Figure  2.  Trigger VR, MIP, MPR, and cross-sectional bitmaps with a threshold of 200 HU

    但是继续提高触发阈值到250 HU,2 s后启动扫描时,虽然对比剂能够充分的灌注到肺动脉5~6级分支,上腔静脉潴留更少,但是肺静脉显影浓度增加,会干扰到肺动脉的观察。本C组右肺静脉的CT值达到(252.05±71.64),对肺动脉末梢细小血栓的判断有干扰(图3)。

    图  3  触发阈值为250 HU的VR、MIP、MPR及横断位图
    (a)~(d)患者体重64 kg,使用对比剂为45 mL,注射流速为4.5 mL/s,触发阈值为250 HU,肺动脉显示清晰肺动脉显示清晰,与周围组织界限清晰可辨,可以显示肺动脉第5~6级血管属支,5分;肺静脉显影密度高(红箭头所指),密度接近肺动脉,干扰肺动脉,-1分。综合评分:4分。
    Figure  3.  Trigger VR, MIP, MPR, and cross-sectional bitmaps with a threshold of 250 HU

    综上所述,肺动脉CTA触发阈值在200 HU,2 s后启动扫描时,既能保证肺动脉主干CT值的浓度,满足临床诊断,又保证对比剂充分地灌注到5~6级分支,上腔静脉潴留少,肺静脉显影淡,肺动脉CTA的图像质量最佳。

    本组研究尚存在一定的不足:①本组研究样本量较小并且为单中心研究,在主观评分上可能存在个体偏好,因此仍需进一步扩大样本量和组织多中心研究,以进一步研究证实本研究的有效性;②由于在本次研究设计中,分组间的触发阈值的跨度较大,在接下来的研究中逐步缩小阈值间的间距,如间距为20 HU,可进一步精确触发阈值的具体数字;③同时本次缺乏个性化设计,特别在肺主干栓塞、右心功能不全者、肺动脉高压等情况下,可借助心脏超声检查肺动脉压力和右心射血分数,来满足不同患者的个性化要求。

  • 期刊类型引用(1)

    1. 杨谨旭,陈颖,雷玉洁,黄云超. 囊腔型肺癌的诊疗进展. 中国肺癌杂志. 2023(10): 774-781 . 百度学术

    其他类型引用(2)

计量
  • 文章访问数:  641
  • HTML全文浏览量:  23
  • PDF下载量:  87
  • 被引次数: 3
出版历程
  • 收稿日期:  2021-01-09
  • 网络出版日期:  2021-11-05

目录

/

返回文章
返回
x 关闭 永久关闭