一种超声图像中小脑蚓部分割方法技术

本发明专利技术公开了一种超声图像中小脑蚓部分割方法，包括以下步骤：S1:获取小脑蚓部CT图像并将其转换为灰度图；S2:读取所述灰度图的边缘坐标，以灰度图的中心点为圆心，求出以该圆心0‑360°内每个角度上到图像边缘点形成的线段上的所有像素点坐标及其在图像上对应的像素值；S3:以圆心向外每个角度形成的射线上寻找此射线上的初始分割点，S4:采用误差点判断的方式对粗调整分割点进行误差判断：查找单峰误差点，对单峰误差点周围的分割点进行误差判断，S5:对待精确调整数组中的所有点进行调整得到精准分割点；S6:将调整后的精准分割点顺次连接构成小脑蚓部边缘图像。

【技术实现步骤摘要】
一种超声图像中小脑蚓部分割方法
本专利技术涉及医学图像处理领域，尤其涉及超声图像中小脑蚓部分割方法。
技术介绍
随着计算机技术在医学上的广泛应用，使得早先的人工分割图像逐渐被计算机自动分割所替代。目前对于如何检测胎儿的小脑蚓部是否发育正常，通常采用人工分割CT图像中小脑蚓部区域，并计算其面积来检测出是否发育正常。如何简化图像分割，并使分割效果更加精细成为关键的问题。目前对于CT图像分割小脑蚓部问题，没有已知的有效的计算机自动处理办法，都是由人工来人为分割。这需要专业的知识储备，也花费了不少时间，并且效果不一定精细，所以准确性不能保证。
技术实现思路
根据现有技术中没有存在针对小脑蚓部分割的现状，本专利技术公开了一种超声图像中小脑蚓部分割方法，具体包括以下步骤：获取小脑蚓部CT图像并将其转换为灰度图；读取所述灰度图的边缘坐标，以灰度图的中心点为圆心，计算以该圆心0-360°内每个角度上到图像边缘点形成的线段上的所有像素点坐标及其在图像上对应的像素值；以圆心向外每个角度形成的射线上寻找此射线上的初始分割点，对于获取的初始分割点、以每隔区间角度值形成的多圈分割点采用误差点判断的方式进行误差粗判断并修改误差点，剔除杂点得到粗调整分割点；采用误差点判断的方式对粗调整分割点进行误差判断，查找单峰误差点，对单峰误差点周围的分割点进行误差判断从而重新确定应该删除的分割点，对所有要删除的分割点保存到待精确调整数组内；对待精确调整数组中的所有点进行调整得到精准分割点；将调整后的精准分割点顺次连接构成小脑蚓部边缘图像。进一步的，所述寻找此射线上的分割点时：设置区间角度值，分别以每个整数角度对应的线段为起始、每隔所述区间角度值寻找该角度线段上的分割点；其中在每个角度线段上选择分割点时采用如下方式：首先获取该角度射线上所有点的个数length，累加该射线[length/3,length/2]区域上背景点的个数，如果个数小于length/6并且length/2处像素值大于设定阈值A时，则采用判断区间起始值length/2进行分割点寻找，否则采用length/3进行分割点寻找；在线段上从判断区间起始位置到边缘点中，若按顺序找到一个点的像素值小于设定阈值A则认为其之后的像素为背景像素,将该点设为初始分割点，否则选择像素值最小的点作为初始分割点；对得到的初始分割点采用误差点判断的方式进行调整：以区间角度值为区域，分别以每个整数角度对应的分割点为起始、将每隔区间角度值形成的一圈分割点定义为分割点圆周，并对该分割点圆周进行误差判断：如果当前点(xj,yj)到中心点O(x0,y0)的距离同时小于分割点圆周上前一分割点和后一分割点到中心点的距离、并且距离差大于设定阈值B则认为当前点(xj,yj)是误差点，则改变当前点(xj,yj)的位置，其中改变方式如下：将分割点圆周上前一分割点到中心点距离的一半和后一分割点到中心点距离的一半的加和作为当前点(xj,yj)到中心点的距离M；如果当前点(xj,yj)到中心点O(x0,y0)的距离都大于分割点圆周上前一分割点和后一分割点到中心点的距离，并且距离差大于设定阈值B，则认为当前点(xj,yj)是误差点，则改变当前点(xj,yj)的位置，其中改变方式如下：将分割点圆周上前一分割点到中心点距离的一半和后一分割点到中心点距离的一半的加和作为当前点(xj,yj)到中心点的距离M，根据获得的距离M调整(xj,yj)的新点坐标作为粗调整分割点。进一步的，对粗调整分割点再次进行误差调整:采用误差点判断的方式对分割点圆周上360个分割点依次进行误差判断并修改误差点；对修改后的得到的分割点进行单峰误差点剔除：将所有前后相邻三个粗调整分割点组成的矢量夹角的补角的平均值定义为角度阈值Tθ，将所有粗调整分割点到中心点的距离的平均值的Q倍定义为弦高阈值Ts；对于分割点圆周上360个粗调整分割点中的每个点计算其到中心点的距离d，如果每点到中心点的距离大于弦高阈值Ts并且角度大于角度阈值Tθ，则将该点存放在待删除列表中；在待删除列表中进行判断，并在待删除列表的基础上重新确定需要删除的分割点：如果某分割点和它的下一个需要删除分割点到中心点的距离差大于设定阈值B，则将某分割点和它的下一个需要删除分割点以及二者的中间分割点都删除；如果某分割点和前一个需要删除的分割点到中心点的距离差大于设定阈值B，并且某分割点之后的分割点连续都为需要删除的分割点，此时这些分割点都要删掉，将所有需要删除的分割点保存至待精确调整数组中。进一步的，将待精确调整数组中重复的分割点删除，对分割点按照角度的先后进行排序，如果当前分割点(xr,yr)在原始数据中的前一个相邻分割点和待精确调整数组中的前一个分割点不一致，则将当前分割点(xr,yr)到中心点的距离tj改为前一分割点到中心点的距离tj-1，根据新获取的距离tj-1计算当前分割点(xr,yr)的新坐标，将当前分割点(xr,yr)从待精确调整数组中移除；采用上述方式对待精确调整数组中的分割点循环处理直到数组中没有成员为止。由于采用了上述技术方案，本专利技术提供的超声图像中小脑蚓部分割方法，本方法只需要用户给出期望分割的小脑蚓部的CT图，即可得到所需的分割出的小脑蚓部边缘，减少了日常生活中需要医生手工来标记边缘的操作，从而减少一些由于人工失误导致的误差，使得准确性和效率提高，并且不需要有相应的医学技术就可以操作，并且占用内存较小，运行时间快。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术方法实现的流程图；图2为本专利技术中将图像转换为灰度图的效果图；图3为本专利技术中求图像边缘点的示意图；图4为本专利技术中单峰误差点剔除的示意图；图5为本专利技术中删除列表中添加删除点的示意图；图6为本专利技术中小脑蚓部CT图像分割后的效果图。具体实施方式为使本专利技术的技术方案和优点更加清楚，下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚完整的描述：如图1所示的一种超声图像中小脑蚓部分割方法，在实施过程中，将图像转换为灰度图的效果图如图2所示，然后计算出图像一圈360度的边缘点，在图像中心和各个边缘点之间找到该角度对应的分割点，然后对误差点进行不断调整，使之更加贴近真实情况，最后将分割点显示出来。本专利技术公开的方法具体步骤如下：S1：输入小脑蚓部CT图像，将图像转换为灰度图，具体采用如下方式：遍历图像的每个像素点的像素值(rk,gk,bk)(k＝0...num，num为像素点的个数)，计算出的灰度值为：grayk＝0.299*rk+0.587*gk+0.114*bkrk＝gray本文档来自技高网...

