Browse Source

Code and output for a global tSNR map across all scans

Michael Hanke 8 years ago
parent
commit
b4fd1274dd

+ 26 - 0
code/calc_global_tsnr_map.sh

@@ -0,0 +1,26 @@
+#!/bin/bash
+
+set -e
+set -u
+
+odir="qa/tsnr/bold3Tp2"
+mkdir -p $odir
+
+for s in sub-*; do
+  echo "$s"
+  $FSLDIR/bin/fslmerge -t $odir/${s}_alltsnr $s/in_bold3Tp2/${s}_*_bold_tsnr.nii.gz
+  $FSLDIR/bin/fslmaths $odir/${s}_alltsnr -Tmean $odir/${s}_avgtsnr
+  $FSLDIR/bin/applywarp \
+    -i $odir/${s}_avgtsnr \
+    -r src/tnt/templates/grpbold3Tp2/brain.nii.gz \
+    -o $odir/${s}_avgtsnr_grpbold3Tp2 \
+    -w src/tnt/${s}/bold3Tp2/in_grpbold3Tp2/subj2tmpl_warp.nii.gz \
+    --interp=trilinear
+done
+
+$FSLDIR/bin/fslmerge -t $odir/alltsnr_grpbold3Tp2 $odir/*_avgtsnr_grpbold3Tp2.nii.gz
+$FSLDIR/bin/fslmaths $odir/alltsnr_grpbold3Tp2 -Tmean $odir/avgtsnr_grpbold3Tp2
+$FSLDIR/bin/fslmaths $odir/alltsnr_grpbold3Tp2 -Tstd $odir/stdtsnr_grpbold3Tp2
+$FSLDIR/bin/fslmaths $odir/avgtsnr_grpbold3Tp2 -uthr 60 -bin $odir/avgtsnr_grpbold3Tp2_lt60
+$FSLDIR/bin/fslmaths $odir/avgtsnr_grpbold3Tp2 -uthr 50 -bin $odir/avgtsnr_grpbold3Tp2_lt50
+

+ 1 - 0
qa/tsnr/bold3Tp2/avgtsnr_grpbold3Tp2.nii.gz

@@ -0,0 +1 @@
+../../../.git/annex/objects/G3/jG/MD5E-s1103203--6a3339e5f6754f16a8a7d863e4c32d60.nii.gz/MD5E-s1103203--6a3339e5f6754f16a8a7d863e4c32d60.nii.gz

+ 1 - 0
qa/tsnr/bold3Tp2/stdtsnr_grpbold3Tp2.nii.gz

@@ -0,0 +1 @@
+../../../.git/annex/objects/mv/jZ/MD5E-s1087175--c97bd8112787b2d9b6718f6b0a4fb594.nii.gz/MD5E-s1087175--c97bd8112787b2d9b6718f6b0a4fb594.nii.gz