Retauracion de una imagen en matlab

Mejía Romero, Sergio.
Instituto Nacional de Astrofísica, Óptica y Electrónica, (INAOE)
Apartado Postal 23, Tonantzintla, Puebla.
e-mail: smejia@inaoep.mx,
Entregado el 3 de Julio 2010%--------------------------------------------------------------------------
% Examen de Transformada Rapida de Fourier
% Dr. Luis Raul Berriel Valdos
%Alumno: Sergio Mejía Romero
%--------------------------------------------------------------------------
clc;clear all;
%--------------------------------------------------------------------------
%1). Obtener 20 fotografías (color) de un mismo objeto-inmóvil y con la
% cámara fija.
% 2). Pasarlas a tonos de grises una vez que reduzcan el formato de la
% fotografía a 256X256 píxeles.
% 3).Promediar las 20 fotografías en su nuevo formato 256X256 píxeles.
% 4). Quitarle el ruido a una sola de estas fotografías y comparar ésta con el
% promedio de las 20 fotografías.
% 5). Desenfocarnuméricamente la fotografía sin ruido o la fotografía
% promedio.
% 6). Restaurar la fotografía del punto 5 (de esta forma se debe obtener una
% fotografía muy parecida a la fotografía promedio).
%%--------------------------------------------------------------------------
%1) y 2) Leyendo las 19 imagenes.
% Se cambian a escala de gris y se recorta la imagen a un formato de
% 256 X 256pixels.
% De igual forma se cambia el formato de la matrix-imagen de uint8 que solo
% tiene valores por pixel de 0 a 255 a un formato double que puede tener
% valores arbitrarios por pixel , parapoder realizar la suma de pixel por
% pixel y el promedio de ella.
%--------------------------------------------------------------------------
I1=imread('DSCN3148.jpg');I1= rgb2gray(I1);
I1=imcrop(I1,[930,690,255,255]);I1 = double(I1);
I2=imread('DSCN3149.jpg');I2= rgb2gray(I2);
I2= imcrop(I2,[930,690,255,255]);I2 = double(I2);
I3=imread('DSCN3150.jpg');I3= rgb2gray(I3);
I3=...