ISSN 1004-4140
CN 11-3017/P

虚拟单能图像与常规CT图像的图像质量对比

陈琳钰, 方姝, 陈勇, 董海鹏, 王兰, 常蕊

陈琳钰, 方姝, 陈勇, 等. 虚拟单能图像与常规CT图像的图像质量对比[J]. CT理论与应用研究, 2022, 31(2): 219-226. DOI: 10.15953/j.ctta.2021.041.
引用本文: 陈琳钰, 方姝, 陈勇, 等. 虚拟单能图像与常规CT图像的图像质量对比[J]. CT理论与应用研究, 2022, 31(2): 219-226. DOI: 10.15953/j.ctta.2021.041.
CHEN L Y, FANG S, CHEN Y, et al. Comparison of image quality between virtual monochromatic images and conventional CT images[J]. CT Theory and Applications, 2022, 31(2): 219-226. DOI: 10.15953/j.ctta.2021.041. (in Chinese).
Citation: CHEN L Y, FANG S, CHEN Y, et al. Comparison of image quality between virtual monochromatic images and conventional CT images[J]. CT Theory and Applications, 2022, 31(2): 219-226. DOI: 10.15953/j.ctta.2021.041. (in Chinese).

虚拟单能图像与常规CT图像的图像质量对比

详细信息
    作者简介:

    陈琳钰: 女,上海交通大学附属瑞金医院放射科技师,主要从事MRI、CT、X-ray等工作,E-mail:cly04164@rjh.com.cn

  • 中图分类号: R  814

