《2022年数字图像程序文件》由会员分享，可在线阅读，更多相关《2022年数字图像程序文件（5页珍藏版）》请在金锄头文库上搜索。
1、1%计算图象经过膨胀运算后图象面积的变化% BW1=imread(Boat.bmp); subplot(121),imshow(BW1) title( 原始图象 ); SE=ones(5); BW2=imdilate(BW1,SE); subplot(122),imshow(BW2) title( 膨胀后的图象); zengjia=(bwarea(BW2)-bwarea(BW1)/bwarea(BW1) 实验结果： zengjia = 3.3379e-006 2%计算矩阵的面积和质心坐标% A=0 1 1 1 0 1 1 0;. 0 1 1 1 1 1 1 0;. 0 1 1 1 0 0 0
2、0;. 1 1 1 1 1 1 0 0;. 1 1 1 1 1 1 1 1;. 1 1 1 1 1 0 0 0;. 0 1 1 1 1 1 1 0;. 0 0 0 1 1 1 1 1 regionprops(A,area) regionprops(A,centroid) 实验结果：A = 0 1 1 1 0 1 1 0 0 1 1 1 1 1 1 0 0 1 1 1 0 0 0 0 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - -
3、 - 第 1 页，共 5 页 - - - - - - - - - 1%计算图象经过膨胀运算后图象面积的变化% BW1=imread(Boat.bmp); subplot(121),imshow(BW1) title( 原始图象 ); SE=ones(5); BW2=imdilate(BW1,SE); subplot(122),imshow(BW2) title( 膨胀后的图象); zengjia=(bwarea(BW2)-bwarea(BW1)/bwarea(BW1) 实验结果： zengjia = 3.3379e-006 2%计算矩阵的面积和质心坐标% A=0 1 1 1 0 1 1 0;.
4、 0 1 1 1 1 1 1 0;. 0 1 1 1 0 0 0 0;. 1 1 1 1 1 1 0 0;. 1 1 1 1 1 1 1 1;. 1 1 1 1 1 0 0 0;. 0 1 1 1 1 1 1 0;. 0 0 0 1 1 1 1 1 regionprops(A,area) regionprops(A,centroid) 实验结果：A = 0 1 1 1 0 1 1 0 0 1 1 1 1 1 1 0 0 1 1 1 0 0 0 0 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - -
5、 - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 5 页 - - - - - - - - - 1%计算图象经过膨胀运算后图象面积的变化% BW1=imread(Boat.bmp); subplot(121),imshow(BW1) title( 原始图象 ); SE=ones(5); BW2=imdilate(BW1,SE); subplot(122),imshow(BW2) title( 膨胀后的图象); zengjia=(bwarea(BW2)-bwarea(BW1)/bwarea(BW1) 实验结果： zengjia = 3.3379e-006 2%
6、计算矩阵的面积和质心坐标% A=0 1 1 1 0 1 1 0;. 0 1 1 1 1 1 1 0;. 0 1 1 1 0 0 0 0;. 1 1 1 1 1 1 0 0;. 1 1 1 1 1 1 1 1;. 1 1 1 1 1 0 0 0;. 0 1 1 1 1 1 1 0;. 0 0 0 1 1 1 1 1 regionprops(A,area) regionprops(A,centroid) 实验结果：A = 0 1 1 1 0 1 1 0 0 1 1 1 1 1 1 0 0 1 1 1 0 0 0 0 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 名师资料总结 -
7、- -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 5 页 - - - - - - - - - MFN=imadd(MF,im2uint8(noise); %生成运动模糊 +高斯噪声图象figure(2);imshow(MFN,); %显示运动模糊叠加高斯噪声的图象NSR=sum(noise(:).2)/sum(MFN(:).2); %计算信噪比figure(3); imshow(deconvwnr(MFN,PSF),); %显示逆滤波复原图象figure(4); imshow(deconvw
8、nr(MFN,PSF,NSR),); %显示维纳滤波复原图象NP=0.02*prod(size(I); % 噪声强度reg1,LAGRA=deconvreg(MFN,PSF,NP/3.0); % 有噪声强度约束最小二乘滤波复原figure(5);imshow(reg1); edged=edgetaper(MFN,PSF); % 对图象边缘信息提取reg2=deconvreg(edged,PSF,LAGRA/300); %有拉格郎日乘法算子最小乘滤波复原figure(6);imshow(reg2); 6.%利用 deconvwnr 函数对有噪声模糊图象进行复原重建% clear I=imread
9、(1.jpg);% 读入已经模糊的图像I=I(5+1:256,40+1:256,:); subplot(3,3,1),imshow(I); LEN=31;% 指定运动位移为31 像素，运动角度为11 度THETA=11; PSF=fspecial(motion,LEN,THETA);%产生运动模糊的PSF Blurred=imfilter(I,PSF ,circular,conv); subplot(3,3,2),imshow(Blurred); wnr1=deconvwnr(Blurred,PSF); subplot(3,3,3),imshow(wnr1);%逆滤波复原noise=0.1.*
10、rand(size(Blurred);%添加随机噪声J=im2uint8(noise); BlurredNoise=imadd(Blurred,J); subplot(3,3,4),imshow(BlurredNoise);%显示添加随机噪声的图像wnr2=deconvwnr(Blurred,PSF);% 调用 deconvwer 函数直接对一副有噪声的图象进行复原名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 5 页 - - - - - - - - - subplot(
11、3,3,5),imshow(wnr2); %显示对有噪声的图象复原NSR=sum(noise(:).2)/sum(im2double(I(:).2);以信噪比NSR 为噪声参数进行图象恢复wnr3=deconvwnr(BlurredNoise,PSF,NSR); subplot(3,3,6);imshow(wnr3); wnr31=deconvwnr(BlurredNoise,PSF,NSR/2); subplot(3,3,7);imshow(wnr3); NP=abs(fftn(noise).2; NCORR=fftshift(real(ifftn(NP); %计算噪声的自相关函数IP=abs(fftn(im2double(Blurred).2; ICORR=fftshift(real(ifftn(IP); wnr4=deconvwnr(BlurredNoise,PSF,NCORR,ICORR); %计算图象自相关函数subplot(3,3,8);imshow(wnr4); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 5 页 - - - - - - - - -
