医学图像分割.ppt

《医学图像分割.ppt》由会员分享，可在线阅读，更多相关《医学图像分割.ppt（54页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、医学图像分割 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望讨论内容图像分割概述阈值分割1、图像分割概述将不同区域区分开来，这些区域是互不相交的，每一个区域都满足特定区域的一致性。其分割的目的是为了将感兴趣区域提取出来，从而为定量、定性分析提供基础，同时它也是三维可视化的基础。1、图像分割概述 P(gk(x,y)U gj(x,y)=FALSE.任意相邻部分的合并都会破坏这种一致性。任意相邻部分的合并都会破坏这种一致性。1、图像分割概述 如果连通性的约束被取消，那么

2、对像素集合如果连通性的约束被取消，那么对像素集合的划分就称为分类的划分就称为分类(Classification),(Classification),每一个像素每一个像素集称为类集称为类(Class)(Class)。经典的分割和像素分类通称为分割。经典的分割和像素分类通称为分割。基于区域的分割方法 基于边缘的分割方法 结合区域与边界信息的方法 基于模糊集理论的方法 基于神经网络的方法 基于数学形态学的方法图谱引导（Atlas-guided）方法1、图像分割概述1、图像分割概述基于区域的分割方法 利用区域内的相似性（一致性）阈值分割区域生长和分裂合并分类器和聚类基于随机场的方法其它基于统计学的方法

3、1、图像分割概述基于边缘的分割方法 利用区域之间差异性 并行微分算子曲面拟合法基于边界曲线拟合的方法串行边界查找医学图像特点:模糊、不均匀、个体差异、复杂多样灰度不均匀:不均匀的组织器官、磁场等伪影和噪声:成像设备局限性、组织的蠕动边缘模糊 :局部体效应边缘不明确:病变组织1、图像分割概述局部体效应(partial volume effects)1、图像分割概述Ideal ImageAcquired Image医学图像分割方法的公共特点：分割算法面向具体的分割任务，没有通用的方法更加重视多种分割算法的有效结合需要利用医学中的大量领域知识交互式分割方法受到日益重视 医学图像分割是一项十分困难的任

4、务，至今仍然没有获得圆满的解决。1、图像分割概述2、阈值分割阈值分割是最常见的一种分割方法。它基于对灰度图像的一种假设：目标或背景内的相邻象素间的灰度值是相似的，但不同目标或背景的象素在灰度上有差异，反映在图像的直方图上，不同目标和背景则对应不同的峰。选取的阈值应位于两个峰之间的谷，从而将各个峰分开CTCT图像图像中皮肤中皮肤骨骼的骨骼的分割分割2、阈值分割阈值分割的三种技术方案直接门限法间接门限法 对图像进行预处理后再运用门限法。拉氏或梯度运算，邻域平均多门限法2、阈值分割多门限法2、阈值分割乳腺钼靶图像乳腺钼靶图像单门限分割单门限分割多门限分割多门限分割 门限的确定方法 根据直方图确定门限

5、最小误判概率准则下的最佳门限最大类间距准则下的最佳门限最大类间类内距离比准则下的最佳门限最大熵准则下的最佳门限根据二维直方图确定图像分割门限边缘灰度作为分割门限分水岭方法2、阈值分割阈值分割的优点 简单，常作为预处理方法阈值分割的缺点不适用于多通道图像不适用于特征值相差不大的图像不适用于各物体灰度值有较大重叠的图像对噪声和灰度不均匀敏感2、阈值分割ThresholdingThe simplest and most efficient image segmentation method is thresholding.Thresholding is to segment the image in

6、to two regions according to the gray level of image pixels.If the gray level is higher than the given threshold T,the output at this pixel is set to 1,otherwise it is set to 0.Image ThresholdingOriginal image Segmented image(T=128,145)Determination of ThresholdIn thresholding method,the most difficu

7、lt is to determine a proper value of the threshold.There are different types of the threshold:Global threshold(constant threshold)Adaptive thresholdDetermination of Global thresholdIf the object and background have different distributions,the value of the global threshold can be determined by calcul

8、ating the histogram of the image.The global threshold can also be determined interactively.The threshold can also be determined by optimization.Determination of the globalthreshold from histogramT=150The Otsu Algorithml If t is chosen as a threshold,and p(i)is the normalized histogram0K-1N bits mean

9、s K=2NtThe Otsu AlgorithmmeansvariancesMeans and variance for each classThe Otsu AlgorithmlStatistical discrimination measure based on variance between classes:lRun through all possible values of t,and pick the one that maximizes the discrimination measure:Chosen Threshold The Otsu AlgorithmFor each

10、 potential threshold T,1.Separate the pixels into two clusters according to the threshold.2.Find the mean of each cluster.3.Square the difference between the means.4.Calculate the object function of .5.Find the optimal threshold T*that maximizes the value of .Determination of Otsus thresholdAutomati

11、c Threshold based on mean and standard deviationwAutomatic threshold based on mean and standard deviation:where are the automatic threshold at the point(i,j),the mean and standard deviation of the neighbors of(i,j),i.e.,a local window,k is the weight and can be a real number.Determination of thresho

12、ld by maximum entropylWhat is an entropy?lEntropy is the measurement of the information content in a probability distributionlMaximum entropy segmentation is to select such a threshold that the entropies in both object and background areas have maximum distributions.根据二维直方图确定图像分割门限 灰度平均灰度直方图 平均灰度局部方

13、差直方图 最大熵 灰度梯度直方图 采用聚类的方法，分三类 平均灰度局部方差直方图 最大熵Determination of threshold by 2-D HistogramlDefinition of 2D histogram:Suppose f(x,y)to be an image of NxN pixels.Its gray level is from 0 to L-1.Segment the image by using the following equation:where lFor the 2D thresholding method,it considers the avera

14、ge gray level of the point(x,y)simultaneously as follows.Determination of threshold by 2-D HistogramlThe average gray level at the point(x,y)of its nxn neighbors is:where lFor the 2D thresholding method,it considers the average gray level of the point(x,y)simultaneously,i.e.,use(f(x,y),g(x,y)to repr

15、esent an image and to segment the image with 2D vector threshold(S,T):Determination of threshold by 2-D Histogram where lFor one image,let rij to be the occurrence number of gray level i and the average gray level j,we can define the joint probability as:lP is called the 2D histogram of the image f(

16、x,y)Determination of threshold by 2-D Histogram If the threshold vector is(S,T),the 2D histogram will be divided into 4 parts:In Part 0 and Part 1,i.e.,the object or background,the gray level and the average is close,while in Part 2 and part 3,the difference between the gray level and the average is

17、 big,which is corresponding to the boundary points.2D histogram of imageDetermination of threshold by 2-D HistogramlThe maximum entropy for the 2D histogram is to determine a threshold vector(S,T)such that we can divide the image into object(A)and background(B)with the probability of where Determina

18、tion of threshold by 2-D HistogramlThe goal of segmentation is to let the entropies in the object and background areas as big as possible,lThe maximum entropies of the object and background will correspond to the optimal threshold vector(S,T).Determination of threshold by 2-D Histogram-ExperimentDet

19、ermination of threshold by Fuzzy EntropylThe BlockB and BlockW are defined in Fig.1(a)and(b).Four fuzzy sets,BrightX,DarkX,BrightY,DarkY,are defined based on the S-function and the corresponding Z-functions as follows:(Z()=1-s()Determination of Threshold by Fuzzy EntropyDetermination of threshold by

20、 Fuzzy EntropylThe fuzzy relation Bright is a subset of the full Cartesian product space XY lSimilarly,Definition of Fuzzy EntropylLet A be a fuzzy set with membership function ,where are the possible outputs from source A with the probability.The fuzzy entropy set A is defined as:lThe total image e

21、ntropy is defined as:Determination of threshold by Fuzzy EntropylAs shown in Fig.1(a),the dark block BlockB can be divided into a nonfuzzy region RB and a fuzzy region R1lSimilarly,the bright block BlockW is composed of a nonfuzzy region RW and a fuzzy region R2,as shown in Fig.1(b)Determination of

22、threshold by Fuzzy EntropylThe following four entropies can be calculated:where nxy is the element in the 2-D histogram which represents the number of occurences of the pair(x,y)lTo find the best set of a,b,and c is an optimization problem which can be solved by different optimization methods.For ex

23、ample,we can use genetic algorithm to search for the optimal solution.The proposed method consists of the following three major steps:1)find the 2-D histogram of the image;2)perform fuzzy partition on the 2-D histogram;3)compute the fuzzy entropy.lStep 1)needs to be execute only once while Steps 2)a