Comparison of Image Quality between Virtual Monochromatic Images and Conventional CT Images

  • 摘要: 目的:对比双源双能CT的虚拟单能图像与常规CT图像的图像质量。方法:采用西门子SOMATOM Force双源双能CT扫描图像质量控制体模Catphan 500的高对比度分辨率CTP 528组件和低对比度分辨率CTP 515组件,分别采用双球管扫描方案重建低、中、高3档虚拟单能能谱图像(A组:40 keV,B组:70 keV,C组:100 keV)和一种单球管扫描方案图像(D组:120 kVp)。两名放射科医师分别测量CTP 528的高对比分辨率(最佳线对数)和CTP 515的低对比度分辨率(1% 对比度浓度下的最小孔径);随后测量各组图像的平均噪声、信噪比(SNR)及对比噪声比(CNR)。结果:两名医师的图像一致性评价良好(高/低对比度分辨率:Kappa系数为0.667;噪声参数:Spearman相关系数r=0.920。各组图像中,D组高对比度分辨率最佳(线对6 lp/cm);A、B组图像低对比度分辨率最佳(孔径均为4 mm)。C组图像平均噪声最低(2.333±0.289),SNR表现最佳(37.944±3.949);噪声和SNR在各组图像间差异有统计学意义。B组图像CNR最高(3.361±0.875),但各组图像间CNR差异不具有统计学意义。结论:双能CT虚拟单能图像相比常规CT图像可以在保持低对比度分辨率的基础上,改善图像噪声;随着虚拟单能能级的增加,图像噪声减小,信噪比增加。
    Abstract: Objective: To compare the image quality between virtual monochromatic images (VMI) of dual-source dual-energy CT (DECT) and conventional CT images. Siemens SOMATOM Force dual source dual energy CT was adopted in the scanning of the high contrast resolution CTP 528 component and low contrast resolution CTP 515 component of image quality control phantom catphan 500. Dual-tube dual-energy scanning scheme was used to reconstruct low, medium and high-level VMI (group A: 40 keV, group B: 70 keV, group C: 100 keV) and a single tube scanning scheme was used to reconstruct conventional images (Group D: 120 kVp). Two radiologiss carried out the measurement of high resolution by the smallest line pair one could distinguish and low-contrast sensitivity was measured by the smallest diameter under 1% contrast. Mean noise, signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were also evaluated for each group. Results: The image consistency of the two radiologists turned out well (high reolution and low-contrast sensitivity: Kappa coefficient=0.667; noise parameter: Spearman r = 0.920). Among all images, D group performed the best in high resolution (6 lp/cm) and A and B groups showed the best low-contrast sensitivity (diameter of 4 mm). C group had the lowest noise (2.333±0.289) and highest SNR (37.944±3.949). Mean noise and SNR were found significantly difference among the four groups. B group showed the highest CNR (3.361±0.875), but no statistical significance was found among the groups. Conclusion: Compared with conventional CT, VMI images from DECT could improve image quality while retain low-contrast sensitivity; image noise decreases while SNR increases as VMI elevates.
  • CT是目前临床上诊治疾病最常用的影像学手段之一,但由于X线的辐射特性,剂量问题一直受到公众的担忧。CT检查以每年10% 以上的速度增长,而辐射剂量也由上世纪80年代的0.54 mSv增加到了3 mSv以上[1-2]。CT检查虽只占整个放射检查的很小比例,但造成的辐射却占到了所有放射检查的一半左右[1]。以胸部检查举例,CT造成的辐射剂量是平片的100倍以上[3]。调查显示,每年大约有1.2%~2% 的新发肿瘤由CT扫描引起[4-6]

    如何在保持图像质量的同时降低辐射剂量,是目前CT临床研究的主要改进方向和研究热点[7-8]。管电压除与剂量明显相关以外,还是图像质量的重要影响因素,但目前针对不同CT机器的管电压对比的相关研究较少。研究表明,由于厂商调校的差异,不同机型对CT值等定量参数的准确性方面存在一定的差异[9-12],对于拥有多机型的医院来说统一的诊断标准非常重要。因此,本研究旨在通过体模研究对比两种机型两组管电压下的图像质量。

    选用Catphan 500图像质量控制模体的CTP 528模块和CTP 515模块进行图像质量评价。其中CTP 528模块为高对比度检测组件,直径15 cm,厚度4 cm,有21组高密度线对结构呈放射状分布,分辨率为1~21 lp/cm;CTP 515的直径15 cm,厚4 cm,内外两组低密度孔经结构呈放射状分布,内层孔阵浓度分别为0.3%、0.5% 和1.0%,内层孔径直径分别对应为3、5、7和9 mm;外层孔阵浓度分别为0.3%、0.5% 和1.0%,外层孔径直径分别对应为2、3、4、5、6、7、8、9和15 mm。

    采用西门子SOMATOM Force CT和飞利浦IQon Spectral CT分别扫描CTP 528和CTP 515两个组件。每台机器分别采取120 kVp和140 kVp两种单能电压扫描方案共产生4组图像,分别为A组Force120kVp、B组Force140kVp、C组IQon120kVp和D组IQon140kVp图像。4组图像扫描剂量均为20 mGy,重建层厚均为5 mm。

    具体扫描参数。A组管电流299 mAs,B组图像205 mAs,其余参数两组一致:螺距1.0,旋转时间1 s,准直宽度192×0.6 mm,迭代重建算法ADMIRE(level 2,0~5),标准滤波函数Br40。C组管电流221 mAs,D组管电流154 mAs,其余参数两组一致:螺距1.014,旋转时间0.75 s,准直宽度64×0.625 mm,迭代重建算法iDose4(level 3,0~6),标准滤波函数Standard(B)。

    由两名放射科医师分别记录CTP 528高对比度分辨率组件的最佳的线对(1~2 lp/cm,图1(a))以及CTP 515低对比度分辨率组件的1% 浓度下的最小孔径(直径2~15 mm,要求分辨80% 以上有效面积,图1(b)),比较两名医师的主观评价结果,如有分歧则由另一名高年资医生做最终评判。两个指标均在组件中心及上下共3个层面进行测量以记录各组图像最优结果。

    图  1  Force140 kVp图像的Catphan 500体模CTP 528高对比度分辨率和CTP 515低对比度分辨率测量
    Figure  1.  Measurement of CTP 528 high contrast resolution and CTP 515 low contrast resolution of Catphan 500 phantom of Force 140 kVp image

    将Froce和IQon的图像分别传至西门子Syngo Via(版本VB10)和飞利浦Intellispace Portal(版本10.1)后处理工作站上,由以上两名放射科医师对CTP 515模块图像进行测量。采用两个面积均为100 mm2的圆形感兴趣区域(ROI)分别测量1% 浓度组的15 mm孔径和邻近固态水背景的CT值(CTm和CTw,单位HU)。将15 mm孔径ROI标准差(SD,单位HU)定义为噪声,随后计算各组图像的信噪比(signal noise ratio,SNR=CTm/SD)和对比噪声比[contrast noise ratio,CNR= (CTm - CTw)/SD]。

    所有主观及客观评价指标均在各组件中心的连续3个层面上进行测量;取两名医师测量的SD、CTm和CTw的均值计算各组图像的SNR和CNR。

    运用SPSS 22.0版本软件进行数据分析,连续性指标采用平均值±标准差表示。4组图像的噪声值、SNR及CNR的比较采用单因素方差分析。两名医生对高/低对比度分辨率的主观一致性评价采用Cohen's Kappa检验,对图像噪声、CTm和CTw的客观一致性评价采用Spearman相关性分析。两台机器内部的噪声参数的比较采用独立样本t检验;以P<0.05认为差异具有统计学意义。

    Force120kVp组和IQon120kVp组图像具有最好的高对比度分辨率,线对均为6 lp/cm,其余两组的线对均为5 lp/cm。4组图像低对比度分辨率相同,孔径均为5 mm。

    各组图像的噪声值和SNR间差异均具有统计学意义,其中IQon120kVp组图像噪声最大(3.400±0.050),SNR最低(17.040±0.410);Force120kVp组图像噪声最小(2.500±0.000),SNR最高(30.806±1.398)。Force组图像SNR均高于IQon组。IQon140kVp组图像的CNR最高(3.325±0.300),但CNR在各组图像间差异无统计学意义。

    组内比较方面,Force两组图像噪声差异具有统计学意义,但SNR和CNR差异均无统计学意义;IQon两组图像之间噪声、SNR及CNR之间差异均无统计学意义(具体结果见表1)。

    表  1  四组图像噪声参数组间对比
    Table  1.  Comparison of image noise parameters among the four groups
    噪声参数Force120kVpForce140kVpIQon120kVpIQon140kVp$F $$P $
    SD2.500±0.0002.833±0.2893.400±0.0502.783±0.3518.1530.008
    SNR30.806±1.39829.472±4.42717.040±0.41023.844±2.82315.9050.001
    CNR2.778±0.2412.556±0.4192.604±0.2083.325±0.3004.0650.050
    下载: 导出CSV 
    | 显示表格

    两名医师对各组图像的主观指标的一致性评价良好,均为0.750,差异具有统计学意义;同样,两人对各组图像的客观噪声参数一致性较好(0.630~0.857),差异均具有统计学意义,具体结果见表2

    表  2  两名放射科医师的一致性分析结果
    Table  2.  Reliability analysis between the two radiologists
    主客观指标一致性指数P  主客观指标一致性指数P
    高对比度分辨率10.7500.007  CTm20.8570.000
    低对比度分辨率10.7500.007  CTw20.7680.004
    噪声20.6300.028 
    注:1-使用Cohen's Kappa检验;2-使用Spearman相关性分析。
    下载: 导出CSV 
    | 显示表格

    在CT临床应用日益普及和快速增长的情况下,如何在达到诊断标准要求的图像质量的基础上进一步降低辐射剂量是目前努力的主要方向。辐射剂量可能受扫描次数、管电压、管电流、旋转时间、螺距及患者体型等多因素的影响。降低辐射剂量的常用手段包括管电流调节、自动管电压选择及图像重建优化算法等,但这些技术目前主要应用在头、胸部CT扫描方面,其他部位的研究很多还处在探索阶段[13]。不同研究表明,管电压在同等扫描条件下相比其他扫描参数对图像质量和剂量的影响更为显著[14-16]

    对于拥有多CT机型的医疗机构来说,由于辐射和伦理的问题,很难对同一个患者进行重复扫描,这样,由CT机器或患者本身导致的图像差异便无法验证。体模研究在这种情况下就显得尤为重要。在相同的剂量条件下,不同物质在不同管电压下的衰减系数存在区别;低电压(如80 kVp)其X线光子能量接近43.7 keV,接近碘原子的k-edge临界值(33 keV),因此图像的对比度会升高[17]

    本研究采用固定剂量,每台机器两组管电压的扫描方案,结果显示两台机器120 kVp图像的高对比度分辨率均略好于140 kVp的图像,但4组图像的低对比度分辨率相同,说明120 kVp的图像相较于140 kVp空间分辨率更高,能显示更多细节,但是软组织分辨率相近。客观指标图像噪声参数部分,两台机器之间的差异较为明显,尤其对于西门子Force来说,120 kVp图像的3个噪声参数均好于140 kVp图像。而飞利浦IQon的140 kVp图像略均好于120 kVp图像,这可能是由于IQon光谱CT在以上两种管电压下可以同时实行单能和双能扫描,双层探测器结构加上与之相对应的物质分离算法使得高能量图像具有更广的物质分离光谱,图像质量有一定提升,这与Wellenberg等[18]的发现相类似。值得注意的是,虽然本研究中两者噪声有区别,但差异无统计学意义。

    由于厂商对机器校准和算法的区别,要求采用同一套扫描标准以达到相同的图像质量仍较为困难。因此,因地制宜地针对不同CT机器调整适用于影像诊断要求图像质量的扫描参数显得尤为重要。

    本研究不足,①两台机型的 CT扫描参数较单一,并未使用低电压低剂量扫描作对比,加上是体模研究,其临床适用性还需进一步验证;②迭代重建算法可能是本研究一个潜在影响因素,特别是本研究中Force使用了全模迭代算法ADMIRE,而IQon使用的是部分迭代算法iDose4;③不同厂商迭代重建算法对图像质量的影响还需要进一步研究。

    综上所述,相同条件下,不同厂商不同机型之间图像质量存在差异;两台机器120 kVp图像的高对比度分辨率好于140 kVp图像;Force120kVp的图像噪声值低于其140 kV图像,而IQon两组图像的噪声参数无显著区别。

  • 图  1   Capthan 500体模组件构造图

    Figure  1.   Schematic components of the Capthan 500 phantom

    图  2   三档虚拟单能图像的示例图

    (a)~(c)分别代表40、70和100 keV图像;(a1)、(b1)和(c1)分别代表相同窗宽窗位下(400,0)40、70和100 keV的图像对比;(a2)和(a3)、(b2)和(b3)、(c2)和(c3)分别代表相应虚拟单能图像的CTP 515低对比度分辨率和CTP 528高对比度分辨率模块图像。

    Figure  2.   Schematic diagram for the three series of virtual monochromatic images

    表  1   四组图像高对比度分辨率和低对比度分辨率比较

    Table  1   Comparison of high-contrast resolution and low-contrast resolution of four groups of images

    对比度分辨率A组
    (40 keV)
    B组
    (70 keV)
    C组
    (100 keV)
    D组
    (120 kVp)
    高对比对分辨率(线对/(lp/cm))5556
    低对比度分辨率(孔径直径/mm)4455
    下载: 导出CSV

    表  2   四组图像图噪声参数评价结果比较

    Table  2   Comparison of the evaluation results of the noise parameters of the five groups of images

    组别SD/HUSNRCNR
    A组(40 keV)4.167±0.28910.942±0.4652.383±0.336
    B组(70 keV)2.500±0.50024.028±5.0003.361±0.875
    C组(100 keV)2.333±0.28937.944±3.9493.139±0.337
    D组(120 kVp)2.500±0.00030.806±1.3982.778±0.241
    $P $<0.001<0.0010.179
    下载: 导出CSV
  • [1]

    SUN H, QIU S, LOU S, et al. A correction method for nonlinear artifacts in CT imaging[C]//IEEE Engineering in Medicine and Biology Society. Conference Proceedings 2004, 2004: 1290-1293.

    [2]

    RODRÍGUEZ-GRANILLO G A, ROSALES M A, DEGROSSI E, et al. Signal density of left ventricular myocardial segments and impact of beam hardening artifact: Implications for myocardial perfusion assessment by multidetector CT coronary angiography[J]. International Journal of Cardiovascular Imaging, 2010, 26(3): 345−354. doi: 10.1007/s10554-009-9531-5

    [3] 任国庆, 滑炎卿, 李剑颖. CT能谱成像的基本原理及临床应用[J]. 国际医学放射学杂志, 2011,34(6): 559−563. doi: 10.3784/j.issn.1674-1897.2011.06.Z0612

    REN G Q, HUA Y Q, LI J Y. The basic principle and clinical application of CT spectrum. imaging[J]. International Journal of Medical Radiology, 2011, 34(6): 559−563. (in Chinese). doi: 10.3784/j.issn.1674-1897.2011.06.Z0612

    [4]

    Mc COLLOUGH C H, LENG S, YU L, et al. Dual- and multi-energy CT: Principles, technical approaches, and clinical applications[J]. Radiology, 2015, 276(3): 637−653. doi: 10.1148/radiol.2015142631

    [5] 王文杰, 陈平, 潘晋孝, 等. 基于参考组分的双能CT成像方法[J]. CT理论与应用研究, 2021,30(1): 61−69. DOI: 10.15953/j.1004-4140.2021.30.01.06.

    WANG W J, CHEN P, PAN J X, et al. Dual-energy CT imaging method based on reference components[J]. CT Theory and Applications, 2021, 30(1): 61−69. DOI: 10.15953/j.1004-4140.2021.30.01.06. (in Chinese).

    [6] 许伟雄, 黄旭璇, 林顺发, 等. 腹部动脉能谱CT混合能量与单能量图像的比较研究[J]. 中国CT和MRI杂志, 2014,12(8): 91−94. doi: 10.3969/j.issn.1672-5131.2014.08.30

    XU W Q, HUANG X X, LIN S F, et al. Comparative study of abdominal artery energy spectrum CT mixed energy and single energy image[J]. Chinese Journal of CT and MRI, 2014, 12(8): 91−94. (in Chinese). doi: 10.3969/j.issn.1672-5131.2014.08.30

    [7] 朱晓红, 刘斌, 周勇, 等. 冠状动脉能谱CT单能量成像与混合能量成像质量的比较−静态模型的实验研究[J]. 临床放射学杂志, 2012,31(5): 738−742.

    ZHU X H, LIU B, ZHOU Y, et al. Comparison of the quality of coronary artery energy spectrum CT single-energy imaging and mixed-energy imaging−An experimental study of a static model[J]. Journal of Clinical Radiology, 2012, 31(5): 738−742. (in Chinese).

    [8] 韦炜, 邓克学, 赵英明, 等. 腹部双能量能谱CT成像中混合能量模式与单能量模式重建图像的比较研究[J]. 安徽医科大学学报, 2016,51(11): 1650−1653.

    WEI W, DENG K X, ZHAO Y M, et al. A comparative study of mixed energy mode and single energy mode in dual-energy spectrum CT imaging of the abdomen[J]. Journal of Anhui Medical University, 2016, 51(11): 1650−1653. (in Chinese).

    [9] 李林坤, 王秋燕, 丛湘华, 等. 双下肢单能量60 keV与混合能量120 kVp CTA应用价值比较[J]. 临床放射学杂志, 2019,38(11): 2180−2184.

    LI L K, WANG Q Y, CONG X H, et al. Comparison of the application value of single energy 60 keV and mixed energy 120 kVp CTA for both lower limbs[J]. Journal of Clinical Radiology, 2019, 38(11): 2180−2184. (in Chinese).

    [10] 王伟新, 张秋奂, 郭鹏德, 等. 能谱CT单能量结合ASiR技术对腹部静脉成像质量的优化研究[J]. CT理论与应用研究, 2019,28(1): 61−72. DOI: 10.15953/j.1004-4140.2019.28.01.07.

    WANG W X, ZHANG Q H, GUO P D, et al. Study on the optimization of energy spectrum CT single energy combined with ASiR technology on the quality of abdominal vein imaging[J]. CT Theory and Applications, 2019, 28(1): 61−72. DOI: 10.15953/j.1004-4140.2019.28.01.07. (in Chinese).

    [11]

    ALBRECHT M H, VOGL T J, MARTIN S S, et al. Review of clinical applications for virtual monoenergetic dual-energy CT[J]. Radiology, 2019, 293(2): 260−271. doi: 10.1148/radiol.2019182297

    [12]

    LEITHNER D, MAHMOUDI S, WICHMANN J L, et al. Evaluation of virtual monoenergetic imaging algorithms for dual-energy carotid and intracerebral CT angiography: Effects on image quality, artefacts and diagnostic performance for the detection of stenosis[J]. European Journal Radiology, 2018, 99: 111−117. doi: 10.1016/j.ejrad.2017.12.024

    [13]

    KAUP M, SCHOLTZ J E, ENGLER A, et al. Dual-energy computed tomography virtual monoenergetic imaging of lung cancer: Assessment of optimal energy levels[J]. Journal of Computer Assisted Tomography, 2016, 40(1): 80−85. doi: 10.1097/RCT.0000000000000319

    [14]

    MARTIN S S, WICHMANN J L, PFEIFER S, et al. Impact of noise-optimized virtual monoenergetic dual-energy computed tomography on image quality in patients with renal cell carcinoma[J]. European Journal Radiology, 2017, 97: 1−7. doi: 10.1016/j.ejrad.2017.10.008

    [15]

    BODANAPALLY U K, ARCHER-ARROYO K, DREIZIN D, et al. Dual-energy computed tomography imaging of head: Virtual high-energy monochromatic (190 keV) images are more reliable than standard 120 kV images for detecting trau-matic intracranial hemorrhages[J]. Journal of Neurotrauma, 2019, 36(8): 1375−1381. doi: 10.1089/neu.2018.5985

    [16]

    LEE S M, KIM S H, AHN S J, et al. Virtual monoenergetic dual-layer, dual-energy CT enterography: Optimization of keV settings and its added value for Crohn’s disease[J]. European Radiology, 2018, 28(6): 2525−2534. doi: 10.1007/s00330-017-5215-z

    [17]

    SYMONS R, CHOI Y, CORK T E, et al. Optimized energy of spectral coronary CT angiography for coronary plaque detection and quantification[J]. Journal of Cardiovascular Computed Tomography, 2018, 12(2): 108−114. doi: 10.1016/j.jcct.2018.01.006

    [18]

    APFALTRER P, SUDARSKI S, SCHNEIDER D, et al. Value of monoenergetic low-kV dual energy CT datasets for improved image quality of CT pulmonary angiography[J]. European Journal Radiology, 2014, 83(2): 322−328. doi: 10.1016/j.ejrad.2013.11.005

    [19]

    Mc NAMARA M M, LITTLE M D, ALEXANDER L F, et al. Multireader evaluation of lesion conspicuity in small pancreatic adenocarcinomas: Complimentary value of iodine material density and low keV simulated monoenergetic images using multiphasic rapid kVp-switching dual energy CT[J]. Abdominal Imaging, 2015, 40(5): 1230−1240. doi: 10.1007/s00261-014-0274-y

    [20]

    DARRAS K E, Mc LAUGHLIN P D, KANG H, et al. Virtual monoenergetic reconstruction of contrast-enhanced dual energy CT at 70 keV maximizes mural enhance-ment in acute small bowel obstruction[J]. European Journal Radiology, 2016, 85(5): 950−956. doi: 10.1016/j.ejrad.2016.02.019

    [21]

    MATSUMOTO K, JINZAKI M, TANAMI Y, et al. Virtual monochromatic spectral imaging with fast kilovoltage switching: Improved image quality as compared with that obtained with conventional 120-kVp CT[J]. Radiology, 2011, 259(1): 257−62. doi: 10.1148/radiol.11100978

    [22] 周旸, 曾勇明, 周蜜, 等. 双能CT单能谱70 keV图像与常规CT 120 kVp图像的等效性[J]. 中国医学影像技术, 2015,31(7): 1100−1105.

    ZHOU Y, CENG Y M, ZHOU M, et al. The equivalence of dual-energy CT single energy spectrum 70 keV image and conventional CT 120 kVp image[J]. China Medical Imaging Technology, 2015, 31(7): 1100−1105. (in Chinese).

  • 期刊类型引用(3)

    1. 刘雨欣,陈凡秀,孙洁,王远,王潇,于洋,顾焱吉. 管电压管电流和开机时间对CT图像质量的影响. 实验力学. 2024(01): 17-26 . 百度学术
    2. 朱丽娟,马瑞,沈云,汪芳,杨彦兵,曹永佩,杨利莉,吴小红. 探讨不同管电压下DLIR重建算法对冠状动脉CTA图像质量的影响. 宁夏医学杂志. 2023(09): 787-791+865 . 百度学术
    3. 王振宇,黄蕾. 分析CT图像的质量控制与设备维护保养. 中国设备工程. 2023(18): 85-87 . 百度学术

    其他类型引用(0)

图(2)  /  表(2)
计量
  • 文章访问数:  424
  • HTML全文浏览量:  544
  • PDF下载量:  75
  • 被引次数: 3
出版历程
  • 收稿日期:  2021-11-20
  • 录用日期:  2021-12-22
  • 网络出版日期:  2021-12-27
  • 发布日期:  2022-03-31

目录

/

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