/* # This file is part of the Astrometry.net suite. # Licensed under a 3-clause BSD style license - see LICENSE */ #include #include #include #include #include "cutest.h" #include "anwcs.h" #include "sip.h" #include "sip_qfits.h" #include "errors.h" #include "ioutils.h" #include "log.h" struct walk_token2 { dl* radecs; }; void test_mercator_2(CuTest* tc) { double ra,dec; double ra2,dec2; anwcs_t* wcs = anwcs_create_mercator_2(180., 0., 128.5, 128.5, 1., 256, 256, TRUE); printf("RA,Dec in corners:\n"); anwcs_pixelxy2radec(wcs, 1., 1., &ra, &dec); CuAssertDblEquals(tc, 85.0, dec, 0.1); CuAssertDblEquals(tc, 360.0, ra, 1.0); printf(" 1,1 -> %.1f, %.1f\n", ra, dec); anwcs_pixelxy2radec(wcs, 256., 256., &ra, &dec); printf(" W,H -> %.1f, %.1f\n", ra, dec); CuAssertDblEquals(tc, -85.0, dec, 0.1); CuAssertDblEquals(tc, 0.0, ra, 1.0); printf("Delta near the center:\n"); anwcs_pixelxy2radec(wcs, 128., 128., &ra, &dec); anwcs_pixelxy2radec(wcs, 129., 128., &ra2, &dec2); printf("dx = 1 -> dRA, dDec %.1f, %.1f\n", ra2-ra, dec2-dec); anwcs_pixelxy2radec(wcs, 128., 129., &ra2, &dec2); printf("dy = 1 -> dRA, dDec %.1f, %.1f\n", ra2-ra, dec2-dec); } /* static void walk_callback2(const anwcs_t* wcs, double ix, double iy, double ra, double dec, void* token) { struct walk_token2* walk = token; dl_append(walk->radecs, ra); dl_append(walk->radecs, dec); } */ void test_walk_outline(CuTest* tc) { anwcs_t* plotwcs = anwcs_create_allsky_hammer_aitoff(180., 0., 400, 200); #if 0 tan_t tanwcs; tanwcs.crval[0] = 10.; tanwcs.crval[1] = 0.; tanwcs.crpix[0] = 25.; tanwcs.crpix[1] = 25.; tanwcs.cd[0][0] = 1.; tanwcs.cd[0][1] = 0.; tanwcs.cd[1][0] = 0.; tanwcs.cd[1][1] = 1.; tanwcs.imagew = 50.; tanwcs.imageh = 50.; tanwcs.sin = 0; #endif dl* rd; /* anwcs_t* wcs = anwcs_new_tan(&tanwcs); struct walk_token2 token2; double stepsize = 1000; */ log_init(LOG_ALL); /* token2.radecs = dl_new(256); anwcs_walk_image_boundary(wcs, stepsize, walk_callback2, &token2); rd = token2.radecs; logverb("Outline: walked in %i steps\n", dl_size(rd)/2); // close dl_append(rd, dl_get(rd, 0)); dl_append(rd, dl_get(rd, 1)); int i; for (i=0; i crpix. x = tan->crpix[0]; y = tan->crpix[1]; ra = 1; ra2 = 2; dec = 3; dec2 = 4; tan_pixelxy2radec(tan, x, y, &ra, &dec); CuAssertDblEquals(tc, tan->crval[0], ra, 1e-6); CuAssertDblEquals(tc, tan->crval[1], dec, 1e-6); ok = anwcs_pixelxy2radec(anwcs, x, y, &ra2, &dec2); CuAssertIntEquals(tc, 0, ok); CuAssertDblEquals(tc, tan->crval[0], ra2, 1e-6); CuAssertDblEquals(tc, tan->crval[1], dec2, 1e-6); ra = tan->crval[0]; dec = tan->crval[1]; x = 1; x2 = 2; y = 3; y2 = 4; ok = tan_radec2pixelxy(tan, ra, dec, &x, &y); CuAssertIntEquals(tc, TRUE, ok); CuAssertDblEquals(tc, tan->crpix[0], x, 1e-6); CuAssertDblEquals(tc, tan->crpix[1], y, 1e-6); ok2 = anwcs_radec2pixelxy(anwcs, ra, dec, &x2, &y2); CuAssertIntEquals(tc, 0, ok2); CuAssertDblEquals(tc, tan->crpix[0], x2, 1e-6); CuAssertDblEquals(tc, tan->crpix[1], y2, 1e-6); // check pixel (0,0). x = y = 0.0; ra = 1; ra2 = 2; dec = 3; dec2 = 4; tan_pixelxy2radec(tan, x, y, &ra, &dec); ok = anwcs_pixelxy2radec(anwcs, x, y, &ra2, &dec2); CuAssertIntEquals(tc, 0, ok); CuAssertDblEquals(tc, ra, ra2, 1e-6); CuAssertDblEquals(tc, dec, dec2, 1e-6); // check RA,Dec (85, -4) ra = 85; dec = -4; x = 1; x2 = 2; y = 3; y2 = 4; ok = tan_radec2pixelxy(tan, ra, dec, &x, &y); CuAssertIntEquals(tc, TRUE, ok); ok2 = anwcs_radec2pixelxy(anwcs, ra, dec, &x2, &y2); CuAssertIntEquals(tc, 0, ok2); printf("x,y (%g,%g) vs (%g,%g)\n", x, y, x2, y2); CuAssertDblEquals(tc, x, x2, 1e-6); CuAssertDblEquals(tc, y, y2, 1e-6); anwcs_free(anwcs); } free(tan); }