ISSN 1004-4140
CN 11-3017/P

基于非局部广义总变分正则化的稀疏角度CT重建算法

蒋敏, 陶红伟, 程凯

蒋敏, 陶红伟, 程凯. 基于非局部广义总变分正则化的稀疏角度CT重建算法[J]. CT理论与应用研究(中英文), 2025, 34(1): 129-139. DOI: 10.15953/j.ctta.2023.170.
引用本文: 蒋敏, 陶红伟, 程凯. 基于非局部广义总变分正则化的稀疏角度CT重建算法[J]. CT理论与应用研究(中英文), 2025, 34(1): 129-139. DOI: 10.15953/j.ctta.2023.170.
JIANG M, TAO H W, Cheng K. Sparse View CT Reconstruction Algorithm Based on Non-Local Generalized Total Variation Regularization[J]. CT Theory and Applications, 2025, 34(1): 129-139. DOI: 10.15953/j.ctta.2023.170. (in Chinese).
Citation: JIANG M, TAO H W, Cheng K. Sparse View CT Reconstruction Algorithm Based on Non-Local Generalized Total Variation Regularization[J]. CT Theory and Applications, 2025, 34(1): 129-139. DOI: 10.15953/j.ctta.2023.170. (in Chinese).

基于非局部广义总变分正则化的稀疏角度CT重建算法

基金项目: 河南省科技攻关项目(222102210032);河南省科技攻关项目(232102210014)。
详细信息
    作者简介:

    蒋敏,女,讲师,主要从事CT图像重建与CT应用研究,E-mail:jm_zqu@163.com

    通讯作者:

    陶红伟✉,男,副教授、硕士生导师,主要从事计算机网络安全理论与应用研究,E-mail:tthhww_811@163.com

  • 中图分类号: P 631.3;O 242;TP 391.41

