From: Romulo Pinho <pinho@lyon.fnclcc.fr>
Date: Mon, 17 Oct 2011 15:20:46 +0000 (+0200)
Subject: calculation of motion amplitude itself...
X-Git-Tag: v1.3.0~184
X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=commitdiff_plain;h=e2f16c750e87742d77560833f6ac394f3d7aecfc;p=clitk.git

calculation of motion amplitude itself...

- that is, (max - min) over time
---

diff --git a/scripts/calculate_motion_amplitude.sh b/scripts/calculate_motion_amplitude.sh
index 9fa0aa9..b40c766 100755
--- a/scripts/calculate_motion_amplitude.sh
+++ b/scripts/calculate_motion_amplitude.sh
@@ -43,10 +43,16 @@ clitkResampleImage -i ${roi_mask} -o ${roi_mask2} --like ${vector_file}
 dir=( sagittal coronal axial )
 for i in 0 1 2; do
   echo Motion along ${dir[$i]} direction
-  clitkImageStatistics -i ${vector_file} -m ${roi_mask2} -c $i --verbose 2> /dev/null | tail -n 8 | head -n 6
+  clitkImageStatistics -i ${vector_file} -m ${roi_mask2} -c $i --verbose 2> /dev/null | tail -n 8 | head -n 6 > /tmp/res.txt
+  min=`tail -n 2 /tmp/res.txt | head -n 1 | cut -f 2 -d ':'`
+  max=`tail -n 1 /tmp/res.txt | cut -f 2 -d ':'`
+  amp=`echo $max-$min | bc`
+  cat /tmp/res.txt
+  echo "Amplitude: $amp"
   echo
 done
 
 rm `basename $roi_mask .mhd`.{mhd,raw}
 rm `basename $roi_mask2 .mhd`.{mhd,raw}
 rm /tmp/err.txt
+rm /tmp/res.txt