1
+ readarray arr < $1
2
+ newfile=${1} new
3
+ oldiso=0
4
+ idx=0
5
+ echo -n $( echo ${arr[0]} | awk -F " 0," ' {print $1;}' ) > $newfile
6
+ for i in " ${arr[@]} " ; do
7
+ curiso=$( echo $i | cut -d" ," -f 4 | uniq | sed ' s/^ *//' )
8
+ if [[ $curiso -gt $oldiso ]]
9
+ then
10
+ if [[ $idx -gt 0 ]]
11
+ then
12
+ beg=$( echo $i | awk -F " 0," ' {print $1;}' )
13
+ ar=$( echo $ar | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
14
+ ag=$( echo $ag | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
15
+ ab=$( echo $ab | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
16
+ br=$( echo $br | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
17
+ bg=$( echo $bg | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
18
+ bb=$( echo $bb | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
19
+ echo " 0, {$ar ,$ag ,$ab },{$br ,$bg ,$bb }}," >> $newfile
20
+ echo -n $beg >> $newfile
21
+ fi
22
+ oldiso=$curiso
23
+ ar=0
24
+ ag=0
25
+ ab=0
26
+ br=0
27
+ bg=0
28
+ bb=0
29
+ nr=0
30
+ else
31
+ # sum them
32
+ res=$( echo $i | awk -F " 0," ' {print $2;}' )
33
+ ar=$( echo ${res: 2} | cut -d" ," -f 1 | awk -v a=$ar ' {printf "%.14g",a+$0;}' )
34
+ ag=$( echo $res | cut -d" ," -f 2 | awk -v a=$ag ' {printf "%.14g",a+$0;}' )
35
+ x=$( echo $res | cut -d" ," -f 3)
36
+ ab=$( echo ${x: 1: ${# x} -2} | awk -v a=$ab ' {printf "%.14g",a+$0;}' )
37
+ res=$( echo $i | awk -F " }, {" ' {print $2;}' )
38
+ br=$( echo $res | cut -d" ," -f 1 | awk -v a=$br ' {printf "%.14g",a+$0;}' )
39
+ bg=$( echo $res | cut -d" ," -f 2 | awk -v a=$bg ' {printf "%.14g",a+$0;}' )
40
+ x=$( echo $res | cut -d" ," -f 3)
41
+ bb=$( echo ${x: 1: ${# x} -3} | awk -v a=$bb ' {printf "%.14g",a+$0;}' )
42
+ # echo $ar,$ag,$ab $br,$bg,$bb
43
+ nr=$(( $nr + 1 ))
44
+ fi
45
+ idx=$(( $idx + 1 ))
46
+ done ;
47
+ ar=$( echo $ar | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
48
+ ag=$( echo $ag | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
49
+ ab=$( echo $ab | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
50
+ br=$( echo $br | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
51
+ bg=$( echo $bg | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
52
+ bb=$( echo $bb | awk -v n=$nr ' {printf "%.14g",$0/n;}' )
53
+ echo " 0, {$ar ,$ag ,$ab },{$br ,$bg ,$bb }}," >> $newfile
54
+ echo " " >> $newfile
0 commit comments