[发明专利]X射线图像分割方法及终端

权利要求书

1.一种X射线图像分割方法，其特征在于，包括：

获取X射线空气图；

计算所述X射线空气图的噪声层级比例系数；

获取X射线临床图像；

根据所述X射线临床图像和噪声层级比例系数计算得到CNR图像；

对所述CNR图像进行二值化处理，得到CNR图像掩模；

根据所述CNR图像计算得到边界抑制掩模；

根据所述X射线临床图像分别计算得到空气区域掩模和束光器边缘抑制掩模；

对所述CNR图像掩模、边界抑制掩模、空气区域掩模和束光器边缘抑制掩模作与计算，得到组织区域掩模；

所述计算所述X射线空气图的噪声层级比例系数具体包括：

对所述X射线空气图进行AT变换；

采用图像金字塔分解AT变换后的X射线空气图，得到低频信号层G₀、G₃，以及细节信号层D₀、D₃；

在第一预设阈值范围内获取X射线空气图的G₃层中的最大灰度值maxV和最小灰度值minV；

根据公式计算得到G₀层的掩模，其中，T表示第一预设比例；

在G₀层的掩模范围内，根据公式计算得到标准差图像S0，其中m表示窗口大小；

计算标准差图像S₀直方图的最大值，得到X射线空气图的噪声水平nl₀；

根据公式计算得到标准差图像S₃；

计算标准差图像S₃直方图的最大值，得到X射线空气图的噪声水平nl₃；

根据所述X射线空气图的噪声水平nl₀和噪声水平nl₃计算得到噪声层级比例系数；

所述根据所述X射线临床图像和噪声层级比例系数计算得到CNR图像具体包括：

对所述X射线临床图像进行AT变换；

采用图像金字塔分解AT变换后的X射线临床图像，得到低频信号层G'₀、G'₃以及细节信号层D'₀和D'₃；

在第二预设阈值范围内获取X射线临床图像的G₃'层中的最大灰度值maxV'和最小灰度值minV'；

根据公式计算得到G'₀层的掩模，其中，T'表示第二预设比例；

在G'₀层的掩模范围内，根据公式计算得到标准差图像S'₀，其中m'表示窗口大小；

计算标准差图像S'₀直方图的最大值，得到X射线临床图像的噪声水平nl'₀；

根据所述噪声层级比例系数和X射线临床图像的噪声水平nl'₀计算得到X射线临床图像第3层的噪声水平nl'₃；

根据公式计算得到图像标准差S'₃；

将图像标准差S'₃除以噪声水平nl'₃计算得到CNR图像；

所述对所述CNR图像进行二值化处理，得到CNR图像掩模具体为：

在[a,b]的范围内，根据公式计算得到CNR图像掩模mask₀，其中，a、b均为常数；

所述根据所述CNR图像计算得到边界抑制掩模具体为：

获取CNR图像的宽度w和高度h；

根据公式计算得到边界抑制掩模mask₁，其中，bor＝floor(w*c)，c为常数，bor为边界抑制掩膜mask₁的宽度，floor为向下取整函数；

根据所述X射线临床图像分别计算得到空气区域掩模和束光器边缘抑制掩模具体为：

根据公式airThr＝minV'+airRatio·(maxV'-minV')计算得到空气区域阈值airThr，其中，airRatio为常数；

根据公式计算得到空气区域掩模mask₂；

根据公式colThr＝colRatio·(maxV'-minV')计算得到束光器区域阈值colThr，其中colRatio为常数；

计算X射线临床图像G'₃层中任一点周围像素的灰度差异；

若所述灰度差异大于colThr，则将所述任一点的灰度值设为0；若所述灰度差异小于colThr，则将所述任一点的灰度值设为1，得到束光器边缘抑制掩模mask₃。