加减乘除
加
I1=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\加\Noisy_image1.jpg’);
I2=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\加\Noisy_image2.jpg’);
I3=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\加\Noisy_image3.jpg’);
I4=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\加\Noisy_image4.jpg’);
I5=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\加\Noisy_image5.jpg’);
I6=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\加\Noisy_image6.jpg’);
I7=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\加\Noisy_image7.jpg’);
I8=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\加\Noisy_image8.jpg’);
A1=im2double(I1);
B1=im2double(I2);
C1=im2double(I3);
D1=im2double(I4);
E1=im2double(I5);
F1=im2double(I6);
G1=im2double(I7);
H1=im2double(I8);

A=imadd(A1,B1);
B=imadd(C1,D1);
C=imadd(E1,F1);
D=imadd(G1,H1);
F=imadd(A,B);
E=imadd(C,D);
G=imadd(E,F);

imshow(G,[]);

减
I1=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\减\Fig1.tif’);
I2=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\减\Fig2.tif’);
A=imsubtract(I1,I2);
B=im2double(I1);
C=im2double(I2);
D=imsubtract(B,C);
E=imsubtract(C,B);
subplot(121),imshow(D,[]);
subplot(122),imshow(E,[]);
乘
I1=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\乘\Fig1.tif’);
I2=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\乘\Fig2.tif’);
A=immultiply(I2,I1);
B=im2double(I2);
C=im2double(I1);
D=immultiply(B,C);

imshow(D,[]);
除
I1=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\除\Fig1.tif’);
I2=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\加减乘除\除\Fig2.tif’);
A=im2double(I1);
B=im2double(I2);
C=imdivide(A,B);
imshow(C,[]);
显示直方图及直方图匹配
I=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\直方图与直方图匹配\lena.jpg’);
Imatch=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\直方图与直方图匹配\im1.jpg’);
Imatch2=imread(‘C:\Users\1\Desktop\要做的作业\任务1\Part2\直方图与直方图匹配\im2.jpg’);
R=I(:,:,1);
G=I(:,:,2);
B=I(:,:,3);%获取原图像rgb通道
Rmatch=Imatch(:,:,1);
Gmatch=Imatch(:,:,2);
Bmatch=Imatch(:,:,3);
Rmatch2=Imatch2(:,:,1);
Gmatch2=Imatch2(:,:,2);
Bmatch2=Imatch2(:,:,3);%获取匹配图像的rgb通道

Rmatchhist=imhist(Rmatch);
Gmatchhist=imhist(Gmatch);
Bmatchhist=imhist(Bmatch);
Rmatchhist2=imhist(Rmatch2);
Gmatchhist2=imhist(Gmatch2);
Bmatchhist2=imhist(Bmatch2);%获取匹配图像rgb通道直方图

Rout=histeq(R,Rmatchhist);
Gout=histeq(G,Gmatchhist);
Bout=histeq(B,Bmatchhist);
Rout2=histeq(R,Rmatchhist2);
Gout2=histeq(G,Gmatchhist2);
Bout2=histeq(B,Bmatchhist2);%rgb通道直方图匹配

J(:,:,1)=Rout;
J(:,:,2)=Gout;
J(:,:,3)=Bout;
L(:,:,1)=Rout2;
L(:,:,2)=Gout2;
L(:,:,3)=Bout2;
figure;
subplot(131),imshow(J),title(‘匹配im2’);
subplot(132),imshow(L),title(‘匹配im3’);
figure;
subplot(231),imhist(R,256);title(‘R通道直方图’);
subplot(232),imhist(G,256);title(‘G通道直方图’);
subplot(233),imhist(B,256);title(‘B通道直方图’);