%funkcia slúžiaca na výpočet T1 a M0 pomocou funkcie lsqcurvefit
function [aaa1,aaa2,aa]=T1mapy(Z,TI)

%nepotrebný cyklus,slúži na uloženie dát značených a neznačených do konkrétnej matice
for i=1:15
    mjp(:,:,i,1)=cell2mat(Z(1,i));
    mjp(:,:,i,2)=cell2mat(Z(2,i));
end

%cyklus slúži na výber aktuálnej mapy, ktorá bude vstupovať do lsqcurvefit
for k=1:2
       akt_mapa=cell2mat(Z(k,1)); %mapa saprevedie z cell na matrix
       [i_lim j_lim]=size(akt_mapa);%určí sa veľkosť mapy
       %stanovia sa parametre vstupujúce do lsqcurvefit
       TIa=TI(k,:);
       x0=[1000 1000]; %predstavuje T1 a M0
       x=x0;
       lb=[1e2 1e2];
       ub=[1e10 1e10];
       for i=1:i_lim
           for j=1:j_lim
               for u=1:15
            akt_pix(u)=mjp(i,j,u,k);
               end
    [aa]=lsqcurvefit(@rovnica,x0,TI(k,:),akt_pix,lb,ub);
            aaa1(i,j,k)=aa(1);   
            aaa2(i,j,k)=aa(2);
           end
       end
    
end
end
%end