24、nd 3)are performed iteratively for each set of(a,b,c).The optimum(a,b,c)determines the fuzzy region(i.e.,interval a,c).The threshold is selected as the crossover point of the membership function which has membership 0.5 implying the largest fuzziness.Determination of threshold by Fuzzy EntropyDeterm

25、ination of threshold by Fuzzy EntropyDetermination of threshold by Fuzzy Entropy-Experiment1Comparison of global and local threshold segmentationDetermination of threshold by Fuzzy Entropy-Experiment2HP DCE/9000K-means clusteringK-means follow a simple and easy way to classify a given data set throu

26、gh a certain number of clusters(assume k clusters)fixed a priori.The main idea is to define k centroids,one for each cluster.These centroids shoud be placed in a cunning way because of different location causes different result.So,the better choice is to place them as much as possible far away from

27、each other.The next step is to take each point belonging to a given data set and associate it to the nearest centroid.When no point is pending,the first step is completed and an early groupage is done.At this point we need to re-calculate k new centroids as barycenters of the clusters resulting from

28、 the previous step.After we have these k new centroids,a new binding has to be done between the same data set points and the nearest new centroid.A loop has been generated.As a result of this loop we may notice that the k centroids change their location step by step until no more changes are done.In

29、 other words centroids do not move any more.K-means clustering Finally,this algorithm aims at minimizing an objective function,in this case a squared error function.The objective function where is a chosen distance measure between a data point xji and the cluster centre cj,is an indicator of the dis

30、tance of the n data points from their respective cluster centroids.K-means clustering AlgorithmThe algorithm is composed of the following steps:1.Place K points into the space represented by the objects that are being clustered.These points represent initial group centroids.2.Assign each object to t

31、he group that has the closest centroid.3.When all objects have been assigned,recalculate the positions of the K centroids.4.Repeat Steps 2 and 3 until the centroids no longer move.This produces a separation of the objects into groups from which the metric to be minimized can be calculated.Fuzzy K-me

32、ans clustering Fuzzy K-means clustering algorithm aims at minimizing the following objective function with respect to the membership function ij and the centroids cj:(1)where K is a number of clusters or classes,n is the total number of feature points or vectors and is a weighting exponent.And we ha

33、ve:(2)Fuzzy K-means clustering wFor each input vector xji,using Lagrangian multiplier method,for m1,local minimum solutions of equation(1)was demonstrated if and only if:阈值分割的改进2、阈值分割利用像素邻域的局部信息：基于过渡区的方法利用像素邻域的局部信息：基于过渡区的方法利用像素点空间位置：变化阈值法利用像素点空间位置：变化阈值法结合局部灰度结合局部灰度结合连通信息结合连通信息基于最大熵原则的阈值选择方法基于最大熵原则的阈值选择方法