【正文】 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==2 clc。 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]。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。203。189。207。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。184。196。189。248。181。196。sobel39。)。 end%======================================================================== if chos==11 clc。 close all。 img=imnoise(img,39。gaussian39。,)。 % imshow(img)。 endend 51 沈陽大學(xué)畢業(yè)設(shè)計(jì)（論文） No