【技术保护点】
1.一种超声图像中小脑蚓部分割方法，其特征在于包括：/n获取小脑蚓部CT图像并将其转换为灰度图；/n读取所述灰度图的边缘坐标，以灰度图的中心点为圆心，计算以该圆心0-360°内每个角度上到图像边缘点形成的线段上的所有像素点坐标及其在图像上对应的像素值；/n以圆心向外每个角度形成的射线上寻找此射线上的初始分割点，对于获取的初始分割点、以每隔区间角度值形成的多圈分割点采用误差点判断的方式进行误差粗判断并修改误差点，剔除杂点得到粗调整分割点；/n采用误差点判断的方式对粗调整分割点进行误差判断，查找单峰误差点，对单峰误差点周围的分割点进行误差判断从而重新确定应该删除的分割点，对所有要删除的分割点保存到待精确调整数组内；/n对待精确调整数组中的所有点进行调整得到精准分割点；/n将调整后的精准分割点顺次连接构成小脑蚓部边缘图像。/n

【技术特征摘要】
1.一种超声图像中小脑蚓部分割方法，其特征在于包括：
获取小脑蚓部CT图像并将其转换为灰度图；
读取所述灰度图的边缘坐标，以灰度图的中心点为圆心，计算以该圆心0-360°内每个角度上到图像边缘点形成的线段上的所有像素点坐标及其在图像上对应的像素值；
以圆心向外每个角度形成的射线上寻找此射线上的初始分割点，对于获取的初始分割点、以每隔区间角度值形成的多圈分割点采用误差点判断的方式进行误差粗判断并修改误差点，剔除杂点得到粗调整分割点；
采用误差点判断的方式对粗调整分割点进行误差判断，查找单峰误差点，对单峰误差点周围的分割点进行误差判断从而重新确定应该删除的分割点，对所有要删除的分割点保存到待精确调整数组内；
对待精确调整数组中的所有点进行调整得到精准分割点；
将调整后的精准分割点顺次连接构成小脑蚓部边缘图像。


2.根据权利要求1所述的一种超声图像中小脑蚓部分割方法，其特征还在于：所述寻找此射线上的分割点时：设置区间角度值，分别以每个整数角度对应的线段为起始、每隔所述区间角度值寻找该角度线段上的分割点；
其中在每个角度线段上选择分割点时采用如下方式：首先获取该角度射线上所有点的个数length，累加该射线[length/3,length/2]区域上背景点的个数，如果个数小于length/6并且length/2处像素值大于设定阈值A时，则采用判断区间起始值length/2进行分割点寻找，否则采用length/3进行分割点寻找；在线段上从判断区间起始位置到边缘点中，若按顺序找到一个点的像素值小于设定阈值A则认为其之后的像素为背景像素,将该点设为初始分割点，否则选择像素值最小的点作为初始分割点；
对得到的初始分割点采用误差点判断的方式进行调整：以区间角度值为区域，分别以每个整数角度对应的分割点为起始、将每隔区间角度值形成的一圈分割点定义为分割点圆周，并对该分割点圆周进行误差判断：
如果当前点(xj,yj)到中心点O(x0,y0)的距离同时小于分割点圆周上前一分割点和后一分割点到中心点的距离、并且距离差大于设定阈值B则认为当前点(xj,yj)是误差点，则改变当前点(xj,yj)的位置，其中改变方式如下：将分割点圆周上前一分割点到中心点距离的一半和后一分割点到中心点距离的...

【专利技术属性】
技术研发人员：刘斌，陈彦洁，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21