基于最大流策略的三维体数据分割方法技术

本发明专利技术公开了一种基于最大流策略的三维体数据分割方法，包括以下步骤：S1：读取虚拟人序列化切片图像中的虚拟人切片信息，并将读取的数据信息保存在体数据集合中；S2：对虚拟人部分切片进行前景像素点和背景像素点的标记，将欲提取的部位设置为前景像素点，不提取的部分设置为背景像素点；S3：对读取到的体数据按照颜色划分到不同的颜色区间，将体数据中所有体素和每个颜色区间对应的辅助节点建立有向连通图，采用最大流算法对有向连通图中每个节点划分前景体素节点或背景体素节点；S4：对三维体数据的角度、阴影和光线参数进行设置只显示三维体数据的前景体素不显示背景体素。

3-D volume data segmentation method based on maximum flow strategy

【技术实现步骤摘要】
基于最大流策略的三维体数据分割方法
本专利技术涉及医学三维图像处理领域，尤其涉及基于最大流策略的三维体数据分割方法。
技术介绍
随着三维技术在医学上的广泛应用，使得早先的二维医学图像逐渐被三维成像所替代，如何简化直接提取三维图像中有效的数据的方法成为关键的问题。目前已知的方法是基于区域生长、分裂合并或者虚拟人切片技术的三维体数据分割算法，其中区域生长技术需要用大量的内存保存待生长的体数据点，时间复杂度较大。其中分裂合并技术需要依据用户需求的某种规则划分子区域，子区域划分的结果会直接影响分割效果，因此该算法具有一定的限定性。其中阈值分割技术是根据灰度阈值将三维体数据分割成多阈值体数据，从而进行分割的，由于阈值是在给定的某一范围内进行选择，所以精确性不能保证。
技术实现思路
根据现有技术中关于三维体数据分割方法精度低的问题，本专利技术公开了一种基于最大流策略的三维体数据分割方法，具体包括以下步骤：S1：读取虚拟人序列化切片图像中的虚拟人切片信息，并将读取的数据信息保存在体数据集合中；S2：对虚拟人部分切片进行前景像素点和背景像素点的标记，将欲提取的部位设置为前景像素点，不提取的部分设置为背景像素点；S3：对读取到的体数据按照颜色划分到不同的颜色区间，将体数据中所有体素和每个颜色区间对应的辅助节点建立有向连通图，采用最大流算法对有向连通图中每个节点划分前景体素节点或背景体素节点；S4：对三维体数据的角度、阴影和光线参数进行设置只显示三维体数据的前景体素不显示背景体素。进一步的，S3中具体采用如下方式：S31：将每个体素根据红、绿、蓝通道的值划分到不同的颜色区间，每个体素对应的颜色区间编号保存在一个数组中，将每个颜色区间包括的体素的编号保存到一个队列中；设体素f的颜色区间已包括其他的体素h，计算体素f和体素h的相似度，当两个体素相似时，体素f对应的颜色区间设置为体素h所属的颜色区间，把体素f加入到该颜色区间对应体素的队列中；当两个体素不相似时，继续判断体素f和其余有相同区间值的体素的相似性，直到该区间对应的所有体素判断完毕或找到相似的体素，如果没有相似体素，则增加一个新的颜色区间，最终所有的相似体素属于同一个颜色区间；S32：计算体数据集合中所有体素和其邻居体素的权重的方差；S33：为所有体素建立有向连通图，每个节点对应一个体素，称为体素节点，体素节点的值为对应体素的RGB颜色值，并为每个颜色区间内的所有体素建立一个辅助节点；S34：将S2中的前景像素点对应的体素和背景像素点对应的体素标记在有向连通图中，分别定义为前景体素节点和背景体素节点；S35：采用最大流算法对有向连通图进行分割：将所有已经标记过的体素节点加入到队列中，对队列中每个体素节点进行遍历扩展，即遍历其所有子节点，若子节点此时没有设置前景或背景类型，则该子节点需要被扩展，其中被扩展的子节点类型设置为与当前体素节点类型相同，把被扩展的子节点加入到队列中，并且判断该子节点是否通过当前体素节点到达源节点的距离最短，若不是则从该子节点到达源节点的路径中找到某节点M，该子节点通过节点M到达源节点距离最短，把节点M设为该子节点的父节点；当找到一条可行流后调整饱和边，其中可行流为从前景类型的体素节点经过一些节点后到达背景类型的体素节点的路径；即给饱和边指向的节点p找个新的父节点q，父节点q需要满足如下条件：到p的边的权重大于0，父节点q是已经被扩展过的节点，通过父节点q到达源节点的距离最近；如果找不到满足上述条件的节点，则重新寻找节点p的父节点、节点p的所有子节点的父节点和节点p的邻居节点的父节点。由于采用了上述技术方案，本专利技术提供的基于最大流策略的三维体数据分割方法，本方法只需要用户标记部分提取目标，简化了图像分割步骤，减少了人工操作，从而减少了一些由于人工失误导致的误差，使得准确性增加，并且占用内存较小，运行时间少。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术方法实现的流程图图2为本专利技术中体数据建立效果图图3为本专利技术中计算方差相邻点示意图图4为本专利技术中建立图myGraph的示意图图5为本专利技术中计算权重相邻点示意图图6为本专利技术中三维体数据分割后的效果图具体实施方式为使本专利技术的技术方案和优点更加清楚，下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚完整的描述：如图1所示的一种基于最大流策略的三维体数据分割方法，在实施过程中，体数据建立效果图如图2所示，图中体数据和坐标存在对应关系，可以根据坐标来获取每个数据，然后对体数据中的一些层的数据点进行标记，而后把每个数据点划分到不同颜色区间，并建立包括所有数据点和代表不同颜色区间的辅助节点的图，对图进行最大流处理，最后对数据点进行体绘制。本专利技术公开的方法具体步骤如下：S1：输入虚拟人序列化切片图像并存储体数据：对虚拟人切片进行读取，并把读取到的数据保存在体数据reader中。具体采用如下方式：由于体数据为长方体，可以设置体数据位于xyz坐标系中，每个数据间隔为1，体数据的建立如图2所示，所以每个数据和坐标有个对应关系，可以根据坐标来获取每个数据。例点A坐标为(ax,ay,az)，其对应的体数据编号为id＝ax*image.length+ay*image.width+az*height，其中image.length和image.width分别为图片的长和宽，height为图像的层数。S2：标记一些图像的前景和背景：对虚拟人部分切片进行前景像素点和背景像素点的标记。将要提取的部位设置为前景像素点，用白色标记；其余为背景像素点，用黑色标记。S3：三维体数据分割：对读取到的体数据按照颜色划分到不同的颜色区间，并将所有体数据的体素和每个颜色区间对应的辅助节点建立有向连通图采用最大流算法对有向连通图中每个节点划分前景体素节点或背景体素节点。具体采用如下方式：S31.将每个体素根据rgb值划分到不同的颜色区间，数组id_bin中保存每个体素对应的颜色区间编号。先计算每个体素的rgb值的区间其中，ri、gi、bi分别为该体素对应的rgb值，numBinsPerChannel为切片图像的位深度。并用变量newBinIdx记录此时已有颜色区间的数量，作为每个颜色区间的编号。将每个颜色区间包括的体素的编号保存到一个队列中。数组occupiedBins来标记该颜色区间是否已经有体素存在。若没有存在，则标记该颜色区间已有体素存在，编号为newBinIdx，并保存此时颜色区间对应体素的rgb值，该体素所属的颜色区间为newBinIdx，newBinIdx自加；否则要判断下该体素f＝(rf,gf,bf)'和该颜色区间已包括其他的体素h＝(rh,gh,bh)'是否应该划分到同一个颜色区间中，要计算下两个体素的相似度。相似度C(f,h)计算方式如下所示：C(f,h)＝λRR(f,h)+λBB(f,h)，其中设置λR＝0.85，λB＝0.15。若C(f,h)<0.995，则认为两个体素不相似，不属本文档来自技高网...

