大三下all dip实验图像2017.pptxVIP

数字图像处理实验实验十、图像分割1

3.约束最小二乘滤波恢复（提示：按照书P268－272操作）模糊和加噪；利用真实PSF函数和不同的噪声强度作为参数，对模糊加噪图像进行约束最小二乘滤波恢复；对噪声估计过小引起的振铃现象进行进一步处理；利用不同的拉各朗日乘法算子进行恢复。

一、图像阈值分割（P204页）读入一幅灰度图像，先做直方图(imhist)，然后对其进行阈值分割(im2bw)。由直方图确定阈值Inew=im2bw(I,level);%level为归一化阈值。

二、图像边缘检测1.读入一幅灰度图像，分别对图像进行Robert、Sobel、Prewitt三种算子边缘检测(edge).(P208页)edge(I,’method’)%method是可以Sobel、Prewitt、Roberts等

2.画LOG算子(拉普拉斯-高斯算子)的三维图形，并用该算子对灰度图像进行边缘检测。（P209-210页）LOG算子：原图像：高斯运算：拉普拉斯变换：

(1)LOG算子三维图形—P210页中间部分程序段:r=1/(2*pi*sigma^4)*((x.^2+y.^2)/(sigma^2)-2).*...exp(-(x.^2+y.^2)/(sigma^2));mesh(x,y,r)(2)LOG算子边缘检测:edge(I,’log’);LOG算子Matlab实现:

三、图像分割的应用1.通过图像分割检测细胞(P217页)(1)读入图像；I=imread(‘cell.tif’)图像处理工具包内置的图像

8E=edge(I,’sobel’,(graythresh(I)*0.1);%graythresh()用最大类间方差法获得合适的阈值。(2)细胞边缘检测边缘不连续

(3)确定结构单元(用于填补边缘缝隙)se90=strel(‘line’,3,90);se0=strel(‘line’,3,0);%strel()构造结构元素，line为线性结构长度角度

10(4)用结构单元进行膨胀操作Ed=imdilate(E,[se90,se0]);%用结构元素进行膨胀操作；imdilate()对图像进行膨胀操作存在孔洞

(5)填充空洞Edf=imfill(Ed,’holes’);%imfill()对图像进行填充操作。存在与图像边界相连的对象

12(6)移除与边界连通的目标Edfb=imclearborder(Edf,4);%imclearborder()抑制与图像边界相连的亮对象。边缘不平滑，有突起或下凹

(7)腐蚀操作进行平滑seD=strel(‘diamond’,1)C=imerode(Edfb,seD);%imerode()对图像进行平滑处理(8)标示结果F=bwperim(C);I(F)=255;%bwperim查找二值图像的边缘F保存了细胞边缘的坐标值可以多次操作边缘点置白

2.测试图像的微小结构（P219页）(1)读入图像I；(2)对图像求补(3)对图像进行阈值分割，得到图像中的微小结构。Ic=plement(I);BW=im2bw(Ic,graythresh(Ic));%用最大类间方差法获得的阈值

(4)利用闭、开操作，得到大结构图像；Se=strel(‘disk’,6);%确定结构元素BWc=imclose(BW,Se);%形态闭合操作BWo=imopen(BWc,Se);%形态开操作(5)小结构图像和大结构图像进行“与”操作，得到原始图像中的所有小对象。Mask=BWBWo;

图像开操作：使图像的轮廓变得光滑，断开狭窄的间断和消除细的突出物；图像闭操作：使图像的轮廓变得光滑，消除狭窄的间断和长细的鸿沟，消除小的空洞，并填补轮廓线中的裂痕。闭操作结果

期末考试时间：5月25日形式：上机实验考试内容：与实验1-10内容差不多；要求：可带资料（书、实验报告等），禁止带移动存储设备（U盘等），考试中不能上网。