%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%           !!!!!!!!!!!!"dicomwork" erforderlich!!!!!!!!!!!!          %%%
%%%            !!!!!!!!!!!!"roiread" erforderlich!!!!!!!!!!!!           %%%
%%%             !!!!!!!!!!!"timeread" erforderlich!!!!!!!!!!            %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Signalstärke über Zeit auslesen 
% Signalstärke über Ort auslesen
%
% 
% 
% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Signalstärke auslesen
% posall = Koordinaten der ROI's (aus "roiread")
% alld = alle Dicombilder für bestimmte SD und SL (aus "dicomwork")
% allp = Namen aller Dicombilder in Folder ("dicomwork")
% infoall = Position aller Dicombilder für SD und SL in allp ("dicomwork")
% timeread = function "timeread"
% roir = ROI-Art z.B. Zeit oder Ort

infoalls = size(infoall);
infoalls = infoalls(2);

clear signd; 
clear X; 
clear Sg; 
clear signdall;
clear c1; 
clear c2; 
clear c3; 
clear c4; 
clear c5;
Sg = 0;
c5 = 1;
c2 = 1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% über Zeit auslesen
switch(roir)
    case 'z'
    
      for c1 = 1:1:infoalls                    % Bild auswählen
            infoall(c1);
            Pic = allp((infoall(c1)),:);
            infopic = dicominfo (Pic);
            testsd = infopic.SeriesDescription;
            testsl = infopic.SliceLocation;
            X(c1) = timeread(Pic);               % function "timeread" 
            Picsign = dicomread(Pic);

        for c2 = 1:1:roim                    % ROI auswählen
            for c3 = posall(c2,1):1:(posall(c2,1)+posall(c2,3))     % Pixel auswählen (Spalte)
                for c4 = posall(c2,2):1:(posall(c2,2)+posall(c2,4)) % Pixel auswählen (Reihe)
                    S = Picsign(c4, c3);                      % Pixel einlesen
                    Sg = S+Sg;                                % Pixel aufsummieren
                end
            end
            signd(c2) = Sg/(posall(c2,3) * posall(c2,4)); % Signalldurschnitt in ROI berechnen
       end
       Sg = 0;
       signdall(c1,:) = signd;
      end

      for c5 = 1:1:c1         % absolute Zeit auf dt umstellen
          Xt(c5) = X(c5)-X(1);
      end

      ALL(:,1) = Xt;          % alle Daten in Matrix einlesen
      for c6 = 2:1:roim+1
          c7 = c6-1;
          ALL(:,c6) = signdall(:,c7);
      end

      All = sortrows(ALL);    % Signalwerte nach Zeit Sortieren
                        % da Spooler manchmal Reihenfolge der Bilder
                        % vertauscht