[发明专利]一种全自动的检测腹盆腔CT体积内ASIS的方法

说明书

本发明公开了一种全自动的检测腹盆腔CT体积内ASIS的方法，包括以下步骤：S1、将全身CT图像中的中轴平面以上部分裁剪掉，得到下半身图像；S2、设定前景像素、背景像素；S3、建立坐标系并将下半身图像置于坐标系中；S4、从上到下对下半身图像进行水平切片，并从下半身图像中最顶端的切片开始对切片进行搜索；S5、根据前景像素和背景像素，在当前切片中找到L_ASIS候选位置和R_ASIS候选位置；S6、对L_ASIS候选位置、R_ASIS候选位置进行选择；S7、最终确定，切片中i坐标值最低的L_ASIS候选位置、R_ASIS候选位置即为最终检测位置。

技术领域

本发明涉及医学影像领域，尤其涉及一种全自动的检测腹盆腔CT体积内ASIS的方法。

背景技术

ASIS(髂前上棘)是人类前骨盆中易于检测(通过简单触诊)的骨结构，通常用作不同医学应用的参考标志。但是现有技术中通常采用人工分辨的方式在CT影像中寻找ASIS的位置，其操作方式较为复杂，严重影响了医学影像操作的工作效率，因此，有必要研究一种在CT影像中自动测量ASIS位置的方法来解决上述问题。

发明内容

本发明目的是针对上述问题，提供一种提高操作效率的全自动的检测腹盆腔CT体积内ASIS的方法。

为了实现上述目的，本发明的技术方案是：

一种全自动的检测腹盆腔CT体积内ASIS的方法，包括以下步骤：

S1、将全身CT图像中的中轴平面以上部分裁剪掉，得到下半身图像；

S2、设定下阈值LT和上阈值UT，将下半身图像中位于区间[LT,UT]内的体素值设置为前景像素，将下半身图像中位于区间[LT,UT]外的体素值设置为背景像素；

S3、建立坐标系并将下半身图像置于坐标系中；

S4、从上到下对下半身图像进行水平切片，并从下半身图像中最顶端的切片开始对切片进行搜索；

S5、根据前景像素和背景像素，在当前切片中找到L_ASIS候选位置和R_ASIS候选位置；

S6、如果当前切片中L_ASIS候选位置、R_ASIS候选位置的i坐标值低于上一级切片中L_ASIS候选位置、R_ASIS候选位置的i坐标值，则接受当前切片中的L_ASIS候选位置、R_ASIS候选位置并选择下一级切片继续重复步骤S6；如果当前切片中L_ASIS候选位置、R_ASIS候选位置的i坐标值不低于上一级切片中L_ASIS候选位置、R_ASIS候选位置的i坐标值，则选择结束；

S7、最终确定，切片中i坐标值最低的L_ASIS候选位置、R_ASIS候选位置即为最终检测位置。

进一步的，所述步骤S2中，下阈值LT＝WL-WW/2，上阈值UT＝WL+WW/2；其中，WL和WW是仅保留骨组织的统计估计常数。

进一步的，所述步骤S5中，在当前切片中找到L_ASIS候选位置和R_ASIS候选位置包括以下步骤：

S51、在切片中设定L_ROI区域和R_ROI区域，在L_ROI区域中搜索L_ASIS候选位置，在R_ROI区域中搜索R_ASIS候选位置；

S52、在切片中将坐标值jJcutoff的部分划分出来作为L_ROI区域；

S53、在切片中将坐标值jN-Jcutoff的部分划分出来作为R_ROI区域；

S54、在L_ROI区域中找出前景像素，则L_ASIS候选位置为(i0，j0)，i0≤L_ROI区域中所有的i坐标值；

S55、在L_ROI区域中找出找出前景像素，则R_ASIS候选位置为(i1，j1)，i1≤R_ROI区域中所有的i坐标值。