| OBJECTID | Id | depth | Elevation | Year | RASTERVALU | method | |
1 | 1 | 0 | 11.6 | -11.6 | 1984 | -11.6 | a |
2 | 2 | 0 | 12.2 | -12.2 | 1984 | -12.2001 | a |
3 | 3 | 0 | 12.4 | -12.4 | 1984 | -12.4003 | a |
4 | 4 | 0 | 13 | -13 | 1984 | -13.0005 | a |
5 | 5 | 0 | 13.8 | -13.8 | 1984 | -13.7989 | a |
6 | 6 | 0 | 13.4 | -13.4 | 1984 | -13.3996 | a |
7 | 7 | 0 | 12 | -12 | 1984 | -12.0007 | a |
8 | 8 | 0 | 12 | -12 | 1984 | -12.0012 | a |
9 | 9 | 0 | 12.6 | -12.6 | 1984 | -12.6004 | a |
10 | 10 | 0 | 11.4 | -11.4 | 1984 | -11.3999 | a |
11 | 11 | 0 | 11.6 | -11.6 | 1984 | -11.6006 | a |
12 | 12 | 0 | 12.6 | -12.6 | 1984 | -12.5993 | a |
1281 | 128 | 0 | 5 | -5 | 1984 | -4.74908 | b |
1291 | 129 | 0 | 5 | -5 | 1984 | -5.02803 | b |
1301 | 130 | 0 | 5 | -5 | 1984 | -4.96928 | b |
1311 | 131 | 0 | 5 | -5 | 1984 | -4.81961 | b |
1321 | 132 | 0 | 5 | -5 | 1984 | -4.55515 | b |
1331 | 133 | 0 | 5 | -5 | 1984 | -5.07164 | b |
1341 | 134 | 0 | 5 | -5 | 1984 | -5.05858 | b |
1351 | 135 | 0 | 5 | -5 | 1984 | -5.05777 | b |
1361 | 136 | 0 | 5 | -5 | 1984 | -4.98249 | b |
1371 | 137 | 0 | 0 | 0 | 1984 | -0.24859 | b |
1381 | 138 | 0 | 0 | 0 | 1984 | -0.30513 | b |
1391 | 139 | 0 | 0 | 0 | 1984 | -0.34433 | b |
1401 | 140 | 0 | 0 | 0 | 1984 | -0.40141 | b |
1411 | 141 | 0 | 0 | 0 | 1984 | -0.2499 | b |
1421 | 142 | 0 | 0 | 0 | 1984 | -0.3198 | b |
zwc <-function(Elevation, RASTERVALU){
1/length(Elevation)*sum((Elevation-RASTERVALU)^2)
}
如上所示,有这么一些数据,我写了这么一个函数,我想自动根据method列计算a,b的zwc值,求教大神如何应用这个公式呢?
我想的是应该是按照上面那个流程图那样通过判断实现吧?但是不只是ab两种方法啊,好多种,就是这样示意一下,不会写程序啊,只会简单的函数,求大神教啊


雷达卡




京公网安备 11010802022788号







