如何在Shell下实现这样的数值统计计算?

Osirix

原始数据格式如下,有很多个SECTION,共n个;
表中每个SECTION都有两列,(X,Y)可以看成各个点的坐标,每个点自成一行.
现在要对各个点在所有SECTION中的坐标作统计平均.
另外,每个SECTION之间都由点的个数(现等于9)和"SECTION"标记隔开.

1. 
   
2. ---------------------------------------------
   
3.      X         Y
   
4. ---------------------------------------------
   
5. ....
   
6. SECTION: 1             #<-这是第k行,从这里开始读取#
   
7. 1.848619   1.582796    #<-想把它俩分别装入数组XA[1]和YA[1],并在每个SECTION中累加
   
8. 20.17222   20.90597    #<-想把它俩分别装入数组XA[2]和YA[2],并在每个SECTION中累加
   
9. 300.1233   300.8589    #<-想把它俩分别装入数组XA[3]和YA[3],并在每个SECTION中累加
   
10. 0.081323   0.962245    #<-想把它俩分别装入数组XA[4]和YA[4],并在每个SECTION中累加
    
11. -10.0806   10.46079    #<-想把它俩分别装入数组XA[5]和YA[5],并在每个SECTION中累加
    
12. 20.15859   20.55888    #<-想把它俩分别装入数组XA[6]和YA[6],并在每个SECTION中累加
    
13. 30.04000   30.39677    #<-想把它俩分别装入数组XA[7]和YA[7],并在每个SECTION中累加
    
14. -40.2340   40.84182    #<-想把它俩分别装入数组XA[8]和YA[8],并在每个SECTION中累加
    
15. -50.0834   50.96503    #<-想把它俩分别装入数组XA[9]和YA[9],并在每个SECTION中累加
    
16. 9                      #<-表示总共有9个点
    
17. SECTION: 2              
    
18. 1.039361   1.083161   
    
19. 2.062854   2.139038
    
20. -3.00659   3.216112
    
21. 0.243749   0.543359
    
22. 10.03200   10.73950
    
23. 20.44733   20.76969
    
24. 30.41897   30.42323
    
25. 40.04939   40.24543
    
26. 50.28102   50.28565
    
27. 9
    
28. SECTION: 3
    
29. 1.501221   -1.17865
    
30. 2.591375   2.032659
    
31. 3.547986   -3.01310
    
32. 0.623798   0.667415
    
33. 10.50620   10.50388
    
34. 20.43490   20.72506
    
35. 30.80939   30.81383
    
36. 40.63154   40.05711
    
37. 50.42599   50.92859
    
38.   9
    
39. SECTION: 4
    
40. .... ........ ......
    
41. .... ........ ......
    
42.   9
    
43. SECTION: n
    
44. .... ........ ......
    
45. .... ........ ......
    
46. ------------------------------------------
    

复制代码

只晓得用sed可以控制处理哪些行,但现在似乎对行号也要设一个变量进行计数,而且还要同时进行循环累加,就不晓得该怎么处理了?

哪位大侠可否给个方子?

many 3ks!

mingfal

这种问题用C多方便呢？
ps:我也不知道怎么用shell来实现:-(

HuiCai

利用awk实现：
awk -f f.awk YourFile，其中f.awk内容为：

1. BEGIN {
   
2.     sumx=0.0
   
3.     sumy=0.0
   
4. }
   
5. {
   
6. if ( $1 ~ "SECTION:" ) {
   
7.     sec=$2
   
8.     p=1
   
9.     sumx=0.0
   
10.     sumy=0.0
    
11.     } else if ( NF == 2 ) {
    
12.         sumx=sumx + $1
    
13.         sumy=sumy + $2
    
14.     } else if ( NF == 1 && p == 1 ) {
    
15.         print "SECTION",sec,": ", sumx/$1 , sumy/$1
    
16.     }
    
17. }

复制代码

这里#<-及其后面的内容删除掉才成