beskrivning och utvärdering av diffusions mr - Örebro universitet
beskrivning och utvärdering av diffusions mr - Örebro universitet
beskrivning och utvärdering av diffusions mr - Örebro universitet
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Examensarbete 10p <strong>Örebro</strong> Universitet VT 2006<br />
%b-värde 1000 Isotropisk bild<br />
im4_info = dicominfo(char(d(cnt9+16)));<br />
%Motsvarande bild i ADC map<br />
im5_info = dicominfo(char(g(cnt13)));<br />
%Referens med b-värde 0 omskalad bild<br />
im0_bild = double(dico<strong>mr</strong>ead(im0_info)*im0_info.RescaleSlope+...<br />
im0_info.RescaleIntercept);<br />
%Bilderna omskalas <strong>och</strong> läses in<br />
im1_bild = double(dico<strong>mr</strong>ead(im1_info)*im1_info.RescaleSlope+...<br />
im1_info.RescaleIntercept);<br />
im2_bild = double(dico<strong>mr</strong>ead(im2_info)*im2_info.RescaleSlope+...<br />
im2_info.RescaleIntercept);<br />
im3_bild = double(dico<strong>mr</strong>ead(im3_info)*im3_info.RescaleSlope+...<br />
im3_info.RescaleIntercept);<br />
im4_bild = double(dico<strong>mr</strong>ead(im4_info)*im4_info.RescaleSlope+...<br />
im4_info.RescaleIntercept);<br />
%ADC map bilden omskalas <strong>och</strong> läses in<br />
im5_bild = double(dico<strong>mr</strong>ead(im5_info)*im5_info.RescaleSlope+...<br />
im5_info.RescaleIntercept);<br />
mask=logical(im5_bild>0); %mask för att begränsa beräkningsytan<br />
[M,N] = size(im1_bild);<br />
ADC_value = zeros(M,N);<br />
%b-värden läggs i en vektor<br />
bvalues=[im0_info.Private_2001_1003...<br />
im1_info.Private_2001_1003 im2_info.Private_2001_1003...<br />
im3_info.Private_2001_1003 im4_info.Private_2001_1003];<br />
index = find(mask==1);<br />
%Loopar igenom pixlarna i bilderna för att jämföra<br />
%signalintensiteterna<br />
cnt4 = 1;<br />
for cnt3 = 1:length(index)<br />
S0 = im0_bild(index(cnt3));<br />
S1 = im1_bild(index(cnt3));<br />
S2 = im2_bild(index(cnt3));<br />
S3 = im3_bild(index(cnt3));<br />
S4 = im4_bild(index(cnt3));<br />
% Signalvärdes vektor<br />
SI_matrix = [S1 S2 S3 S4];<br />
index1 = find(SI_matrix==0);<br />
if S0 ~= 0 % Undvik att dividera med 0<br />
if length(index1)==4 %Hitta pixlar som saknar värde<br />
PositionToFill(cnt4) = index(cnt3);<br />
cnt4 = cnt4 + 1;<br />
else<br />
% Hitta nollor<br />
46