pkg load image

clear
clc

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

I3 = uint8(zeros(256));
I4 = uint8(zeros(256));

for i = 1 : 256
  for j = 1 : 256
    if (I1(i,j) <= I2(i,j))
      I3(i,j) = I1(i,j);
    else 
      I3(i,j) = I2(i,j);
    endif
    if (I1(i,j) >= I2(i,j))
      I4(i,j) = I1(i,j);
    else
      I4(i,j) = I2(i,j);
    endif
    
  endfor
endfor

I1min = min(min(I1));
I1max = max(max(I1));
I1avg = mean(mean(I1));
I1median = median(median(I1));

I2min = min(min(I2));
I2max = max(max(I2));
I2avg = mean(mean(I2));
I2median = median(median(I2));

I3min = min(min(I3));
I3max = max(max(I3));
I3avg = mean(mean(I3));
I3median = median(median(I3));

I4min = min(min(I4));
I4max = max(max(I4));
I4avg = mean(mean(I4));
I4median = median(median(I4));

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