%% Imagem
close all
figure

imagem = imread('Foto.jpg');
subplot(3,1,1)
imshow(imagem)
title('Imagem original')

imagem_exp = uint8(zeros(2*size(imagem,1),2*size(imagem,2),3));
imagem_exp(1:2:end,1:2:end,:) = imagem;

subplot(3,2,3)
imshow(imagem_exp)
title('Imagem expandida')

imagem_interp = imagem_exp;

for k = 1:3
    for i = 2:2:size(imagem_exp,1)-1
        for j = 2:2:size(imagem_exp,2)-1
            
            pixels = [imagem_exp(i-1,j-1,k), imagem_exp(i-1,j+1,k);
                      imagem_exp(i-1,j-1,k), imagem_exp(i+1,j+1,k)];
            
            pixel = (1/4)*sum(sum(pixels));
        
            imagem_interp(i,j,k) = uint8(pixel);
            
        end
    end
end

subplot(3,2,4)
imshow(imagem_interp)
title('Interpolação dos pixels de linhas e colunas pares')

for k = 1:3
    for i = 3:2:size(imagem_exp,1)-1
        for j = 2:2:size(imagem_exp,2)-1
            
            pixels = [imagem_interp(i-1,j,k), imagem_interp(i+1,j,k);
             imagem_interp(i,j-1,k), imagem_interp(i,j+1,k)];
            
            pixel = (1/4)*sum(sum(pixels));
        
            imagem_interp(i,j,k) = uint8(pixel);
            
        end
    end
end

subplot(3,2,5)
imshow(imagem_interp)
title('Interpolação dos pixels de linhas ímpares e colunas pares')

for k = 1:3
    for i = 2:2:size(imagem_exp,1)-1
        for j = 3:2:size(imagem_exp,2)-1
            
            pixels = [imagem_interp(i-1,j,k), imagem_interp(i+1,j,k);
             imagem_interp(i,j-1,k), imagem_interp(i,j+1,k)];
            
            pixel = (1/4)*sum(sum(pixels));
        
            imagem_interp(i,j,k) = uint8(pixel);
            
        end
    end
end

subplot(3,2,6)
imshow(imagem_interp)
title('Interpolação dos pixels de linhas pares e colunas ímpares')