1 #include "global.h" 2 get_aver_elev(struct Ortho_Control_Points * cpz,double * aver_z)3int get_aver_elev(struct Ortho_Control_Points *cpz, double *aver_z) 4 { 5 double meanz; 6 double *cp = cpz->z2; 7 int i, n; 8 9 /* Need 1 control points */ 10 if (cpz->count <= 0) { 11 return (-1); 12 } 13 14 /* set average elevation from mean values of CONZ points */ 15 meanz = 0; 16 n = 0; 17 for (i = 0; i < cpz->count; i++) { 18 if (cpz->status[i] <= 0) 19 continue; 20 21 n++; 22 meanz += *(cp++); 23 G_debug(3, "In ortho meanz = %f", meanz); 24 } 25 26 *aver_z = meanz / n; 27 28 G_debug(1, "In ortho aver_z = %f", *aver_z); 29 30 return 0; 31 } 32