Changeset 10875 for trunk/psModules/test/astrom/tap_pmAstrometryWCS.c
- Timestamp:
- Jan 1, 2007, 11:06:12 AM (20 years ago)
- File:
-
- 1 edited
-
trunk/psModules/test/astrom/tap_pmAstrometryWCS.c (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/test/astrom/tap_pmAstrometryWCS.c
r10824 r10875 46 46 pmChip *chip = pmChipAlloc (NULL, "test"); 47 47 48 // toFPA carries pixel scale ( pixels per micron)49 // to TPA carries plate scale (microns per arcsecond)50 bool status = pmAstromReadWCS (fpa, chip, header, 25.0*PM_RAD_DEG/3600.0);48 // toFPA carries pixel scale (microns per pixel) 49 // toSky carries plate scale (radians per micron) 50 bool status = pmAstromReadWCS (fpa, chip, header, 10.0); 51 51 ok (status, "converted WCS keywords to FPA astrometry"); 52 52 skip_start (!status, 1, "*** WCS Conversion FAILS *** : skipping related tests"); … … 58 58 ok_float(fpa->toSky->R*PM_DEG_RAD, 0.0, "projection center DEC %f", fpa->toSky->R*PM_DEG_RAD); 59 59 60 ok_float(fpa->toSky->Xs, 25.0*PM_RAD_DEG/3600.0, "projection X scale %f", fpa->toSky->Xs);61 ok_float(fpa->toSky->Ys, 25.0*PM_RAD_DEG/3600.0, "projection X scale %f", fpa->toSky->Ys);60 ok_float(fpa->toSky->Xs, PM_RAD_DEG/3600.0/10.0, "projection X scale %f", fpa->toSky->Xs); 61 ok_float(fpa->toSky->Ys, PM_RAD_DEG/3600.0/10.0, "projection X scale %f", fpa->toSky->Ys); 62 62 63 ok_float(fpa->toTPA->x->coeff[1][0], 0.04, "TP scale (mm per pixel): %f", fpa->toTPA->x->coeff[1][0]);64 ok_float(fpa->toTPA->y->coeff[0][1], 0.04, "TP scale (mm per pixel): %f", fpa->toTPA->x->coeff[1][0]);63 ok_float(fpa->toTPA->x->coeff[1][0], 1.0, "TP scale (unity): %f", fpa->toTPA->x->coeff[1][0]); 64 ok_float(fpa->toTPA->y->coeff[0][1], 1.0, "TP scale (unity): %f", fpa->toTPA->x->coeff[1][0]); 65 65 66 66 ok_float(chip->toFPA->x->coeff[0][0], 0.0, "ref pixel X: %f", chip->toFPA->x->coeff[0][0]); 67 67 ok_float(chip->toFPA->y->coeff[0][0], 0.0, "ref pixel Y: %f", chip->toFPA->y->coeff[0][0]); 68 68 69 ok_float(chip->toFPA->x->coeff[1][0], 1 .0, "CD1_1: %f", chip->toFPA->x->coeff[1][0]);69 ok_float(chip->toFPA->x->coeff[1][0], 10.0, "CD1_1: %f", chip->toFPA->x->coeff[1][0]); 70 70 ok_float(chip->toFPA->x->coeff[0][1], 0.0, "CD1_2: %f", chip->toFPA->x->coeff[0][1]); 71 71 ok_float(chip->toFPA->y->coeff[1][0], 0.0, "CD2_1: %f", chip->toFPA->y->coeff[1][0]); 72 ok_float(chip->toFPA->y->coeff[0][1], 1 .0, "CD2_2: %f", chip->toFPA->y->coeff[0][1]);72 ok_float(chip->toFPA->y->coeff[0][1], 10.0, "CD2_2: %f", chip->toFPA->y->coeff[0][1]); 73 73 74 74 psFree (fpa); … … 115 115 // toFPA carries pixel scale (pixels per micron) 116 116 // toTPA carries plate scale (microns per arcsecond) 117 bool status = pmAstromReadWCS (fpa, chip, header, 25.0*PM_RAD_DEG/3600.0);117 bool status = pmAstromReadWCS (fpa, chip, header, 10.0); 118 118 ok (status, "converted WCS keywords to FPA astrometry"); 119 119 skip_start (!status, 1, "*** WCS Conversion FAILS *** : skipping related tests"); … … 125 125 ok_float(fpa->toSky->R*PM_DEG_RAD, 0.0, "projection center DEC %f", fpa->toSky->R*PM_DEG_RAD); 126 126 127 ok_float(fpa->toSky->Xs, 25.0*PM_RAD_DEG/3600.0, "projection X scale %f", fpa->toSky->Xs);128 ok_float(fpa->toSky->Ys, 25.0*PM_RAD_DEG/3600.0, "projection X scale %f", fpa->toSky->Ys);127 ok_float(fpa->toSky->Xs, PM_RAD_DEG/3600.0/10.0, "projection X scale %f", fpa->toSky->Xs); 128 ok_float(fpa->toSky->Ys, PM_RAD_DEG/3600.0/10.0, "projection X scale %f", fpa->toSky->Ys); 129 129 130 ok_float(fpa->toTPA->x->coeff[1][0], 0.04, "TP scale (mm per pixel): %f", fpa->toTPA->x->coeff[1][0]);131 ok_float(fpa->toTPA->y->coeff[0][1], 0.04, "TP scale (mm per pixel): %f", fpa->toTPA->x->coeff[1][0]);130 ok_float(fpa->toTPA->x->coeff[1][0], 1.0, "TP scale (unity): %f", fpa->toTPA->x->coeff[1][0]); 131 ok_float(fpa->toTPA->y->coeff[0][1], 1.0, "TP scale (unity): %f", fpa->toTPA->x->coeff[1][0]); 132 132 133 ok_float(chip->toFPA->x->coeff[0][0], -10 .0, "ref pixel X: %f", chip->toFPA->x->coeff[0][0]);134 ok_float(chip->toFPA->y->coeff[0][0], -10 .0, "ref pixel Y: %f", chip->toFPA->y->coeff[0][0]);133 ok_float(chip->toFPA->x->coeff[0][0], -100.0, "ref pixel X: %f", chip->toFPA->x->coeff[0][0]); 134 ok_float(chip->toFPA->y->coeff[0][0], -100.0, "ref pixel Y: %f", chip->toFPA->y->coeff[0][0]); 135 135 136 ok_float(chip->toFPA->x->coeff[1][0], 1 .0, "CD1_1: %f", chip->toFPA->x->coeff[1][0]);136 ok_float(chip->toFPA->x->coeff[1][0], 10.0, "CD1_1: %f", chip->toFPA->x->coeff[1][0]); 137 137 ok_float(chip->toFPA->x->coeff[0][1], 0.0, "CD1_2: %f", chip->toFPA->x->coeff[0][1]); 138 138 ok_float(chip->toFPA->y->coeff[1][0], 0.0, "CD2_1: %f", chip->toFPA->y->coeff[1][0]); 139 ok_float(chip->toFPA->y->coeff[0][1], 1 .0, "CD2_2: %f", chip->toFPA->y->coeff[0][1]);139 ok_float(chip->toFPA->y->coeff[0][1], 10.0, "CD2_2: %f", chip->toFPA->y->coeff[0][1]); 140 140 141 141 psFree (fpa);
Note:
See TracChangeset
for help on using the changeset viewer.
