From aacc8194e42edb7a3bb1fbbd91ab3e954cdf823e Mon Sep 17 00:00:00 2001
From: Saad Jbabdi <saad@fmrib.ox.ac.uk>
Date: Fri, 26 Oct 2012 08:23:26 +0000
Subject: [PATCH] rotate the bvecs finally in FDT
---
fdt_rotate_bvecs | 70 ++++++++++++++++++++++++++++++++++++++++++++++++
1 file changed, 70 insertions(+)
create mode 100644 fdt_rotate_bvecs
diff --git a/fdt_rotate_bvecs b/fdt_rotate_bvecs
new file mode 100644
index 0000000..e96a86a
--- /dev/null
+++ b/fdt_rotate_bvecs
@@ -0,0 +1,70 @@
+#!/bin/sh
+
+if [ "$3" == "" ] ; then
+ echo "Usage: <original bvecs> <rotated bvecs> <ecclog>"
+ echo ""
+ echo "<ecclog> is the output log file from ecc"
+ echo ""
+ exit 1;
+fi
+
+i=$1
+o=$2
+ecclog=$3
+
+if [ ! -e $1 ] ; then
+ echo "Source bvecs $1 does not exist!"
+ exit 1
+fi
+if [ ! -e $ecclog ]; then
+ echo "Ecc log file $3 does not exist!"
+ exit 1
+fi
+
+ii=1
+rm -f $o
+tmpo=${o}$$
+cat ${ecclog} | while read line; do
+ echo $ii
+ if [ "$line" == "" ];then break;fi
+ read line;
+ read line;
+ read line;
+
+ echo $line > $tmpo
+ read line
+ echo $line >> $tmpo
+ read line
+ echo $line >> $tmpo
+ read line
+ echo $line >> $tmpo
+ read line
+
+ m11=`avscale $tmpo | grep Rotation -A 1 | tail -n 1| awk '{print $1}'`
+ m12=`avscale $tmpo | grep Rotation -A 1 | tail -n 1| awk '{print $2}'`
+ m13=`avscale $tmpo | grep Rotation -A 1 | tail -n 1| awk '{print $3}'`
+ m21=`avscale $tmpo | grep Rotation -A 2 | tail -n 1| awk '{print $1}'`
+ m22=`avscale $tmpo | grep Rotation -A 2 | tail -n 1| awk '{print $2}'`
+ m23=`avscale $tmpo | grep Rotation -A 2 | tail -n 1| awk '{print $3}'`
+ m31=`avscale $tmpo | grep Rotation -A 3 | tail -n 1| awk '{print $1}'`
+ m32=`avscale $tmpo | grep Rotation -A 3 | tail -n 1| awk '{print $2}'`
+ m33=`avscale $tmpo | grep Rotation -A 3 | tail -n 1| awk '{print $3}'`
+
+ X=`cat $i | awk -v x=$ii '{print $x}' | head -n 1 | tail -n 1 | awk -F"E" 'BEGIN{OFMT="%10.10f"} {print $1 * (10 ^ $2)}' `
+ Y=`cat $i | awk -v x=$ii '{print $x}' | head -n 2 | tail -n 1 | awk -F"E" 'BEGIN{OFMT="%10.10f"} {print $1 * (10 ^ $2)}' `
+ Z=`cat $i | awk -v x=$ii '{print $x}' | head -n 3 | tail -n 1 | awk -F"E" 'BEGIN{OFMT="%10.10f"} {print $1 * (10 ^ $2)}' `
+ rX=`echo "scale=7; ($m11 * $X) + ($m12 * $Y) + ($m13 * $Z)" | bc -l`
+ rY=`echo "scale=7; ($m21 * $X) + ($m22 * $Y) + ($m23 * $Z)" | bc -l`
+ rZ=`echo "scale=7; ($m31 * $X) + ($m32 * $Y) + ($m33 * $Z)" | bc -l`
+
+ if [ "$ii" -eq 1 ];then
+ echo $rX > $o;echo $rY >> $o;echo $rZ >> $o
+ else
+ cp $o $tmpo
+ (echo $rX;echo $rY;echo $rZ) | paste $tmpo - > $o
+ fi
+
+ let "ii+=1"
+
+done
+rm -f $tmpo
--
GitLab