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 -iname "*[0-99].[0-9].mhd" -o -iname "*[0-99]*\]*.mhd" | wc -l`
32 if [ $nb_phase_file = 0 ]
34 echo "Error: no phase found"
38 ########## CT #########
41 list_phase_file=`find $1 -iname "*[0-99].[0-9].mhd"`
42 for phase_file in $list_phase_file
44 phase_file_name=`basename $phase_file`
45 if [[ ! -z `echo "$phase_file_name" | grep ","` ]]
47 prefix=`echo $phase_file_name | sed "s/,_.*/,_/"`
51 if [[ -z `echo "$list_pattern" | grep "$prefix"` ]]
53 list_pattern="$list_pattern $prefix"
58 for pattern in $list_pattern
61 if [ "$pattern" = "NONE" ]
66 nbph0=`find $1 -iname "${pattern}[0-9].[0-9].mhd" | wc -l`
67 orig0=`find $1 -iname "${pattern}[0-9].[0-9].mhd" | sort | head -n 1`
68 listph0=`find $1 -iname "${pattern}[0-9].[0-9].raw" | sort`
70 nbph1=`find $1 -iname "${pattern}[0-9][0-9].[0-9].mhd" | wc -l`
71 orig1=`find $1 -iname "${pattern}[0-9][0-9].[0-9].mhd" | sort | head -n 1`
72 listph1=`find $1 -iname "${pattern}[0-9][0-9].[0-9].raw" | sort`
74 nbph=$(( nbph0+nbph1 ))
75 orig=`echo -e "$orig0\n$orig1" | head -n 1`
76 listph="$listph0 $listph1"
78 file_name_4D="${pattern}_4D.mhd"
85 ############ PET ###########
88 list_phase_file=`find $1 -iname "*[0-99]*\]*.mhd"`
89 for phase_file in $list_phase_file
91 phase_file_name=`basename $phase_file`
92 if [[ ! -z `echo "$phase_file_name" | grep "[0-99]-.*\]"` ]]
94 prefix=`echo $phase_file_name | sed "s/.*[0-99]-/-/;s/\]_.*//"`
98 if [[ -z `echo "$list_pattern" | grep -- "$prefix"` ]]
100 list_pattern="$list_pattern $prefix"
105 for pattern in $list_pattern
108 if [ "$pattern" = "NONE" ]
113 nbph=`find $1 -iname "*[0-99]${pattern}\]*.mhd" | wc -l`
114 orig=`find $1 -iname "*[0-99]${pattern}\]*.mhd" | sort | head -n 1`
115 listph=`find $1 -iname "*[0-99]${pattern}\]*.raw" | sort`
117 file_name_4D=`basename "$orig" | sed "s/[0-99]${pattern}\]/${pattern}\]/;s/_.mhd/_4D.mhd/"`