Correlation between Emphysema Volume and Small Shadow Density on CT and the Lung Function in Pneumoconiosis
-
摘要: 目的:比较尘肺患者肺气肿程度及肺部小阴影对肺功能的影响。方法:选择2015年至2020年间十堰市太和医院收治的明确诊断尘肺且合并肺气肿的患者96名,进行胸部CT检查及肺功能测定;采用64排螺旋CT(GE Optima CT 680)进行三维容量测量系统定量分析全肺、左右肺、各肺上中下肺野低衰减区容积百分比(LAV%),并根据全肺LAV% 将其分为A组:LAV%<15%,B组:LAV%15%~30%,C组:LAV%>30% 三组,分析肺气肿分布特点及与肺功能的相关性;将胸片评定小阴影密集度的标准应用于CT中得出CT图像的全肺、左右肺、上中下肺野小阴影密集度(SSD),并分析与肺功能的相关性。结果:全肺SSD与FEV1、TLC、PEF、MEF75%、MEF25% 呈负相关,与RV/TLC呈正相关;全肺、左肺、右肺、及上中下肺野LAV% 与FEV1、FEV1%、DLCO、PEF、MEF75%、MEF50%、MEF25% 呈负相关,与RV、RV/TLC% 呈正相关;当LAV%>30% 时,即C组LAV% 与肺功能FEV1、FEV1%、DLCO、RV/TLC、PEF、MEF75%、MEF50%、MEF25% 间的相关性较A组和B组好;肺部小阴影密集度与肺通气功能、容积参数有一定相关性,但与弥散功能相关性较弱。结论:尘肺患者肺气肿程度及小阴影密集度与肺功能均有一定相关性,但肺气肿较小阴影密集度与肺功能相关性更好,且能在一定程度上评估尘肺患者肺损伤程度,值得临床推广应用。Abstract: Objective: To compare the influence of emphysema and lung small shadow on the lung function of patients with pneumoconiosis. Methods: We selected 96 patients who were admitted to the Taihe Hospital (Shiyan city) for pneumoconiosis complicated by emphysema between 2015 and 2020. Chest computed tomography (CT) and lung function assessments were performed. The low attenuation volumes (LAV%) of the whole lung; left and right lungs; and upper, middle, and lower lobes of each lung were quantified using 64-slice spiral CT (GE Optima CT680). Based on the LAV% of the whole lung, the patients were divided into groups A (LAV%<15%), B (LAV%: 15%–30%), and C (LAV% >30%). The small shadow density (SSD) of the whole lung; left and right lungs; and upper, middle, and lower lung fields was calculated by applying the criteria for evaluating SSD on chest radiographs to CT images. The correlation between SSD and lung function was also analyzed. Results: The SSD of the whole lung was negatively correlated with the forced expiratory volume in the first second (FEV1), total lung capacity (TLC), peak expiratory flow (PEF), maximum expiratory flow after 75% of the forced vital capacity (FVC) has not been exhaled (MEF75%), and maximum expiratory flow after 25% of the FVC has not been exhaled (MEF25%), and it was positively correlated with residual volume (RV)/TLC. LAV% was negatively correlated with FEV1, ratio of FEV to FEV1 (FEV1%), carbon monoxide diffusion capacity (DLCO), PEF, MEF75%, maximum expiratory flow after 50% of the FVC has not been exhaled (MEF50%), and MEF25% in the whole lung; left and right lungs; and upper, middle, and lower lung lobes; however, it was positively correlated with RV and RV/TLC%. When LAV% was >30%, the correlation between LAV% and FEV1, FEV1%, DLCO, RV/TLC, PEF, MEF75%, MEF50%, and MEF25% was stronger in group C than in groups A and B. A certain degree of correlation was observed between the SSD in the lungs and the ventilation function and volume parameters; however, the correlation between SSD and the diffusion function was weak. Conclusion: The degree of emphysema and SSD in patients with pneumoconiosis have a certain degree of correlation with the lung function; however, the SSD of emphysema has a relatively better correlation with the lung function. It can be used to evaluate the degree of lung injury in patients with pneumoconiosis to a certain extent and is worthy of clinical application.
-
Keywords:
- CT /
- pneumoconiosis /
- emphysema /
- lung function
-
X射线计算机断层成像(computed tomography,CT)因其独特的三维层析能力与密度分辨力优势,在医学诊断、手术导航、精准医疗等方面得到了广泛的应用[1]。然而,CT成像中X射线产生的电离辐射可能诱发基因突变,增大了癌症等疾病的风险[2]。自CT成像技术问世以来,研究人员一直致力于降低成像所需射线的剂量,主要方法包括高灵敏度探测技术、高效率扫描机制、新型信号处理技术等[3-5]。得益于上述技术的进步,目前主流的胸部CT的辐射剂量约为7 mSv,腹盆CT增强检查的辐射剂量约为20 mSv,均低于100 mSv/年的安全阈值[1]。另一方面,X射线辐射的随机性效应决定了癌症的患病概率会随着辐射剂量的增加而增加,因此,进一步降低成像所需的最小剂量仍然是CT成像技术的重点问题[6]。
近年来,随着新型光电成像技术与信号处理理论的发展,Brady小组[7]在2009年首先提出了编码孔径CT成像的概念,通过放置于X射线源前的编码掩膜板降低CT成像投影之间的结构性冗余,结合信号稀疏测量理论,可以有效提高采样效率,降低辐射剂量。该文分析4种不完全采样条件下的X射线CT成像策略,实验结果表明,在X射线辐射剂量相同的条件下,随机分布探测器的采样策略在稀疏测量与重建的理论框架内性能最优。相比于传统扇形扫描和锥形束扫描CT成像,在重建质量相同的条件下,可降低90% 以上的辐射剂量;而相比于目前辐射剂量较低的CT成像方法,例如基于稀疏约束的有限(稀疏)角度CT成像方法,在重建质量相同的条件下,可降低30%~50% 的辐射剂量[8-12]。
2017年Mojica等[13]提出了超分辨率编码孔径X射线CT成像,利用低分辨率的X射线探测器测量经高分辨率编码掩膜板调制的X射线投影,结合压缩感知重建算法,从低分辨率、高信噪比的测量中恢复出高分辨率、高信噪比的三维CT图像。
2019年Cuadros等[14]提出了编码孔径X射线光谱CT成像,该编码孔径CT成像系统利用X射线K边缘能谱特性,通过设计合适的编码掩膜板产生特定的空间-光谱X射线投影,实现对四维目标(三维空间+光谱)的线性多路复用测量,并通过交替方向乘子法重建出四维断层图像(三维空间+光谱)。2021年,Cuadros等[15-16]首次通过实验验证了编码孔径X射线CT成像,在3D打印形成的0-1随机分布的二值图案中,使用钨粉末填充非闭合元素,并通过环氧树脂密封,最终形成X射线波段的编码掩膜板。实验表明,放置于X射线源之前的0-1分布的编码掩膜板并不会导致金属伪影,但X射线经过编码掩膜板会产生一定的散射伪影,可以通过掩膜校正抑制该伪影[17]。
此外,现阶段编码孔径CT成像重建算法的研究主要集中于通过优化编码掩膜板中非闭合元素的分布,结合稀疏测量与重建理论,通过迭代运算,提高重建质量。其主要方法包括基于有限等距性质(restricted isometry property,RIP)的随机梯度下降算法[18],基于均匀采样(uniform sensing)策略的优化算法[19-20]和基于稀疏主成分分析(sparse principle component analysis,SPCA)的快速搜索算法[21]。
虽然上述工作大幅度降低了CT成像所需的辐射剂量,但极度欠定采样的策略与目前的主流的解析类CT重建框架并不兼容,重建与优化均是建立在代数类重建框架内,算法复杂度高,重建时间长,难以满足临床需求。本文提出基于U型生成对抗网络的编码孔径CT成像方法,本方法通过构建编码孔径CT成像的非连续性投影数据动态博弈模型,利用医学图像的高维稀疏特性,预测欠定采样策略下非连续性CT投影的缺失数据,使其与解析类CT重建框架兼容,实现高质量的快速CT图像重建。
1. 原理与方法
1.1 编码孔径CT成像
编码孔径CT成像的光学结构与传统CT成像类似,这里我们以扇形X射线束照明,平板型探测器测量的扇形编码孔径CT成像结构为例,其光学几何结构如图1所示。编码掩膜板
$ {\lambda _i} \in {\mathbb{R}^M} $ ,$ i = 1,2, \cdots ,M $ ,放置于扇形X射线光源与目标之间,并且编码掩膜板上的像素与线性探测器上的像素一一对应。其中,M是编码掩膜板中元素的数量。经过编码掩膜板调制与目标吸收衰减的 X射线投影被平板型探测器测量,探测器每个像素所测得的X射线强度是其与X射线光源所构成直线上所有点的衰减的积分;X射线光源和对应的探测器绕着待测目标旋转P个角度,$ {S_1} $ ,$ {S_2} $ $, \cdots,$ $ {S_P} $ ,测量得到不同角度下的目标X射线投影。根据Beer-Lamber定理[22],X射线在路径l的衰减可以描述为:$$ {I^l} = I_0^l\exp \left( { - \int\limits_0^{\; \propto } {{\lambda _i}\alpha (l)} {\rm{d}}l} \right) + I_n^l , $$ (1) 这里,
$ {I^l} $ 是探测器的测量值,$ I_0^l $ 是单色X射线光源的强度,$ I_n^l $ 是X射线CT成像的系统噪声,α是沿路径l的线性衰减系数。式(1)通过对数运算可以变形为线性方程组,以矩阵形式可以表示为:
$$ {\boldsymbol{y}} = {\boldsymbol{CWf}} 。 $$ (2) 详细推导过程可以参考本课题组前期工作[23]。其中,
${\boldsymbol{y}} \in {\mathbb{R}^{MP}}$ 是向量化的通过对数运算得到的测量值,其可以通过探测器的测量值和X射线的强度计算得到;$ W \in {\mathbb{R}^{MP \times N}} $ 是结构矩阵,${\boldsymbol{f}} \in {\mathbb{R}^N}$ 是向量化的线性衰减系数。具体的,编码矩阵可以表示为:$$ {\boldsymbol{C}} = \left( {\begin{array}{*{20}{c}} {{\text{diag(}}{{\boldsymbol{c}}_1}{\text{)}}}&0& \cdots &0 \\ 0&{{\text{diag(}}{{\boldsymbol{c}}_2}{\text{)}}}& \cdots &0 \\ \vdots & \vdots & \ddots & \vdots \\ 0&0& \cdots &{{\text{diag(}}{{\boldsymbol{c}}_p}{\text{)}}} \end{array}} \right) \text{,} $$ (3) 其中,编码掩膜板
${{\boldsymbol{c}}_i} \in {\mathbb{R}^M},i = 1,2, \cdots ,P$ 。从式(3)可以看出,X射线编码孔径CT成像是对X射线CT成像在空域的0-1随机采样,根据傅里叶切片定理,即物体
$ r(x,y) $ 在方向$ (t,\theta ) $ 上的投影$ P(t,\theta ) $ 的一维傅里叶变换变为:$$ F\left\{ {c(t,\theta )r(t,\theta )} \right\} = F\left\{ {c(t,\theta )} \right\} * F\left\{ {r(t,\theta )} \right\} \text{,} $$ (4) 其中,
$ c(t,\theta ) $ 是方向$ (t,\theta ) $ 上的编码掩膜板分布。从式(4)可以看出,经过编码掩膜板,投影变换的频域混叠严重,并不满足解析类重建算法的基本要求。同时,X射线编码孔径CT成像中的编码矩阵C的对角向量是一个 0-1向量,其采样频率远小于传统X射线CT成像乃至有限角度X射线CT成像。因此,FBP算法等解析重建算法无法直接应用于X射线编码孔径CT成像。
1.2. 基于U型生成对抗网络的正弦图重建
如图2所示,与传统CT成像不同,编码孔径的引入使得其正弦图具有两大特征:非连续性和稀疏性。本文提出一种基于U型生成对抗网络的非连续正弦图重建模型,以解决在不连续、稀疏条件下编码孔径正弦图的反演问题。该方法的核心思想是生成对抗网络[24],包括训练模型和测试模型两个模块。在训练模型中,构建U型编解码全卷积网络作为训练模型的生成器,将编码孔径CT成像系统的随机投影数据作为输入图像对生成器网络进行训练。
如图3所示,生成器网络的主要作用是对输入的非连续投影数数据进行预测,随后用预测的投影数据结果和完整投影数据对鉴别网络进行训练,鉴别器网络的作用是估计生成器网络所生成的投影数据是否为真实的概率。同时在整个训练过程加入联合损失函数以保证编码孔径投影数据和完整投影数据的一致性。在测试阶段,主要是将新获得的编码孔径随机投影数据输入至测试阶段训练好的生成器中得到预测的非连续的投影数据,从而得到预测的完整投影数据。
1.3 U型生成对抗网络
本文拟采用如图4所示的U型生成对抗网络,包括生成网络(图4(a)所示)和鉴别网络(图4(b))所示。生成网络分为主干特征提取网络、加强特征提取网络和预测网络3个部分。生成网络的输入为一系列大小为512×512大小的投影数据;主干特征提取网络分为5层,由1层输入块C1和4层下采样块C2~C5组成。输入块C1对输入图像进行2次3×3的卷积、归一化和整流线性单位,进而生成特征图谱。每个下采样块C2~C5由两组3×3卷积层-批标准化层-ReLU激活函数层和一组2×2的最大池化层组成,增加特征通道的数量,以便获得4个初步有效特征层,在图中依次表示为S1至S4。同时,下采样模块可以使得网络更具有鲁棒性,抑制过拟合发生的概率。
加强特征提取网络由4个上采样块组成。每个上采样块由两组3×3卷积层-批标准化层-ReLU激活函数层和一组最近邻采样层组成,其作用是将下采样降级的图像特征还原至原始图像大小。在加强特征提取部分中,利用主干特征提取获得的初步有效特征层S1至S4,进行4次上采样块操作,并进行特征融合。最后,在预测网络部分利用一个1×1的卷积进行通道调整,将获得的最后一个有效特征层对每一个特征点进行分类,获得预测结果。
鉴别器网络结构与生成器结构相近,该结构由5个卷积层C6~C10以及1个输出层C11组成。其中,每个卷积层包括卷积运算-批标准化层-LeakyReLU激活函数层。为了避免棋盘效应以及过拟合现象,每层卷积核的大小设定为4×4,横轴与纵轴方向步长均设置为2个像素值,填充数为1个像素。鉴别器将生成器预测出的结果及对应的真实数据作为输入,其输出是预测出来的结果被鉴别为真实数据的概率,值为0到1之间,值越接近1,代表预测出的数据越接近为真实数据。
1.4 联合损失函数
在训练过程中,对于生成器输入的编码正弦图,其重建损失函数Lrec可以为:
$$ {L_{{\rm{rec}}}}\Big( {G(p \odot M),\,p} \Big) = {\Big\| {\varphi \Big( {G(p \odot M)} \Big) - \varphi (p)} \Big\|_2} \text{,} $$ (5) 其中,G表示训练时使用的生成器,p表示完整的正弦图,M表示编码掩膜板,⊙表示哈达玛积,φ为使用基于 Image-Net训练权重的VGG16网络第4组卷积网络的第3层输出特征图。同时,本文采取Wasserstein距离作为其的对抗性损失函数
$$ {L_G}(p \odot M,\,p) = \mathop {\max }\limits_G \bigg( {{{E}}\Big( {{{D}}\big( {{G}(p \odot M)} \big)} \Big)} \bigg) \text{,} $$ (6) $$ {L_D}(p \odot M,\,p) = \mathop {\max }\limits_D \bigg( {E\Big( {D(p)} \Big) - E\Big( {D\big(G(p \otimes M)\big)} \Big)} \bigg) + \delta {\Big( {{{\left\| {{\nabla _{\hat p}}D(\hat p)} \right\|}_2} - 1} \Big)^2} \text{,} $$ (7) 其中,
$ \hat p $ 表示随机噪声,D表示鉴别器,δ表示梯度惩罚。损失函数LG的主要目的是使鉴别器对合成数据判定概率尽可能的接近 1,提高生成器输出图像的质量。LD的主要目的是提高鉴别器对真实数据判定概率并降低预测数据的判定概率,以此来鉴别真实数据和预测数据。通过生成器与鉴别器训练过程中的互相博弈,使得预测图像的数据分布接近于真实图像的数据分布。2. 实验结果
本文采用课题组前期工作的编码孔径CT成像的光学几何结构[20,22],CT图像分辨率为256×256,探测器分辨率为512×1,扫描次数为512。因此,正弦图的空间分辨率为512×512。本文选用北美放射协会(Radiological Society of North America,RSNA)的2500张颅内出血CT图像作为数据集并将其空间分辨率缩放至256×256。其中,选取2000张用于测试CT图像及其对应的编码正弦图进行训练,500张CT图像及其对应的编码正弦图。
训练参数上的设置:Batch Size大小4,训练epoch数115,学习率0.0001。训练平台intel i9 9900 k,GeForce RTX 2080 Ti,128 GB RAM。由于在神经网络的训练过程中,当预测图像的准确度出现震荡或loss不再下降时,进行适当的学习率衰减可以有效提高准确率。因此,本文在训练过程中加入学习率衰减函数且在训练15个epoch后执行学习率衰减函数,可以使训练结果达到最佳。由于经过U型生成对抗网络的编码孔径CT成像的正弦图在空间域上是连续的,因此,本文可以采用滤波反投影算法进行重建。相比于代数类重建算法和稀疏约束重建算法,滤波反投影算法在内存需求上可以降低2~3个数量级,重建算法复杂度可以降低2~3个数量级。
图5表示编码掩膜板透过率为20%,训练次数为115个epoch的训练结果图。其中,图5(a)~图5(c)表示对抗性损失函数曲线图,在经过80 k的迭代之后,损失函数逐渐收敛,趋于稳定。图5(d)~图5(g)表示10% 透过率下每个epoch训练和测试的峰值信噪比(peak signal-to-noise ratio,PSNR)和结构相似性(structural similarity,SSIM),从图中可以看出,测试集重建图像质量略有下降,但整体质量较好(PSNR>40 dB,SSIM>0.98)。其中,利用训练好的网络从编码正弦图反演完整的正弦图,运算时间为1~2 s,滤波反投影算法运算时间为1~2 s。而同等分辨率条件下,编码孔径CT成像方法迭代反演方法需要十多分钟。因此,相比于目前已有的编码孔径CT成像方法,本文所提出的方法可以将重建时间降低约两个数量级。
图6表示编码掩膜板透过率为25% 时,基于U型生成对抗网络的编码孔径CT成像的结果图。其中(a)为原始图像,(b)编码掩膜板,(c)为编码孔径CT成像的正弦图,(d)是经过U型生成对抗网络反演预测的正弦图,(e)为预测的正弦图(d)经FBP重建算法反演的CT图像。这里,经网络预测的正弦图与真实正弦图的PSNR分别为35.2 、46.5 和46.3 dB,SSIM分别为0.98、0.98和0.99;而根据预测的正弦图经FBP重建算法反演的CT图像与原始图像的PSNR分别为41.6 、38.3 和36.1 dB,SSIM分别为0.96、0.94和0.90。实验结果表明,基于U型生成对抗网络的编码孔径CT成像可以在兼顾大幅度降低辐射剂量(仅25% 的透过率)和保证重建效率(预测正弦图<2 s)的情况下,高质量的重建出CT图像(PSNR>35 dB,SSIM>0.9)。
图7表示编码掩膜板透过率为5%、15% 以及25% 时,基于U型生成对抗网络的编码孔径CT成像的结果图。其中(a)为原始图像,(b)是透过率为5%、15% 和25% 的编码掩膜板,(c)为编码孔径CT成像的正弦图,(d)是经过U型生成对抗网络反演预测的正弦图,(e)为预测的正弦图(d)经FBP重建算法反演的CT图像。这里,经网络预测的正弦图与真实正弦图的PSNR分别为40.8 、41.9 和45.2 dB,SSIM分别为0.99、0.98和0.98;而根据预测的正弦图经FBP重建算法反演的CT图像与原始图像的PSNR分别为38.7 、41.3 和41.6 dB,SSIM分别为0.95、0.96和0.96。实验结果表明,基于U型生成对抗网络的编码孔径CT成像方法鲁棒性较强,在编码掩膜板透过率5% 时仍可以较高质量的重建CT图像。
表1表示在不同编码掩膜板透过率下,测试集中预测的正弦图和CT重建图像的平均值。其中,编码掩膜板透过率5%~25% 的CT重建图像平均PSNR仅相差3.3 dB @ PSNR>30 dB,平均SSIM仅相差0.03 @ SSIM>0.87。因此,基于U型生成对抗网络的编码孔径CT成像方法对于编码掩膜板透过率的变化鲁棒性较好,在实际应用中可以尽可能的降低编码掩膜板的透过率,进一步降低辐射剂量。
表 1 不同编码掩膜板透过率下,测试集中预测的正弦图和CT重建图像的平均值Table 1. The average results of sonograms and CT reconstructions at different transmission rates测试数据集 图像质量评价指标 编码掩膜板透过率 5% 10% 15% 20% 25% 预测的正弦图 PSNR/dB 39.1 41.6 40.9 42.6 43.8 SSIM 0.98 0.97 0.98 0.99 0.98 CT重建图像 PSNR/dB 33.7 32.0 33.5 35.3 33.6 SSIM 0.88 0.88 0.89 0.91 0.90 3. 结论
本文提出一种基于U型生成对抗网络的编码孔径CT成像方法。通过构建基于U型生成对抗网络的非连续稀疏投影的动态博弈模型,结合联合损失函数,预测正弦图的结构性缺失,实现编码孔径CT成像分析类(非迭代)快速重建。实验结果表明,相比于编码孔径CT成像目前使用的迭代类算法,基于U型生成对抗网络的编码孔径CT成像方法可以将重建时间降低约两个数量级。同时,在辐射剂量降低95% 的条件下,仍然可以获得峰值信噪比大于30 dB @ 256×256的高质量重建结果。该方法对于下一代低剂量CT成像设备的研究具有一定的理论基础和参考依据。
-
图 1 分别为A、B、C组部分尘肺患者肺部CT冠状位图像及轴位图像,冠状位图像由气管纵向水平上用等分点的水平线将左右肺各分为上、中、下3个肺野,蓝色区域为肺气肿低衰减区三维图像
Figure 1. Coronal and axial CT images of some patients with pneumoconiosis in groups A, B, and C. The left and right lungs were divided into upper, middle, and lower lung fields from the vertical level of the trachea using equidistant horizontal lines; blue areas present emphysema as three-dimensional images of low-attenuation areas
表 1 各肺野LAV% 数据表
Table 1 LAV% data for each lung field
数据范围 右肺 左肺 右肺 上野 中野 下野 左肺 上野 中野 下野 -1024~-950/% 24.030 24.120 22.428 25.237 23.080 23.053 19.105 26.269 -950~3071/% 75.970 75.881 77.572 74.763 76.920 76.947 80.895 73.731 Total/L 4.300 0.984 1.456 1.860 3.384 0.843 1.128 1.413 表 2 A、B、C组组间总体SSD间差异性分析
Table 2 Comparison of the total SSD among groups A, B, and C
项目 SSD A B B C A C M(Q25,Q75) 2(1,3) 2(2,3) 2(2,3) 2(1,3) 2(1,3) 2(1,3) Z -0.454 -0.373 -0.537 P 0.650 0.767 0.591 注:M(Q25,Q75)为A、B、C三组总体小阴影密集度中位数及四分位数,组间对比不符合正态分布采用非参数检验得到P>0.05,
3组总体小阴影密集度对比差异无统计学意义。表 3 A组各肺野SSD、LAV% 与肺功能指标的相关性分析
Table 3 Analysis of correlation between SSD and LAV% of the lung field and lung function indexes in group A
A组 右上 右中 右下 左上 左中 左下 SSD LAV% SSD LAV% SSD LAV% SSD LAV% SSD LAV% SSD LAV% FEV1 -0.099 0.019 -0.157 -0.093 -0.330* -0.073 -0.228 0.132 -0.317* -0.203 -0.172 0.065 FEV1% 0.159 -0.068 0.037 0.060 -0.124 -0.159 -0.037 0.005 -0.085 -0.337* -0.076 0.082 FVC -0.258 0.020 -0.374* -0.034 -0.388* -0.047 -0.296 0.084 -0.501** -0.231 -0.385* 0.123 DLCO 0.035 -0.022 -0.007 -0.055 0.053 0.119 -0.077 -0.091 0.011 -0.233 -0.079 0.268 RV -0.079 0.192 -0.124 -0.048 0.027 0.070 -0.036 0.194 -0.082 0.082 -0.019 -0.060 TLC -0.004 0.256 -0.030 -0.030 -0.152 0.080 -0.010 0.386* -0.082 0.262 -0.002 0.008 RV/TLC 0.140 0.089 0.058 -0.199 0.187 -0.162 0.176 -0.025 0.164 0.179 0.171 -0.148 PEF 0.020 -0.034 -0.144 -0.145 -0.267 -0.083 -0.091 0.134 -0.250 -0.253 -0.122 0.057 MEF75% -0.055 -0.059 -0.055 -0.127 -0.262 -0.250 -0.061 0.048 -0.225 -0.322* -0.106 0.097 MEF50% -0.057 -0.018 -0.057 -0.032 -0.249 -0.189 -0.110 0.048 -0.213 -0.310* -0.165 0.114 MEF25% -0.116 0.040 -0.116 0.218 -0.251 0.107 -0.095 0.121 -0.227 -0.253 -0.286 0.106 注:* 为P<0.05,**为P<0.01。 表 4 B组各肺野SSD、LAV% 与肺功能指标的相关性分析
Table 4 Analysis of correlation between SSD and LAV% of the lung field and lung function indexes in group B
B组 右上 右中 右下 左上 左中 左下 SSD LAV% SSD LAV% SSD LAV% SSD LAV% SSD LAV% SSD LAV% FEV1 -0.053 0.247 -0.149 -0.036 -0.096 -0.314 -0.123 0.102 -0.009 -0.211 -0.080 -0.521* FEV1% -0.024 0.314 -0.190 0.025 -0.075 -0.304 -0.074 0.128 -0.104 -0.133 -0.147 -0.615* FVC -0.179 -0.078 -0.263 -0.177 -0.079 -0.133 -0.209 -0.038 -0.009 -0.111 -0.161 -0.407* DLCO 0.062 0.157 -0.365 -0.269 -0.154 0.020 0.041 -0.025 0.239 -0.227 -0.288 -0.224 RV -0.170 -0.095 -0.024 -0.112 0.048 0.035 -0.227 -0.025 -0.069 -0.138 0.009 0.009 TLC -0.149 0.155 -0.056 -0.020 -0.268 -0.270 -0.219 0.196 0.147 -0.038 -0.041 -0.113 RV/TLC 0.170 -0.092 0.187 -0.179 0.191 0.450* 0.118 -0.148 0.039 0.051 0.120 0.321 PEF 0.085 0.215 -0.153 -0.015 -0.235 -0.152 -0.044 0.060 0.028 -0.121 -0.140 -0.280 MEF75% 0.072 0.266 -0.055 0.084 -0.002 -0.155 -0.033 0.060 0.021 -0.133 0.013 -0.478* MEF50% 0.048 0.201 -0.067 -0.112 0.015 -0.243 -0.018 0.037 -0.004 -0.259 -0.043 -0.533* MEF25% -0.152 0.199 -0.214 0.024 -0.162 -0.180 -0.238 0.012 -0.176 -0.201 -0.174 -0.479* 注:* 为P<0.05。 表 5 C组各肺野SSD、LAV% 与肺功能指标的相关性分析
Table 5 Analysis of correlation between SSD and LAV% of the lung field and lung function indexes in group C
C组 右上 右中 右下 左上 左中 左下 SSD LAV% SSD LAV% SSD LAV% SSD LAV% SSD LAV% SSD LAV% FEV1 0.324 -0.453* 0.242 -0.435* -0.126 -0.555* 0.030 -0.575* 0.053 -0.554* 0.149 -0.519* FEV1% 0.302 -0.513* 0.209 -0.420* -0.007 -0.532* 0.076 -0.552* 0.068 -0.581* 0.205 -0.389* FVC -0.138 0.010 -0.181 -0.086 -0.284 -0.205 0.036 -0.137 -0.042 -0.156 -0.102 -0.291 DLCO 0.338 -0.282 0.511* -0.506* 0.420* -0.503* 0.395 -0.326 0.492* -0.523* 0.648* -0.410* RV -0.358 0.065 -0.315 0.221 0.260 0.347 -0.191 0.268 -0.199 0.490* -0.563* 0.325 TLC -0.033 -0.285 0.043 0.085 -0.362 -0.110 -0.248 -0.414* -0.136 -0.203 -0.136 -0.233 RV/TLC 0.150 0.042 0.048 0.130 -0.134 0.211 0.242 0.151 0.197 0.299 0.115 0.332 PEF 0.250 -0.535* 0.069 -0.324 -0.211 -0.469* -0.111 -0.575* -0.069 -0.476* 0.003 -0.530* MEF75% 0.161 -0.498* -0.036 -0.297 -0.266 -0.423* -0.199 -0.529* -0.178 -0.395 -0.053 -0.496* MEF50% 0.136 -0.376 -0.077 -0.229 -0.338 -0.329 -0.201 -0.367 -0.153 -0.286 -0.139 -0.431* MEF25% 0.040 -0.312 -0.031 -0.202 -0.377 -0.298 -0.221 -0.306 -0.114 -0.239 -0.141 -0.404* 注:* 为P<0.05。 -
[1] 韩欢, 徐绍德. 胸部多排螺旋CT扫描诊断职业性尘肺病的效果及准确率分析[J]. 影像研究与医学应用, 2021,5(8): 64−65. doi: 10.3969/j.issn.2096-3807.2021.08.030 HAN H, XU S D. Analysis of effect and accuracy of multi-slice spiral CT in the diagnosis of occupational pneumoconiosis[J]. Imaging Research and Medical Applications, 2021, 5(8): 64−65. (in Chinese). doi: 10.3969/j.issn.2096-3807.2021.08.030
[2] JONES C M, PASRICHA S S, HEINZE S B, et al. Silicosis in artificial stone workers: Spectrum of radiological high-resolution CT chest findings[J]. Journal Medical Imaging and Radiation Oncology, 2020, 64(2): 241−249. doi: 10.1111/1754-9485.13015
[3] PEREZ T, GARCIA G, ROCHE N, et al. Société de pneumologie de langue française. Recommandation pour la pratique clinique. Prise en charge de la BPCO. Mise à jour 2012. Exploration fonctionnelle respiratoire. Texte long [French Pulmonary Medicine Society. Guidelines for clinical practice. Management of COPD. Update 2012: Pulmonary function tests: Full length text][J]. Revue Des Maladies Respiratoires, 2014, 31(3): 263−294. doi: 10.1016/j.rmr.2013.11.007
[4] 刘阿茹, 魏华, 邓永红. HRCT肺气肿定量分析与COPD患者疾病相关性分析[J]. 中国CT和MRI杂志, 2021,19(8): 74−76. DOI: 10.3969/j.issn.1672-5131.2021.08.024. LIU A R, WEI H, DENG Y H. Analysis of correlation between quantitative analysisi of HRCT emphysema and COPD patients[J]. Chinese Journal of CT and MRI, 2021, 19(8): 74−76. DOI: 10.3969/j.issn.1672-5131.2021.08.024. (in Chinese).
[5] 卞明敏, 胡茂能. 职业性尘肺病肺部影像表现与肺功能改变关系的研究[J]. 职业卫生与应急救援, 2021,39(5): 519−523. DOI: 10.16369/j.oher.issn.1007-1326.2021.05.008. BIAN M M, HU M N. Study on the relationship pulmonary imaging manifestations and pulmonary function changes in occupational pnuemoconiosis[J]. Occupational Health and Emergency Rescue, 2021, 39(5): 519−523. DOI: 10.16369/j.oher.issn.1007-1326.2021.05.008. (in Chinese).
[6] ZAPOROZHAN J, LEY S, EBERHARDT R, et al. Paired inspiratory/expiratory volumetric thin-slice CT scan for emphysema analysis: Comparison of different quantitative evaluations and pulmonary function test[J]. Chest, 2005, 128(5): 3213−3220.
[7] 张雨洁, 俞同福. COPD患者肺气肿程度对气道径线与肺功能相关性的影响[J]. 实用放射学杂志, 2016,32(7): 1028−1032. DOI: 10.3969/j.issn.1002-1671.2016.07.009. ZHANG Y J, YU T F. Infiuence of the degree of emphysema on the correlation between airway line and lung function in COPD patients[J]. Journal of Applied Radiology, 2016, 32(7): 1028−1032. DOI: 10.3969/j.issn.1002-1671.2016.07.009. (in Chinese).
[8] VOGELMEIER C F, CRINER G J, MARTINEZ F J, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive lung disease 2017 report: GOLD executive summary[J]. American Journal of Respiratory and Critical Care Medicine, 2017, 49(6): 557−582.
[9] GBZ70-2015, 职业性尘肺病的诊断[S]. [10] 翟荣存, 李年春, 刘晓东, 等. 尘肺病小阴影密集度的CT分级方法及参考片探讨[J]. 中华劳动卫生职业病杂志, 2021,39(6): 453−457. DOI: 10.3760/cma.j.cn121094-20200917-00537. ZHAI R C, LI N C, LIU X D, et al. Discussion on CT grading method of small shadow density in pneumoconiosis and refernce film[J]. Chinese Journal of Labor Health and Occupational Diseases, 2021, 39(6): 453−457. DOI: 10.3760/cma.j.cn121094-20200917-00537. (in Chinese).
[11] MASHIMO S, CHUBACHI S, TSUTSUMI A, et al. Relationship between diminution of small pulmonary vessels and emphysema in chronic obstructive pulmonary disease[J]. Clinical Imaging, 2017, 46: 85−90. doi: 10.1016/j.clinimag.2017.07.008
[12] 余先超, 王朗, 温晓玲, 等. CT肺容积及肺密度与尘肺肺功能的相关性[J]. 中国医学影像学杂志, 2020,28(11): 846−848, 851. DOI: 10.3969/j.issn.1005-5185.2020.11.012. YU X C, WANG L, WEN X L, et al. Correlation of lung volume and lung density with pneumoconiosis[J]. Chinese Journal of Medical Imaging, 2020, 28(11): 846−848, 851. DOI: 10.3969/j.issn.1005-5185.2020.11.012. (in Chinese).
[13] 王雷, 沈聪, 邹常咏, 等. 基于CT定量评估肺气肿肺叶分布特点与肺功能的相关性[J]. 西安交通大学学报(医学版), 2018,39(6): 902−906. DOI: 10.7652/jdyxb201806026. WANG L, SHEN C, ZOU C Y, et al. Quantitative evealuation of the correlation between pulmonary lobe distribution characteristics and lung function in emphysema based on CT[J]. Journal of Xi'an Jiaotong University (Medical Edition), 2018, 39(6): 902−906. DOI: 10.7652/jdyxb201806026. (in Chinese).
[14] 兰长青, 王洁, 黄梅萍, 等. CT肺气肿指数及空间分布对慢性阻塞性肺疾病肺损伤的评估价值[J]. 中国临床医学影像杂志, 2021,32(5): 320−325. DOI: 10.12117/jccmi.2021.05.004. LAN C Q, WANG J, HUANG M P, et al. Evaluation of CT emphysema index and its spatial distribution in patients with chronic obstructive pulmonary disease[J]. Chinese Journal of Clinical Medical Imaging, 2021, 32(5): 320−325. DOI: 10.12117/jccmi.2021.05.004. (in Chinese).
[15] 王岩. 探析多层螺旋CT及后处理功能对早期尘肺诊断的应用价值[J]. 中国医疗器械信息, 2020,26(8): 151−152. DOI: 10.15971/j.cnki.cmdi.2020.08.076. WANG Y. Application of multislice spiral CT and post-processing in early diagnosis of pneumoconiosisi[J]. China Medical Device Information, 2020, 26(8): 151−152. DOI: 10.15971/j.cnki.cmdi.2020.08.076. (in Chinese).
[16] 胡必锋, 朱胜康, 翟荣存, 等. CT小阴影密集度判定方法及参考片在矽肺诊断中的应用[J]. 中华放射学杂志, 2021,55(11): 1172−1177. DOI: 10.3760/cma.j.cn112149-20210201-00085. HU B F, ZHU S K, ZHAI R C, et al. Determination of small shadow intensity in CT and application of reference film in diagnosis of silicosis[J]. Chinese Journal of Radiology, 2021, 55(11): 1172−1177. DOI: 10.3760/cma.j.cn112149-20210201-00085. (in Chinese).
[17] 夏养萱, 张健杰, 朱德香, 等. 尘肺病患者胸部螺旋CT影像与肺功能变化研究[J]. 职业卫生与应急救援, 2016,30(4): 267−269, 274. DOI: 10.16369/j.oher.issn.1007-1326.2016.04.001. XIA Y X, ZHANG J J, ZHU D X, et al. Study on spiral CT imaging of chest and pulmonary function in patients with pneumoconiosis[J]. Occupational Health and Emergency Rescue, 2016, 30(4): 267−269, 274. DOI: 10.16369/j.oher.issn.1007-1326.2016.04.001. (in Chinese).
[18] CHAMPLIN J, EDWARDS R, PIPAVATH S. Imaging of occupational lung disease[J]. Radiologic Clinics of North America, 2016, 54(6): 1077−1096. doi: 10.1016/j.rcl.2016.05.015