Spectral CT Multi-material Decomposition in Evaluation of Secondary Lower Extremity Lymphedema: A Prospective Study
-
摘要: 目的：探讨能谱CT多物质解析算法在继发性下肢淋巴水肿分级评价中的应用价值。材料及方法：前瞻性收集2019年6月至2020年6月因继发性下肢淋巴水肿入我院就诊并行能谱CT断层扫描的患者，由两名具有5年以上CT诊断经验的医师共同对能谱CT图像进行分析，采用阈值法测量患侧与健侧下肢体积，计算（患侧下肢体积 - 健侧下肢体积）/健侧下肢体积的百分比，根据2016版国际淋巴学会共识，将入组患者分为轻度、中度、重度；使用能谱CT多物质解析算法分别测量患侧与健侧下肢脂肪含量，得出脂肪体积分数，采用配对<i<t</i<检验分别分析轻、中、重度患者患侧与健侧下肢脂肪体积分数之间的差异；采用秩和检验分别分析轻、中、重度患者两两之间患侧下肢脂肪体积分数的差异。结果：最终入组患者 40例，女/男（36/4），年龄范围32～71岁，中位年龄（53±10）岁，病程（3±5）年，轻度8例，中度14例，重度18例。轻、中、重度患者的患侧与健侧下肢脂肪体积分数之间差异均有统计学意义；轻、中、重度患者两两之间患侧下肢脂肪体积分数差异无统计学意义。结论：能谱CT多物质解析算法可以量化继发性下肢淋巴水肿的脂肪含量比，但不能单独用于下肢病变的分级评价中。Abstract: Objective: To explore the application value of spectral CT multi-material decomposition (MMD) algorithm in the grading of secondary lower extremity Lymphedema. Materials and methods: Prospective collection of patients who underwent spectral CT admitted to our hospital from June 2019 to June 2020 were prospectively collected. The spectral CT images were analyzed by two radiologists with more than five years of experience in CT diagnosis. The volume of the both lower limbs were measured using the threshold value software. Based on the percentage increase in the volume of the affected side compared to the healthy side, the enrolled patients were divided into three grades: mild, moderate, and severe. Percentage fat volume fraction (FVF) images were generated by using the MMD algorithm on spectral CT to measure the fat content of the affected and healthy lower limbs respectively. Paired <i<t</i< test were used to analyze the difference of the percentage FVF between the affected and healthy side in mild, moderate, severe patients. Results: Forty patients (female/male, 36/4 cases, with age range of 32-71 years old, median age of (53±10) years) were finally enrolled, and the course of disease was (3±5) years. There were 8 cases in mild group, 14 cases in moderate group and 18 cases in severe group. There were statistical differences in the percentage FVF between the affected and healthy side in mild, moderate and severe patients. There was no statistical difference in the percentage FVF of the affected side among mild, moderate and severe patients.Conclusion: The spectral CT MMD algorithm can quantify the percentage FVF of Secondary Lower Extremity Lymphedema, but it cannot be used alone in the grading of secondary lower extremity lymphedema.
-
Keywords:
- spectrum CT /
- MMD /
- lower extremity lymphedema /
- grade
-
近些年，机器学习方法特别是深度学习方法在地震学中的应用愈发广泛。深度学习方法在地震数据处理中表现出了接近人工分析的精度，并且处理效率高于人工处理[1-2]。典型的应用包括地震震相拾取[3-8]、地震检测[9-10]、P波初动分析[11-14]、震源参数分析[15-16]、地震波形去噪[17-20]、震相关联[21-24]、地震定位[25-30]和地震分类[31-34]等。
目前地震学中多数神经网络模型着眼于单一任务，比如用于震相拾取的神经网络模型以地震波形为输入并以P、S震相作为输出。在一些研究中，尝试将更多信息作为输入，比如Münchmeyer等[35]尝试将地震定位和震级估计同时进行处理，SeisCLIP[36]尝试将震相到时、震中距、方位角等信息输入到神经网络中，并构建了一个能够处理更多信息的神经网络模型。
2023年，中国地震台网中心发布了CSNCD数据集[37]，其中包含了超过
4500 万条人工标注震相到时、震相初动、地震类型、地震震级等信息。在此基础上，我们构建一个用于单台数据分析的双向神经网络预训练模型，称为先导数据大模型。作为预训练模型，我们加入了Pg、Sg、Pn、Sn 4种震相、P波初动方向、地震类型和原始波形4种输出，其中原始波形输出为自监督训练，其他为监督训练。测试表明，我们的模型不仅可以完成常规的Pg、Sg、Pn震相检测、P波初动方向判定和事件类型判断工作，还可以通过迁移学习将模型用于其他波形数据分析工作中。
1. 预训练模型和迁移学习方法简介
通过大量数据训练供不同机器学习任务使用的模型称为预训练模型，目前广泛使用的预训练模型包括单向模型和双向模型。其中单向模型以GPT系列为代表，在处理数据中通常只考虑前文，其通常用于文本生成、文本图像音频生成等多模态工作中。双向模型则以BERT系列为代表，其在处理文本过程中会考虑前后文数据。生成式模型在目前可训练参数数量已经可以达到千亿量级，在文本处理中表现出了良好的性能。SeisCLIP使用生成式模型来处理多台地震频谱并生成所需信息。但是地震学研究中通常需要对原始波形进行分析并从中提取更加准确的信息，在此需求中双向模型通常能够获得更高的精度[38]。
迁移学习指的是通过海量数据训练后的模型在处理其他任务中可以通过少量样本训练即可迁移至新的工作中。迁移学习目前在图像处理、文本处理中均表现出良好的性能。地震学数据分析中需要处理大量的数据，而分析数据的需求各不相同。如果针对每个需求训练深度学习模型，那么在每个场景中都需要海量的人工标注数据，但是在诸如滑坡微震等场景人工标注数据是偏少的。因此，在地震学中同样有研究者使用迁移学习方法以减少对训练数据量的需求。
研究表明，由于预训练模型使用海量数据进行训练，在只使用少量数据进行迁移学习后，神经网络模型在新的场景中表现出良好的效果，比如地震事件分类和远震震相检测。
2. 训练数据
预模型训练需要海量的各种类型数据。因此，我们使用中国地震台网中心的CSNCD数据集来对模型进行训练。CSNCD数据集是一个大规模的全球地震信号标记数据集，其特点在于覆盖了多样化的地理区域和各种地质环境下的地震事件。我们使用的事件和台站如图1所示，这些事件包含了丰富的震相类型和大量的初动数量。
![]()
图 1 使用的台站（蓝色三角形）和事件（红色点）的分布Figure 1. Distribution of stations (blue triangles) and events (red points)除此之外，这些数据还覆盖了从 − 2.0到9.0级不同震级的地震事件。数据集中还包含的广泛的信噪比（SNR）范围，这一特点使得模型能够在不同质量的数据上进行有效的学习，从而提高模型的适应性和泛化能力。
CSNCD数据集的多样性和复杂性为先导数据大模型提供了坚实的基础。
3. 模型结构
在处理更多地震波波形信息时需要较大的模型容量，即需要更多可训练参数。因此我们在构建神经网络过程中以Tansformer模型为核心，同时为了处理波形输入和输出，我们加入卷积神网络，具体网络结构如图2所示。
因为地震波形数据中包含较多高频信息，因此我们使用卷积神经网络构建编码器来处理原始波形，得到波形特征。但波形特征中缺少位置相关的信息，为此我们加入双向循环神经网络构建的位置编码器来对特征位置信息进行编码。在此之上我们加入多层双向Transformer，这可以综合考虑全部波形的特征，并输出特征向量。由于所得特征需要输出事件类型等信息，我们加入了多任务编码（Task ID），这可以提供额外特征供后续分析。
4. 先导大模型的应用
4.1 震相检测和P波初动检测
目前大多数地震学震相检测神经网络模型只能检测Pg和Sg震相，如PhaseNet和EQTransformer等[3-8]。除了Pg、Sg震相拾取模型外，一些研究者还构建了其他的地震学神经网络模型。例如，Yuan等[39]提出的远震震相拾取方案，可以检测PcP震相和PKiKP震相；Ding等[40]构建的PmPNet可以检测PmP震相；Ross等[11]和Zhao等[14]分别构建了用于确定P波初动极性的神经网络模型。
基于CSNCD数据的特点，我们构建的先导大模型通过多个解码器同时实现了震相检测和P波初动检测工作。而且在震相检测方面，不仅可以检测常规的Pg和Sg震相，还可以检测传统上认为难以检测的Pn和Sn震相。图3展示了一个模型震相检测的一个例子。
![]()
图 3 模型拾取结果示例，实线为人工标注震相到时，虚线为模型预测震相到时Figure 3. A example of picks in the test dataset. The solid vertical lines represent the arrival times of phases picked by analysts. The dashed lines represent the probability distributions of phases predicted by the model4.2 迁移学习应用
我们分别使用内蒙地区数据和远震数据测试模型的迁移学习能力。其中，内蒙地区的数据用于事件分类测试，在预训练过程中也有不同事件类型的标签，但是由于预训练中样本极度不均衡（天然地震占比99.5%），这使得原始的预训练模型对事件类型的预测均为天然地震，因此需要进行迁移学习。
事件分类测试数量总计419个，其中251个地震事件用于训练，168个地震事件用于测试。主要包含3种事件类型：天然地震（EQ，103个事件），爆破（EP，244个事件）和塌陷（SS，72个事件）。在经过迁移学习后，模型的分类准确度达到了92.2%（图4）。
我们整理全球的远震数据，挑选
1000 条样本进行远震迁移学习训练和测试。测试中，原始的预训练模型的查准率较高，但查全率偏低，这代表大量样本未被检测到，这说明远震波形特征和近震特征之间存在差异。经过迁移学习后，模型对远震的检测精度明显提高，查准率由原来的0.270提高到了0.708。5. 结论
先导数据大模型在地震学中具有广泛的应用前景，通过利用先导数据大模型，地震台站可以实时处理和分析地震波形数据，快速识别地震事件、震相和初动方向，这可以显著提高地震监测的效率和准确性，尤其在地震多发地区和紧急情况下，能够提供更及时的预警信息。此外，先导数据大模型可以通过迁移学习，适应不同地理区域和不同类型的地震数据。这样可以减少数据集的依赖和节约训练时间成本，实现在不同的地区和任务中快速部署。
先导数据大模型在地震监测、预警、科学研究和数据管理等方面具有重要的应用前景，能够显著提升地震数据处理的自动化和智能化水平，为地震学的发展和地震灾害的防范提供有力支持。
感兴趣的专家读者可以通过数据大模型网址https://github.com/cangyeone/seismological-ai-tools来更好的了解和使用模型。
-
![]()
图 1 能谱CT平扫下肢轴面图像，分别选取患侧与健侧同层面为感兴趣区，该截面患侧和健侧肢体脂肪体积（cm3）分别为26.1±152.4和22.0±169.1
Figure 1. The axial image of lower extremity by spectral CT plain scan. The same section of the affected and healthy side were selected as the ROIs, the fat volume of the affected and healthy side on this section (cm3) is 26.1±152.4 and 22.0±169.1, respectively
![]()
图 2 下肢轴面脂肪体积分数CT图像，分别选取患侧与健侧同层面为感兴趣区，该截面患侧和健侧肢体脂肪体积分数（％）分别为：66.42±37.83和64.83±40.79
Figure 2. The axial FVF image of the lower extremity. The same section of the affected and healthy side were selected as the ROIs, the percentage FVF of the affected and healthy side on this section is 66.42±37.83 and 64.83±40.79, respectively
表 1 能谱CT阈值法测量患侧与健侧下肢体积
Table 1 The volume of the affected and healthy side by spectral CT threshold value software
分级\部位 体积/cm3 统计检验 患侧 健侧 t P 全组（n=40） 9962±2261 6762±1422 11.302 0.000 轻度（n=8） 8358±2213 7379±1862 8.126 0.000 中度（n=14） 8705±1701 6547±1314 15.460 0.000 重度（n=18） 10684±1864 6927±1376 15.487 0.000 
下载: 导出CSV
表 2 能谱CT多物质解析算法测量患侧与健侧下肢脂肪体积分数
Table 2 The percentage FVF of the affected and healthy side by spectral CT MMD algorithm
分级\部位 脂肪体积分数/％ 统计检验 患侧 健侧 t P 全组（n=40） 58.00±6.98 54.96±8.71 5.195 0.000 轻度（n=8） 57.98±5.96 55.37±7.53 2.430 0.045 中度（n=14） 58.20±9.08 54.36±10.19 3.987 0.002 重度（n=18） 57.86±5.84 55.25±8.42 2.693 0.015
下载: 导出CSV
-
[1] 常鲲, 夏松, 孙宇光, 等. 联合应用抽吸减容术与淋巴静脉吻合术治疗下肢继发性淋巴水肿的临床效果[J]. 中华外科杂志, 2017,55(4): 274−278. DOI: 10.3760/cma.j.issn.0529-5815.2017.04.008. CHANG K, XIA S, SUN Y G, et al. Liposuction combined with lymphatico-venous anastomosis for treatment of secondary lymphedema of the lower limbs: A report of 49 cases[J]. Chinese Journal of Surgery, 2017, 55(4): 274−278. DOI: 10.3760/cma.j.issn.0529-5815.2017.04.008. (in Chinese).
[2] 中国妇幼保健协会妇科肿瘤防治专业委员会. 妇科肿瘤治疗后下肢淋巴水肿专家共识[J]. 中国临床医生杂志, 2021,49(2): 149−155. DOI: 10.3969/j.issn.2095-8552.2021.02.006. Professional Committee of Gynecological Cancer Prevention and Treatment of China Maternal and Child Health Association. Expert consensus on lower limb lymphedema after treatment of gynecological tumors[J]. Chinese Journal for Clinicians, 2021, 49(2): 149−155. DOI: 10.3969/j.issn.2095-8552.2021.02.006. (in Chinese).
[3] FORTE A J, KHAN N, HUAYLLANI M T, et al. Lymphaticovenous anastomosis for lower extremity lymphedema: A systematic review[J]. Indian Journal of Plastic Surgery, 2020, 53(1): 17−24. DOI: 10.1055/s-0040-1709372.
[4] Executive Committee. The diagnosis and treatment of peripheral lymphedema: 2016 consensus document of the international society of lymphology[J]. Lymphology, 2016, 49(4): 170−184.
[5] MACLELLAN R A, GREENE A K. Lymphedema[J]. Seminars in Pediatric Surgery, 2014, 23(4): 191−197. DOI: 10.1053/j.sempedsurg.2014.07.004.
[6] 信建峰, 孙宇光, 夏松, 等. 淋巴脂肪抽吸减容术在下肢原发性淋巴水肿中的治疗及分析[J]. 中华整形外科杂志, 2019,35(2): 142−147. DOI: 10.3760/cma.j.issn.1009-4598.2019.02.009. XIN J F, SUN Y G, XIA S, et al. Application of liposuction in treating the primary end-stage lymphedema of lower extremities[J]. Chinese Journal of Plastic Surgery, 2019, 35(2): 142−147. DOI: 10.3760/cma.j.issn.1009-4598.2019.02.009. (in Chinese).
[7] HYODO T, YADA N, HORI M, et al. Multimaterial decomposition algorithm for the quantification of liver fat content by using fast-kilovolt-peak switching dual-energy CT: Clinical evaluation[J]. Radiology, 2017, 283(1): 108−118. DOI: 10.1148/radiol.2017160130.
[8] YAMADA K, SHINAOKA A, KIMATA Y. Three-dimensional imaging of lymphatic system in lymphedema legs using interstitial computed tomography-lymphography[J]. Acta Medica Okayama, 2017, 71(2): 171−177. DOI: 10.18926/AMO/54986.
[9] OLSZEWSKI W L. Lymphovenous microsurgical shunts in treatment of lymphedema of lower limbs: A 45-year experience of one surgeon/one center[J]. European Journal of Vascular and Endovascular Surgery, 2013, 45(3): 282−290. DOI: 10.1016/j.ejvs.2012.11.025.
[10] DELTOMBE T, JAMART J, RECLOUX S, et al. Reliability and limits of agreement of circumferential, water displacement, and optoelectronic volumetry in the measurement of upper limb lymphedema[J]. Lymphology, 2007, 40(1): 26−34.
[11] TSUKADA A, UCHIDA K, AIKAWA J, et al. Unilateral-dominant reduction in muscle volume in female knee osteoarthritis patients: Computed tomography-based analysis of bilateral sides[J]. Journal of Orthopaedic Surgery and Research, 2020, 15(1): 543. DOI: 10.1186/s13018-020-02074-x.
[12] 杨献峰, 杨尚文, 胡安宁, 等. 基于CT值阈值法测量骨骼肌体积[J]. 中国医学计算机成像杂志, 2017,23(3): 287−291. DOI: 10.3969/j.issn.1006-5741.2017.03.019. YANG X F, YANG S W, HU A N, et al. Measurement of skeletal muscle volume based on CT value threshold method[J]. Chinese Journal of Medical Computer Imaging, 2017, 23(3): 287−291. DOI: 10.3969/j.issn.1006-5741.2017.03.019. (in Chinese).
[13] 张国来, 包强, 陈光辉, 等. CT阈值法测量乳突气房体积与面神经管垂直段位置的相关性[J]. 中国医学影像学杂志, 2016,24(3): 175−178. DOI: 10.3969/j.issn.1005-5185.2016.03.004. ZHANG G L, BAO Q, CHEN G H, et al. Correlation between the volume of the mastoid air chamber and the position of the vertical segment of the facial nerve canal measured by the CT threshold method[J]. Chinese Journal of Medical Imaging, 2016, 24(3): 175−178. DOI: 10.3969/j.issn.1005-5185.2016.03.004. (in Chinese).
[14] 郭佳, 沈文彬, 信建峰, 等. CT评价继发性下肢淋巴水肿分级的价值研究[J]. CT理论与应用研究, 2022,31(2): 237−244. DOI: 10.15953/j.ctta.2021.039. GUO J, SHEN W B, XIN J F, et al. CT imaging in secondary lower extremity lymphedema: A prospective study[J]. CT Theory and Applications, 2022, 31(2): 237−244. DOI: 10.15953/j.ctta.2021.039. (in Chinese).
[15] CUCCHI F, ROSSMEISLOVA L, SIMONSEN L, et al. A vicious circle in chronic lymphedema pathophysiology? An adipocentric view[J]. Obesity Reviews, 2017, 18(10): 1159−1169. DOI: 10.1111/obr.12565.
[16] CAO Q, YAN C, HAN X, et al. Quantitative evaluation of nonalcoholic fatty liver in rat by material decomposition techniques using rapid-switching dual energy CT[J]. Academic Radiology, 2021, S1076-6332(21): 00365−2. DOI: 10.1016/j.acra.2021.07.027.
[17] MENDONCA P R, LAMB P, SAHANI D V, et al. A flexible method for multi-material decomposition of dual-energy CT images[J]. IEEE IEEE Transactions on Medical Imaging, 2014, 33(1): 99−116. DOI: 10.1109/TMI.2013.2281719.
[18] WU X, HE P, LONG Z, et al. Multi-material decomposition of spectral CT images via fully convolutional densenets[J]. Journal of X-ray Science and Technology, 2019, 27(3): 461−471. DOI: 10.3233/XST-190500.
计量
- 文章访问数: 292
- HTML全文浏览量: 89
- PDF下载量: 17
