Rss Feed
Tweeter button
Facebook button
Technorati button
Reddit button
Myspace button
Linkedin button
Webonews button
Delicious button
Digg button
Flickr button
Stumbleupon button
Newsvine button
Youtube button

Notes:

  • The site will be morphing over the next little while.
  • I am having some issues with tabs/spaces in the Python code. Sometimes they are correct, sometimes they get eaten. I am trying to figure it out.

Please feel free to leave a comment on a post if it interests you or if you have questions

Quantitiative Diffusion Measures

Add comments

Given the eigenvalues \lambda_1, \lambda_2, \lambda_3 and eigenvectors v_1, v_2, v_3 one can define many metrics.

Fractional anisotropy

Given a set of three eigenvalues \lambda_1, \lambda_2, \lambda_3  (variable ‘evalues’ in Matlab) the fractional anisotropy (FA) is calculated as:

FA = \frac{ \left(\lambda_1 - \lambda_2\right)^2+ \left( \lambda_1 - \lambda_3\right)^2 + \left( \lambda_2 - \lambda_3\right)^2}{2\left( \lambda_1^2 + \lambda_2^2 + \lambda_3^2 \right)}

which can be calculated using the following Matlab code:

1
2
3
4
sum_sq_evalues = sum(evalues.^2);

sde = (evalues(1)-evalues(2))^2 + (evalues(2)-evalues(3))^2 + (evalues(1)-evalues(3))^2;
fa =  sqrt( sde / (2*sum_sq_evalues) );

For this, it doesn’t matter that the eigenvalues are sorted, though typically they are sorted from highest to lowest.

Apparent Diffusion Coefficient

Relative Anisotropy

Lambda Parallel

Lambda Perpendicular

Posted by craig at 11:49 am

Leave a Reply

(required)

(required)

Spam Protection by WP-SpamFree

© 2010 Math, Computing and Research Please leave a comment or contact me craig@mri.brechmos.org if you have any questions. Suffusion WordPress theme by Sayontan Sinha