[发明专利]一种基于形状信息的腕骨区域分割方法

权利要求书

1.一种基于形状信息的腕骨区域分割方法，其特征在于，包括以下步骤：

步骤1：采用星型中值滤波，灰度图像二值化方法，依次对X光图像进行预处理；对二值化图像进行8连通域的标记，使得同一连通区域内的所有像素相同，不同连通区域像素不同，并统计所有连通区域内白色像素点S_n，n代表标记为n的连通域；寻找白色像素点最多的一个连通区域m，即为手掌区域，满足S_m≥S_n，N为连通区域的数量；其余的连通区域当作背景区域，使用0像素值填充：

其中，X表示经过预处理后的图像；

步骤2：根据腕骨区域分布特点，使用双向迭代扫描算法，确定腕骨区域上边界的纵坐标位置；具体包括：

21)求出连通区域S_m的最小外接矩阵，高和宽分别为BH，BW；预处理后的图像高和宽分别为XH，XW；

22)定义一个函数F(R_t)，接收图像中某一行t的所有像素值集合R_t作为参数，R_t＝{X(t,j)|j＝1,2,…,XW}，t∈{1,2,…,XH}，用于判断第t行是否满足以下条件：

D＝5

p_i＝len(G(R_i)) (2-2)

X(i,a-1)＝0,X(i,b+1)＝0}

其中，D值表示判定取样间隙，初始值为5，函数G(R_i)表示找出R_i集合中所有像素大于0的连续子序列段，p_i表示G(R_i)中序列段的数量；

23)从下往上，每隔50像素点对图像X进行扫描，起始位置为XH-BH/3，结束位置为XH-BH，当某行t，满足F(R_t)为TRUE时，XH-BH/3≤t≤XH-BH，跳转到24)，否则继续扫描；

24)从上往下，每隔5像素点对图像X进行扫描，起始位置为t，结束位置为t+50，当某行u，满足F(R_u)为TRUE时，t≤u≤t+50，结束扫描，否则继续扫描；

25)确定腕骨区域上边界的纵坐标位置top，公式为：

步骤3：根据手腕关节形状信息，计算手腕两侧边缘梯度变化，确定下边界左右候选点纵坐标位置；

步骤4：综合上边界、左右候选点和腕骨区域信息，对坐标位置进行优化矫正，并根据得到的左上角、右下角坐标进行腕骨区域裁剪。