4 #################################################
5 # create_mhd_4D argument : repertoire #
6 #################################################
9 echo "Usage: create_mhd_4D.sh DIRECTORY"
13 nb_phase_file=`find $1 -iname "*[0-99].[0-9].mhd" | wc -l`
14 if [ $nb_phase_file = 0 ]
16 echo "Error: no phase found"
21 list_phase_file=`find $1 -iname "*[0-99].[0-9].mhd"`
22 for phase_file in $list_phase_file
24 phase_file_name=`basename $phase_file`
25 if [[ ! -z `echo "$phase_file_name" | grep ","` ]]
27 preffix=`echo $phase_file_name | sed "s/,_.*/,_/"`
31 if [[ -z `echo "$list_pattern" | grep "$preffix"` ]]
33 list_pattern="$list_pattern $preffix"
38 for pattern in $list_pattern
41 if [ "$pattern" = "NONE" ]
46 nbph0=`find $1 -iname "${pattern}[0-9].[0-9].mhd" | wc -l`
47 orig0=`find $1 -iname "${pattern}[0-9].[0-9].mhd" | sort | head -n 1`
48 listph0=`find $1 -iname "${pattern}[0-9].[0-9].raw" | sort`
50 nbph1=`find $1 -iname "${pattern}[0-9][0-9].[0-9].mhd" | wc -l`
51 orig1=`find $1 -iname "${pattern}[0-9][0-9].[0-9].mhd" | sort | head -n 1`
52 listph1=`find $1 -iname "${pattern}[0-9][0-9].[0-9].raw" | sort`
54 nbph=$(( nbph0+nbph1 ))
55 orig=`echo -e "$orig0\n$orig1" | head -n 1`
56 listph="$listph0 $listph1"
58 cat $orig | sed "s/NDims = .*/NDims = 4/
59 s/TransformMatrix = .*/TransformMatrix = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1/
61 /CenterOfRotation/ s/.*/& 0/
62 s/AnatomicalOrientation = .*/AnatomicalOrientation = ????/
63 /ElementSpacing/ s/.*/& 1/
64 /DimSize/ s/.*/& $nbph/
65 s/ElementDataFile = .*/ElementDataFile = LIST/" > "$1/${pattern}_4D.mhd"
70 echo "$ph" >> "$1/${pattern}_4D.mhd"