[发明专利]一种X光图像处理方法、系统及X光图像处理设备

权利要求书

1.一种X光图像处理方法，其特征在于，所述方法包括：

获取待处理的X光图像；

预处理所述X光图像，以去除所述X光图像中的噪声；

计算所述X光图像的梯度图像；

根据所述梯度图像的梯度直方图计算梯度阈值，将梯度图像中梯度值小于所述梯度阈值的所有像素点标记为种子点；

将所述梯度图像和所述种子点作为分水岭算法的输入进行分水岭分割，得到X光图像的标记图像，所述标记图像包括多个标记区域，其中，每个标记代表一个区域；

对所述标记图像进行区域合并；

二值化处理所述区域合并后的X光图像；所述预处理X光图像，以去除所述X光图像中的噪声具体为：

计算所述X光图像的灰度直方图；

根据所述灰度直方图，由高向低依次将每一灰度值对应的像素点个数进行累加，并计算累加像素点个数与所述X光图像像素点个数的比值，当所述比值大于预设噪声像素点个数与所述X光图像像素点个数的比值时，将当前的灰度值作为灰度阈值；

遍历所述X光图像，将所述X光图像中所有大于所述灰度阈值的像素点的灰度值设置为所述灰度阈值。

2.如权利要求1所述的方法，其特征在于，所述计算所述X光图像的梯度图像具体为：计算所述X光图像所有像素点的梯度值，其中，计算每一个像素点的梯度的过程具体为：

在水平方向计算当前像素点右面两个像素和与左边两个像素和之差的绝对值Dx；

在垂直方向计算当前像素下面两个像素和与上面两个像素和之差的绝对值Dy；

取Dx和Dy中的较大值作为当前像素点的梯度值。

3.如权利要求1所述的方法，其特征在于，所述根据所述梯度图像的梯度直方图计算梯度阈值，将梯度图像中梯度值小于所述梯度阈值的所有像素点标记为种子点具体为：

计算所述梯度图像的梯度直方图；

根据所述梯度直方图，由低向高依次将每一灰度值对应的像素点个数进行累加，并计算累加像素点个数与所述X光图像像素点个数的比值，当所述比值大于预设值时，将当前的梯度值作为梯度阈值；

遍历所述梯度图像，将所述梯度图像中所有小于所述梯度阈值的像素点作为种子点。

4.如权利要求1所述的方法，其特征在于，所述对所述标记图像进行区域合并具体为：

判断当前区域面积是否大于预设面积阈值；

若大于，则当前区域不参与区域合并；

若小于，则判断当前区域平均灰度值与所有区域的最大平均灰度值的比值是否大于预设值；

若大于，则当前区域不参与区域合并；

若小于，则计算当前区域的所有相邻区域的平均灰度，并将所述相邻区域中平均灰度最小的区域与当前区域进行合并。

5.如权利要求1所述的方法，其特征在于，所述二值化处理所述区域合并后的X光图像包括：

对所述区域合并后的X光图像进行背景预分割，获取背景区域和前景区域的灰度比例s_p，其中，所述s_p满足开口方向向下的二次曲线y_p＝a*s_p²+b*s_p+c；

通过训练获取所述二次曲线的参数值a，b，c；

根据所述s_p，a，b，c，计算y_p；

根据所述y_p计算所述X光图像的二次分割阈值为：t＝X光图像背景区域的平均灰度值*y_p；

根据所述二次分割阈值t对所述X光图像进行二次分割，得到所述X光图像的二值图像，具体分割过程为：如果当区域的灰度平均值小于t，则当前区域为前景，如果当区域的灰度平均值大于t，则当前区域为背景。