【正文】 termins a class. This set should include a number of images for each person, with some variations in expression and in the lighting.39。 while chos~=possibility, chos=menu(39。邊緣檢測,39。Roberts39。,39。Prewitt39。,39。Kirsch39。,39。Laplacian39。,39。LoG39。,39。Canny39。,39。Sobel39。,39。八方向Sobel39。, 39。改進的Sobel39。,39。加入噪聲39。,39。退出39。)。 %========================================================================== if chos==1 clc。 close all。 selezionato=0。 while selezionato==0 [namefile,pathname]=uigetfile({39。*.bmp。*.tif。*.tiff。*.jpg。*.jpeg。*.png。*.gif39。,39。IMAGE Files (*.bmp,*.tif,*.tiff,*.jpg,*.jpeg,*.png,*.gif)39。},39。Chose GrayScale Image39。)。 if namefile~=0 [img,map]=imread(strcat(pathname,namefile))。 selezionato=1。 else disp(39。Select a grayscale image39。)。 end if (any(namefile~=0) amp。amp。 (~isgray(img))) disp(39。Select a grayscale image39。)。 selezionato=0。 end end figure(39。Name39。,39。Selected image39。)。 imshow(img)。 if isa(img,39。uint839。) graylevmax=2^81。 end if isa(img,39。uint1639。) graylevmax=2^161。 end if isa(img,39。uint3239。) graylevmax=2^321。 end end %========================================================================== if chos==2clc。 close all。 A = img。 x_mask = [1 0。0 1]。 y_mask = rot90(x_mask)。 I = im2double(A)。 dx = imfilter(I, x_mask)。 dy = imfilter(I, y_mask)。 grad = sqrt(dx.*dx + dy.*dy)。 grad = mat2gray(grad)。 level = graythresh(grad)。 BW = im2bw(grad,level)。 figure, imshow(BW)。 title(39。Roberts39。) end %========================================================================== if chos==3 clc。 close all。 A = img。 y_mask = [1 1 1。0 0 0。1 1 1]。 x_mask = y_mask39。 I = im2double(A)。 dx = imfilter(I, x_mask)。 dy = imfilter(I, y_mask)。 grad = sqrt(dx.*dx + dy.*dy)。 grad = mat2gray(grad)。 level = graythresh(grad)。 BW = im2bw(grad,level)。 figure, imshow(BW)。 title(39。Prewitt39。) end %========================================================================= if chos==4 clc。 close all。 A = img。 mask1=[3,3,3。3,0,5。3,5,5]。 mask2=[3,3,5。3,0,5。3,3,5]。 mask3=[3,5,5。3,0,5。3,3,3]。 mask4=[3,3,3。3,0,3。5,5,5]。 mask5=[5,5,5。3,0,3。3,3,3]。 mask6=[3,3,3。5,0,3。5,5,3]。 mask7=[5,3,3。5,0,3。5,3,3]。 沈陽大學(xué)畢業(yè)設(shè)計(論文) No 46 mask8=[5,5,3。5,0,3。3,3,3]。 I = im2double(A)。 d1 = imfilter(I, mask1)。 d2 = imfilter(I, mask2)。 d3 = imfilter(I, mask3)。 d4 = imfilter(I, mask4)。 d5 = imfilter(I, mask5)。 d6 = imfilter(I, mask6)。 d7 = imfilter(I, mask7)。 d8 = imfilter(I, mask8)。 dd = max(abs(d1),abs(d2))。 dd = max(dd,abs(d3))。 dd = max(dd,abs(d4))。 dd = max(dd,abs(d5))。 dd = max(dd,abs(d6))。 dd = max(dd,abs(d7))。 dd = max(dd,abs(d8))。 grad = mat2gray(dd)。 level = graythresh(grad)。 BW = im2bw(grad,level)。 figure, imshow(BW)。 title(39。Kirsch39。) end %========================================================================== if chos==5 clc。 close all。 A=img。 mask=[0,1,0。1,4,1。0,1,0]。 I = im2double(A)。 dx = imfilter(I, mask)。 grad = mat2gray(dx)。 BW = im2bw(grad,)。 figure, imshow(BW)。 title(39。Laplacian39。) end %========================================================================== if chos==6 clc。 close all。 A=img。 mask=[0,0,1,0,0。0,1,2,1,0。1,2,16,2,1。0,1,2,1,0。0,0,1,0,0]。 I = im2double(A)。 dx = imfilter(I, mask)。 grad = mat2gray(dx)。 BW = im2bw(grad,)。 figure, imshow(BW)。 title(39。log39。) end %========================================================================== if chos==7 clc。 close all。 A=img。 BW1 = edge(I,39。canny39。)。 figure,imshow(BW1)。 title(39。Canny39。) end %========================================================================= if chos==8 clc。 close all。 A=img。 y_mask = [1 2 1。0 0 0。1 2 1]。 x_mask = y_mask39。 I = im2double(A)。 dx = imfilter(I, x_mask)。 dy = imfilter(I, y_mask)。 grad = sqrt(dx.*dx + dy.*dy)。 grad = mat2gray(grad)。 level = graythresh(grad)。 BW = im2bw(grad,level)。 figure, imshow(BW)。 title(39。sobel39。) end %========================================================================== if chos==9 clc。 close all。 A=img。 mask1=[1 2 1。0 0 0。1 2 1]。 mask2=[2 1 0。1 0 1。0 1 2]。 mask3=[1 0 1。2 0 2。1 0 1]。 mask4=[0 1 2。1 0 1。2 1 0]。 mask5=[1 2 1。0 0 0。1 2 1]。 mask6=[2 1 0。1 0 1。0 1 2]。 mask7=[1 0 1。2 0 2。1 0 1]。 mask8=[0 1 2。1 0 1。2 1 0]。 I = im2double(A)。 d1 = imfilter(I, mask1)。 d2 = imfilter(I, mask2)。 d3 = imfilter(I, mask3)。 d4 = imfilter(I, mask4)。 d5 = imfilter(I, mask5)。 d6 = imfilter(I, mask6)。 d7 = imfilter(I, mask7)。 d8 = imfilter(I, mask8)。 dd = max(abs(d1),abs(d2))。% dd = max(dd,abs(d3))。 dd = max(dd,abs(d4))。 dd = max(dd,abs(d5))。 dd = max(dd,abs(d6))。 dd = max(dd,abs(d7))。 dd = max(dd,abs(d8))。 grad = mat2gray(dd)。 level = graythresh(grad)。 BW = im2bw(grad,level)。 figure, imshow(BW)。 title(39。sobel176。242。39。) end %========================================================================== if chos==10 clc。 close all。 A=img。 mask1=[1 2 1。0 0 0。1 2 1]。 mask2=[2 1 0。1 0 1。0 1 2]。 mask3=[1 0 1。2 0 2。1 0 1]。 mask4=[0 1 2。1 0 1。2 1 0]。 mask5=[1 2 1。0 0 0。1 2 1]。 mask6=[2 1 0。1 0 1。0 1 2]。 mask7=[1 0 1。2 0 2。1 0 1]。 mask8=[0 1 2。1 0 1。2 1 0]。 I = im2double(A)。 d1 = imfilter(I, mask1)。 d2 = imfilter(I, mask2)。 d3 = imfilter(I, mask3)。 d4 = imfilter(I, mask4)。 d5 = imfilter(I, mask5)。 d6 = imfilter(I, mask6)。 d7 = imfilter(I, mask7)。 d8 = imfilter(I, mask8)。 dd = max(abs(d1),abs(d2))。 dd = max(dd,abs(d3))。 dd = max(dd,abs(d4))。 dd = max(dd,abs(d5))。 dd = max(dd,abs(d6))。 dd = max(dd,abs(d7))。 dd = max(dd,abs(d8))。 grad = mat2gray(dd)。 level = graythresh(grad)。 BW = im2bw(grad,level)。 BB = grad。 FW=median(BB(:))/。 B = BB.*BB。 B = sum(B(:))。 [M,N] = size(A)。 FX= sqrt(B/(M*N))。 FS=sqrt(max(abs(FX^2FW^2)))。 T=sqrt(2)*FW^2/FS。 grad = mat2gray(BB)。 BW2=im2bw(grad,T)。 figure, imshow(BW2)。 title(39。sobel39。)。 end %======================================================================== if chos==11 clc。 close all。 img=imnoise(img,39。gaussian39。,)。 % imshow(img)。 end end