|
|
@@ -57,12 +57,12 @@ int main(int argc, char *argv[]) |
|
|
|
|
|
|
|
/* Compute the standard error */ |
|
|
|
tmp = pipi_gaussian_blur(dither, sigma); |
|
|
|
e0 = pipi_measure_rmsd(gauss, tmp); |
|
|
|
e0 = pipi_measure_msd(gauss, tmp); |
|
|
|
pipi_free(tmp); |
|
|
|
|
|
|
|
/* Compute the fast error */ |
|
|
|
tmp = pipi_gaussian_blur_ext(dither, sigma, sigma, 0.16, 0.26); |
|
|
|
e1 = pipi_measure_rmsd(gauss, tmp); |
|
|
|
e1 = pipi_measure_msd(gauss, tmp); |
|
|
|
pipi_free(tmp); |
|
|
|
|
|
|
|
/* Compute displacement */ |
|
|
@@ -74,7 +74,7 @@ int main(int argc, char *argv[]) |
|
|
|
tmp = pipi_gaussian_blur_ext(dither, sigma, sigma, |
|
|
|
fx + step * dx / Z, |
|
|
|
fy + step * dy / Z); |
|
|
|
e = pipi_measure_rmsd(gauss, tmp); |
|
|
|
e = pipi_measure_msd(gauss, tmp); |
|
|
|
pipi_free(tmp); |
|
|
|
if(e < best) |
|
|
|
{ |
|
|
|