clear all; close all; clc

%-----------------------------------------------------------
% SEGUNDA PARTE

% Lista de nombres de las imagenes del agujero de ozono de cada aņo
path='E:\Docencia\CTE 1\2017\Practicas 2017\PRACTICA 4\Ozono14\Aņos';
type='*.png'; PNG=dir(fullfile(path,type));
display('IMAGENES: ')
for j=1:length(PNG)
    Nombres{j}=[PNG(j).name];
    display([num2str(j) '. ' Nombres{j}])
end

% Vector tiempo (sin el aņo 1995)
Tiempo=1978:2012;
Tiempo(find(Tiempo==1995))=[];

% Creamos una estructura 'imagenes' que contenga las mismas imagenes que 'nombres' pero ordenadas por aņos
n1=1;n2=1;n3=1;n4=1;
for j=1:length(Nombres)
    figure(1)
    A=imread([Nombres{j}]);
    imagesc(A)
    axis equal
    caxis([0 255])
    
    if ((j >= 13)&&(j <= 27)) % n7t
        load n7t.mat
        colormap(n7t)
        [xN7T,yN7T]=ginput;
        AreaN7T(n1)=polyarea(xN7T,yN7T);
        n1=n1+1;
        
    elseif ((j >= 11)&&(j <= 12)) % mt3
        load mt3.mat
        colormap(mt3)
        [xMT3,yMT3]=ginput;
        AreaMT3(n2)=polyarea(xMT3,yMT3);
        n2=n2+1;
       
    elseif ((j >= 1)&&(j <= 10)) % ept
        load ept.mat
        colormap(ept)
        [xEPT,yEPT]=ginput;
        AreaEPT(n3)=polyarea(xEPT,yEPT);
        n3=n3+1;
        
    elseif ((j >= 30)&&(j <= 36)) % omi (descartamos imagenes de 2004 y 2005
        load omi.mat
        colormap(omi)
        [xOMI,yOMI]=ginput;
        AreaOMI(n4)=polyarea(xOMI,yOMI);
        n4=n4+1;   
    end
end



% Definimos el vector de Area que incluye a todo los vectores areas definidos en el loop
Areas=[AreaN7T AreaMT3 AreaEPT AreaOMI];

figure(2)
plot(Tiempo,Areas,'*-r')
title('Evolucion del tamaņo del agujero de Ozono en la Antartida con el paso de los aņos')
xlabel('Aņos')
ylabel('Area del agujero (en px^2)')