CN113450399A - 一种正位胸片心胸比测量方法及装置 - Google Patents

一种正位胸片心胸比测量方法及装置 Download PDF

Info

Publication number
CN113450399A
CN113450399A CN202110590001.6A CN202110590001A CN113450399A CN 113450399 A CN113450399 A CN 113450399A CN 202110590001 A CN202110590001 A CN 202110590001A CN 113450399 A CN113450399 A CN 113450399A
Authority
CN
China
Prior art keywords
image
chest
heart
transverse diameter
lung field
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110590001.6A
Other languages
English (en)
Other versions
CN113450399B (zh
Inventor
贾潇
王子腾
代玉婷
王立威
丁佳
吕晨翀
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Yizhun Intelligent Technology Co.,Ltd.
Original Assignee
Beijing Yizhun Medical AI Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Yizhun Medical AI Co Ltd filed Critical Beijing Yizhun Medical AI Co Ltd
Priority to CN202110590001.6A priority Critical patent/CN113450399B/zh
Publication of CN113450399A publication Critical patent/CN113450399A/zh
Application granted granted Critical
Publication of CN113450399B publication Critical patent/CN113450399B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/12—Edge-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30048—Heart; Cardiac
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30061—Lung

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

本发明提供一种正位胸片心胸比测量方法及装置。所述方法包括：从输入胸片图像中分割出心脏轮廓、左肺野轮廓和右肺野轮廓；分别求出心脏轮廓横坐标的最大值和最小值，计算所述最大值和最小值的差，得到心脏横径；在右肺野轮廓的下部找到斜率绝对值最小的点，过此点作一水平线，计算所述水平线与左肺野轮廓、右肺野轮廓的外侧交点之间的线段长度，得到胸廓横径；用所述心脏横径除以所述胸廓横径得到心胸比。本发明通过求解斜率绝对值最小点找到右膈顶，过右膈顶作水平线计算胸廓横径，解决了现有技术利用包围左、右肺野轮廓的矩形竖边之间的距离得到的胸廓横径大于实际值，从而导致心胸比偏小的问题，提高了心胸比的计算精度。

Description

一种正位胸片心胸比测量方法及装置
技术领域
本发明属于医学影像技术领域，具体涉及一种正位胸片心胸比测量方法及装置。
背景技术
胸部X线检查是诊断肺部和心脏疾病的关键技术，是体检筛查肺部疾病的主要手段。临床中，放射科医生通常将心胸比作为评估心脏肥大的重要量化指标。心胸比是指正位胸片上心脏横径与胸廓横径的比值。心脏横径等于左、右心缘至体中线的最大距离之和，胸廓横径即右膈顶水平胸廓的内径。
目前的心胸比测量方法可以分为三类：第一类是手动测量，即在胸片上通过视觉定位左、右心缘上离体中线最远的点，然后手通测量这两个点之间的水平距离从而得到心脏横径；通过视觉定位右膈顶水平线，然后沿该水平线测量胸廓内径从而得到胸廓横径。最后再计算心胸比。该方法操作繁琐，误差大，效率低；第二类是利用传统的图像处理技术进行心肺区域的分割，然后利用得到的心肺区域的边缘来确定几个关键点，来计算心脏横径和胸廓横径，从而计算心胸比。这类算法由于采用了传统的图像处理技术，如阈值分割和区域生长等，往往十分依赖像素阈值的选择和计算，容易受到心肺的邻近组织的干扰，同时对图像的噪声十分敏感。因此，此类算法的适用性不高；第三类是利用深度学习技术先对心脏和肺部进行定位，然后再在相应的感兴趣区域内进行心肺区域的分割，最后利用分割出来的心肺区域的矩形包围盒来计算心脏横径和胸廓横径，从而计算心胸比。这类算法虽然能较好地分割出心肺区域，但是其算法设计复杂，且利用包围左、右肺野轮廓的矩形竖边之间的距离得到的胸廓横径大于实际值，从而导致心胸比偏小，尤其是在胸片图像有倾斜角或心肺区域的边界与相应的包围盒的边界有较大夹角的情况下，误差更大。
发明内容
为了解决现有技术中存在的上述问题，本发明提供一种正位胸片心胸比测量方法及装置。
为了实现上述目的，本发明采用以下技术方案。
第一方面，本发明提供一种正位胸片心胸比测量方法，包括以下步骤：
从输入胸片图像中分割出心脏轮廓、左肺野轮廓和右肺野轮廓；
分别求出心脏轮廓横坐标的最大值和最小值，计算所述最大值和最小值的差，得到心脏横径；
在右肺野轮廓的下部找到斜率绝对值最小的点，过此点作一水平线，计算所述水平线与左肺野轮廓、右肺野轮廓的外侧交点之间的线段长度，得到胸廓横径；
用所述心脏横径除以所述胸廓横径得到心胸比。
进一步地，所述方法还包括将输入胸片图像缩放到800*1333像素*像素，灰度值范围为[0,255]。
进一步地，所述方法还包括输入胸片图像倾斜校正步骤：将输入的胸片图像输入一个空间变换网络STN，基于STN的输出与模板图像的差异计算损失函数，利用所述损失函数对STN进行训练，通过优化网络参数(α,sx,sy,tx,ty)对输入胸片图像进行旋转、伸缩和平移变换，使其与模板图像对正；其中，α为旋转角度，sx、sy分别为x、y方向的缩放倍数，tx、ty分别x、y方向的平移量；所述模板图像通过对多个体中线为竖直线的标准正位胸片图像求均值获得。
更进一步地，从输入胸片图像到模板图像的坐标变换方程为：
Figure BDA0003088999650000021
即：
Figure BDA0003088999650000031
式中，T为坐标变换矩阵，(x,y)为输入胸片图像的像素点坐标，(x′,y′)为模板图像的像素点坐标。
进一步地，所述斜率绝对值的计算公式为：
Figure BDA0003088999650000032
式中，ki为右肺野轮廓上第i个像素点(xi,yi)处的斜率的绝对值，(xi-δ,yi-δ)、(xi+δ,yi+δ)分别为(xi,yi)前面和后面第δ个像素点的坐标，yi＞(ymax+ymin)/2，ymax、ymin分别为右肺野轮廓上纵坐标的最大值和最小值，坐标原点在图像的左上角，纵坐标y向下为正方向。
第二方面，本发明提供一种正位胸片心胸比测量装置，包括：
轮廓分割模块，用于从输入胸片图像中分割出心脏轮廓、左肺野轮廓和右肺野轮廓；
心脏横径计算模块，用于分别求出心脏轮廓横坐标的最大值和最小值，计算所述最大值和最小值的差，得到心脏横径；
胸廓横径计算模块，用于在右肺野轮廓的下部找到斜率绝对值最小的点，过此点作一水平线，计算所述水平线与左肺野轮廓、右肺野轮廓的外侧交点之间的线段长度，得到胸廓横径；
心胸比计算模块，用于用所述心脏横径除以所述胸廓横径得到心胸比。
进一步地，所述装置还包括将图像标准化模块，用于将输入胸片图像缩放到800*1333像素*像素，灰度值范围为[0,255]。
进一步地，所述装置还包括倾斜校正模块，用于将输入的胸片图像输入一个空间变换网络STN，基于STN的输出与模板图像的差异计算损失函数，利用所述损失函数对STN进行训练，通过优化网络参数(α,sx,sy,tx,ty)对输入胸片图像进行旋转、伸缩和平移变换，使其与模板图像对正；其中，α为旋转角度，sx、sy分别为x、y方向的缩放倍数，tx、ty分别x、y方向的平移量；所述模板图像通过对多个体中线为竖直线的标准正位胸片图像求均值获得。
更进一步地，从输入胸片图像到模板图像的坐标变换方程为：
Figure BDA0003088999650000041
即：
Figure BDA0003088999650000042
式中，T为坐标变换矩阵，(x,y)为输入胸片图像的像素点坐标，(x′,y′)为模板图像的像素点坐标。
进一步地，所述斜率绝对值的计算公式为：
Figure BDA0003088999650000043
式中，ki为右肺野轮廓上第i个像素点(xi,yi)处的斜率的绝对值，(xi-δ,yi-δ)、(xi+δ,yi+δ)分别为(xi,yi)前面和后面第δ个像素点的坐标，yi＞(ymax+ymin)/2，ymax、ymin分别为右肺野轮廓上纵坐标的最大值和最小值，坐标原点在图像的左上角，纵坐标y向下为正方向。
与现有技术相比，本发明具有以下有益效果。
本发明通过分别求出心脏轮廓横坐标的最大值和最小值，计算所述最大值和最小值的差得到心脏横径，在右肺野轮廓的下部找到斜率绝对值最小的点，过此点作一水平线，计算所述水平线与左肺野轮廓、右肺野轮廓的外侧交点之间的线段长度得到胸廓横径，用所述心脏横径除以所述胸廓横径得到心胸比，实现了心胸比的自动计算。本发明通过求解斜率绝对值最小点找到右膈顶，过右膈顶作水平线计算胸廓横径，解决了现有技术利用包围左、右肺野轮廓的矩形竖边之间的距离得到的胸廓横径大于实际值，从而导致心胸比偏小的问题，提高了心胸比的计算精度。
附图说明
图1为本发明的实施例一种正位胸片心胸比测量方法的流程图。
图2为心脏轮廓、左肺野轮廓和右肺野轮廓的示意图。
图3为心脏横径示意图。
图4为胸廓横径示意图。
图5为对正前后图像的示意图，(a)为对正前图像，(b)为对正后图像。
图6为STN网络的结构示意图。
图7为本发明的实施例一种正位胸片心胸比测量方法装置的方框图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚、明白，以下结合附图及具体实施方式对本发明作进一步说明。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例一种正位胸片心胸比测量方法的流程图，包括以下步骤：
步骤101，从输入胸片图像中分割出心脏轮廓、左肺野轮廓和右肺野轮廓；
步骤102，分别求出心脏轮廓横坐标的最大值和最小值，计算所述最大值和最小值的差，得到心脏横径；
步骤103，在右肺野轮廓的下部找到斜率绝对值最小的点，过此点作一水平线，计算所述水平线与左肺野轮廓、右肺野轮廓的外侧交点之间的线段长度，得到胸廓横径；
步骤104，用所述心脏横径除以所述胸廓横径得到心胸比。
本实施例给出一种正位胸片心胸比测量方法。心胸比是放射科医生用于评估心脏肥大的重要量化指标。心胸比是指正位胸片上心脏横径与胸廓横径的比值。心脏横径与胸廓横径是两个常用的影像名词，心脏横径等于左、右心缘至体中线的最大距离之和，胸廓横径即右膈顶水平胸廓的内径。本实施例通过分别计算心脏横径与胸廓横径，然后求二者的比值得到心胸比。
本实施例中，步骤101主要用于从输入胸片图像中分割出心脏轮廓、左肺野轮廓和右肺野轮廓。本实施例采用深度学习中的语义分割技术从胸片图像上自动分割出心脏和左、右肺野的轮廓。深度卷积神经网络能自动提取整合图像的纹理、形状等特征，避免了传统方法中出现的对像素灰度阈值的选择和依赖的情况。同时，由语义分割技术得到的物体轮廓边界光滑，很好地保持了物体边界的几何性质，尤其是对像心脏、肺野这类较大的物体。通常的语义分割网络，如Unet及其各种变体、FCN、Deeplab等，都可以用于本任务。经对比试验，本实施例采用Deeplab_v3网络，一次性分割出左、右肺野和心脏的轮廓。在训练分割网络时，采用了结合DICE损失和交叉熵(cross entropy)损失作为分割结果的损失函数。并且在计算交叉熵损失时，还采用了Hard pixel mining技术，将损失集中在预测误差较大的像素上(通常是边界处的像素)，这样可以让网络在物体边界上的预测效果更好。
本实施例中，步骤102主要用于计算心脏横径。心脏横径等于左、右心缘至体中线的最大距离之和，如图3所示。体中线正常情况下是竖直线，平行于纵轴y，因此心脏横径等于心脏轮廓横坐标的最大值和最小值的差。如果输入的胸片图像发生了倾斜，体中线将不再是竖直线，按照上述方法计算心脏横径将会产生偏差。因此，需要对胸片图像进行倾斜校正即对正。
本实施例中，步骤103主要用于计算胸廓横径。胸廓横径也就是右膈顶水平胸廓的内径，如图4所示。由于正位胸片上人体两侧胸廓的边界并不总是竖直的，因此使用包围左、右肺野轮廓的矩形的竖边之间的最大水平距离来计算胸廓横径，通常会比右膈顶水平胸廓内径要大，从而使计算出来的心胸比偏小。为了解决这个问题，本实施例通过找到右膈顶，过右膈顶作水平线计算右膈顶水平胸廓的内径。由于右膈顶是右肺野下方穹窿形边界的顶部，也就是轮廓线的极大值点，此点的导数为0，过该点的切线的斜率为0或斜率绝对值最小。因此，可以在右肺野下半部分的边界上搜索切线斜率绝对值最小的点，此点即为右膈顶点(图4中的p点)；过该点的水平线与左、右肺野轮廓的外侧交点之间的线段长度(图4中线段p1p2的长度)就是右膈顶水平胸廓的内径，即胸廓横径。值得注意的是，右肺野上方边界也有一个极大值点，因此，但该点不是右膈顶，通常称为肺尖，计算时需要避开肺尖部分，否则会出现错误。
本实施例中，步骤104主要用于计算心胸比。有了心脏横径和胸廓横径，计算二者的比值就得到了心胸比。
作为一可选实施例，所述方法还包括将输入胸片图像缩放到800*1333像素*像素，灰度值范围为[0,255]。
本实施例给出了对输入胸片图像进行预处理的一种技术方案。由于市面上的DR(直接数字化X射线摄影装置)机型多样，拍摄后处理参数各不相同，成像质量差异较大，造成不同DR影像具有不同的灰度值范围。如果直接使用DR影像原始的灰度值作为输入，将不利于用于轮廓分割的深度学习网络的训练和预测。因此，需要对DR影像的灰度值进行统一化处理，使不同的DR影像具有基本一致的灰度值范围。具体的灰度值统一化方法为：将DR影像的灰度值截断在以窗位为中心，以窗宽为长度的区间(称之为窗口区间)里，然后线性地将窗口区间映射到目标区间，如[0,255]。对于有些缺失窗宽或窗位信息的图像，直接对其图像的灰度值进行线性映射，将灰度值的最小和最大值分别映射到目标区间的最小和最大值。这样，DR影像的灰度值就全部被映射到了一个统一的范围里。同时，为了应对不同DR影像尺寸差异较大的问题，还对DR影像进行了尺寸标准化处理，将胸片图像缩放到800*1333像素*像素。具体做法是：将图像保持宽高比按较小边缩放到800，此时如果其较大边超过了1333，就保持宽高比按较大边缩放到1333。
作为一可选实施例，所述方法还包括输入胸片图像倾斜校正步骤：将输入的胸片图像输入一个空间变换网络STN(Spatial Transformer Network)，基于STN的输出与模板图像的差异计算损失函数，利用所述损失函数对STN进行训练，通过优化网络参数(α,sx,sy,tx,ty)对输入胸片图像进行旋转、伸缩和平移变换，使其与模板图像对正；其中，α为旋转角度，sx、sy分别为x、y方向的缩放倍数，tx、ty分别x、y方向的平移量；所述模板图像通过对多个体中线为竖直线的标准正位胸片图像求均值获得。
本实施例给出了对输入胸片图像进行倾斜校正的一种技术方案。标准的正位胸片的体中线应为竖直线，即与纵轴平行。但实际上很多胸片图像都有一定程度的倾斜，使其体中线偏离竖直方向，如图5所示。由前面实施例给出的心脏横径和胸廓横径的求解方法可知，如果胸片图像的体中线偏离竖直方向，得到的心脏横径和胸廓横径都将有一定的偏差，从而影响心胸比的计算精度。因此需要对倾斜的胸片图像进行校正。
本实施例利用一个STN网络实现胸片图像的校正。STN网络的结构如图6所示，待校正的胸片图像作为STN网络的输入，STN网络的Localization net部分用来回归出空间变换的参数。在训练STN网络时，将STN网络的输出(即图中的V)与模板图像进行对比，计算它们之间的感知损失(perceptual loss)，并以此感知损失为损失函数优化网络参数使得感知损失足够小。当感知损失足够小时，输出V和模板图像在整体上足够接近，从而使待校正的胸片图像的体中线也接近竖直。网络参数(α,sx,sy,tx,ty)分别是旋转角度、缩放倍数和平移量，因此，网络参数的优化过程是实际上是对输入胸片图像逐步进行旋转、伸缩和平移的变换过程；也是输入胸片图像与模板图像的对正过程，因此对胸片图像的倾斜校正过程通常也称为对正或对齐。STN网络中的Localization net可以是深度学习中常用的分类/回归网络，如VGG，ResNet以及DenseNet等。
本实施例通过利用STN网络对胸片图像进行倾斜校正，使胸片图像的体中线成为竖直线，可提高心脏横径和胸廓横径的计算精度，从而提高心胸比的计算精度。
作为一可选实施例，从输入胸片图像到模板图像的坐标变换方程为：
Figure BDA0003088999650000081
即：
Figure BDA0003088999650000091
式中，T为坐标变换矩阵，(x,y)为输入胸片图像的像素点坐标，(x′,y′)为模板图像的像素点坐标。
本实施例给出了倾斜校正的等效变换方程。倾斜校正实际上是对一个坐标系的旋转、伸缩和平移的变换过程，即从输入胸片图像坐标系到模板图像坐标系的变换过程。变换矩阵为T，矩阵形式和方程组形式的变换方程分别如上式所示。
作为一可选实施例，所述斜率绝对值的计算公式为：
Figure BDA0003088999650000092
式中，ki为右肺野轮廓上第i个像素点(xi,yi)处的斜率的绝对值，(xi-δ,yi-δ)、(xi+δ,yi+δ)分别为(xi,yi)前面和后面第δ个像素点的坐标，yi＞(ymax+ymin)/2，ymax、ymin分别为右肺野轮廓上纵坐标的最大值和最小值，坐标原点在图像的左上角，纵坐标y向下为正方向。
本实施例给出了斜率绝对值的计算公式。本实施例采用相邻点的差商作为切线斜率的近似。同时，考虑到可能出现“拐点”，如图4下侧最左边和最右边的点，采用双边差商的绝对值的平均值来近似代替切线斜率的绝对值，如上式所示。如果采用通常的单边差商计算斜率，那么在图4下侧“拐点”会出现斜率等于0或斜率绝对值最小的情况，从而被误判为是右膈顶。而采用双边差商后，只有一边可能为0，即yi+δ-yi＝0或yi-δ-yi＝0，却不可能出现双边同时为0即yi-δ＝yi＝yi+δ的情况，因为在“拐点”上方作水平线与轮廓只能有两个交点而不可能有三个交点。δ为一较小的整数，一般取1～3，取1时为最邻近点，取3时最稳定。为了避开右肺野上方边界那个极大值点即“肺尖”的干扰，这里限定了yi＞(ymax+ymin)/2。由于坐标原点在图像的左上角，纵坐标y向下为正，如此限定后使求解范围限制在右肺野的下方边界。
图7为本发明实施例一种正位胸片心胸比测量装置的组成示意图，所述装置包括：
轮廓分割模块11，用于从输入胸片图像中分割出心脏轮廓、左肺野轮廓和右肺野轮廓；
心脏横径计算模块12，用于分别求出心脏轮廓横坐标的最大值和最小值，计算所述最大值和最小值的差，得到心脏横径；
胸廓横径计算模块13，用于在右肺野轮廓的下部找到斜率绝对值最小的点，过此点作一水平线，计算所述水平线与左肺野轮廓、右肺野轮廓的外侧交点之间的线段长度，得到胸廓横径；
心胸比计算模块14，用于用所述心脏横径除以所述胸廓横径得到心胸比。
本实施例的装置，可以用于执行图1所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。后面的实施例也是如此，均不再展开说明。
作为一可选实施例，所述装置还包括将图像标准化模块，用于将输入胸片图像缩放到800*1333像素*像素，灰度值范围为[0,255]。
作为一可选实施例，所述装置还包括倾斜校正模块，用于将输入的胸片图像输入一个空间变换网络STN，基于STN的输出与模板图像的差异计算损失函数，利用所述损失函数对STN进行训练，通过优化网络参数(α,sx,sy,tx,ty)对输入胸片图像进行旋转、伸缩和平移变换，使其与模板图像对正；其中，α为旋转角度，sx、sy分别为x、y方向的缩放倍数，tx、ty分别x、y方向的平移量；所述模板图像通过对多个体中线为竖直线的标准正位胸片图像求均值获得。
作为一可选实施例，从输入胸片图像到模板图像的变换方程为：
Figure BDA0003088999650000101
即：
Figure BDA0003088999650000102
式中，T为坐标变换矩阵，(x,y)为输入胸片图像的像素点坐标，(x′,y′)为模板图像的像素点坐标。
作为一可选实施例，所述斜率绝对值的计算公式为：
Figure BDA0003088999650000111
式中，ki为右肺野轮廓上第i个像素点(xi,yi)处的斜率的绝对值，(xi-δ,yi-δ)、(xi+δ,yi+δ)分别为(xi,yi)前面和后面第δ个像素点的坐标，yi＞(ymax+ymin)/2，ymax、ymin分别为右肺野轮廓上纵坐标的最大值和最小值，坐标原点在图像的左上角，纵坐标y向下为正方向。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

Claims (10)

1.一种正位胸片心胸比测量方法，其特征在于，包括以下步骤：
从输入胸片图像中分割出心脏轮廓、左肺野轮廓和右肺野轮廓；
分别求出心脏轮廓横坐标的最大值和最小值，计算所述最大值和最小值的差，得到心脏横径；
在右肺野轮廓的下部找到斜率绝对值最小的点，过此点作一水平线，计算所述水平线与左肺野轮廓、右肺野轮廓的外侧交点之间的线段长度，得到胸廓横径；
用所述心脏横径除以所述胸廓横径得到心胸比。
2.根据权利要求1所述的正位胸片心胸比测量方法，其特征在于，所述方法还包括将输入胸片图像缩放到800*1333像素*像素，灰度值范围为[0,255]。
3.根据权利要求1所述的正位胸片心胸比测量方法，其特征在于，所述方法还包括输入胸片图像倾斜校正步骤：将输入的胸片图像输入一个空间变换网络STN，基于STN的输出与模板图像的差异计算损失函数，利用所述损失函数对STN进行训练，通过优化网络参数(α,sx,sy,tx,ty)对输入胸片图像进行旋转、伸缩和平移变换，使其与模板图像对正；其中，α为旋转角度，sx、sy分别为x、y方向的缩放倍数，tx、ty分别x、y方向的平移量；所述模板图像通过对多个体中线为竖直线的标准正位胸片图像求均值获得。
4.根据权利要求3所述的正位胸片心胸比测量方法，其特征在于，从输入胸片图像到模板图像的坐标变换方程为：
Figure FDA0003088999640000011
即：
Figure FDA0003088999640000012
式中，T为坐标变换矩阵，(x,y)为输入胸片图像的像素点坐标，(x′,y′)为模板图像的像素点坐标。
5.根据权利要求1所述的正位胸片心胸比测量方法，其特征在于，所述斜率绝对值的计算公式为：
Figure FDA0003088999640000021
式中，ki为右肺野轮廓上第i个像素点(xi,yi)处的斜率的绝对值，(xi-δ,yi-δ)、(xi+δ,yi+δ)分别为(xi,yi)前面和后面第δ个像素点的坐标，yi＞(ymax+ymin)/2，ymax、ymin分别为右肺野轮廓上纵坐标的最大值和最小值，坐标原点在图像的左上角，纵坐标y向下为正方向。
6.一种正位胸片心胸比测量装置，其特征在于，包括：
轮廓分割模块，用于从输入胸片图像中分割出心脏轮廓、左肺野轮廓和右肺野轮廓；
心脏横径计算模块，用于分别求出心脏轮廓横坐标的最大值和最小值，计算所述最大值和最小值的差，得到心脏横径；
胸廓横径计算模块，用于在右肺野轮廓的下部找到斜率绝对值最小的点，过此点作一水平线，计算所述水平线与左肺野轮廓、右肺野轮廓的外侧交点之间的线段长度，得到胸廓横径；
心胸比计算模块，用于用所述心脏横径除以所述胸廓横径得到心胸比。
7.根据权利要求6所述的正位胸片心胸比测量装置，其特征在于，所述装置还包括将图像标准化模块，用于将输入胸片图像缩放到800*1333像素*像素，灰度值范围为[0,255]。
8.根据权利要求6所述的正位胸片心胸比测量装置，其特征在于，所述装置还包括倾斜校正模块，用于将输入的胸片图像输入一个空间变换网络STN，基于STN的输出与模板图像的差异计算损失函数，利用所述损失函数对STN进行训练，通过优化网络参数(α,sx,sy,tx,ty)对输入胸片图像进行旋转、伸缩和平移变换，使其与模板图像对正；其中，α为旋转角度，sx、sy分别为x、y方向的缩放倍数，tx、ty分别x、y方向的平移量；所述模板图像通过对多个体中线为竖直线的标准正位胸片图像求均值获得。
9.根据权利要求8所述的正位胸片心胸比测量装置，其特征在于，从输入胸片图像到模板图像的坐标变换方程为：
Figure FDA0003088999640000031
即：
Figure FDA0003088999640000032
式中，T为坐标变换矩阵，(x,y)为输入胸片图像的像素点坐标，(x′,y′)为模板图像的像素点坐标。
10.根据权利要求6所述的正位胸片心胸比测量装置，其特征在于，所述斜率绝对值的计算公式为：
Figure FDA0003088999640000033
式中，ki为右肺野轮廓上第i个像素点(xi,yi)处的斜率的绝对值，(xi-δ,yi-δ)、(xi+δ,yi+δ)分别为(xi,yi)前面和后面第δ个像素点的坐标，yi＞(ymax+ymin)/2，ymax、ymin分别为右肺野轮廓上纵坐标的最大值和最小值，坐标原点在图像的左上角，纵坐标y向下为正方向。
CN202110590001.6A 2021-05-28 2021-05-28 一种正位胸片心胸比测量方法及装置 Active CN113450399B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110590001.6A CN113450399B (zh) 2021-05-28 2021-05-28 一种正位胸片心胸比测量方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110590001.6A CN113450399B (zh) 2021-05-28 2021-05-28 一种正位胸片心胸比测量方法及装置

Publications (2)

Publication Number Publication Date
CN113450399A true CN113450399A (zh) 2021-09-28
CN113450399B CN113450399B (zh) 2022-02-25

Family

ID=77810378

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110590001.6A Active CN113450399B (zh) 2021-05-28 2021-05-28 一种正位胸片心胸比测量方法及装置

Country Status (1)

Country Link
CN (1) CN113450399B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114445377A (zh) * 2022-01-28 2022-05-06 北京大学第一医院 一种心胸比计算方法及心脏肥大辨别方法、装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521476A (zh) * 2011-11-03 2012-06-27 丁国成 基于成人x线正位胸片的异常心脏大血管计算机筛选系统及筛选方法
CN106991697A (zh) * 2017-04-17 2017-07-28 心医国际数字医疗系统(大连)有限公司 胸部数字影像的心胸比测量方法
CN107665497A (zh) * 2016-07-29 2018-02-06 上海联影医疗科技有限公司 一种医学图像中计算心胸比的方法
CN108918526A (zh) * 2018-04-10 2018-11-30 华南理工大学 一种柔性ic封装基板线路的缺口缺陷检测方法
CN111462020A (zh) * 2020-04-24 2020-07-28 上海联影医疗科技有限公司 心脏图像的运动伪影校正方法、系统、存储介质和设备
CN111784616A (zh) * 2020-07-29 2020-10-16 中科汇金数字科技(北京)有限公司 一种基于图像处理的老唱片数字化音频处理方法
CN112508918A (zh) * 2020-12-11 2021-03-16 上海商汤智能科技有限公司 图像处理方法及装置、电子设备和存储介质

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521476A (zh) * 2011-11-03 2012-06-27 丁国成 基于成人x线正位胸片的异常心脏大血管计算机筛选系统及筛选方法
CN107665497A (zh) * 2016-07-29 2018-02-06 上海联影医疗科技有限公司 一种医学图像中计算心胸比的方法
CN106991697A (zh) * 2017-04-17 2017-07-28 心医国际数字医疗系统(大连)有限公司 胸部数字影像的心胸比测量方法
CN108918526A (zh) * 2018-04-10 2018-11-30 华南理工大学 一种柔性ic封装基板线路的缺口缺陷检测方法
CN111462020A (zh) * 2020-04-24 2020-07-28 上海联影医疗科技有限公司 心脏图像的运动伪影校正方法、系统、存储介质和设备
CN111784616A (zh) * 2020-07-29 2020-10-16 中科汇金数字科技(北京)有限公司 一种基于图像处理的老唱片数字化音频处理方法
CN112508918A (zh) * 2020-12-11 2021-03-16 上海商汤智能科技有限公司 图像处理方法及装置、电子设备和存储介质

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
武杰等: "心胸比率的计算机辅助测量和判断", 《中国医学物理学杂志》 *
王亚丽等: "探讨研究CT判断心脏增大的测量方法", 《中国临床医学影像杂志》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114445377A (zh) * 2022-01-28 2022-05-06 北京大学第一医院 一种心胸比计算方法及心脏肥大辨别方法、装置

Also Published As

Publication number Publication date
CN113450399B (zh) 2022-02-25

Similar Documents

Publication Publication Date Title
CN106485695B (zh) 基于统计形状模型的医学图像Graph Cut分割方法
CN107665497B (zh) 一种医学图像中计算心胸比的方法
US20180005372A1 (en) Methods and systems for detecting a centerline of a vessel
US20180114313A1 (en) Medical Image Segmentation Method and Apparatus
CN111105424A (zh) 淋巴结自动勾画方法及装置
WO2022088581A1 (zh) 图像检测模型的训练方法及相关装置、设备、存储介质
CN106611416B (zh) 一种医学图像中肺分割的方法及装置
CN106056596B (zh) 基于局部先验信息和凸优化的全自动三维肝脏分割方法
CN106296613B (zh) 一种基于dr机器的双能量减影方法
CN112785632B (zh) 基于epid的图像引导放疗中dr和drr影像跨模态自动配准方法
Kurazume et al. 3D reconstruction of a femoral shape using a parametric model and two 2D fluoroscopic images
Yim et al. Correction of segmented lung boundary for inclusion of pleural nodules and pulmonary vessels in chest CT images
CN111951215A (zh) 一种图像检测方法、装置及计算机可读存储介质
CN113450399B (zh) 一种正位胸片心胸比测量方法及装置
Liu et al. An enhanced neural network based on deep metric learning for skin lesion segmentation
CN117274216B (zh) 基于水平集分割的超声下颈动脉斑块检测方法及系统
CN113012127A (zh) 基于胸部医学影像的心胸比测量方法
CN110533667B (zh) 基于图像金字塔融合的肺部肿瘤ct影像3d分割方法
CN112529900A (zh) 匹配乳腺图像中roi的方法、装置、终端与存储介质
CN111513743A (zh) 一种骨折检测方法及装置
CN115294401B (zh) 一种基于ct图像的分类方法、装置、设备及存储介质
CN107451610B (zh) 一种提高特征匹配精度的图像检测方法
Zhuo et al. Local adaptive segmentation algorithm for 3-D medical image based on robust feature statistics
CN109300138B (zh) 一种基于两步纳什均衡改进c-v模型的医学图像分割方法
CN112258536A (zh) 一种胼胝体及小脑蚓部的一体化定位及分割方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A method and device for measuring cardiothoracic ratio of chest radiograph

Effective date of registration: 20230111

Granted publication date: 20220225

Pledgee: Haidian Beijing science and technology enterprise financing Company limited by guarantee

Pledgor: Beijing Yizhun Intelligent Technology Co.,Ltd.

Registration number: Y2023110000013

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: Room 3011, 2nd Floor, Building A, No. 1092 Jiangnan Road, Nanmingshan Street, Liandu District, Lishui City, Zhejiang Province, 323000

Patentee after: Zhejiang Yizhun Intelligent Technology Co.,Ltd.

Address before: Room 1106, Weishi building, Haidian District, Beijing 100083

Patentee before: Beijing Yizhun Intelligent Technology Co.,Ltd.

PC01 Cancellation of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Granted publication date: 20220225

Pledgee: Haidian Beijing science and technology enterprise financing Company limited by guarantee

Pledgor: Beijing Yizhun Intelligent Technology Co.,Ltd.

Registration number: Y2023110000013