31 G_debug(3,
"I_cluster_means(nbands=%d,nclasses=%d)", C->nbands,
34 for (band = 0; band < C->nbands; band++) {
35 s = C->band_sum[band];
37 v = C->band_sum2[band] - s * m;
38 v = sqrt(v / (C->npoints - 1));
39 for (
class = 0;
class < C->nclasses;
class ++)
40 C->mean[band][class] = m;
42 for (
class = 0;
class < C->nclasses;
class ++)
43 C->mean[band][class] +=
44 ((2.0 *
class) / (C->nclasses - 1) - 1.0) * v;