1、根据实验程序一2幅灰度图像相加的程序,用matlab实现2幅3通道rgb图像的相加。
a=imread(''
b=imread(''
a1=a(:,1);
a2=a(:,2);
a3=a(:,3);
b1=b(:,1);
b2=b(:,2);
b3=b(:,3);
c(:,1)=imadd(a1,b1,'uint16');
c(:,2)=imadd(a2,b2,'uint16');
c(:,3)=imadd(a3,b3,'uint16');
subplot(2,2,1);
imshow(a,
subplot(2,2,2);
imshow(b,
subplot(2,2,3);
imshow(uint8(c));
2、针对程序十一,实现对彩色图像的均值滤波处理。
i=imread(''读入源图像。
a=i(:,1)
b=i(:,2)
c=i(:,3)
figure;
subplot(3,2,1);
imshow(i);
title('原图');
aa=im2double(a);%将原图转换为double类型,值在0 1之间,imnoise要求。
bb=im2double(b)
cc=im2double(c)
j=imnoise(aa,'gaussian',0,0.005);%人为给图像添加高斯噪声,噪声均值为0,方差为0.005
k=imnoise(bb,'gaussian',0,0.005)
l=imnoise(cc,'gaussian',0,0.005)
e(:,1)=j
e(:,2)=k
e(:,3)=l
subplot(3,2,2);
imshow(e);%显示噪声图像。
title('加入高斯噪声后的图像');
用filter2函数对加入高斯噪声后的图像进行均值滤波。
f1=filter2(fspecial('**erage',3),j);%对图像进行模板大小为3的均值滤波。
f2=filter2(fspecial('**erage',3),k)
f3=filter2(fspecial('**erage',3),l)
d(:,1)=f1
d(:,2)=f2
d(:,3)=f3
subplot(3,2,3);
imshow(d); 显示模板大小为3的均值滤波图像。
title('模板为3的均值滤波');
f4=filter2(fspecial('**erage',5),j);%对图像进行模板大小为3的均值滤波。
f5=filter2(fspecial('**erage',5),k)
f6=filter2(fspecial('**erage',5),l)
g(:,1)=f4
g(:,2)=f5
g(:,3)=f6
subplot(3,2,4);
imshow(g); 显示模板大小为3的均值滤波图像。
title('模板为5的均值滤波');
3、根据程序五的定位程序,换张**(text_重新定位。
w=imread('d:\ccx\photo\text_')
将字母a从图象中切割出来。
a= w(33:45,88:98);;
figure; imshow(w);
figure;imshow(a);
将字母a和图象进行快速傅立叶变换,然后计算字母a和图象的卷积。
c=real(ifft2(fft2(w).*fft2(rot90(a,2),220,220)))
figure,imshow(c,
thresh=60;
显示象素值超过阈值的象素。
figure,imshow(c>thresh)
4、图像空间域增强。空域增强包括空域变换增强(线性对比度增强,非线性对比度增强等),选择一个非线性对比度增强的方法实现,把增强公式体现出来。
a=imread(''
subplot(1,2,1);
imshow(a);
a=double(a);
i=log(1+a);
subplot(1,2,2);
imshow(mat2gray(i));
5、将一幅图像的dct系数大于20的系数进行缩小2倍,然后进行dct反变换得到处理后的图像。
a=imread( '
subplot(2,2,1);
imshow(a);
i=rgb2gray(a);
subplot(2,2,2);
imshow(i);
j=dct2(i);
subplot(2,2,3);
imshow(log(abs(j)),
i=abs(j)>20;
i=abs(i)*0.5;
if abs(j)>20 ;
j=i;k=idct2(j);
subplot(2,2,4);
imshow(k,[0,255]);
6、将图像进行fft变换,将高频系数放大倍数,然后反变换,得到图像锐化效果。
a=imread(''
b=rgb2gray(a);
g=fftshift(fft2(b));
n1 n2]=size(b);
result2=zeros(n1,n2);
d0=5;n=2;l=1.2;
n1=fix(n1/2);
n2=fix(n2/2);
for i=1:n1
for j=2:n2
d=sqrt((i-n1)^2+(j-n2)^2);
h=1/((1+d0/d)^(2*n));
if g(i,j)>d0
result2(i,j)=g(i,j)*l;
elseresult2(i,j)=g(i,j);
endend
endr2=ifftshift(result2);
r2=ifft2(r2);
figure,subplot(2,2,1),imshow(a);
subplot(2,2,2),imshow(b);
subplot(2,2,3),imshow(g);
subplot(2,2,4),imshow(uint8(real(r2)),
7、实现一种图像分割方法。
function thresh
a=rgb2gray(imread(''
subplot(221);imshow(a);
imgbw=im2bw(a);
subplot(222);imshow(imgbw);
imgbw=im2bw(a,0.25);
subplot(223);imshow(imgbw);
i=graythresh(a);imgbw=im2bw(a,i);
subplot(224);imshow(imgbw);
数字图像作业
数字图像处理大作业,请 要求每个同学完成1题,按照学号进行分配 完成第一题 完成第二题,以此类推。每道题要求有理论说明和实验验证,完成后撰写报告并打印,报告页数不少于5页,包括,理论 实验 和实验结果,报告完成后由班长收齐后交到计算机楼530,时间7月2日截止。1 什么是图像增强,常用的图像增强方法...
数字图像作业
第三次作业完成情况报告。课程 数字图像处理课程编号 21909601 摘要 这次的作业主要任务是学习利用拉布拉斯公式对图像进行处理,并学会调用matlab中的函数对图像进行处理,最后是绘制所需的函数图像。技术说明 task1.利用拉布拉斯公式求矩阵的梯度场 输出梯度场图像以及角度的矩阵。task2....
数字图像处理作业
2017暑假小学期。专业。班级。姓名。学号。摘要 数字图像处理的目的是改善 质量,利用matlab设计几个简单的图像处理程序。该程序具备图像处理的日常功能 图像代数与几何运算 直方图与图像明暗关系分析 直方图均衡化 滤波器的使用与图像分割技术。通过运行 解释结合图像效果,分别介绍了处理原理与过程。关...