d = dir('*.dcm');
Q = zeros (3,numel(d));
Metadata1 = dicominfo(d(1).name);
for p = 1:numel(d)
metadata = dicominfo(d(p).name);
%overskriver data fra en header med data fra en anden header.
metadata.(dicomlookup('0008','0032'))=Metadata1.(dicomlookup('0008','0032'));
%if (p>1)
%metadatatemp =dicominfo(d(p-1).name);
% if ~(metadatatemp.(dicomlookup('0020', '0013')) == metadata.(dicomlookup('0020', '0013'))-1)
% metadata.(dicomlookup('0020','0013'))=metadatatemp.(dicomlookup('0020','0013'))+2;
% end
%end
Y=dicomread(metadata);
dicomwrite(Y,d(1).name,metadata);
Q(1,p)= metadata.(dicomlookup('0020', '1041')); %Slice Location - Har en intern sammenhæng imellem snittykkelse og position - kan bruges til at sammenstykke de to skans med.
Q(2,p)= metadata.(dicomlookup('0020', '0011')); %Series Number - Afhænger kun af skan men er den samme for alle billeder indenfor skanningen
Q(3,p)= metadata.(dicomlookup('0020', '0013')); %Instance Number - Stiger med 1 igennem hver skan internt for hver skan
%disp( metadata.(dicomlookup('0020','0032'))); %ImagePositionPatient - kun den tredje variable (long) ændrer sig med 5 altså snittykkelsen. afhængig af skan.
%disp( metadata.(dicomlookup('0008', '0016'))); %SOP Class UID - Er den samme for alle uanset skan eller billede
%disp( metadata.(dicomlookup('0008', '0018'))); %SOP Instance UID - Afhænger af billede og skan
%disp( metadata.(dicomlookup('0020', '000D'))); %Study Instance UID - Afhænger kun af skan men er den samme for alle billeder indenfor skanningen
disp( metadata.(dicomlookup('0020', '000E'))); %Series Instance UID - Afhænger kun af skan men er den samme for alle billeder indenfor skanningen
disp( metadata.(dicomlookup('0020', '0013'))); %Instance Number - Stiger med 1 igennem hver skan internt for hver skan
disp( metadata.(dicomlookup('0020','0037'))); %Image Orientation
disp( metadata.(dicomlookup('0020','0032'))); %Image Position
disp( metadata.(dicomlookup('0018','0050'))); %SliceThickness
% disp( metadata.(dicomlookup('0018','0088'))); %DistanceBetweenSlices
%disp( metadata.(dicomlookup('0008', '0030'))); %Study Time - Afhænger kun af skan men er den samme for alle billeder indenfor skanningen
disp( metadata.(dicomlookup('0008', '0032'))); %Acquisition Time - Afhænger kun af skan men er den samme for alle billeder indenfor skanningen
end
Q=Q';
Q