[发明专利]一种自动提取磁共振图像胼胝体的方法及系统

权利要求书

1.一种自动提取磁共振图像胼胝体的方法，其特征在于，包括：

A、根据正中矢状面T1加权核磁颅脑图像，确定颅脑轮廓阈值与胼胝体阈值；

B、根据所述颅脑轮廓阈值与所述胼胝体阈值，对正中矢状面T1加权核磁颅脑图像进行阈值分割，分别得到颅脑轮廓二值图像与胼胝体二值图像；

C、对所述颅脑轮廓二值图像与所述胼胝体二值图像分别进行边缘跟踪，获取颅脑轮廓位置与胼胝体边缘族；

D、根据胼胝体的解剖特性、几何特性，创建决策树，运用决策树对胼胝体边缘族进行筛选，确定胼胝体区域；

E、根据所述胼胝体区域确定颅脑扫描线组的倾斜角度及中心位置。

2.根据权利要求1所述自动提取磁共振图像胼胝体的方法，其特征在于，所述步骤A具体包括：

A1、统计正中矢状面T1加权核磁颅脑图像各灰度级的像素数目，并计算累积直方图；

A2、计算累计直方图导数值，根据导数值区间分类，获取累计直方图的转折点；

A3、根据图像特性，从转折点中分别获取颅脑轮廓阈值和胼胝体阈值。

3.根据权利要求1所述自动提取磁共振图像胼胝体的方法，其特征在于，所述步骤B具体包括：

B1、根据胼胝体阈值对T1加权核磁颅脑图像进行阈值分割，当第一灰度值大于或等于所述胼胝体阈值则为1，否则为0，根据所述第一灰度值获得胼胝体二值图像；

B2、根据颅脑轮廓阈值对T1加权核磁颅脑图像进行阈值分割，当第二灰度值大于或等于所述颅脑轮廓阈值则为1，否则为0，根据所述第二灰度值获得颅脑轮廓二值图像。

4.根据权利要求1所述自动提取磁共振图像胼胝体的方法，其特征在于，所述步骤D中胼胝体的解剖特性、几何特性包括胼胝体处于颅脑的中央位置、最小包裹矩形的长宽比、周长、面积信息。

5.根据权利要求1所述自动提取磁共振图像胼胝体的方法，其特征在于，所述步骤E具体包括：

E1、计算颅脑轮廓的区域中心与胼胝体的区域中心，并根据颅脑轮廓区域中心及胼胝体区域中心位置计算颅脑扫描线中心位置；

E2、从胼胝体水平中心向两边扫描计算胼胝体的左右底脚位置，并根据胼胝体区域的左右底脚位置，计算扫描线的倾斜角度，并确定颅脑扫描线组的角度。

6.一种自动提取磁共振图像胼胝体的系统，其特征在于，包括：

阈值确定模块，用于根据正中矢状面T1加权核磁颅脑图像，确定颅脑轮廓阈值与胼胝体阈值；

二值图像获取模块，用于根据所述颅脑轮廓阈值与所述胼胝体阈值，对正中矢状面T1加权核磁颅脑图像进行阈值分割，分别得到颅脑轮廓二值图像与胼胝体二值图像；

图像跟踪模块，用于对所述颅脑轮廓二值图像与所述胼胝体二值图像分别进行边缘跟踪，获取颅脑轮廓位置与胼胝体边缘族；

胼胝体区域获取模块，用于根据胼胝体的解剖特性、几何特性，创建决策树，运用决策树对胼胝体边缘族进行筛选，确定胼胝体区域；

扫描线组调节模块，用于根据所述胼胝体区域确定颅脑扫描线组的倾斜角度及中心位置。

7.根据权利要求6所述自动提取磁共振图像胼胝体的系统，其特征在于，所述阈值确定模块具体包括：

累积直方图获取单元，用于统计正中矢状面T1加权核磁颅脑图像各灰度级的像素数目，并计算累积直方图；

转折点获取单元，用于计算累计直方图导数值，根据导数值区间分类，获取累计直方图的转折点；

阈值获取单元，用于根据图像特性，从转折点中分别获取颅脑轮廓阈值和胼胝体阈值。

8.根据权利要求6所述自动提取磁共振图像胼胝体的系统，其特征在于，所述二值图像获取模块具体包括：

胼胝体二值图像获取单元，用于根据胼胝体阈值对T1加权核磁颅脑图像进行阈值分割，当第一灰度值大于或等于所述胼胝体阈值则为1，否则为0，根据所述第一灰度值获得胼胝体二值图像；

颅脑轮廓二值图像获取单元，用于根据颅脑轮廓阈值对T1加权核磁颅脑图像进行阈值分割，当第二灰度值大于或等于所述颅脑轮廓阈值则为1，否则为0，根据所述第二灰度值获得颅脑轮廓二值图像。

9.根据权利要求6所述自动提取磁共振图像胼胝体的系统，其特征在于，所述胼胝体区域获取模块中胼胝体的解剖特性、几何特性包括胼胝体处于颅脑的中央位置、最小包裹矩形的长宽比、周长、面积信息。

10.根据权利要求6所述自动提取磁共振图像胼胝体的系统，其特征在于，所述扫描线组调节模块具体包括：

扫描线中心位置获取单元，用于计算颅脑轮廓的区域中心与胼胝体的区域中心，并根据颅脑轮廓区域中心及胼胝体区域中心位置计算颅脑扫描线中心位置；

颅脑扫描线组的角度获取单元，用于从胼胝体水平中心向两边扫描计算胼胝体的左右底脚位置，并根据胼胝体区域的左右底脚位置，计算扫描线的倾斜角度，并确定颅脑扫描线组的角度。