#!/bin/tcsh

### THIS IS THE CODE TO WHICH THE PREVIOUS POST REFERRED.

rm tmp__*
rm {$1}.MSK*

# ARG1: Input file
# ARG2: Bucket number from the input data
# ARG3: Cluster t-value
# ARG4: Cluster vmul (from AlphaSim)
# ARG5: Group

set group = $5
set subjects = `cat ~/afnidata/DDpilot/{$group}/SubList `

3dclust -1clip {$3} -1Dformat 5 {$4} {$1}+tlrc'['{$2}']' > tmp__clust

if (-e {$1}.MSK+tlrc.BRIK ) then
else
3dmerge -1clip {$3} -1clust_order 5 {$4} -prefix {$1}.MSK+tlrc {$1}+tlrc'['{$2}']'
endif

3dresample -dxyz 1.0 1.0 1.0 -prefix {$1}.order+tlrc -inset {$1}.MSK+tlrc

set numclusts = `1dcat tmp__clust'[0]'`
set COGx = `1dcat tmp__clust'[1]'`
set COGy = `1dcat tmp__clust'[2]'`
set COGz = `1dcat tmp__clust'[3]'`

echo FileVersion\: 3 > tmp__bvqx.voi
echo CoordsType\: TAL >> tmp__bvqx.voi
echo SubjectVOINamingConvention\: \<SUBJ\>_\<VOI\> >> tmp__bvqx.voi


echo NrOfVOIs\: $#numclusts >> tmp__bvqx.voi

foreach clust ( `count -digits 1 1 $#numclusts` )
set COGxx = $COGx[$clust]
set COGyy = $COGy[$clust]
set COGzz = $COGz[$clust]
whereami $COGxx $COGyy $COGzz 1
echo NameOfVOI\: Cluster{$clust} >> tmp__bvqx.voi
echo ColorOfVOI\: 255 75 0 >> tmp__bvqx.voi
3dmaskdump -mask {$1}.MSK.order+tlrc -mrange {$clust} {$clust} -xyz {$1}.MSK.order+tlrc > tmp__ClRep.1D
set numvox = `1dcat tmp__ClRep.1D'[0]'`
echo NrOfVoxels\: $#numvox >> tmp__bvqx.voi
1dcat tmp__ClRep.1D'[3]' tmp__ClRep.1D'[4]' tmp__ClRep.1D'[5]' > tmp__coords
1deval -expr 'int(0-a)' -a tmp__coords'[0]' > tmp__LPIX.1D
1deval -expr 'int(0-a)' -a tmp__coords'[1]' > tmp__LPIY.1D
1deval -expr 'int(a)' -a tmp__coords'[2]' > tmp__LPIZ.1D
1dcat tmp__LPIX.1D tmp__LPIY.1D tmp__LPIZ.1D >> tmp__bvqx.voi
end

cat tmp__bvqx.voi > {$1}.BVQX.voi

rm tmp__*