Soeren Gebbert
3a539ed71b
Fix: return the allocated buffer
|
15 years ago | |
|---|---|---|
| .. | ||
| Makefile | 15 years ago | |
| README | 17 years ago | |
| alloc.c | 17 years ago | |
| eol.c | 17 years ago | |
| find.c | 17 years ago | |
| fopen.c | 17 years ago | |
| georef.c | 17 years ago | |
| group.c | 15 years ago | |
| imagerylib.dox | 16 years ago | |
| list_gp.c | 16 years ago | |
| list_subgp.c | 16 years ago | |
| loc_info.c | 15 years ago | |
| points.c | 15 years ago | |
| ref.c | 17 years ago | |
| sig.c | 17 years ago | |
| sigfile.c | 16 years ago | |
| sigset.c | 17 years ago | |
| sigsetfile.c | 17 years ago | |
| target.c | 17 years ago | |
| title.c | 15 years ago | |
| var.c | 17 years ago | |
README
To merge 2 signatures
n1 number of points in sig 1
n2 number of points in sig 2
mean1[nbands] means per band for sig 1
mean2[nbands] means per band for sig 2
var1[b1][b2] covariance band 1 with band 2 for sig 1
var2[b1][b2] covariance band 1 with band 2 for sig 2
the meger is
n = n1+n2
mean[b] = (mean1[b]*n1 + mean2[b]*n2)/n
sum1 = var1[b1][b2] * (n1-1) + n1 * mean[b1] * mean[b2];
sum2 = var2[b1][b2] * (n2-1) + n2 * mean[b1] * mean[b2];
var[b1][b2] = (sum1+sum2 - n*mean[b1]*mean[b2) / (n-1)
n1 number of points in sig 1
n2 number of points in sig 2
mean1[nbands] means per band for sig 1
mean2[nbands] means per band for sig 2
var1[b1][b2] covariance band 1 with band 2 for sig 1
var2[b1][b2] covariance band 1 with band 2 for sig 2
the meger is
n = n1+n2
mean[b] = (mean1[b]*n1 + mean2[b]*n2)/n
sum1 = var1[b1][b2] * (n1-1) + n1 * mean[b1] * mean[b2];
sum2 = var2[b1][b2] * (n2-1) + n2 * mean[b1] * mean[b2];
var[b1][b2] = (sum1+sum2 - n*mean[b1]*mean[b2) / (n-1)