2021-05-16

PART2

1.

（1）.图片加法

a=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\加\Noisy_image1.jpg')
b=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\加\Noisy_image2.jpg')
c=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\加\Noisy_image3.jpg')
d=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\加\Noisy_image4.jpg')
e=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\加\Noisy_image5.jpg')
f=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\加\Noisy_image6.jpg')
g=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\加\Noisy_image7.jpg')
h=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\加\Noisy_image8.jpg')
i=imadd(a,b)
j=imadd(c,d)
k=imadd(e,f)
l=imadd(g,h)
m=imadd(i,j)
n=imadd(k,l)
o=imadd(m,n)
imshow(o)

（2）.图片减法

A=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\减\Fig1.tif')
B=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\减\Fig2.tif')
C=imsubtract(A,B)figure;%显示原图像、匹配图像和匹配后的图像
subplot(1,3,1),imshow(A);title('原图像');
subplot(1,3,2),imshow(B);title('匹配图像');
subplot(1,3,3),imshow(C);title('匹配之后图像');

 

（3）.图片乘法

I=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\乘\Fig1.tif')
J=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\乘\Fig2.tif')
I2=imdivide(I,J)

figure;%显示原图像、匹配图像和匹配后的图像
subplot(1,3,1),imshow(I);title('原图像');
subplot(1,3,2),imshow(J);title('匹配图像');
subplot(1,3,3),imshow(I2);title('匹配之后图像');

（4）.图片除法

I=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\除\Fig1.tif')
J=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\加减乘除\除\Fig2.tif')
I1=immultiply(I,J) 

figure;%显示原图像、匹配图像和匹配后的图像
subplot(1,3,1),imshow(I);title('原图像');
subplot(1,3,2),imshow(J);title('匹配图像');
subplot(1,3,3),imshow(I1);title('匹配之后图像');

2

（1）.输出图片RGB直方图

img=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\直方图与直方图匹配\lena.jpg')
R=img(:,:,1)
G=img(:,:,2)
B=img(:,:,3)

figure(1)
imhist(R);
figure(2)
imhist(G);
figure(3)
imhist(B);

（2）.图片匹配

I=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\直方图与直方图匹配\lena.jpg');%读取图像
Imatch=imread('D:\study\WeChat Files\wxid_sovkik65v05322\FileStorage\File\2021-05\任务1\Part2\直方图与直方图匹配\im1.jpg');%读取匹配图像
R=I(:,:,1);%获取原图像R通道
G=I(:,:,2);%获取原图像G通道
B=I(:,:,3);%获取原图像B通道
Rmatch=Imatch(:,:,1);%获取匹配图像R通道
Gmatch=Imatch(:,:,2);%获取匹配图像G通道
Bmatch=Imatch(:,:,3);%获取匹配图像B通道
Rmatch_hist=imhist(Rmatch);%获取匹配图像R通道直方图
Gmatch_hist=imhist(Gmatch);%获取匹配图像G通道直方图
Bmatch_hist=imhist(Bmatch);%获取匹配图像B通道直方图
Rout=histeq(R,Rmatch_hist);%R通道直方图匹配
Gout=histeq(G,Gmatch_hist);%G通道直方图匹配
Bout=histeq(B,Bmatch_hist);%B通道直方图匹配
J(:,:,1)=Rout;
J(:,:,2)=Gout;
J(:,:,3)=Bout;
figure;%显示原图像、匹配图像和匹配后的图像
subplot(1,3,1),imshow(I);title('原图像');
subplot(1,3,2),imshow(Imatch);title('匹配图像');
subplot(1,3,3),imshow(J);title('匹配之后图像');