Changeset 10786 for trunk/psModules/test/astrom/tap_pmAstrometryWCS.c
- Timestamp:
- Dec 16, 2006, 11:48:15 PM (20 years ago)
- Location:
- trunk/psModules/test/astrom
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
tap_pmAstrometryWCS.c (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psModules/test/astrom
- Property svn:ignore
-
old new 9 9 tap_pmAstrometryWCS_DVO2 10 10 tap_pmAstrometryWCS_DVO3 11 tap_pmAstrometryWCS_DVO
-
- Property svn:ignore
-
trunk/psModules/test/astrom/tap_pmAstrometryWCS.c
r10615 r10786 9 9 int main (void) 10 10 { 11 plan_tests(33); 11 12 12 pmModelGroupInit (); 13 14 plan_tests(240); 15 16 diag("pmAstrometryWCS tests"); 13 note("pmAstrometryWCS tests"); 17 14 18 15 { 19 diag("test pmAstromReadWCS");16 note("test pmAstromReadWCS"); 20 17 psMemId id = psMemGetId(); 21 18 … … 51 48 // toFPA carries pixel scale (pixels per micron) 52 49 // toTPA carries plate scale (microns per arcsecond) 53 bool status = pmAstromReadWCS (fpa, chip, header, 25.0*P S_RAD_DEG/3600.0, false);50 bool status = pmAstromReadWCS (fpa, chip, header, 25.0*PM_RAD_DEG/3600.0, false); 54 51 ok (status, "converted WCS keywords to FPA astrometry"); 55 52 skip_start (!status, 1, "*** WCS Conversion FAILS *** : skipping related tests"); … … 58 55 59 56 // make these tests double 60 ok_float(fpa->toSky->R*P S_DEG_RAD, 0.0, "projection center RA %f", fpa->toSky->R*PS_DEG_RAD);61 ok_float(fpa->toSky->R*P S_DEG_RAD, 0.0, "projection center DEC %f", fpa->toSky->R*PS_DEG_RAD);57 ok_float(fpa->toSky->R*PM_DEG_RAD, 0.0, "projection center RA %f", fpa->toSky->R*PM_DEG_RAD); 58 ok_float(fpa->toSky->R*PM_DEG_RAD, 0.0, "projection center DEC %f", fpa->toSky->R*PM_DEG_RAD); 62 59 63 ok_float(fpa->toSky->Xs, 25.0*P S_RAD_DEG/3600.0, "projection X scale %f", fpa->toSky->Xs);64 ok_float(fpa->toSky->Ys, 25.0*P S_RAD_DEG/3600.0, "projection X scale %f", fpa->toSky->Ys);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); 65 62 66 63 ok_float(fpa->toTPA->x->coeff[1][0][0][0], 0.04, "TP scale (mm per pixel): %f", fpa->toTPA->x->coeff[1][0][0][0]); … … 84 81 85 82 { 86 diag("test pmAstromReadWCS");83 note("test pmAstromReadWCS"); 87 84 psMemId id = psMemGetId(); 88 89 psPlaneTransform *wcsTrans = psPlaneTransformAlloc (1, 1);90 wcsTrans->x->coeff[0][0] = 0.0;91 wcsTrans->y->coeff[0][0] = 0.0;92 wcsTrans->x->coeff[1][0] = 1.0;93 wcsTrans->x->coeff[0][1] = 0.0;94 wcsTrans->y->coeff[1][0] = 0.0;95 wcsTrans->y->coeff[0][1] = 1.0;96 97 98 85 99 86 // construct a header with a simple set of WCS values … … 128 115 // toFPA carries pixel scale (pixels per micron) 129 116 // toTPA carries plate scale (microns per arcsecond) 130 bool status = pmAstromReadWCS (fpa, chip, header, 25.0*P S_RAD_DEG/3600.0, false);117 bool status = pmAstromReadWCS (fpa, chip, header, 25.0*PM_RAD_DEG/3600.0, false); 131 118 ok (status, "converted WCS keywords to FPA astrometry"); 132 119 skip_start (!status, 1, "*** WCS Conversion FAILS *** : skipping related tests"); … … 135 122 136 123 // make these tests double 137 ok_float(fpa->toSky->R*P S_DEG_RAD, 0.0, "projection center RA %f", fpa->toSky->R*PS_DEG_RAD);138 ok_float(fpa->toSky->R*P S_DEG_RAD, 0.0, "projection center DEC %f", fpa->toSky->R*PS_DEG_RAD);124 ok_float(fpa->toSky->R*PM_DEG_RAD, 0.0, "projection center RA %f", fpa->toSky->R*PM_DEG_RAD); 125 ok_float(fpa->toSky->R*PM_DEG_RAD, 0.0, "projection center DEC %f", fpa->toSky->R*PM_DEG_RAD); 139 126 140 ok_float(fpa->toSky->Xs, 25.0*P S_RAD_DEG/3600.0, "projection X scale %f", fpa->toSky->Xs);141 ok_float(fpa->toSky->Ys, 25.0*P S_RAD_DEG/3600.0, "projection X scale %f", fpa->toSky->Ys);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); 142 129 143 130 ok_float(fpa->toTPA->x->coeff[1][0][0][0], 0.04, "TP scale (mm per pixel): %f", fpa->toTPA->x->coeff[1][0][0][0]); … … 152 139 ok_float(chip->toFPA->y->coeff[0][1], 1.0, "CD2_2: %f", chip->toFPA->y->coeff[0][1]); 153 140 154 // apply both systems to real data:155 156 157 141 psFree (fpa); 158 142 psFree (chip);
Note:
See TracChangeset
for help on using the changeset viewer.
