From 2cbcfe688ed379473b4e85ea4f15e6fa88b4c819 Mon Sep 17 00:00:00 2001
From: Saad Jbabdi <saad@fmrib.ox.ac.uk>
Date: Wed, 1 Aug 2007 17:10:16 +0000
Subject: [PATCH] renamed bedpost_datacheck, and removed checking for nodif
---
bedpostx_datacheck | 58 ++++++++++++++++++++++++++++++++++++++++++++++
1 file changed, 58 insertions(+)
create mode 100755 bedpostx_datacheck
diff --git a/bedpostx_datacheck b/bedpostx_datacheck
new file mode 100755
index 0000000..2679ea7
--- /dev/null
+++ b/bedpostx_datacheck
@@ -0,0 +1,58 @@
+#!/bin/sh
+
+if [ $# -lt 1 ]; then
+ echo "Usage: $0 data_dir"
+ exit 1;
+fi
+
+DIR=$1;
+for img in ${DIR}/data ${DIR}/nodif_brain_mask ;do
+ if [ `${FSLDIR}/bin/imtest $img` -eq 1 ] ; then
+ echo $img
+ ${FSLDIR}/bin/fslinfo $img;
+ echo ""
+ else
+ echo $img does not exist;
+ fi
+done
+
+for bv in ${DIR}/bvals ${DIR}/bvecs;do
+ echo " num lines in $bv "
+ cat $bv |grep -v "^[ ]*$"|wc -l
+ echo " num words in $bv "
+ cat $bv|wc -w
+done
+
+if [ `${FSLDIR}/bin/imtest ${DIR}/data` -eq 1 -a `${FSLDIR}/bin/imtest ${DIR}/nodif_brain_mask` -eq 1 ];then
+
+ dx=`${FSLDIR}/bin/fslval ${DIR}/data dim1`
+ dy=`${FSLDIR}/bin/fslval ${DIR}/data dim2`
+ dz=`${FSLDIR}/bin/fslval ${DIR}/data dim3`
+ dt=`${FSLDIR}/bin/fslval ${DIR}/data dim4`
+
+ nbmx=`${FSLDIR}/bin/fslval ${DIR}/nodif_brain_mask dim1`
+ nbmy=`${FSLDIR}/bin/fslval ${DIR}/nodif_brain_mask dim2`
+ nbmz=`${FSLDIR}/bin/fslval ${DIR}/nodif_brain_mask dim3`
+
+
+ if [ $dx -ne $nbmx -o $dy -ne $nbmy -o $dz -ne $nbmz ];then
+ echo "data dimensions do not match mask dimensions"
+ fi
+
+ bvallen=`cat ${DIR}/bvals|grep -v "^[ ]*$"|wc -w`
+ bveclenw=`cat ${DIR}/bvecs|wc -w`
+ bveclenl=`cat ${DIR}/bvecs|grep -v "^[ ]*$"|wc -l`
+ bveclen=`echo "$bveclenw / $bveclenl"|bc`
+
+
+ if [ $bveclen -ne $bvallen ];then
+ echo "bvecs and bvals are incompatible"
+ fi
+ if [ $bvallen -ne $dt ];then
+ echo "number of elements in bvals is not equal to number of vols in data"
+ fi
+ if [ $bveclen -ne $dt ];then
+ echo "number of elements per line in bvecs is not equal to number of vols in data"
+ fi
+
+fi
--
GitLab