END_TEST
START_TEST(test_old_pixel_to_metric)
{
/* input */
double xc = 0.0; // [mm]
double yc = 0.0; // [mm]
int imx = 1024; // standard image size
int imy = 1008;
double pix_x = 0.010; // 10 micron pixel size
double pix_y = 0.010;
int field = 0; // simple image format, not interlaced
/* output */
double xp, yp;
/* compare the xc, yc to the original */
double xc1, yc1;
old_metric_to_pixel (&xp, &yp, xc, yc, imx, imy, pix_x, pix_y, field);
old_pixel_to_metric (&xc1, &yc1,xp, yp, imx, imy, pix_x, pix_y, field);
ck_assert_msg( fabs(xc1 - xc) < EPS &&
fabs(yc1 - yc) < EPS,
"Expected %f, %f but got %f %f\n", xc, yc, xc1, yc1);
xc = 1.0;
yc = 0.0;
old_metric_to_pixel (&xp, &yp, xc, yc, imx, imy, pix_x, pix_y, field);
old_pixel_to_metric (&xc1, &yc1,xp, yp, imx, imy, pix_x, pix_y, field);
ck_assert_msg( fabs(xc1 - xc) < EPS &&
fabs(yc1 - yc) < EPS,
"Expected %f, %f, but got %f %f\n", xc,yc,xc1, yc1);
xc = 0.0;
yc = -1.0;
old_metric_to_pixel (&xp, &yp, xc, yc, imx, imy, pix_x, pix_y, field);
old_pixel_to_metric (&xc1, &yc1,xp, yp, imx, imy, pix_x, pix_y, field);
ck_assert_msg( fabs(xc1 - xc) < EPS &&
fabs(yc1 - yc) < EPS,
"Expected %f, %f, but got %f %f\n", xc,yc,xc1, yc1);
}
/* wrap metric_to_pixel */
void metric_to_pixel(double * x_pixel
, double * y_pixel
, double x_metric
, double y_metric
, control_par* parameters
){
old_metric_to_pixel(x_pixel
, y_pixel
, x_metric
, y_metric
, parameters->imx
, parameters->imy
, parameters->pix_x
, parameters->pix_y
, parameters->chfield);
}
