Skip to content
Snippets Groups Projects
Commit 74161965 authored by Stephen Smith's avatar Stephen Smith
Browse files

CVS: ----------------------------------------------------------------------

parent 3e508845
No related branches found
No related tags found
No related merge requests found
...@@ -156,6 +156,8 @@ else ...@@ -156,6 +156,8 @@ else
echo "" >> ${I}.sienax echo "" >> ${I}.sienax
echo "---------- convert brain volume into normalised volume --------------" >> ${I}.sienax echo "---------- convert brain volume into normalised volume --------------" >> ${I}.sienax
echo "" >> ${I}.sienax
echo "tissue volume unnormalised-volume" >> ${I}.sienax
if [ $regional = 1 ] ; then if [ $regional = 1 ] ; then
${FSLDIR}/bin/avwmaths_32R ${I}_talmaskbrain_pve_1 -mas ${I}_talmask_segperiph ${I}_talmaskbrain_pve_1_segperiph ${FSLDIR}/bin/avwmaths_32R ${I}_talmaskbrain_pve_1 -mas ${I}_talmask_segperiph ${I}_talmaskbrain_pve_1_segperiph
S=`${FSLDIR}/bin/avwstats ${I}_talmaskbrain_pve_1_segperiph -m -v` S=`${FSLDIR}/bin/avwstats ${I}_talmaskbrain_pve_1_segperiph -m -v`
...@@ -163,9 +165,7 @@ else ...@@ -163,9 +165,7 @@ else
xb=`echo $S | awk '{print $3}'` xb=`echo $S | awk '{print $3}'`
uxg=`echo "2 k $xa $xb * 1 / p" | dc -` uxg=`echo "2 k $xa $xb * 1 / p" | dc -`
xg=`echo "2 k $xa $xb * $vscale * 1 / p" | dc -` xg=`echo "2 k $xa $xb * $vscale * 1 / p" | dc -`
echo "" >> ${I}.sienax echo "pgrey $xg $uxg (peripheral grey)" >> ${I}.sienax
echo "unnormalised-pgrey $uxg" >> ${I}.sienax
echo "pgrey $xg" >> ${I}.sienax
${FSLDIR}/bin/avwmaths_32R ${I}_talmaskbrain_pve_0 -mas ${I}_talmask_segvent ${I}_talmaskbrain_pve_0_segvent ${FSLDIR}/bin/avwmaths_32R ${I}_talmaskbrain_pve_0 -mas ${I}_talmask_segvent ${I}_talmaskbrain_pve_0_segvent
S=`${FSLDIR}/bin/avwstats ${I}_talmaskbrain_pve_0_segvent -m -v` S=`${FSLDIR}/bin/avwstats ${I}_talmaskbrain_pve_0_segvent -m -v`
...@@ -173,34 +173,25 @@ else ...@@ -173,34 +173,25 @@ else
xb=`echo $S | awk '{print $3}'` xb=`echo $S | awk '{print $3}'`
uxg=`echo "2 k $xa $xb * 1 / p" | dc -` uxg=`echo "2 k $xa $xb * 1 / p" | dc -`
xg=`echo "2 k $xa $xb * $vscale * 1 / p" | dc -` xg=`echo "2 k $xa $xb * $vscale * 1 / p" | dc -`
echo "" >> ${I}.sienax echo "vcsf $xg $uxg (ventricular CSF)" >> ${I}.sienax
echo "unnormalised-vcsf $uxg" >> ${I}.sienax
echo "vcsf $xg" >> ${I}.sienax
fi fi
S=`${FSLDIR}/bin/avwstats ${I}_talmaskbrain_pve_1 -m -v` S=`${FSLDIR}/bin/avwstats ${I}_talmaskbrain_pve_1 -m -v`
xa=`echo $S | awk '{print $1}'` xa=`echo $S | awk '{print $1}'`
xb=`echo $S | awk '{print $3}'` xb=`echo $S | awk '{print $3}'`
ugrey=`echo "2 k $xa $xb * 1 / p" | dc -` ugrey=`echo "2 k $xa $xb * 1 / p" | dc -`
ngrey=`echo "2 k $xa $xb * $vscale * 1 / p" | dc -` ngrey=`echo "2 k $xa $xb * $vscale * 1 / p" | dc -`
echo "" >> ${I}.sienax echo "GREY $ngrey $ugrey" >> ${I}.sienax
echo "unnormalised-GREY $ugrey" >> ${I}.sienax
echo "GREY $ngrey" >> ${I}.sienax
S=`${FSLDIR}/bin/avwstats ${I}_talmaskbrain_pve_2 -m -v` S=`${FSLDIR}/bin/avwstats ${I}_talmaskbrain_pve_2 -m -v`
xa=`echo $S | awk '{print $1}'` xa=`echo $S | awk '{print $1}'`
xb=`echo $S | awk '{print $3}'` xb=`echo $S | awk '{print $3}'`
uwhite=`echo "2 k $xa $xb * 1 / p" | dc -` uwhite=`echo "2 k $xa $xb * 1 / p" | dc -`
nwhite=`echo "2 k $xa $xb * $vscale * 1 / p" | dc -` nwhite=`echo "2 k $xa $xb * $vscale * 1 / p" | dc -`
echo "" >> ${I}.sienax echo "WHITE $nwhite $uwhite" >> ${I}.sienax
echo "unnormalised-WHITE $uwhite" >> ${I}.sienax
echo "WHITE $nwhite" >> ${I}.sienax
fi fi
ubrain=`echo "2 k $uwhite $ugrey + 1 / p" | dc -` ubrain=`echo "2 k $uwhite $ugrey + 1 / p" | dc -`
nbrain=`echo "2 k $nwhite $ngrey + 1 / p" | dc -` nbrain=`echo "2 k $nwhite $ngrey + 1 / p" | dc -`
echo "" >> ${I}.sienax echo "BRAIN $nbrain $ubrain" >> ${I}.sienax
echo "unnormalised-BRAIN $ubrain" >> ${I}.sienax
echo "BRAIN $nbrain" >> ${I}.sienax
echo "" >> ${I}.sienax
${FSLDIR}/bin/overlay 1 1 -c ${I} -a ${I}_talmaskbrain_seg 1.9 5 ${I}_render ${FSLDIR}/bin/overlay 1 1 -c ${I} -a ${I}_talmaskbrain_seg 1.9 5 ${I}_render
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment