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 />
Nästa steg är att plocka ut de sex unika elementen, d = [ Dxx<br />
Dyy<br />
Dzz<br />
Dxy<br />
Dxz<br />
Dyz<br />
] ,ur<br />
tensorn D. Detta görs med hjälp <strong>av</strong> att vi från början vet ADC värdena för de sex olika gradienterna,<br />
y = [ ADC1<br />
ADC2<br />
ADC3<br />
ADC4<br />
ADC5<br />
T<br />
ADC6]<br />
. Detta sätts i relation till den konstruerade<br />
matrisen X (5.5) som beror på riktningsvektorerna ur ekvationerna 5.3 ovan.<br />
⎡ g<br />
⎢<br />
X = ⎢ ...<br />
⎢<br />
⎣ g<br />
2<br />
1x<br />
2<br />
nx<br />
g<br />
g<br />
2<br />
1y<br />
2<br />
ny<br />
g<br />
g<br />
2<br />
1z<br />
2<br />
nz<br />
2g<br />
2g<br />
1x<br />
Sammantaget ger detta ekvationen 5.6 nedan<br />
nx<br />
g<br />
g<br />
1y<br />
ny<br />
2g<br />
2g<br />
1x<br />
nx<br />
g<br />
g<br />
1z<br />
nz<br />
2g<br />
2g<br />
1y<br />
ny<br />
g<br />
g<br />
1z<br />
nz<br />
⎤<br />
⎥<br />
⎥<br />
⎥<br />
⎦<br />
(5.5)<br />
y = Xd<br />
(5.6)<br />
För att komma åt den eftertraktade d krävs lite omvandling <strong>av</strong> formeln 5.6 så att nedanstående<br />
ekvation ger oss d.<br />
d = X −1 y (5.7)<br />
Utifrån denna matris kan man plocka ut de unika tensorelementen <strong>och</strong> representera dem som en<br />
vektor.<br />
Ett räkneexempel på hur det kan se ut, när vi använder de gradienter som är angivna i 5.3 ovan. En<br />
start med att sätta in dessa i matris X enligt 5.5 ovan ger uttrycket<br />
X<br />
=<br />
⎡ 1<br />
⎢ 2<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
⎢<br />
0<br />
1<br />
2<br />
0<br />
1<br />
2<br />
1<br />
⎣ 2<br />
1<br />
2<br />
1<br />
2<br />
0<br />
1<br />
2<br />
1<br />
2<br />
0<br />
0<br />
1<br />
2<br />
1<br />
2<br />
1<br />
2<br />
0<br />
1<br />
2<br />
1<br />
0<br />
0<br />
0<br />
− 1<br />
0<br />
0<br />
0<br />
1<br />
0<br />
0<br />
− 1<br />
⎤<br />
0<br />
⎥<br />
⎥<br />
1 ⎥ ⎡ D<br />
⎥ ⎢<br />
⎢<br />
D<br />
0<br />
⎥<br />
⎥ ⎢ D<br />
⎥ ⋅ ⎢<br />
⎥ ⎢ D<br />
− 1<br />
⎥ ⎢ D<br />
⎥ ⎢<br />
0<br />
⎥ ⎢⎣<br />
D<br />
⎥<br />
0 ⎥<br />
⎦<br />
36<br />
xx<br />
yy<br />
zz<br />
xy<br />
xz<br />
yz<br />
⎡ 1 ⎛ S ⎞ ⎤ 1<br />
⎢ − ln<br />
⎜<br />
⎟ ⎥<br />
⎢ b ⎝ S0<br />
⎠ ⎥<br />
⎢ 1 ⎛ S ⎞ ⎥<br />
2<br />
⎢ − ln<br />
⎜<br />
⎟<br />
⎤<br />
⎥<br />
⎡ ADC1<br />
⎤<br />
⎢ b ⎝ S0<br />
⎠<br />
⎥<br />
⎥<br />
⎢ ⎥<br />
⎥<br />
⎢<br />
⎥<br />
⎢<br />
ADC2<br />
⎥ 1 ⎛ S ⎞ 3<br />
⎢ − ln<br />
⎜<br />
⎟<br />
⎥<br />
⎥<br />
⎢ ADC ⎥ 3 ⎢ b ⎝ S0<br />
⎠<br />
⎥ =<br />
⎥<br />
⎢ ⎥ =<br />
⎥<br />
⎢ ⎛ ⎞ ⎥<br />
⎢ ADC4<br />
⎥ 1 S4<br />
⎢ − ln<br />
⎜<br />
⎟<br />
⎥<br />
⎥<br />
⎢ ADC ⎥<br />
5<br />
⎥<br />
⎢ b ⎝ S0<br />
⎠ ⎥<br />
⎢ ⎥<br />
⎥ ⎢ ⎥ ⎢ ⎛ ⎞ ⎥<br />
⎦ ⎣ ADC6<br />
⎦ 1 S5<br />
⎢ − ln<br />
⎜<br />
⎟ ⎥<br />
⎢ b ⎝ S0<br />
⎠ ⎥<br />
⎢ 1 ⎛ S ⎞ ⎥<br />
6<br />
⎢ − ln<br />
⎜<br />
⎟ ⎥<br />
⎢⎣<br />
b ⎝ S0<br />
⎠ ⎥⎦