西安电子科技大学。
数字图像处理。
上机大作业。
教师 : 教授。
学院 :通信工程学院。
班级 : 专业 :电子与通信工程。
姓名 : 学号 :
f1与f2的幅度谱相同,因为f1与f2振幅的绝对值是相同的。
f3的幅度谱与f2的幅度谱按顺时针旋转90度后相同。
f4的幅度谱与f1的幅度谱按顺时针旋转90度后相同。
f5的幅度谱是f1与f4的幅度谱相加,f6的幅度谱是f2与f3的幅度谱相加。
程序:f1=zeros(256,256);
f1(65:192,113:144)=100;
figure
subplot(121)
imshow(f1)
title('原图f1')
f=fft2(f1);
subplot(122)
imshow(abs(f))
title('f1幅度谱图')
for m=(1:256);
for n=(1:256);
f2(m,n)=power(-1,m+n)*f1(m,n);
endend
figure
subplot(121)
imshow(f2)
title('f2图')
f1=fft2(f2);
subplot(122)
imshow(abs(f1))
title('f2幅度谱图')
f3=imrotate(f2,-90);
figure
subplot(121)
imshow(f3)
title('f3图')
f2=fft2(f3);
subplot(122)
imshow(abs(f2))
title('f3幅度谱图')
f4=imrotate(f1,-90);
figure
subplot(121)
imshow(f4)
title('f4图')
f5=fft2(f4);
subplot(122)
imshow(abs(f5))
title('f4幅度谱图')
f5=f1+f4;
figure
subplot(121)
imshow(f5)
title('f5图')
f3=fft2(f5)
subplot(122)
imshow(abs(f3))
title('f5幅度谱图')
f6=f2+f3;
figure
subplot(121)
imshow(f6)
title('f6图')
f4=fft2(f6)
subplot(122)
imshow(abs(f4))
title('f6幅度谱图')
256*256图像如下:
从中值滤波和均值滤波后的图像可以看出,中值滤波把黑色方快完全完全隔离开,而白色方块则通过边角的两个多出来的像素点相连,而均值滤波则把图像变模糊了,边缘模糊的比较厉害。中值滤波后图像的直方图没有变化,而均值滤波后的直方图变化了。
程序:a=zeros(256,256);
for i=(1:64:256);
for j=(1:64:256);
a(j:j+32,i:i+32)=1;
endend
for m=(32:64:256);
for n=(32:64:256);
a(n:n+32,m:m+32)=1;
endend
figure
imshow(a)
title('原图像')
figure
k=medfilt2(a,[3,3])
subplot(121)
imshow(k)
title('3*3中值滤波')
m=filter2(fspecial('**erage',3),a)
subplot(122)
imshow(m)
title('3*3均值值滤波')
figure
subplot(131)
imhist(a)
axis([0 2 0 2])
title('原图直方图')
subplot(132)
imhist(k)
axis([0 2 0 2])
title('中值滤波后的直方图')
subplot(133)
imhist(m)
axis([0 2 0 2])
title('均值滤波后的直方图')
整体而言中值滤波的结果比均值滤波的结果要好,均值滤波模糊了轮廓边缘,但是中值滤波对椒盐噪声的滤波破坏了轮廓边缘。中值滤波对高斯噪声的水平滤波不是很好。
程序:i=zeros(256,256);
for i=(32:24:224);
i(23:233,i:i+7)=1;
endfigure
imshow(i)
title('原图像')
a=imnoise(i,'gaussian',0,0.025);
figure
subplot(121)
imshow(a)
title('高斯白噪声')
b=imnoise(i,'salt & pepper',0.025);
subplot(122)
imshow(b)
title('椒盐噪声')
h=[1 1 1;1 1 1;1 1 1];%领域平均法消除噪声模板为1/9[1 1 1;1 1 1;1 1 1]
h=9./h;
ja=conv2(a,h);
figure
subplot(121)
imshow(ja,
title('对高斯白噪声的均值滤波')
jb=conv2(b,h);
subplot(122)
imshow(jb,
title('对椒盐噪声的均值滤波')
ka=medfilt2(a,[3,3])
figure
subplot(121)
imshow(ka,
title('对高斯白噪声的中值滤波')
kb=medfilt2(b,[3,3]);
subplot(122)
imshow(kb,
title('对椒盐噪声的中值波')
roberts自动选择的阈值为:
sobel自动选择的阈值为:
prewitt自动选择的阈值为:
程序:i=imread(''
figure
imshow(i)
title('原始图像')
bw1,thresh1]=edge(i,'roberts');
disp('roberts自动选择的阈值为:')
disp(thresh1)
figure
subplot(131)
imshow(bw1)
title('roberts算子')
bw2,thresh2]=edge(i,'sobel');
disp('sobel自动选择的阈值为:')
disp(thresh2)
subplot(132)
imshow(bw2)
title('sobel算子')
bw3,thresh3]=edge(i,'prewitt');
disp('prewitt自动选择的阈值为:')
disp(thresh3)
subplot(133)
imshow(bw3)
title('prwitte算子')
bw4=filter2(fspecial('prewitt'),i);
bw5=filter2(fspecial('prewitt'),i);
bw6=filter2(fspecial('sobel'),i);
figure
subplot(2,2,1);
imshow(i);
title('灰度图像');
subplot(2,2,2);
imshow(bw4);
title('roberts的锐化结果');
subplot(2,2,3);
imshow(bw5);
title('prewitt的锐化结果');
subplot(2,2,4);
imshow(bw6);
title('sobel的锐化结果');
迭代后的阈值:131
程序:clear all;
i=imread(''
zmax=max(max(i));
zmin=min(min(i));
tk=(zmax+zmin)/2;
bcal=1;
isize=size(i);
while(bcal)
ifground=0;
ibground=0;
fgrounds=0;
bgrounds=0;
for i=1:isize(1)
for j=1:isize(2)
tmp=i(i,j);
if(tmp>=tk)
ifground=ifground+1;
fgrounds= fgrounds+double(tmp);
elseibground=ibground+1;
bgrounds= bgrounds+double(tmp);
endend
endzo=fgrounds/ifground;
zb=bgrounds/ibground;
tktmp=uint8((zo+zb)/2);
if(tktmp==tk)
bcal=0;
elsetk=tktmp;
endend
disp(strcat('迭代后的阈值:',num2str(tk)))
newi=im2bw(i,double(tk)/255);
subplot(121)
imshow(i)
title('原始图像')
subplot(122)
imshow(newi)
title('迭代法分割得到的二值图像')
数字图像处理大作业
1图像变换。1.1实验背景。在数字图像应用领域,图像需要进行分析 变换 压缩或者增强来提高图像的可处理性和视觉效果。其中,图像变换是将图像从空间域变换到频率域,变换的目的是根据图像在变换域的某些性质对其进行处理,而这些性质在空间域难以获取,通常在频率域才能获取,在变换域处理完后再反变换到空间域,恢复...
数字图像处理大作业
图像处理技术。大作业。院 系 物联网工程学院。专业 计算机科学与技术。班级 学号 姓名 摘要。图像分割就是指把图像分成各具特性的区域并提取出感兴趣目标的技术和过程。它是图像处理 模式识别和人工智能等多个领域中的重要课题,也是计算机视觉技术中首要的 重要的关键步骤。图像分割的应用非常广泛,几乎出现在有...
数字图像处理大作业
图像处理 特点 输入是图像,输出也是图像,即图像之间进行的变换。图像分割。特点 输入是图像,输出是数据。图像识别。特点 以客观世界为中心,借助知识 经验等来把握整个客观世界。输入是数据,输出是理解。rgb 红 绿 蓝 模型。cmy 青 品红 黄 模型。hsi 色调 饱和度 亮度 模型。1 采样。采样...