6 cat $orig | sed "s/NDims = .*/NDims = 4/
7 s/TransformMatrix = .*/TransformMatrix = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1/
9 /CenterOfRotation/ s/.*/& 0/
10 s/AnatomicalOrientation = .*/AnatomicalOrientation = ????/
11 /ElementSpacing/ s/.*/& 1/
12 /DimSize/ s/.*/& $nbph/
13 s/ElementDataFile = .*/ElementDataFile = LIST/" > "$1/$file_name_4D"
18 echo "$ph" >> "$1/$file_name_4D"
22 #################################################
23 # create_mhd_4D argument : repertoire #
24 #################################################
27 echo "Usage: create_mhd_4D.sh DIRECTORY"
31 nb_phase_file=`find $1 -maxdepth 1 -iname "*[0-9]*.mhd" | wc -l`
32 if [ $nb_phase_file = 0 ]
34 echo "No phase found in $1"
40 list_phase_file=`find $1 -maxdepth 1 -iname "*[0-9]*.mhd"`
41 for phase_file in $list_phase_file
43 phase_file_name=`basename $phase_file`
44 if [[ ! -z `echo "$phase_file_name" | grep "__[0-9]"` ]]
46 prefix=`echo $phase_file_name | sed "s/__[0-9].*/__/"`
48 if [[ ! -z `echo "$phase_file_name" | grep "[0-9]-.*\]"` ]]
50 if [[ ! -z `echo "$phase_file_name" | grep "[0-9][0-9]-.*\]"` ]]
52 prefix=`echo $phase_file_name | sed "s/[0-9][0-9]-.*//"`
54 prefix=`echo $phase_file_name | sed "s/[0-9]-.*//"`
61 if [[ -z `echo "$list_prefix" | grep "$prefix"` ]]
63 list_prefix="$list_prefix $prefix"
68 for prefix in $list_prefix
70 if [ "$prefix" = "NONE" ]
76 list_phase_file_prefix=`find $1 -maxdepth 1 -iname "${prefix}[0-9]*.mhd"`
77 for phase_file_prefix in $list_phase_file_prefix
79 phase_file_prefix_name=`basename $phase_file_prefix`
80 if [[ ! -z `echo "$phase_file_prefix_name" | grep "__[0-9]"` ]]
84 if [[ ! -z `echo "$phase_file_prefix_name" | grep "[0-9]-.*\]"` ]]
86 suffix=`echo $phase_file_prefix_name | sed "s/.*[0-9]-//;s/_\.mhd//;s/\.mhd//"`
92 if [[ -z `echo "$list_suffix" | grep "$suffix"` ]]
94 list_suffix="$list_suffix $suffix"
98 for suffix in $list_suffix
100 if [ "$suffix" = "NONE" ]
104 nbph=`find $1 -maxdepth 1 -iname "*${prefix}*[0-9]*${suffix}*.mhd" | wc -l`
105 orig=`find $1 -maxdepth 1 -iname "*${prefix}*[0-9]*${suffix}*.mhd" | sort | head -n 1`
106 listph=`find $1 -maxdepth 1 -iname "*${prefix}*[0-9]*${suffix}*.raw" | sort`
108 file_name_4D="${prefix}4D${suffix}.mhd"