beskrivning och utvärdering av diffusions mr - Örebro universitet

Examensarbete 10p Örebro Universitet VT 2006
10. MatLab kod
%Programmet beräknar och lägger samman bilder för ADC-karta, 4 olika
%varianter:
%1) 2 b-värden Isotropiska bilder
%2) 5 b-värden Isotropiska bilder
%3) 2 b-värden och alla spatiala riktningar på diffusion, M, P och S
%4) 5 b-värden och alla spatiala riktningar på diffusion, M, P och S
clear variables;
% Variables;
Column = 60;
%Anger vilka bilder som ska läsas in för beräkning av ADC-karta
im_nr =[1516:1887];
%Anger vika bilder som ska läsas in för jämförelse, färdiga ADC-kartor
im_nr1 =[1920:1934];
%Lägger bildnamn i vektorn d
for cnt7 = 1:length(im_nr)
tmp = sprintf('IM_%.04d', im_nr(cnt7));
d(cnt7) = {tmp};
end
%Lägger bildnamn i vektorn g, jämförelsebilder dvs färdiga ADC-kartor
for cnt8 = 1:length(im_nr1)
tmp1 = sprintf('IM_%.04d', im_nr1(cnt8));
g(cnt8) = {tmp1};
end
%Loopar igenom bilderna och beräknar ADC-kartor
for cnt13 = 1:length(im_nr1)
tmp3 = dicominfo(char(g(cnt13)));
for cnt9 = 1:17:length(im_nr)
d(cnt9);
tmp4 = dicominfo(char(d(cnt9)));
%För att ta samma snitt att jämföra med som man beräknat
if int8(tmp4.SliceLocation) == int8(tmp3.SliceLocation)
int8(tmp3.SliceLocation)
int8(tmp4.SliceLocation)
%Referens med b värde 0
im0_info = dicominfo(char(d(cnt9)));
%b-värde 250 Isotropisk bild
im1_info = dicominfo(char(d(cnt9+4)));
%b-värde 500 Isotropisk bild
im2_info = dicominfo(char(d(cnt9+8)));
%b-värde 750 Isotropisk bild
im3_info = dicominfo(char(d(cnt9+12)));
45