Sparse View CT Reconstruction Algorithm Based on Non-Local Generalized Total Variation Regularization

  • 摘要:

    基于广义总变分(TGV)正则化的CT图像重建算法可以有效克服总变分(TV)正则化的阶梯效应,从而能保护重建图像过渡区域的结构特征。尽管TGV重建方法优于TV重建方法,但它仍然忽略了非局部自相似先验信息在恢复CT图像细节方面的显著作用。为了克服TGV重建方法的上述局限性,本文引入一种非局部广义总变分(NLTGV)正则项,并提出基于NLTGV正则化的稀疏角度CT重建算法。该方法不仅可以利用不同阶的非局部变分信息来保护图像结构特征,而且还可以利用非局部自相似性来恢复重建图像的细节。由于重建模型包含双非光滑项,难以直接求解,因此提出基于凸集投影的优化算法,将其分解为几个简单子问题实现有效求解。仿真和实验结果表明,与其他变分正则化重建方法相比,本文重建方法可以有效提高CT图像重建质量。

    Abstract:

    CT image reconstruction algorithm based on generalized total variation (TGV) can overcome the staircase effect of total variation (TV) regularization, thereby protecting the structural features of the reconstructed image transition region. Although the TGV reconstruction method is superior to the TV reconstruction method, it still ignores the role of non-local self-similar prior information in restoring CT image details. To overcome the aforementioned limitations of TGV reconstruction method, we introduce a non-local TGV (NLTGV) regularization term and propose a sparse view CT reconstruction algorithm based on NLTGV regularization. The proposed method can not only utilize non-local variational information of different orders to protect image structural features but can also utilize non-local self-similarity to restore the details of the reconstructed image. Owing to the inclusion of dual non-smooth terms in the reconstruction model, solving it directly is difficult. Therefore, we proposed an optimization algorithm based on convex set projection, which decomposes the problem into several sub-problems to be solved. The simulation and experimental results show that the proposed NLTGV regularization reconstruction method can effectively improve the quality of reconstructed images compared with other variational reconstruction methods.

  • 滤泡性淋巴瘤(follicular lymphoma,FL)是非霍奇金淋巴瘤(non-Hodgkin’s lymphoma,NHL)的常见类型之一,多发生于淋巴结,也可见于胃肠道、皮肤等结外器官或部位,罕见发生于骨骼肌。结外病灶的诊断对评估FL患者预后具有重要价值。能谱增强CT成像能够通过分析不同组织的能谱曲线确定其同质性,因此能帮助诊断继发性结外淋巴瘤。

    本文报告我院收治的1例以能谱增强CT扫描诊断骨骼肌继发性结外FL的病例。

    患者,男,52岁,因发现左侧腹股沟区一枚硬币大小的肿物2月余于2021年9月于外院就诊。2021年9月6日行腹股沟彩超检查示双侧腹股沟淋巴结肿大,行穿刺活检术,10月11日病理报告示倾向为低级别滤泡性淋巴瘤。同时患者自觉消瘦明显,1月内体重下降5 kg,伴夜间出汗明显,肢软乏力。

    于10月13日行正电子发射计算机断层显像(position emission tomography,PET)/计算机断层扫描(computed tomography,CT)检查发现:①全身多区域淋巴结肿大、脾肿大伴氟代脱氧葡萄糖(fluoro deoxy glucose,FDG)代谢增高,全身多发骨骼FDG代谢增高,均考虑淋巴瘤累及;②肝肿大,肝内稍低密度灶伴边缘FDG代谢轻度增高,淋巴瘤累及不除外。

    10月18日患者于外院住院治疗,并于全麻下行功能性颈部淋巴结清除术,术后予谷胱甘肽、兰索拉唑、盐酸左西替利嗪、氟比洛芬酯综合治疗。术后病理免疫组化:CD3(背景T细胞+),CD20(+),CD19(+),CD10(+),cD21(树突细胞+),CD23(树突细胞+),CD5(背景T细胞+),CyclinD1( − ),Bc12(生发中心41),Bc16(生发中心+),Ki67(生发中心15%),PD-1(散在+),结合苏木精-伊红染色法(hematoxylin-eosin stain, HE stain)可符合为滤泡性淋巴瘤I级。

    2021年11月9日,患者因肢麻乏力于我院住院治疗,行上、下腹部CT平扫示:右侧心膈角、肝门区、腹膜后、腹腔、盆壁多发肿大淋巴结,部分融合;脾脏肿大(图1)。行骨髓穿刺+活检示:骨髓增生极度活跃,骨小梁20%,纤维脂肪5%,造血主质75%,三系均见,以中晚幼成熟为主,B淋巴细胞相对明显增生,呈实性、滤泡状生长,表达BCL-2,全片见巨核细胞10个,网状纤维重度增生(MF-3),结合临床所供病史及实验室检查,倾向滤泡性淋巴瘤累及骨髓。

    图  1  上腹部CT平扫:双侧腰大肌形态、密度无明显差异;腹膜后多发肿大淋巴结;脾肿大
    Figure  1.  Computed tomography (CT) scan of the upper abdomen

    患者临床诊断为滤泡性淋巴瘤Ann Arbor IV期,FL国际预后指数(FL international prognosticindex,FLIPI)3分(高危)。目前IV期FL普遍认为不可治愈,且FL为惰性淋巴瘤,大部分患者病变进展缓慢,相当长时间不接受治疗仍可保持良好的生活质量,故该患者未行放、化疗治疗,在患者要求下于我院行中医保守治疗。

    2022年4月患者因腰痛进行性加重2月,妨碍行走就诊我院。血液检查结果:血红蛋白106 g/L,白细胞计数2.75×109/L(正常范围3.50~9.50×109/L),血小板178×109/L,乳酸脱氢酶205.8 U/L,尿β2-微球蛋白2.32 mg/L(正常范围 0.00~0.30 mg/L)。乙型和丙型肝炎病毒以及人类免疫缺陷病毒血清学检测呈阴性。

    2022年4月19日行腰椎磁共振(magnetic resonance imaging,MRI)平扫(图2),因幽闭恐惧症发作无法耐受长时间扫描,且患者拒绝行PET/CT检查及穿刺活检和手术。为明确诊断和病变范围遂于次日行腰椎能谱增强CT检查。

    图  2  腰椎MRI平扫
    (a)和(b)L1-L2水平右侧腰大肌见一肿块,呈T1WI低、T2WI稍高信号。
    Figure  2.  Lumbar spine MRI plain scan

    采用256排CT扫描仪(Revolution CT,GE healthcare)行平扫及增强扫描。采用宝石能谱成像(gemstone spectral imaging,GSI)扫描模式,扫描参数如表1所示。所用对比剂为非离子型碘对比剂(博莱科信宜,国药准字H20053385),浓度370 mg/mL,用量1 mL/kg,流率2 mL/s,生理盐水20 mL,流率3 mL/s。数据采集完成后,重建带能谱数据信息(GSI DATA)序列,层厚0.625 mm,并将所有序列传入工作站进行后处理分析。

    表  1  能谱增强CT扫描方式及参数
    Table  1.  Spectral computed tomography (CT) scan modes and parameters
    扫描方式及参数名具体扫描方式及参数值扫描方式及参数名具体扫描方式及参数值
      采集体位  仰卧位,头先进   层间距  0.625 mm
      扫描方向、范围  自胸椎下缘向下扫描至骶椎   管电压  GSI模式(瞬时切换80 kVp、140 kVp)
      探测器宽度  64×0.625 mm   管电流  自动毫安模式(200~650 mA,NI:9)
      扫描方法  经验法   矩阵  512×512
      层厚  0.625 mm   机架旋转时间  0.6 s/r
    下载: 导出CSV 
    | 显示表格

    将重建后的薄层图像传送到 GE AW4.7工作站,利用GSI Volume Viewer软件进行能谱的多参数分析。在动脉晚期,右侧腰大肌病变强化程度最适宜观察,且受到的血管影干扰最小,因此选择此期进行后处理分析。于右侧肿大腰大肌、对侧镜面位置正常腰大肌及腹膜后肿大淋巴结区域选取合适的感兴趣区(region of interest,ROI)。ROI面积应尽量一致,且均不大于50 mm2,应避开坏死、钙化、囊变、血管影及边缘区。

    利用能谱曲线工具获取ROI能谱曲线,分别计算动脉期3组能谱曲线斜率K。计算公式:$K=({\rm{CT}} $40 keV$-\,{\rm{CT}}$XkeV$)/(X-40)$,其中X分别为70、100及140 keV。

    腹膜后见多发肿大淋巴结,邻近的右侧L1-2水平腰大肌较对侧正常腰大肌增厚。GSI分析:右侧L1-2水平腰大肌能谱曲线走形与腹膜后肿大淋巴结能谱曲线走形接近,与对侧正常形态腰大肌能谱曲线走形差异较大(图3(a))。根据能谱曲线分别计算3条曲线的斜率K:右侧腰大肌KVP40-70=2.38、KVP40-100=1.5及KVP40-140=0.99;左侧腰大肌KVP40-70=1.41、KVP40-100=0.89及KVP40-140=0.6;腹膜后肿大淋巴结KVP40-70=2.76、KVP40-100=1.76及KVP40-140=1.135(表2)。右侧腰大肌的曲线斜率与腹膜后淋巴结斜率接近,与正常腰大肌斜率差异较大,提示右侧腰大肌被淋巴瘤浸润可能性大。

    图  3  腰椎能谱CT增强图像
    (a)动脉晚期 ROI对应能谱曲线(蓝线代表右侧肿大腰大肌,绿线代表左侧正常腰大肌,紫线代表腹膜后肿大淋巴结)。(b)腰椎 CT平扫:腹膜后多发肿大淋巴结;右侧腰大肌较对侧明显肿胀。(c)和(d)腰椎能谱增强 CT动脉晚期轴位图像:双侧腰大肌及腹膜后肿大淋巴结ROI对应的CT值;右侧肿大腰大肌增强后强化程度略高于左侧正常腰大肌,密度均匀。
    Figure  3.  Lumbar spine spectral computed tomography (CT) image
    表  2  双侧腰大肌及腹膜后肿大淋巴结的KVP40-70KVP40-100KVP40-140
    Table  2.  KVP40-70, KVP40-100 and KVP40-140 of the bilateral psoas major muscle and enlarged retroperitoneal lymph nodes
    项目KVP40-70KVP40-100KVP40-140
    右侧肿大腰大肌 2.381.500.99
    左侧正常腰大肌 1.410.890.60
    后腹膜肿大淋巴结2.761.761.14
    下载: 导出CSV 
    | 显示表格

    考虑到中医保守治疗效果不佳,为进一步治疗,患者目前于外院诊疗,病情稳定。

    FL是NHL中的一种常见类型,是淋巴滤泡生发中心B细胞来源的一种惰性B细胞淋巴增殖性疾病,以弥漫性淋巴结肿大、骨髓受累和脾肿大为特征,常见于中、老年人,女性稍多于男性。WHO统计FL约占NHL的22.1%,但在中国人群中发生率较低,只占NHL的10% 左右。

    FL主要累及淋巴结,结外受累较少见。肿瘤区域个数的增加可提高FL患者的FLIPI指数进而增加患者的危险度分级,因此结外淋巴瘤的诊断对评估FL患者的病情及预后具有重要价值。累及骨骼肌的淋巴瘤非常少见,仅占霍奇金淋巴瘤(Hodgkin lymphoma,HL)病例的 0.3%、NHL病例的1.5%[1],多数病理类型为弥漫大B细胞性淋巴瘤(diffuse large B cell lymphoma,DLBCL),而原发性骨骼肌FL的报道极为罕见,且没有关于其发病率的流行病学数据[2]

    骨骼肌淋巴瘤通常通过3种途径发生:①作为一种原发性结外疾病发生;②通过血源性或淋巴途径传播;③通过邻近器官浸润,如骨骼或淋巴结。原发性骨骼肌淋巴瘤发生在腿部损伤后、靠近针注射部位和同性恋男性直肠已有报道[3]

    由于下肢易受损伤,这可能提示原发性骨骼肌淋巴瘤可能和机械刺激有关。继发性骨骼肌淋巴瘤最常见的原因是从邻近淋巴结或其他原发病灶(如骨)扩散或转移而来[1],这与本病例的影像表现相符。

    FL临床表现多为受累肌肉出现进行性增大的软组织肿块、肿胀、疼痛、发热、出汗和体重减轻,最常见于下肢,占所有报告病例的50%[3]。本病例在确诊半年余后出现进行性加重的腰痛,严重影响生活质量,影像学检查提示右侧腰大肌新发的软组织肿块。

    能谱增强CT已广泛应用于淋巴瘤的诊断中,目前研究多集中于与鼻咽癌、纵膈胸腺瘤、结节病、肠道肿瘤和转移性淋巴结等疾病的鉴别[4-9]。鉴于骨骼肌淋巴瘤发生的罕见性,通过能谱增强CT诊断的经验较缺乏,只在超声、平扫CT、MRI及PET/CT的相关研究中有零星报道[2,10-15]

    CT可通过病灶的形态、大小、边界、密度及强化方式来诊断骨骼肌淋巴瘤,表现为受累肌肉肿胀,边界清晰或不清晰,与正常肌肉相比呈均质的等或高或低密度,增强后病变肌肉呈轻微至明显不同程度的强化。作为混合能量图像,普通CT的CT值不能准确反映病灶实际的CT值大小。能谱增强CT成像通过快速切换能量并进行数据采集,可获得不同组织在40~140 keV各单能量图像下的CT值[16]。通过对比能谱曲线的形态及斜率可确定病灶的来源、鉴别良恶性。

    本例患者患有幽闭恐惧症并拒绝再次穿刺、手术及PET/CT检查,遂行能谱增强CT示右侧腰大肌与腹膜后肿大淋巴结能谱曲线走行及形态相近,虽未取得腰大肌病理结果验证,结合既往病史、病理与影像学检查考虑为骨骼肌淋巴瘤浸润的可能性最大。

    迄今为止,大多数FL患者仍然不可治愈。FL的治疗应充分考虑患者耐受情况、肿瘤负荷及复发风险而采取个体化治疗。目前各指南推荐的一线治疗方案均为CD20单抗+化疗;对于晚期(Ⅲ 期、Ⅳ期)患者根据美国国立综合癌症网络(national comprehensive cancer nelwork,NCCN)指南推荐如果低肿瘤负荷且没有症状,可采取观察和等待策略;若存在肿瘤高负荷表现、肿瘤进展、受累器官功能受损等治疗指征时,建议开始治疗[17]

    本病例为 Ⅳ 期,FLIPI 3分(高危)患者,前期临床症状轻微,要求中医保守治疗;后期出现严重腰痛,考虑淋巴瘤浸润骨骼肌可能大,提示病变进展,遂于外院就诊进一步治疗。

    骨骼肌淋巴瘤的诊断需结合临床病史、实验室检查、影像学检查及病理组织学检查结果。当骨骼肌出现进行性增大的软组织肿块,伴疼痛、发热、短期内消瘦,应怀疑本病,尤其在淋巴瘤患者中,应怀疑继发性骨骼肌淋巴瘤。影像学特征具非特异性,需与其他疾病(如原发性软组织肉瘤、其他原发软组织恶性肿瘤、创伤或肌炎等)鉴别。

    (1)软组织肉瘤:包括平滑肌肉瘤、横纹肌肉瘤、脂肪肉瘤、滑膜肉瘤等,表现为受累肌肉的形态和轮廓改变伴不规则软组织肿块,肿块平扫及增强的密度和信号不均匀,实性部分大多强化明显,易出现坏死、液化,瘤周水肿较淋巴瘤明显。

    (2)外周型原始神经外胚瘤:与骨骼肌淋巴瘤同属小细胞肿瘤,好发部位相似,但外周型原始神经外胚瘤多见于10~20岁儿童或青少年,临床发病率低,相对罕见,恶性程度高,易复发转移,预后相对较差,对放化疗敏感。在影像学上肿瘤表现为大软组织肿块,浸润性、跨越性、快速生长,血供丰富,明显不均匀强化,其形态较淋巴瘤更不规则,边界更不清楚,囊变坏死更显著,强化更明显且不均匀。

    (3)神经纤维瘤病:是一类常染色体显性遗传性疾病,可分为3型:Ⅰ型神经纤维瘤病、Ⅱ型神经纤维瘤病和 Ⅲ 型施万细胞瘤病。此类疾病表型差异性大,以皮肤病变、周围神经系统病变和中枢神经系统肿瘤为主,引起多发的、渐进性的损害。神经纤维瘤多发生于肌间隙,与肌肉分界清,邻近肌肉受压变形但密度、信号大多正常;少部分发生于肌肉内,病灶形态不规则,受累肌肉多失去正常形态,多肌肉受累时,肌间隙模糊。

    (4)炎症:起病急,临床症状明显,表现为局部红肿热痛,部分患者全身症状明显,抗炎治疗有效。CT或MRI表现为受累肌肉明显肿胀,肌间隙模糊,筋膜及皮下水肿较明显。

    本文报道了1例以能谱增强CT成像帮助诊断继发性骨骼肌FL的病例,在MRI、PET/CT及穿刺活检无法进行时,能谱增强CT可作为辅助诊断的潜在新方法。通过重建及比较不同组织的能谱衰减曲线,可初步确定病变来源,为广大医师在临床诊疗过程中提供诊断新思路、新方法。

  • 图  1   搜索非局部相似图像块原理

    Figure  1.   Schematic illustrating the searching of non-local similarity patches

    图  2   仿真数据重建结果

    Figure  2.   Reconstruction results for simulated data

    图  3   仿真数据重建图像ROI放大图

    Figure  3.   Zoom-in display of ROI of the images reconstructed from the simulated data

    图  4   仿真数据重建图像的RMSE随迭代次数变化曲线

    Figure  4.   RMSE of the image reconstructed from simulated data versus iterations

    图  5   实际投影数据重建结果

    Figure  5.   Reconstruction results for real projection data

    图  6   实际数据重建图像ROI放大图

    Figure  6.   Zoom-in display of ROI of the images reconstructed from the real data

    图  7   实际数据重建图像的RMSE随迭代次数变化曲线

    Figure  7.   RMSE of the image reconstructed from real data versus iterations

    表  1   不同变分正则项特性的比较

    Table  1   Comparison of characteristics with different variation regularization

    算法 TV TGV NLTV NLTGV
    邻域   局部 局部 非局部 非局部
    变分阶次 一阶 高阶 一阶  高阶 
    下载: 导出CSV

    表  2   基于NLTGV正则化的稀疏角度CT图像重建算法步骤

    Table  2   Procedures of a sparse-view CT image reconstruction algorithm based on NLTGV regularization

    算法1 基于NLTGV的稀疏角度CT重建算法
    (1) 初始化: 给定初值$ {{{\boldsymbol u}}^0}, $$ {t_{\max }} ,$令初始化: 给定初值
      $ {{{\boldsymbol u}}^0} ,$设置$ {t_{\max }} $,$ {q_{\max }} $,$ {k_{\max }} $,令$ k = 1 $
    (2) while (不满足停止准则) do
    (3) 令$ {{{\boldsymbol f}}^0} = {{{\boldsymbol u}}^{k - 1}} $,t=1
    (4) while(t$ \leq $tmax) do
    (5) for j=1,2,···,N
    (6) $ f_j^t = f_j^{t - 1} + \displaystyle\frac{\lambda }{{{A_{ + ,j}}}}\sum\limits_{i = 1}^M {\frac{{{a_{i,j}}}}{{{A_{i, + }}}}} ({p_i} - {\bar p_i}) $
    (7) $ {A_{i, + }} = \displaystyle\sum\limits_{j = 1}^N {{a_{i,j}}} ,\ i = 1,2,\cdots,M $
    (8) $ {A_{ + ,j}} = \displaystyle\sum\limits_{i = 1}^M {{a_{i,j}}} ,\ j = 1,2,\cdots,N $
    (9)$ {\bar p_i} = {A_i}{f^{t - 1}} $
    (10)$ t = t + 1 $
    (11) end for loop
    (12) end while loop
    (13) 非负约束,得到$ {{{\boldsymbol u}}^{{\text{pos}}}} $
    (14) 令$ q = 1, $ b=0, d=0
    (15) while($ q \le {q^{\max }} $) do
    (16) 共轭梯度法求解u, g
    (17) 根据公式(18)与(19)求解sr子问题
    (18) 根据公式(16)与(17)更新变量bd
    (19) $ q = q + 1 $
    (20) end while loop
    (21) $ {{{\boldsymbol u}}^k} = {{{\boldsymbol u}}^{{q_{\max }}}} $
    (22) $ k = k + 1 $
    (23) end while loop
    (24) 输出: 重建图像$ {{{\boldsymbol u}}^*} $
    下载: 导出CSV

    表  3   仿真数据重建图像的PSNR和SSIM结果

    Table  3   PSNR and SSIM for reconstructing images from simulated data

    采样角度数目 指标 算法
    TV TGV NLTV NLTGV
    ROI1 ROI2 ROI1 ROI2 ROI1 ROI2 ROI1 ROI2
    50 PSNR 26.16 25.98 26.19 26.12 26.15 26.18 26.59 26.83
    SSIM×10-2 80.02 79.97 80.37 80.52 78.73 79.71 80.48 81.78
    70 PSNR 26.80 26.10 26.90 25.84 27.38 26.98 28.28 27.82
    SSIM×10-2 81.96 80.90 79.78 80.16 83.05 83.53 85.10 84.93
    90 PSNR 27.07 26.20 27.19 26.10 27.63 26.51 28.51 28.56
    SSIM×10-2 81.64 82.95 82.37 81.80 82.62 83.34 85.11 87.29
    110 PSNR 27.07 27.11 27.77 27.00 27.86 27.59 28.71 28.78
    SSIM×10-2 84.24 84.74 82.71 84.03 83.20 84.72 85.47 87.36
    下载: 导出CSV

    表  4   真实数据重建图像的PSNR和SSIM结果

    Table  4   PSNR and SSIM for reconstructing images from real data

    采样角度数目 指标 算法
    TV TGV NLTV NLTGV
    ROI1 ROI2 ROI1 ROI2 ROI1 ROI2 ROI1 ROI2
    50 PSNR 25.38 26.17 25.85 26.45 26.18 26.77 26.83 26.98
    SSIM×10-2 79.97 79.98 79.52 81.12 79.71 82.74 81.78 82.70
    70 PSNR 26.10 27.02 26.34 27.02 26.98 27.39 27.82 28.94
    SSIM×10-2 80.90 80.00 80.16 83.67 83.53 86.91 84.93 88.18
    90 PSNR 26.20 27.16 26.30 27.32 26.51 27.82 28.36 29.21
    SSIM×10-2 81.95 84.68 82.80 84.27 83.34 84.21 87.29 89.04
    110 PSNR 27.11 27.51 27.00 27.54 27.19 28.64 28.78 29.67
    SSIM×10-2 83.74 84.42 84.03 84.67 84.72 86.22 87.36 89.56
    下载: 导出CSV
  • [1]

    PARCERO E, FLORES L, SÁNCHEZ M G, et al. Impact of view reduction in CT on radiation dose for patients[J]. Radiation Physics and Chemistry, 2017, 137(8): 173-175.

    [2]

    CANDES E J, ROMBERG J, TAO T. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information[J]. IEEE Transactions on Information Theory, 2006, 52(2): 489-509. DOI: 10.1109/TIT.2005.862083.

    [3]

    SIDKY E Y, KAO C M, PAN X. Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT[J]. Journal of X Ray Science & Technology, 2009, 14(2): 119-139.

    [4]

    SIDKY E Y, PAN X. Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization[J]. Physics in Medicine & Biology, 2008, 53(17): 4777-4807.

    [5]

    RUDIN L I, OSHER S, FATEMI E. Nonlinear total variation based noise removal algorithms[J]. Physica D: Nonlinear Phenomena, 1992, 60(1/4): 259-268.

    [6]

    GORDON R, BENDER R, HERMAN G T. Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and X-ray photography[J]. Journal of Theoretical Biology, 1970, 29(3): 471-481. DOI: 10.1016/0022-5193(70)90109-8.

    [7]

    LANGE K, CARSON R E. EM reconstruction algorithms for emission and transmission tomography[J]. Journal of Computer Assisted Tomography, 1984, 8(2): 306-316.

    [8]

    BREDIES, K, KUNISCH, K, POCK, T. Total genera-lized variation[J]. SIAM Journal on Imaging Sciences, 2010, 3(3): 492-526. DOI: 10.1137/090769521.

    [9]

    NIU S, GAO Y, BIAN Z, et al. Sparse-view X-ray CT reconstruction via total generalized variation regularization[J]. Physics in Medicine & Biology, 2014, 59(12): 2997-3017.

    [10]

    SUN T, SUN N, WANG J, et al. Iterative CBCT reconstruction using Hessian penalty[J]. Physics in Medicine and Biology, 2015, 60(5): 1965-1987. DOI: 10.1088/0031-9155/60/5/1965.

    [11]

    ZHANG H M, WANG L Y, YAN B, et al. Constrained total generalized p-variation minimization for few-view X-ray computed tomography image reconstruction[J]. Plos One, 2016, 11(2): e0149899.

    [12]

    XI Y R, QIAO Z R, WANG W J, et al. Study of CT image reconstruction algorithm based on high order total variation[J]. Optik, 2020, 204: 163814. DOI: 10.1016/j.ijleo.2019.163814.

    [13] 闫慧文, 乔志伟. 基于 ASD-POCS 框架的高阶 TpV 图像重建算法[J]. CT 理论与应用研究, 2021, 30(3): 279-289. DOI: 10.15953/j.1004-4140.2021.30.03.01.

    YAN H W, QIAO Z W. 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. (in Chinese).

    [14]

    XI Y, ZHOU P, YU H, et al. Adaptive-weighted high order TV algorithm for sparse-view CT reconstruction[J]. Medical Physics, 2023, 50: 5568-5584. DOI: 10.1002/mp.16371.

    [15]

    HUANG J, MA J, LIU N, et al. Sparse angular CT reconstruction using non-local means based iterative-correction POCS[J]. Computers in Biology and Medicine, 2011, 41(4): 195-205. DOI: 10.1016/j.compbiomed.2011.01.009.

    [16]

    KIM H, CHEN J, WANG A, et al. Non-local total-variation (NLTV) minimization combined with reweighted L1-norm for compressed sensing CT reconstruction[J]. Physics in Medicine & Biology, 2016, 61(18): 6878.

    [17]

    KIM K, EL F G, LI Q. Low-dose CT reconstruction using spatially encoded nonlocal penalty[J]. Physics in Medicine & Biology, 2018, 63(3): 035045.

    [18]

    GILBOA G, OSHER S. Nonlocal operators with applications to image processing[J]. Multiscale Modeling & Simulation, 2009, 7(3): 1005-1028.

    [19]

    RANFTL R, BREDIES K, POCK T. Non-local total generalized variation for optical flow estimation[C]//ECCV, 2014, (6): 439-454.

    [20]

    TOMASI C, MANDUCHI R. Bilateral filtering for gray and color images[C]//ICCV, 1998: 839-846.

    [21]

    JIANG M, WANG G. Convergence of the simultaneous algebraic reconstruction technique (SART)[J]. IEEE Transactions on Image Processing, 2003, 12(8): 957-961. DOI: 10.1109/TIP.2003.815295.

    [22]

    GOLDSTEIN T, OSHER S. The split Bregman method for L1-regularized problems[J]. SIAM Journal on Imaging Sciences, 2009, 2(2): 323-343. DOI: 10.1137/080725891.

    [23]

    WANG Z, BOVIK A C, SHEIKH H R, et al. Image quality assessment: From error visibility to structural similarity[J]. IEEE Transactions on Image Processing, 2004, 13(4): 600-612. DOI: 10.1109/TIP.2003.819861.

    [24]

    WANG J, LI T, XING L. Iterative image reconstruction for CBCT using edge-preserving prior[J]. Medical Physics, 2009, 36(1): 252-260. DOI: 10.1118/1.3036112.

  • 期刊类型引用(0)

    其他类型引用(1)

图(7)  /  表(4)
计量
  • 文章访问数:  207
  • HTML全文浏览量:  84
  • PDF下载量:  59
  • 被引次数: 1
出版历程
  • 收稿日期:  2023-08-29
  • 修回日期:  2023-12-10
  • 录用日期:  2023-12-12
  • 网络出版日期:  2024-06-03
  • 刊出日期:  2025-01-04

目录

/

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