图像形态学的matlab实现 - 副本

file:///C|/Documents and Settings/Administrator/桌面/代码/第17章/代码17-1.txtclc;clear all;k=input('Enter the file name','s'); %输入彩色图像到程序中待处理im=imread(k); %原始彩色图像im1=rgb2gray(im); %灰度图像，im1=medfilt2(im1,[3 3]); %对图像进行中值滤波以便去除图像噪声BW = edge(im1,'sobel'); %寻找边界[imx,imy]=size(BW);msk=[0 0 0 0 0;0 1 1 1 0;0 1 1 1 0;0 1 1 1 0;0 0 0 0 0;];B=conv2(double(BW),double(msk)); %对图像进行平滑处理以便减少连通体个数L = bwlabel(B,8);% 计算连通体mx=max(max(L))% 将会有mx 个连通体。

这里你可以给出L的值为1到mx 中的一个，或者你可以通过一个循环提取所有的连通体 [r,c] = find(L==17);rc = [r c];[sx sy]=size(rc);n1=zeros(imx,imy);for i=1:sxx1=rc(i,1);y1=rc(i,2);n1(x1,y1)=255;end % 将提取的图像存储在数组阵列中figure,imshow(im);figure,imshow(im1);figure,imshow(B);figure,imshow(n1,[]);I = imread('test1.bmp');figure, imshow(I),BW1 = im2bw(I); %一般输入的图像是彩色的，需要转成黑白的se90 = strel('line', 3, 90);%线，长度为3，角度为90度se0 = strel('line', 3, 0); %线，长度为3，角度为0度BW2= imdilate(BW1, [se90 se0]);figure, imshow(BW2);BW3 = imerode(BW1,[se90 se0]);figure, imshow(BW3);% 这个程序实现了对二值图像的数学形态学操作% 第一部分：输入图像到MATLABfprintf('\n This program performs morphological operations on binary image') fprintf('\n Entering the image for MATLAB...')fprintf('\n Note that if You entered a gray scale or RGB image (not Binary),') fprintf('\n It will be converted to Binary Image by Thresholding')fprintf('\n Save the image or its copy in MATLAB working Directory') imagname = input('\n Enter the name of the image file (filename.ext) : ','s');a = imread(imagname);a = im2bw(a,0.5);siz = size(a);figureimshow(a)title('Input Image after conversion to Binary')fprintf('\n Click (1) to Perform Erosion')fprintf('\n Click (2) to Perform Dilation')fprintf('\n Click (3) to Perform Binary Opening')fprintf('\n Click (4) to Perform Binary Closing')fprintf('\n Click (5) to Subtract the Opening from the Input Image')fprintf('\n Click (6) to Subtract the Input Image from its Closing')fprintf('\n Click (7) to Remove isolated pixels (1''s surrounded by 0''s)') fprintf('\n Click (8) to Fill isolated interior pixels (0''s surrounded by 1''s)') fprintf('\n Click (9) to Leave only boundary pixels')fprintf('\n Click (10) to Shrink objects to points')fprintf('\n Click (11) to Make objects Thicker')fprintf('\n Click (12) to Make objects Thinner')choice = input('\n You select the Choice number : ');switch (choice)case 1a = bwmorph(a,'erode');case 2a = bwmorph(a,'dilate');case 3a = bwmorph(a,'open');case 4a = bwmorph(a,'close');case 5a = bwmorph(a,'tophat');case 6a = bwmorph(a,'bothat');case 7a = bwmorph(a,'clean');case 8a = bwmorph(a,'fill');case 9a = bwmorph(a,'remove');case 10a = bwmorph(a,'shrink',Inf);case 11a = bwmorph(a,'thicken',Inf);case 12a = bwmorph(a,'thin',Inf);otherwisefprintf('\n Sorry! Wrong Choice')endfigureimshow(a)title('Result Image')sav = input('\n Do you like to SAVE Result Image? (y/n) : ','s');if (sav == 'y')fprintf('\n You choose to SAVE the Result Image')naming = input('\n Type the name of the new image file (filename.ext) : ','s');fprintf('\n Saving ...')imwrite(a,naming);fprintf('\n The new file is called %s and it is saved in MATLAB working Directory',naming) elsefprintf('\n You choose NOT to SAVE the Result Image')endclearBW = zeros(400,500); BW(40:300,40:400) = 1; figure, imshow(BW,[]);D = bwdist(BW); figure, imshow(D,[]);I = imread('2.jpg');BW1 = im2bw(I); %一般输入的图像是彩色的，需要转成黑白的D = bwdist(BW1);figure, imshow(D,[]), title('Distance transform of 2.jpg')I = imread('test1.bmp');BW1 = im2bw(I); %一般输入的图像是彩色的，需要转成黑白的X = bwlabel(BW1,4);RGB = label2rgb(X, @jet, 'k');imshow(RGB,'notruesize')f = @(x) sum(x(:)) >= 3; lut = makelut(f,3);I = imread('text.jpg'); BW1 = im2bw(I);BW2 = applylut(BW1,lut); imshow(BW1);% figure, imshow(BW2)function intrec(Z)%INTREC 实现数学形态学字符识别A= imread ('zero.bmp');%将0字符图像读入到MATLAB中B= imread ('one.bmp'); %将1字符图像读入到MATLAB中C= imread ('two.bmp'); %将2字符图像读入到MATLAB中D= imread ('three.bmp'); %将3字符图像读入到MATLAB中E= imread ('four.bmp'); %将4字符图像读入到MATLAB中F= imread ('five.bmp'); %将5字符图像读入到MATLAB中G= imread ('six.bmp'); %将6字符图像读入到MATLAB中H= imread ('seven.bmp'); %将7字符图像读入到MATLAB中I= imread ('eight.bmp'); %将8字符图像读入到MATLAB中J= imread ('nine.bmp'); %将9字符图像读入到MATLAB中SE = strel('square',3); %构造结构元素对象K=imdilate(A,SE);A2=K-A; %矩阵A进行膨胀， 得到矩阵A2L=imdilate(B,SE);B2=L-B; %对字符图像矩阵B进行膨胀运算，得到矩阵B2M=imdilate(C,SE);C2=M-C; %对字符图像矩阵C进行膨胀运算，得到矩阵C2 N=imdilate(D,SE);D2=N-D; %对字符图像矩阵D进行膨胀运算，得到矩阵D2 O=imdilate(E,SE);E2=O-E; %对字符图像矩阵E进行膨胀运算，得到矩阵E2 P=imdilate(F,SE);F2=P-F; %对字符图像矩阵F进行膨胀运算，得到矩阵F2Q=imdilate(G,SE);G2=Q-G; %对字符图像矩阵G进行膨胀运算，得到矩阵G2 R=imdilate(H,SE);H2=R-H; %对字符图像矩阵H进行膨胀运算，得到矩阵H2 S=imdilate(I,SE);I2=S-I; %对字符图像矩阵I进行膨胀运算，得到矩阵I2T=imdilate(J,SE);J2=T-J; %对字符图像矩阵J进行膨胀运算，得到矩阵J2disp('The digits found in the image are:');if ~isempty(nonzeros(bwhitmiss(Z,A,A2))) %击中击不中运算disp('0');endif ~isempty(nonzeros(bwhitmiss(Z,B,B2))) %击中击不中运算disp('1');endif ~isempty(nonzeros(bwhitmiss(Z,C,C2))) %击中击不中运算disp('2');endif ~isempty(nonzeros(bwhitmiss(Z,D,D2))) %击中击不中运算disp('3');endif ~isempty(nonzeros(bwhitmiss(Z,E,E2))) %击中击不中运算disp('4');endif ~isempty(nonzeros(bwhitmiss(Z,F,F2))) %击中击不中运算disp('5');endif ~isempty(nonzeros(bwhitmiss(Z,G,G2))) %击中击不中运算 disp('6');endif ~isempty(nonzeros(bwhitmiss(Z,H,H2))) %击中击不中运算 disp('7');endif ~isempty(nonzeros(bwhitmiss(Z,I,I2))) %击中击不中运算 disp('8');endif ~isempty(nonzeros(bwhitmiss(Z,J,J2))) %击中击不中运算 disp('9');end。