Shallow Profile Data Denoising Method Based on Improved Cycle-consistent Generative Adversarial Network
-
摘要: 为解决浅地层剖面数据噪声多、分辨率低问题,本文将环形生成对抗网络的方法应用于浅地层剖面资料的去噪,实现智能去噪。首先,选择具有特殊对称生成对抗网络循环机制以及“循环一致性”损失的环形生成对抗性网络,并对其进行结构改进,提升网络学习和训练的性能。然后,基于优化的浅地层剖面样本集训练网络,实现对于浅地层剖面数据随机噪声的去除,提升数据的信噪比。通过对实验和实际资料的试算,以及与传统带通滤波方法的对比,验证本文方法对浅地层剖面数据去噪的有效性和适应性。Abstract: This study applied the cycle-consistent generative adversarial network method to the denoising of shallow profile data to realize intelligent denoising. This could help resolve the problem of noise and low resolution of shallow profile data. To do this, the cycle generative adversarial network with special symmetric generation countermeasure network cycle mechanism and "cycle consistency loss" was selected. We improved the performance of the network learning and training by optimizing the network structure. Next, based on the optimized shallow profile sample set training network, random noise was removed from the shallow profile data and the signal-to-noise ratio of the data was improved. The effectiveness and adaptability of this method for denoising shallow profile data were verified by trial calculations of experimental and actual data and by comparison with the traditional band-pass filtering method.
-
Keywords:
- artificial intelligence /
- shallow profile data /
- random noise /
- data denoising
-
肝细胞癌(hepatocellular carcinoma,HCC)是具有较高发病率及死亡率的肝脏恶性肿瘤,具有起病隐匿、病情进展迅速及生存率低等特点,其致死率较高,居所有恶性肿瘤的第3位[1-2]。现阶段经导管动脉化疗栓塞(transcatheter arterial chemoembolization,TACE)是临床治疗HCC的重要手段,可以阻断HCC的供血动脉,引起癌组织缩小、坏死,最终发挥消除肿瘤的作用[3-4]。
有研究报道,部分患者在TACE治疗后病灶未完全坏死,加之肿瘤血管再形成,存在较高的肿瘤复发风险[5]。因此,采取安全、快速的方法准确评估HCC患者TACE的治疗效果对指导临床后续治疗具有重要的价值。血清肿瘤标志物监测、计算机断层扫描(computed tomography,CT)及磁共振是临床评估TACE治疗效果的常用方法,但均存在一定不足[6]。肝脏CT动脉增强分数(arterial enhancement fraction,AEF)彩图是近年来逐渐被了解的一种新型方法,有研究发现,AEF在评估肝脏疾病患者治疗效果等方面具有较大优势[7]。基于此,本研究探讨AEF与肝癌TACE治疗效果的相关性,为临床诊疗提供参考。
1. 资料与方法
1.1 一般资料
对我院2019年1月至2020年1月43例HCC患者临床资料进行回顾性分析。纳入标准:①满足《原发性肝癌诊疗规范》[8]中诊断标准,并经病理学检查确诊;②首次确诊,并接受TACE治疗;③术前及术后 1个月均行CT扫描检查;④肿瘤为巨块型或者结节型;⑤肝功能 Child-Pugh分级为A-B级;⑥无肝外转移病灶;⑦患者签署知情同意书。
排除标准:①肝癌复发者;②近期接受放疗或肝脏手术切除者;③存在其他恶性肿瘤者;④合并凝血功能障碍、感染性疾病及免疫功能异常者;⑤心、肺、肾等脏器严重损伤者;⑥临床资料不全者。
43例HCC患者中男性31例,女性12例;年龄36~76岁,平均(53.98±7.69)岁;肿瘤大小3~12 cm,平均(6.16±1.46)cm;肝功能Child-Pugh分级:A级19例,B级24例。
1.2 方法
1.2.1 TACE治疗
患者取平卧位,常规消毒双侧腹股沟区皮肤后铺巾,采用5 mL 2% 利多卡因进行局部麻醉,随后作一5 mm切口,对右侧股动脉进行穿刺,置入动脉鞘,并将导管置于肝总动脉、腹腔干、肝左/右动脉,行DSA,显示肿瘤部位及大小,按方案化疗栓塞,结束后拨出动脉鞘及导管,加压包扎。药物:超液化碘油、羟基喜数碱、表阿霉素。
1.2.2 CT扫描
仪器选择西门子128层螺旋CT扫描仪(西门子双源Definition Flash CT)。检查前,指导患者进行呼吸训练,嘱患者扫描时按口令,在吸气末屏住呼吸。CT平扫及增强扫描动脉期、门脉期呼吸幅度保持一致,在吸气末屏住气扫描。动脉期、门脉期数据采集时间分别为注药后30 s、注药后60 s。入选患者扫描:仰卧位,双手抱头,头先进。在吸气末从膈顶至肝底行全肝扫描(一定把肝组织包全)。
扫描参数:螺旋采集模式,2×128×0.6 mm。吸气相:管电压100 kV,管电流50 mA;呼气相:管电压100 kV,管电流50 mA。机架旋转时间0.33 s/周,矩阵512×512,重建层厚1 mm,重建间隔0.6 mm。
1.2.3 CT数据处理
分别选取TACE治疗前、治疗后1月CT平扫、动脉期和门静脉期图像导入西门子Syngo.VIA后处理工作站,选取MM肿瘤学卡片,生成AEF彩图。
由两名具有10年以上腹部影像诊断经验的放射科医师取得一致意见后选取多个感兴趣区(region of interest,ROI),选取肝肿瘤实质部分,尽量避开囊变、坏死及周围血管区,取病变区域所在相同层面作为ROI进而计算ROI内AEF分布,尽量保证TACE治疗前后病灶的ROI保持一致。
1.2.4 甲胎蛋白水平检测
在治疗前及治疗后1个月所有患者于清晨空腹采集肘静脉血3 mL置于干燥试管中,室温下静置1 h后离心处理,离心时间为5 min,离心速率为3000 r/min。取上清液保存在 -20℃冰备用,采用电化学发光法检测甲胎蛋白(alpha-fetoprotein,AFP)水平。
1.3 评价标准
根据术后CT影像结果并结合临床指标按照mRECIST标准[9]进行疗效评价:①完全缓解(complete remission,CR):所有目标病灶动脉期增强显影均消失;②部分缓解(partial remission,PR):目标病灶(动脉期增强显影)的直径总和缩小≥30%;③病情进展(progressive Disease,PD):目标病灶(动脉期增强显影)的直径总和增加≥20% 或出现新病灶;④病情稳定(stable disease,SD):缩小未达PR或增加未达PD。本研究根据疗效标准,将CR及PR视为有效组(n=28),将PD及SD视为无效组(n=15)。
1.4 观察指标
①记录并比较两组 AEF、AFP、肿瘤大小;②探究 AEF、AFP与肿瘤大小变化的相关性;③分析肝脏动态增强扫描联合 AEF彩图对TACE疗效的评估价值。
1.5 统计学方法
采用统计学软件SPSS 21.0处理分析数据。计数资料以例(%)表示,行χ2检验;计量资料以均数±标准差
$ (\bar x \pm s) $ 表示,两组间比较采用t检验;采用Pearson相关性系数探究AEF、AFP与肿瘤大小变化的相关性;采用受试者工作特征曲线(receiver operating characteristic curve,ROC)及曲线下面积(area under curve,AUC)分析肝脏动态增强扫描联合AEF彩图对TACE疗效的评估价值。以P<0.05为差异具有统计学意义。2. 结果
2.1 两组治疗AEF、AFP、肿瘤大小比较
治疗前,有效组与无效组相比,AEF、AFP、肿瘤大小三者差异均无统计学意义。治疗后,有效组AEF小于无效组,有效组AFP小于无效组、有效组肿瘤大小小于无效组(表1)。
表 1 两组治疗AEF、AFP、肿瘤大小比较($ \bar x \pm s $ )Table 1. Comparison of the AEF, AFP level and tumor size between the two groups ($ \bar x \pm s $ )组别 AEF AFP/(ng/L) 肿瘤大小/cm 治疗前 治疗后 治疗前 治疗后 治疗前 治疗后 有效组 0.51±0.13 0.34±0.05* 385.12±10.29 84.89±8.39* 6.28±1.05 3.01±0.41* 无效组 0.52±0.09 0.48±0.07 382.87±21.54 373.25±17.82 6.05±1.13 5.96±1.05 t 0.265 7.594 0.420 72.436 0.667 13.209 P 0.792 0.000 0.676 0.000 0.509 0.000 注:与本组治疗前,*-P<0.05。 2.2 AEF、AFP与肿瘤大小的相关性
Pearson相关性分析结果显示(表2),经TACE治疗后AEF与肿瘤大小呈正相关(r=0.537),AFP与肿瘤大小呈正相关(r=0.649)。
表 2 AEF、AFP与肿瘤大小的相关性Table 2. Correlation between the AEF, AFP level and tumor size指标 肿瘤大小 P AEF 0.537 0.002 AFP 0.649 0.000 2.3 肝脏动态增强扫描联合AEF对TACE疗效的评估价值
ROC曲线分析结果显示,肝脏动态增强扫描联合AEF评估TACE疗效的AUC为0.902,明显高于肝脏动态增强扫描及AEF单独评估的0.793和0.771(图1和表3)。
表 3 肝脏动态增强扫描联合AEF对TACE疗效的评估价值Table 3. Evaluation value of dynamic enhanced liver scanning combined with AEF for TACE efficacy项目 AUC P 95% CI 肝脏动态增强扫描 0.793 <0.05 0.642~0.901 AEF 0.771 <0.05 0.618~0.885 肝脏动态增强扫描联合AEF 0.902 <0.05 0.773~0.972 2.4 典型病例
男性,58岁,肝硬化30余年,发现肝占位3天,AFP 487.52 ng/mL。肝右叶巨块型肝癌,A为动脉期图像,呈明显高强化,P为门静脉期图像,原动脉期高强化区域强化退出,与肝实质相比呈低密度影,U为CT平扫图像,与肝实质相比为稍低密度区,境界模糊,AEF彩图肿块大部分为红色区域,能够直观显示肿块、并提示肿块为丰富肝动脉供血(图2)。
3. 讨论
HCC是原发性肝癌最为常见的类型,病死率居于我国恶性肿瘤的第3位[1]。TACE是一种有效的非外科手术姑息性治疗方法,将不同药物颗粒经动脉注入肿瘤血管,使肿瘤缺氧缺血而坏死[3]。但TACE治疗后若病灶并未完全坏死,会导致恶性肿瘤再次复发[10-11]。现阶段临床主要通过CT、磁共振等影像学检查评估TACE治疗后肝脏肿瘤的存活情况,但磁共振具有检查用时较长、费用较高等缺点,且对患者呼吸配合要求较高。因此,临床仍以CT为主要手段评估TACE治疗效果。
临床研究发现,HCC患者接受TACE治疗后血供明显减少,可使AEF明显降低[12]。本研究结果显示,TACE治疗后有效组AEF低于TACE治疗前及无效组,与既往研究结果一致,提示AEF可作为临床评估TACE治疗效果的参考指标。肝脏是门静脉及肝动脉的双重供血脏器,在肝脏功能正常情况下,主要以门静脉供血为主,故正常肝实质在门静脉期密度达最高。但HCC患者肿瘤组织主要以肝动脉供血为主,因此在HCC患者肝实质在门静脉期呈低密度。
AEF=(动脉期CT值-平扫期CT值)/(门静脉期CT值 - 平扫期CT值) (1)
AEF是感兴趣区域在肝动脉期内绝对强化增量与门静脉期内绝对强化增量的比值。
HCC是富血供病变,且由肝动脉供血,在动态增强CT扫描中动脉期迅速强化呈稍高密度,门静脉期和平衡期快速退出呈稍低密度,可见,TACE治疗前HCC患者的AEF相对较高。HCC患者在TACE治疗后肿瘤内部和周围动脉发生狭窄,甚至闭塞,使得肝脏固有动脉流速明显降低,对TACE治疗效果良好的肿瘤而言,肿瘤内部血流减少与肝脏固有动脉流有关,从而导致AEF较TACE治疗前明显减少。
蒋长斌等[13]研究发现,AFP在多数HCC患者血清中呈高水平状态,是诊断HCC的常用肿瘤标志物。本研究中有效组治疗后血清AFP水平明显低于治疗前、无效组,与其研究结果一致,表明AFP可作为临床评估HCC患者TACE治疗效果的参考指标。本研究采用Pearson系数分析AEF、AFP与肿瘤大小的相关性发现,HCC患者TACE治疗后AEF、AFP均与肿瘤大小呈正相关,进一步证实了AEF、AFP可用于评估HCC患者TACE的治疗效果。
王曦等[14]认为,相比于常规CT扫描,肝脏动态增强扫描的速度较快、处理技术强大,可以清楚反映肿瘤部位、大小及内部碘油沉积情况,有利于临床判断TACE治疗后肿瘤残留情况。有报道显示,肝脏动态增强扫描评估TACE治疗后肿瘤残留情况的准确度为85.30%[15]。本研究结果显示,肝脏动态增强扫评估TACE疗效的准确度仅为79.30%,低于上述研究结果,这可能是碘油影响对肿瘤血管的遮挡所致。
本研究中AEF评估TACE治疗效果的AUC为0.771,分析原因为TACE治疗后部分患者病灶大小不超过1 cm,此时AEF彩图可能出现失真,易与小血管混淆,造成漏诊;此外,扫描时患者的呼吸、心脏运动也会产生运动伪影,即使经软件矫正后仍会出现剪影,引发误诊[16],故扫描开始前需指导患者进行呼吸训练。本研究结果显示,肝脏动态增强扫描联合AEF评估TACE疗效的AUC为0.902,明显高于肝脏动态增强扫描及AEF单独评估的0.793和0.771,提示肝脏动态增强扫描联合AEF能提高TACE疗效评估的准确度,可用于指导临床疗效评估。
在临床工作中发现,患者在吸气末屏住气扫描与呼气末屏住气扫描相比能够较好的保持呼吸幅度一致,从而在后处理中不同扫描时期的图像能都较好的进行匹配,得到良好的后处理图像,目前尚无统计学支持。
综上所述,AEF对TACE治疗效果具有良好的评估价值,联合肝脏动态增强扫描能进一步提高准确度,可作为评估TACE疗效的辅助手段。本研究尚有不足之处,研究样本较少且为回顾性分析,后期需开展多中心、大样本量、前瞻性研究,以及较长时间的随访,以期得到更为准确的AEF与HCC患者预后的关系,进而指导临床诊疗。
-
表 1 两种去噪方法PSNR对比
Table 1 PSNR comparison of two denoising methods
去噪方法 PSNR/dB 带通滤波去噪 18.729 CycleGAN去噪 25.247 -
[1] 李一保, 张玉芬, 刘玉兰, 等. 浅地层剖面仪在海洋工程中的应用[J]. 工程地球物理学报, 2007,(1): 4−8. LI Y B, ZHANG Y F, LIU Y L, et al. Application of subbottom profiler to ocean engineering[J]. Chinese Journal of Engineering Geophysics, 2007, (1): 4−8. (in Chinese).
[2] 刘玉萍, 丁龙翔, 杨志成, 等. 利用浅剖资料进行海底底质分析[J]. 物探与化探, 2016,40(1): 66−72. LIU Y P, DING L X, YANG Z C, et al. Seabed sediment analysis using sub-bottom profile data[J]. Geophysical and Geochemical Exploration, 2016, 40(1): 66−72. (in Chinese).
[3] 李平, 杜军. 浅地层剖面探测综述[J]. 海洋通报, 2011,30(3): 344−350. LI P, DU J. Review on the probing of sub-bottom profiler[J]. Marine Science Bulletin, 2011, 30(3): 344−350. (in Chinese).
[4] 顾兆峰, 张志珣, 刘怀山. 南黄海西部地区浅层气地震特征[J]. 海洋地质与第四纪地质, 2006,(3): 65−74. GU Z F, ZHANG Z X, LIU H S. Seismic features of shallow gas in the western area of the Yellow Sea[J]. Marine Geology & Quaternary Geology, 2006, (3): 65−74. (in Chinese).
[5] 颜中辉, 栾锡武, 潘军, 等. 海上浅地层剖面处理的关键去噪技术[J]. 海洋地质前沿, 2016,32(9): 64−70. YAN Z H, LUAN X W, PAN J, et al. Key denoising techniques in marine sub-bottom shallow profiling[J]. Marine Geological Frontiers, 2016, 32(9): 64−70. (in Chinese).
[6] 施凤. 浅地层剖面数据精处理关键技术研究[D]. 武汉: 武汉大学, 2017. [7] 冯强强, 温明明, 吴衡, 等. 海洋浅地层剖面资料的数据处理方法[J]. 海洋地质前沿, 2013,29(11): 49−53, 66. FENG Q Q, WEN M M, WU H, et al. Data processing methods for marine sub-bottom profiles[J]. Marine Geological Frontiers, 2013, 29(11): 49−53, 66. (in Chinese).
[8] JEONG B J, LEE Y H. Design of weighted order statistic filters using the perceptron algorithm[J]. IEEE Transactions on Signal Processing, 1994, 42(11): 3264−3269. doi: 10.1109/78.330393
[9] GONZALEZ R C, WOODS R E. 数字图像处理[M]. 阮秋琦, 阮宇智, 译. 北京: 电子工业出版社, 2003. [10] MA Y Y, LI G F, WANG Y J, et al. Random noise attenuation by f-x spatial projection-based complex empirical mode decomposition predictive filtering[J]. Applied Geophysics, 2015, 12(1): 47−54. doi: 10.1007/s11770-015-0467-3
[11] 潘军, 栾锡武, 孙运宝, 等. SRME技术在海洋浅水高分辨率地震勘探中的应用[J]. 地球物理学进展, 2015,30(1): 429−434. PAN J, LUAN X W, SUN Y B, et al. Application of SRME technology in marine shallow water high resolution seismic exploration[J]. Progress in Geophysics, 2015, 30(1): 429−434. (in Chinese).
[12] 王兆湖, 王建民, 高振山, 等. 叠前自适应f-x域相干噪音衰减技术及应用[J]. 地球物理学进展, 2013,28(5): 2605−2610. WANG Z H, WANG J M, GAO Z S, et al. Pre-stack self-adaptive f-x domain coherent noise attenuation technology and application[J]. Progress in Geophysics, 2013, 28(5): 2605−2610. (in Chinese).
[13] YANG J, LIN N, ZHANG K, et al. Reservoir characterization using multi-component seismic data in a novel hybrid model based on clustering and deep neural network[J]. Natural Resources Research, 2021: 1−26.
[14] 杨晶, 丁仁伟, 林年添, 等. 基于深度学习的地震断层智能识别研究进展[J]. 地球物理学进展, 2022,37(1): 298−311. YANG J, DING R W, LIN N T, et al. Research progress of intelligent identification of seismic fault based on deep learning[J]. Progress in Geophysics, 2022, 37(1): 298−311. (in Chinese).
[15] 刘俊, 曹俊兴, 丁蔚楠, 等. 基于双向长短期记忆神经网络的储层孔隙度预测方法研究[J/OL]. 地球物理学进展: 1-9. (2021-07-23)[2021-12-18]. http://kns.cnki.net/kcms/detail/11.2982.P.20210723.1308.045.html. LIU J, CAO J X, DING W N, et al. Research on reservoir porosity prediction method based on bidirectional long-term and short-term memory neural network[J/OL]. Progress in Geophysics: 1-9. (2021-07-23)[2021-12-18]. http://kns.cnki.net/kcms/detail/11.2982.P.20210723.1308.045.html. (in Chinese).
[16] TUKEY J W. Exploratory data analysis[J]. Journal of the American Statistical Association, 1977, 28(1): 163−182.
[17] SUMAN S. Image denoising using new adaptive based median filter[J]. Signal & Image Processing an International Journal, 2014, 5(4): 1−13.
[18] QUAN H D, ZHANG H C, JIA H T. An improved adaptive median filtering algorithm[C]//International Symposium on Test Automation and Instrumentation, 2006: 269-272. [19] JAIN V, SEUNG H S. Natural image denoising with convolutional networks[C]//International Conference on Neural Information Processing Systems. Piscataway: IEEE Press, 2008: 769-776.
[20] ZHANG K, ZUO W, CHEN Y, et al. Beyond a Gaussian denoiser: Residual learning of deep CNN for image denoising[J]. IEEE Transactions on Image Processing, 2016, 26(7): 3142−3155.
[21] GOODFELLOW I J, POUGET-ABADIE J, MIRZA M, et al. Generative adversarial nets[C]// Proceedings of the 27th International Conference on Neural Information Processing Systems, Montreal, 2014. Cambridge: MIT Press, 2014: 2672-2680.
[22] WANG H, LI Y, DONG X. Generative adversarial network for desert seismic data denoising[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020. DOI:10.1109/ TGRS.2020.3030692.
[23] DONG X T, LI Y. Denoising the optical fiber seismic data by using convolutional adversarial network based on loss balance[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020. DOI: 10.1109/TGRS.2020.3036065.
[24] ZHU J Y, PAKR T, ISOLA P, et al. Unpaired image-to-image translation using cycle-consistent adversarial networks[C]//Proceedings of the 2017 IEEE International Conference on Computer Vision, Venice, 2017. Washington: IEEE Computer Society, 2017: 2242-2251.
[25] 梁俊杰, 韦舰晶, 蒋正锋. 生成对抗网络GAN综述[J]. 计算机科学与探索, 2020,14(1): 1−17. LIANG J J, WEI J J, JIANG Z F. Generative adversarial networks GAN overview[J]. Journal of Frontiers of Computer Science & Technology, 2020, 14(1): 1−17. (in Chinese).
[26] 梁雪灿. 基于生成对抗网络的声学图像超分辨率研究[D]. 哈尔滨: 哈尔滨工程大学, 2019. [27] ISOLA P, ZHU J Y, ZHOU T H, et al. Image-to-image translation with conditional adversarial networks[C]//Proceedings of the 2017 IEEE Conference on Computer Vision and Pattern Recognition, Honolulu, 2017. Washington: IEEE Computer Society, 2017: 5967-5976.
[28] 王照, 陈恩庆. 基于深度残差生成对抗网络的本征图像分解算法[J]. 计算机应用与软件, 2022,39(3): 201−206. WANG Z, CHEN E Q. Eigen image decomposition algorithm for generating countermeasure network based on depth residual[J]. Computer Applications and Software, 2022, 39(3): 201−206. (in Chinese).
[29] NAYEF B H, ABDULLAH S N H S, SULAIMAN R, et al. Optimized leaky ReLU for handwritten Arabic character recognition using convolution neural networks[J]. Multimedia Tools and Applications, 2021, 81(2): 2065−2094.
[30] 杨久强, 林年添, 张凯, 等. 深度神经网络模型超参数选取及评价研究: 以含油气性多波地震响应特征提取为例[J]. 石油物探, 2022,61(2): 236−244. YANG J Q, LIN N T, ZHANG K, et al. Hyperparametric selection and evaluation of deep neural network models: A case study of feature extraction of multi-wave seismic response in an oil-gas reservoir[J]. Geophysical Prospecting for Petroleum, 2022, 61(2): 236−244. (in Chinese).
[31] 王议迎, 丁仁伟, 李建平, 等. 联合改进STA/LTA与MLoG算子的微震P波到时自动拾取方法[J]. 山东科技大学学报(自然科学版), 2021,40(6): 1−10. WANG Y Y, DING R W, LI J P, et al. Automatic pickup of microseismic P-wave arrival based on improved STA/LTA and MloG operators[J]. Journal of Shandong University of Science and Technology (Natural Science), 2021, 40(6): 1−10. (in Chinese).
-
期刊类型引用(4)
1. 席晓旭,陈志晔. 光谱CT定量增强参数对肾细胞癌与肾血管平滑肌脂肪瘤的诊断价值. 分子影像学杂志. 2025(02): 218-222 . 百度学术
2. 刘静垚,樊文萍,刘梦琦,葛文浩,姚慧,刘明波,陈志晔. 基于碘图的光谱CT多模态参数成像在甲状腺乳头状癌中的诊断价值. 分子影像学杂志. 2024(02): 138-142 . 百度学术
3. 王兴龙,冯坤鹏,袁牧,孟令武. CT动脉增强分数、甲胎蛋白异质体3、中性粒细胞与淋巴细胞比值联合检测对肝癌患者介入化疗疗效的评估价值. 癌症进展. 2024(10): 1146-1150 . 百度学术
4. 徐石. MSCT增强扫描对乙型肝炎肝硬化背景下原发性肝癌的诊断价值. 影像研究与医学应用. 2023(13): 132-134 . 百度学术
其他类型引用(0)