【技术保护点】
1.一种基于最大流策略的三维体数据分割方法，其特征在于包括以下步骤：S1：读取虚拟人序列化切片图像中的虚拟人切片信息，并将读取的数据信息保存在体数据集合中；S2：对虚拟人部分切片进行前景像素点和背景像素点的标记，将欲提取的部位设置为前景像素点，不提取的部分设置为背景像素点；S3：对读取到的体数据按照颜色划分到不同的颜色区间，将体数据中所有体素和每个颜色区间对应的辅助节点建立有向连通图，采用最大流算法对有向连通图中每个节点划分前景体素节点或背景体素节点；S4：对三维体数据的角度、阴影和光线参数进行设置只显示三维体数据的前景体素不显示背景体素。

【技术特征摘要】
1.一种基于最大流策略的三维体数据分割方法，其特征在于包括以下步骤：S1：读取虚拟人序列化切片图像中的虚拟人切片信息，并将读取的数据信息保存在体数据集合中；S2：对虚拟人部分切片进行前景像素点和背景像素点的标记，将欲提取的部位设置为前景像素点，不提取的部分设置为背景像素点；S3：对读取到的体数据按照颜色划分到不同的颜色区间，将体数据中所有体素和每个颜色区间对应的辅助节点建立有向连通图，采用最大流算法对有向连通图中每个节点划分前景体素节点或背景体素节点；S4：对三维体数据的角度、阴影和光线参数进行设置只显示三维体数据的前景体素不显示背景体素。2.根据权利要求1所述的一种基于最大流策略的三维体数据分割方法，其特征还在于：S3中具体采用如下方式：S31：将每个体素根据红、绿、蓝通道的值划分到不同的颜色区间，每个体素对应的颜色区间编号保存在一个数组中，将每个颜色区间包括的体素的编号保存到一个队列中；设体素f的颜色区间已包括其他的体素h，计算体素f和体素h的相似度，当两个体素相似时，体素f对应的颜色区间设置为体素h所属的颜色区间，把体素f加入到该颜色区间对应体素的队列中；当两个体素不相似时，继续判断体素f和其余有相同区间值的体素的相似性，直到该区间对应的所有体素判断完毕或找到相似的体素，如果没有相似体素，则增加一个新的颜色区间，最终所有的相似体素属于同一...

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