1、数字图像处理实验数字图像处理实验部门： xxx 时间： xxx整理范文，仅供参考，可下载自行编辑实验四。I=rgb2gray(I1。subplot(2,3,1。imshow(I。BW1=edge(I,roberts。subplot(2,3,2。imshow(BW1。BW2=edge(I,prewitt。subplot(2,3,3。imshow(BW2。BW3=edge(I,sobel。subplot(2,3,4。imshow(BW3。BW4=edge(I,log。subplot(2,3,5。imshow(BW4。BW5=edge(I,canny。subplot(2,3,6。imshow(BW5
2、。边缘检测图像：实验四。subplot(131。imshow(I1。subplot(132。imhist(I1。T=90。S=size(I1。maxI,maxP=max(I1(:,:。minI,minP=max(I1(:,:。for i=1:S(1for j=1:S(2if I1(i,j=T I1(i,j=255。elseI1(i,j=0。endendendsubplot(133。imshow(I1。图像分割图像：实验四。colormap(map。title(原图像cA1,cH1,cV1,cD1=dwt2(X,bior3.7。A1= upcoef2(a,cA1,bior3.7,1。H1= up
3、coef2(h,cH1,bior3.7,1。V1= upcoef2(v,cV1,bior3.7,1。D1= upcoef2(d,cD1,bior3.7,1。figure。colormap(map。subplot(2,2,1。image(wcodemat(A1,180。title(Approximation A1subplot(2,2,2。image(wcodemat(H1,255。title(Horizontal Detail H1subplot(2,2,3。image(wcodemat(V1,255。title(Vertical Detail V1subplot(2,2,4。image(wc
4、odemat(D1,255。title(Diagonal Detail D1Y=2.0*IDWT2(A1,H1,V1,D1, bior3.7。Y=imresize(Y,0.5。figure。image(Y。colormap(map。title(逆变后的图像实验一图像：实验二。subplot(2,3,1。imshow(I,。title(a 原始图像f = imnoise(I,salt & pepper,0.04。subplot(2,3,2。imshow(f。title(b 有噪声图像h0 = 1/9.*1 1 1 1 1 1 1 1 1。h1 = 0.1 0.1 0.1。0.1 0.2 0.1。
5、0.1 0.1 0.1。h2 = 1/16.*1 2 1。2 4 2。1 2 1。h3 = 1/8.*1 1 1。1 0 1。1 1 1。g0 = filter2(h0,f。g1 = filter2(h1,f。g2 = filter2(h2,f。g3 = filter2(h3,f。subplot(2,3,3。imshow(g0,。title(c 用模块0处理后的图像subplot(2,3,4。imshow(g1,。title(d 用模块1处理后的图像subplot(2,3,5。imshow(g2,。title(e 用模块2处理后的图像subplot(2,3,6。imshow(g3,。title
6、(f 用模块3处理后的图像平滑处理图像：实验二。subplot(2,2,1。imshow(I。title(a 原图subplot(2,2,2。imhist(I。title(b 原图的直方图J=imadjust(I,0.3 0.7,。subplot(2,2,3。imshow(J。title(c 输出图像subplot(2,2,4。imhist(J。title(d 输出图像的直方图灰度修正图像：实验二。I=rgb2gray(I1。subplot(3,2,1。imshow(I。title(a原始图像h=0 -1 0。-1 4 -1。0 -1 0。J=imfilter(I,h。subplot(3,2,
7、2。imshow(J。title(b拉普拉斯锐化subplot(3,2,3。imhist(J。title(c拉普拉斯锐化后图像的直方图K=imadjust(J,0.0 0.2,。subplot(3,2,4。imhist(J。title(d对锐化后图像的对比度扩展subplot(3,2,5。imshow(J。title(e对比度扩展后的图像拉普拉斯算子锐化图像：实验三图像复原有关程序：clc。clear。I = imread(flowers.jpg。I = I(10+1:1024,222+1:1024,:。subplot(2,4,1。imshow(I。title(原图像。LEN = 31。THE
8、TA= 11。PSF = fspecial(motion,LEN,THETA。blurred = imfilter(I,PSF,circular,conv。subplot(2,4,2。imshow(blurred。title(模糊图像。%len=28。theta=14。PSF=fspecial(motion,len,theta。wnr1=deconvwnr(blurred,PSF。subplot(2,4,3。imshow(wnr1。title(使用真实PSF维纳滤波复原图像。wnr2=deconvwnr(blurred,fspecial(motion,len*2,theta。b5E2RGbCA
9、Psubplot(2,4,4。imshow(wnr2。title(长PSF。wnr3=deconvwnr(blurred,fspecial(motion,len,2*theta。p1EanqFDPwsubplot(2,4,5。imshow(wnr3。title(陡峭PSF。% noise=0.1*randn(size(I。%模糊化BlurredNoisy=imadd(blurred,im2uint8(noise。wnr4=deconvwnr(BlurredNoisy,PSF。subplot(2,4,6。imshow(wnr4。title(维纳滤波复原含噪声图像。NSR=sum(noise(:.
10、2/sum(im2double(I(:.2。wnr5=deconvwnr(BlurredNoisy,PSF,NSR。subplot(2,4,7。imshow(wnr5。title(设置噪声比参数后。NP=abs(fftn(noise.2。NCORR=fftshift(real(ifftn(NP。IP=abs(fftn(im2double(I.2。ICORR=fftshift(real(ifftn(IP。wnr6=deconvwnr(BlurredNoisy,PSF,NCORR,ICORR。subplot(2,4,8。imshow(wnr6。title(设置噪声和图像自相关函数后。图像复原有关图像：申明：所有资料为本人收集整理，仅限个人学习使用，勿做商业用途。
