细胞边缘的精确检测

1、blood = imread(blood1.BMP);x,y=size(blood); % 求出图象大小b=double(blood); N =sqrt(100) * randn(x,y); % 生成方差为10的白噪声I=b+N; % 噪声干扰图象 for i=1:x % 实际图象的灰度为0255 for j=1:y if (I(i,j)255) I(i,j)=255; end if (I(i,j)=allow) % 迭代最佳阈值分割算法 count=count+1; for i=1:x for j=1:y if (I(i,j)=T) S0=S0+I(i,j); n0=n0+1; end if

2、 (I(i,j)=T) Seg(i,j)=1; % 阈值分割的图象 end endend SI=1-Seg; % 阈值分割后的图象求反，便于用腐蚀算法求边缘se1=strel(square,3); % 定义腐蚀算法的结构SI1=imerode(SI,se1); % 腐蚀算法BW=SI-SI1; % 边缘检测%=传统的边缘检测方法=%I=uint8(I);BW1=edge(I,sobel, 0.11);BW2=edge(I,log, 0.015);%=图象显示=%figure(1);imshow(I);title(Original) % 显示阈值分割的图象figure(2);imshow(BW2

3、);title(Soble) % 显示新算法的边缘图象figure(3)imshow(BW1);title(Gauss-Laplace) figure(4)imshow(BW);title(New algorithm) 程序代码解释：关于randn产生正态分布的随机数，方差为1，期望为0.所以要产生需要的方差和期望，可以这样实现：产生一个随机分布的指定均值和方差的矩阵：将randn产生的结果乘以标准差，然后加上期望均值即可。例如，产生均值为0.6，方差为0.1的一个5*5的随机数方式如下： x = .6 + sqrt(0.1) * randn(5) 关于strel：图像处理函数详解strel功

4、能：用于膨胀腐蚀及开闭运算等操作的结构元素对象(本论坛随即对膨胀腐蚀等操作进行讲解)。用法：SE = strel(shape,parameters)创建由指定形状shape对应的结构元素。其中shape的种类有arbitrarypairdiamondperiodiclinediskrectanglelinesquareoctagon参数parameters一般控制SE的大小。例子：se1 = strel(square,6)% 创建6*6的正方形se2 = strel(line,10,45)% 创建直线长度10，角度45se3 = strel(disk,15)% 创建圆盘半径15se4 = st

5、rel(ball,15,5)% 创建椭圆体，半径15，高度5关于腐蚀：功能:腐蚀图像 用法：IM2 = imerode(IM,SE) 腐蚀灰度、二进制或压缩二进制图像 IM ，返回腐蚀图像 IM2 。参数 SE 是函数 strel 返回的一个结构元素体或是结构元素体阵列。如果 IM 是逻辑图像，结构元素是平坦的， imerode 执行二进制膨胀;否则执行灰度腐蚀。如果SE 是一个结构元素体阵列， imerode 接连使用 SE 中的每个结构元素体对输入图像进行多重腐蚀。IM2 = imerode(IM,NHOOD) 腐蚀图像 IM ，这里的 NHOOD 是一个 0 和 1 组成的阵列，指明了结构元素邻域。这等同于 imerode(IM,strel(NHOOD) 的句法规则。函数 imerode 通过 floor(size(NHOOD)+1)/2) 确定了邻域的中心元素。IM2 = imerode(.,PACKOPT,M)指定 IM 是否为一幅压缩二进制图像，如果是，需要提供原始未压缩图像的行的维数 M 。 PACKOPT 可以取以下值之一。默认值由大括号()标出。 Value(值) Description(描述) ispacked IM 可以作为一幅由bwpack 产生的压缩二进制图像。IM 必须是一个2-D 的