15#include <grass/cluster.h>
17#define FAR ((double) -1.0)
37 if (C->count[class1] < 2)
39 if (C->count[class2] < 2)
41 n1 = (double)C->count[class1];
42 n2 = (
double)C->count[class2];
46 for (band = 0; band < C->nbands; band++) {
47 s1 = C->sum[band][class1];
48 s2 = C->sum[band][class2];
56 var = C->sum2[band][class1] - (s1 * m1);
61 var = C->sum2[band][class2] - (s2 * m2);
72 a1 = sqrt(6 * d / a1);
74 a2 = sqrt(6 * d / a2);
double I_cluster_separation(struct Cluster *C, int class1, int class2)
?