33 #include <casu_fits.h>
34 #include <casu_utils.h>
35 #include <casu_mods.h>
36 #include <casu_wcsutils.h>
39 cpl_propertylist *plist,*plist2;
41 double ra,dec,x,y,xi,eta,ra2,dec2,x2,y2;
46 cpl_test_init(PACKAGE_BUGREPORT,CPL_MSG_WARNING);
50 plist = cpl_propertylist_new();
51 cpl_propertylist_update_string(plist,
"CTYPE1",
"RA---ZPN");
52 cpl_propertylist_update_string(plist,
"CTYPE2",
"DEC--ZPN");
53 cpl_propertylist_update_double(plist,
"CRVAL1",81.0578935624986);
54 cpl_propertylist_update_double(plist,
"CRVAL2",-24.5408992630551);
55 cpl_propertylist_update_double(plist,
"CRPIX1",2159.65203390075);
56 cpl_propertylist_update_double(plist,
"CRPIX2",2187.28057518062);
57 cpl_propertylist_update_double(plist,
"CD1_1",-2.96100210581431E-05);
58 cpl_propertylist_update_double(plist,
"CD2_1",-3.54397484069851E-08);
59 cpl_propertylist_update_double(plist,
"CD1_2",-3.93573000440439E-08);
60 cpl_propertylist_update_double(plist,
"CD2_2",2.96102853649076E-05);
61 cpl_propertylist_update_double(plist,
"PV2_1",1.0);
62 cpl_propertylist_update_double(plist,
"PV2_3",400.0);
63 cpl_propertylist_update_int(plist,
"NAXIS",2);
64 cpl_propertylist_update_int(plist,
"NAXIS1",2048);
65 cpl_propertylist_update_int(plist,
"NAXIS2",2048);
66 wcs = cpl_wcs_new_from_propertylist(plist);
73 cpl_test_rel(ra,81.09569,0.001);
74 cpl_test_rel(dec,-24.57600,0.001);
76 cpl_test_rel(x,1000.0,0.00001);
77 cpl_test_rel(y,1000.0,0.00001);
79 cpl_test_rel(xi,0.000600115,0.000001);
80 cpl_test_rel(eta,-0.000612866,0.00001);
83 cpl_test_eq(status,CASU_OK);
84 cpl_test_rel(ra,81.06153,0.001);
85 cpl_test_rel(ra2,81.12822,0.001);
86 cpl_test_rel(dec,-24.60560,0.001);
87 cpl_test_rel(dec2,-24.54493,0.001);
88 plist2 = cpl_propertylist_duplicate(plist);
89 cpl_propertylist_update_double(plist2,
"CRVAL1",81.1578935624986);
90 cpl_propertylist_update_double(plist2,
"CRVAL2",-24.6408992630551);
91 wcs2 = cpl_wcs_new_from_propertylist(plist2);
93 cpl_test_rel(x2,4068.86290,0.001);
94 cpl_test_rel(y2,4382.93236,0.001);
100 cpl_propertylist_delete(plist);
101 cpl_propertylist_delete(plist2);
102 return(cpl_test_end(0));
void casu_radectoxieta(cpl_wcs *wcs, double ra, double dec, double *xi, double *eta)
void casu_radectoxy(cpl_wcs *wcs, double ra, double dec, double *x, double *y)
void casu_xytoradec(cpl_wcs *wcs, double x, double y, double *ra, double *dec)
void casu_xytoxy_list(cpl_wcs *wcs1, cpl_wcs *wcs2, int nc, double *x_1, double *y_1, double *x_2, double *y_2)
int casu_coverage(cpl_propertylist *plist, int fudge, double *ra1, double *ra2, double *dec1, double *dec2, int *status)