[发明专利]双能减影参数的估算方法及计算机可读存储介质

说明书

本发明公开了一种双能减影参数的估算方法及计算机可读存储介质，方法包括：在边缘上的各边缘点的两侧分别获取预设个数的像素点，得到各边缘点对应的软组织点集和骨组织点集；根据像素点的坐标以及预设的不同的减影参数，得到像素点对应不同减影参数的第三灰度值；分别计算得到各边缘点对应各减影参数的软组织灰度和骨组织灰度；分别计算各边缘点对应同一减影参数的软组织灰度的方差；将最小的软组织灰度方差对应的减影参数作为纯骨组织像减影参数；分别计算各边缘点对应同一减影参数的软组织灰度与骨组织灰度的差值；计算所述灰度差的和；将最小的差值和对应的减影参数作为纯软组织像减影参数。本发明可快速准确地确定出合适的减影参数。

技术领域

本发明涉及医学影像处理技术领域，尤其涉及一种双能减影参数的估算方法及计算机可读存储介质。

背景技术

由于肺部有肋骨遮挡，容易造成肺部较小病灶(如结节等)的漏诊。为了降低漏诊率，一般使用虚拟双能减影和真实双能减影两种方法提取纯骨组织像和纯软组织像。虚拟双能减影只需要曝光一次，即可使用模式识别的方法去除肋骨遮挡，但是可能会把钙化灶当作肋骨一并去除。

真实双能减影利用不同组织对不同KV的X射线的吸收特性，通过两次不同KV(例如60KV和120KV)的曝光分别得到低能图像和高能图像。利用骨组织和软组织在不同KV的X射线下衰减率的不同，可以提取纯骨像和纯软组织像。较为简单的提取算法是使用低能图像和高能图像幂的商，公式为：

I_des(x，y)＝I_l(x，y)/I_h(x，y)^ω

其中，I_l(x,y)和I_h(x,y)分别是低能和高能图像中点的灰度，I_des(x,y)是减影图中点的灰度，ω是幂指数，也即减影参数。

其分辨率和图像稳定性都要优于虚拟双能减影。但是，由于不同病人的体厚不同，采用同一组参数可能会造成减影不完全，造成纯软组织像中有骨组织的伪像存在，如图1所示，影响医生诊断。并且不同部位都要设置不同的减影参数，操作较为繁琐。

发明内容

本发明所要解决的技术问题是：提供一种双能减影参数的估算方法及计算机可读存储介质，可快速准确地估算出合适的减影参数。

为了解决上述技术问题，本发明采用的技术方案为：一种双能减影参数的估算方法，包括：

对X光图像进行边缘检测，得到边缘；

在所述边缘上的各边缘点的两侧分别获取预设个数的像素点，得到各边缘点对应的软组织点集和骨组织点集；

根据所述像素点的坐标，获取所述像素点在低能图像中的第一灰度值以及在高能图像中的第二灰度值；

根据所述第一灰度值、第二灰度值以及预设的不同的减影参数，得到所述像素点对应不同减影参数的第三灰度值；

根据各边缘点对应的软组织点集和骨组织点集中各像素点对应不同减影参数的第三灰度值，分别计算得到各边缘点对应各减影参数的软组织灰度和骨组织灰度；

分别计算各边缘点对应同一减影参数的软组织灰度的方差，得到所述同一减影参数对应的软组织灰度方差；

将最小的软组织灰度方差对应的减影参数作为纯骨组织像减影参数；

分别计算各边缘点对应同一减影参数的软组织灰度与骨组织灰度的差值，得到各边缘点对应所述同一减影参数的灰度差；

计算所述灰度差的和，得到所述同一减影参数对应的差值和；

将最小的差值和对应的减影参数作为纯软组织像减影参数。