pkg load image

clear
clc

I1 = imread('lena.bmp');
I2 = imread('peppers.bmp');

I1_8 = bitget(I1,8);
I1_7 = bitget(I1,7);
I1_6 = bitget(I1,6);
I1_5 = bitget(I1,5);
I1_4 = bitget(I1,4);
I1_3 = bitget(I1,3);
I1_2 = bitget(I1,2);
I1_1 = bitget(I1,1);

I2_8 = bitget(I2,8);
I2_7 = bitget(I2,7);
I2_6 = bitget(I2,6);
I2_5 = bitget(I2,5);
I2_4 = bitget(I2,4);
I2_3 = bitget(I2,3);
I2_2 = bitget(I2,2);
I2_1 = bitget(I2,1);

I3 = 128*I1_8 + 64*I2_7 + 32*I1_6 + 16*I2_5 + 8*I1_4 + 4*I2_3 + 2*I1_2 + 1*I2_1;
I4 = 128*I2_8 + 64*I1_7 + 32*I2_6 + 16*I1_5 + 8*I2_4 + 4*I1_3 + 2*I2_2 + 1*I1_1;

subplot(2,2,1) ,imshow(I1,[0 255]), title('Lena')
subplot(2,2,2) ,imshow(I2,[0 255]), title('Peppers')
subplot(2,2,3) ,imshow(I3,[0 255]), title('Combination')
subplot(2,2,4) ,imshow(I4,[0 255]), title('Combination 2')