【正文】 。u(25)=image(i+2,j+2)。 u(27)=image(i2,j+4)。 u(29)=image(i1,j+4)。 u(32)=image(i+1,j+3)。u(34)=image(i+2,j+3)。 u(36)=image(i+3,j2)。 u(38)=image(i+3,j)。u(40)=image(i+3,j+2)。 u(42)=image(i+3,j+4)。 u(44)=image(i+4,j1)。u(46)=image(i+4,j+1)。 u(48)=image(i+4,j+3)。 for p=1:49, for q=1:49p, if u(q)u(q+1)k=u(q)。 u(q+1)=k。end,end,d=uint8(image)。7*7高斯噪聲中值濾波圖像39。 圖33 正方形7*7中值濾波器對椒鹽噪聲和高斯噪聲的去噪圖 （4） 正方形9*9中值濾波器去噪程序I=imread(39。)。B=rgb2gray(I)。Original Image39。J = imnoise(B,39。 pepper39。 %noise density=K= imnoise(B,39。,)。Salt amp。)。Gaussian Noise39。image=double(J)。for i=5:M4, for j=5:N4, u(1)=image(i,j)。 u(3)=image(i1,j+1)。 u(5)=image(i1,j1)。u(7)=image(i+1,j1)。 u(9)=image(i+1,j+1)。 u(11)=image(i2,j1)。u(13)=image(i2,j+1)。 u(15)=image(i1,j2)。 u(17)=image(i,j2)。u(19)=image(i+1,j2)。 u(21)=image(i+2,j2)。 u(23)=image(i+2,j)。u(25)=image(i+2,j+2)。 u(27)=image(i2,j+4)。 u(29)=image(i1,j+4)。 u(32)=image(i+1,j+3)。u(34)=image(i+2,j+3)。 u(36)=image(i+3,j2)。 u(38)=image(i+3,j)。u(40)=image(i+3,j+2)。 u(42)=image(i+3,j+4)。 u(44)=image(i+4,j1)。u(46)=image(i+4,j+1)。 u(48)=image(i+4,j+3)。 u(50)=image(i2,j3)。u(52)=image(i,j3)。 u(54)=image(i+2,j3)。u(56)=image(i+4,j3)。u(58)=image(i1,j4)。u(60)=image(i+1,j4)。u(62)=image(i+3,j4)。u(64)=image(i3,j2)。 u(66)=image(i3,j)。 u(68)=image(i3,j+2)。u(70)=image(i3,j+4)。 u(72)=image(i3,j4)。 u(74)=image(i4,j1)。u(76)=image(i4,j+1)。 u(78)=image(i4,j+3)。 u(80)=image(i4,i3)。 for p=1:81, for q=1:81p, if u(q)u(q+1)k=u(q)。u(q+1)=k。end,end,d=uint8(image)。9*9椒鹽噪聲中值濾波圖像39。 image=double(K)。for i=5:M4, for j=5:N4, u(1)=image(i,j)。 u(3)=image(i1,j+1)。 u(5)=image(i1,j1)。u(7)=image(i+1,j1)。 u(9)=image(i+1,j+1)。 u(11)=image(i2,j1)。u(13)=image(i2,j+1)。 u(15)=image(i1,j2)。 u(17)=image(i,j2)。u(19)=image(i+1,j2)。 u(21)=image(i+2,j2)。 u(23)=image(i+2,j)。u(25)=image(i+2,j+2)。 u(27)=image(i2,j+4)。 u(29)=image(i1,j+4)。 u(32)=image(i+1,j+3)。u(34)=image(i+2,j+3)。 u(36)=image(i+3,j2)。 u(38)=image(i+3,j)。u(40)=image(i+3,j+2)。 u(42)=image(i+3,j+4)。 u(44)=image(i+4,j1)。u(46)=image(i+4,j+1)。 u(48)=image(i+4,j+3)。 u(50)=image(i2,j3)。u(52)=image(i,j3)。 u(54)=image(i+2,j3)。u(56)=image(i+4,j3)。u(58)=image(i1,j4)。u(60)=image(i+1,j4)。u(62)=image(i+3,j4)。u(64)=image(i3,j2)。 u(66)=image(i3,j)。 u(68)=image(i3,j+2)。u(70)=image(i3,j+4)。 u(72)=image(i3,j4)。 u(74)=image(i4,j1)。u(76)=image(i4,j+1)。 u(78)=image(i4,j+3)。 u(80)=image(i4,i3)。 for p=1:81, for q=1:81p, if u(q)u(q+1)k=u(q)。 u(q+1)=k。end,end,d=uint8(image)。9*9高斯噪聲中值濾波圖像39。正方形9*9中值濾波器對椒鹽噪聲和高斯噪聲去噪 圖34 長方形中值濾波器去噪（1）長方形3*5中值濾波器去噪程序I=imread(39。)。B=rgb2gray(I)。Original Image39。J = imnoise(B,39。 pepper39。 %noise density=K= imnoise(B,39。,)。Salt amp。)。Gaussian Noise39。image=double(J)。for i=5:M4, for j=5:N4, u(1)=image(i,j)。 u(3)=image(i1,j+1)。 u(5)=image(i1,j1)。u(7)=image(i+1,j1)。 u(9)=image(i+1,j+1)。 u(11)=image(i,j2)。u(13)=image(i1,j+2)。 u(15)=image(i+1,j+2)。 u(q)=u(q+1)。end,end,end,image(i,j)=u(8)。 subplot(2,3,5), imshow(d), title(39。)。[M,N]=size(B)。 u(2)=image(i,j+1)。u(4)=image(i1,j)。 u(6)=image(i,j1)。 u(8)=image(i+1,j)。u(10)=image(i1,j2)。 u(12)=image(i+1,j2)。 u(14)=image(i,j+2)。for p=1:15, for q=1:15p, if u(q)u(q+1)k=u(q)。 u(q+1)=k。end,end,d=uint8(image)。3*5高斯噪聲中值濾波圖像39。圖35 長方形3*5中值濾波器對椒鹽噪聲和高斯噪聲的去噪 （2）長方形5*7中值濾波器去噪程序I=imread(39。)。B=rgb2gray(I)。Original Image39。J = imnoise(B,39。 pepper39。 %noise density=K= imnoise(B,39。,)。Salt amp。)。Gaussian Noise39。image=double(J)。for i=5:M4, for j=5:N4, u(1)=image(i,j)。 u(3)=image(i1,j+1)。 u(5)=image(i1,j1)。u(7)=image(i+1,j1)。 u(9)=image(i+1,j+1)。 u(11)=image(i,j2)。u(13)=image(i1,j+2)。 u(15)=image(i+1,j+2)。 u(17)=image(i2,j1)。u(19)=image(i2,j+1)。 u(21)=image(i2,j+3)。 u(23)=image(i1,j+3)。u(25)=image(i,j+3)。 u(27)=image(i+1,j+3)。 u(29)=image(i+2,j2)；u(30)=image(i+2,j1)。 u(32)=image(i+2,j+1)。u(34)=image(i+2,j+3)。for p=1:35, for q=1:35p, if u(q)u(q+1)k=u(q)。 u(q+1)=k。end,end,d=uint8(image)。5*7椒鹽噪聲中值濾波圖像39。 image=double(K)。for i=5:M4, for j=5:N4, u(1)=image(i,j)。 u(3)=image(i1,j+1)。 u(5)=image(i1,j1)。u(7)=image(i+1,j1)。 u(9)=image(i+1,j+1)。 u(11)=image(i,j2)。u(13)=image(i1,j+2)。 u(15)=image(i+1,j+2)。 u(17)=image(i2,j1)。u(19)=image(i2,j+1)。 u(21)=image(i2,j+3)。 u(23)=image(i1,j+3)。u(25)=image(i,j+3)。 u(27)=image(i+1,j+3)。 u(29)=image(i+2,j2)；u(30)=image(i+2,j1)。 u(32)=image(i+2,j+1)。u(34)=image(i+2,j+3)。for p=1:35, for q=1:35p, if u(q)u(q+1)k=u(q)。 u(q+1)=k。end,end,d=uint8(image)。5*7高斯噪聲中值濾波圖像39。圖36 長方形5*7中值濾波器對椒鹽噪聲和高斯噪聲去噪 （3）長方形3*9中值濾波器去噪程序I=imread(39。)。B=rgb2gray(I)。Original Image39。J = imnoise(B,39。 pepper39。 %noise density=K= imnoise(B,39。,)。Salt amp。)。Gaussian Noise39。image=double(J)。for i=5:M4, for j=5:N4, u(1)=image(i,j)。 u(3)=image(i1,j+1)。 u(5)=image(i1,j1)。u(7)=image(i+1,j1)。 u(9)=image(i+1,j+1)。 u(11)=image(i,j2)。u(13)=image(i1,j+2)。 u(15)=image(i+1,j+2)。 u(17)=image(i,j+3)。u(19)=image(i1,j+4)。 u(21)=image(i+1,j+4)。 u(23)=image(i,j3)。u(25)=image(i1,j2)。 u(27)=image(i+1,j2)。u(q)=u(q+1)。end,end,end,image(i,j)=u(14)。 subplot(2,3,5),imshow(d),title(39。)。[M,N]=size(B)。u(2)=image(i,j+1)。u(4)=image(i1,j)。u(6)=image(i,j1)。u(8)=image(i+1,j)。u(10)=image(i1,j2)。u(12)=image(i+1,j2)。u(14)=image(i,j+2)。u(16)=image(i1,j+3)。u(18)=image(i+1,j+3)。u(20)=image(i,j+4)。u(22)=image(i1,j3)。u(24)=image(i+1,j3)。u(26)=image(i,j2)。for p=1:27, for q=1:27p, if u(q)u(q+1)k=u(q)。u(q+1)=k。end,end,d=uint8(image)。3*9高斯噪聲中值濾波圖像39。圖37 長方形3*9中值濾波器對椒鹽噪聲和高斯噪聲去噪（4）長方形7*9中值濾波器去噪程序I=imread(39。)。B=rgb2gray(I)。Original Image39。J = imnoise(B,39。 pepper39。 %noise density=K= im