- cat $orig | sed "s/NDims = .*/NDims = 4/
- s/TransformMatrix = .*/TransformMatrix = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1/
- /Offset/ s/.*/& 0/
- /CenterOfRotation/ s/.*/& 0/
- s/AnatomicalOrientation = .*/AnatomicalOrientation = ????/
- /ElementSpacing/ s/.*/& 1/
- /DimSize/ s/.*/& $nbph/
- s/ElementDataFile = .*/ElementDataFile = LIST/" > "$1/${pattern}_4D.mhd"
+ file_name_4D="${pattern}_4D.mhd"
+
+ write_mhd_4D $1
+
+done
+
+
+############ PET ###########
+
+list_pattern=""
+list_phase_file=`find $1 -iname "*[0-99]*\]*.mhd"`
+for phase_file in $list_phase_file
+do
+ phase_file_name=`basename $phase_file`
+ if [[ ! -z `echo "$phase_file_name" | grep "[0-99]-.*\]"` ]]
+ then
+ prefix=`echo $phase_file_name | sed "s/.*[0-99]-/-/;s/\]_.*//"`
+ else
+ prefix="NONE"
+ fi
+ if [[ -z `echo "$list_pattern" | grep -- "$prefix"` ]]
+ then
+ list_pattern="$list_pattern $prefix"
+ fi
+done
+