34#include <casu_utils.h>
39 cpl_image *im,*outim,*outimv;
40 casu_fits *infits[10];
41 unsigned char *rejmask,*rejplus;
42 cpl_propertylist *drs,*phu,*ehu;
48 cpl_test_init(PACKAGE_BUGREPORT,CPL_MSG_WARNING);
53 retval =
casu_imcombine(NULL,NULL,100,2,1,0,100.0,
"EXPTIME",&outim,
54 &outimv,&rejmask,&rejplus,&drs,&status);
55 cpl_test_eq(status,CASU_FATAL);
56 cpl_test_eq(status,retval);
58 cpl_test_null(outimv);
59 cpl_test_null(rejmask);
60 cpl_test_null(rejplus);
66 for (i = 0; i < 10; i++) {
67 im = cpl_image_new(100,100,CPL_TYPE_FLOAT);
75 cpl_image_fill_noise_uniform(im,val1,val2);
76 phu = cpl_propertylist_new();
77 ehu = cpl_propertylist_new();
78 cpl_propertylist_update_float(phu,
"EXPTIME",ave);
80 cpl_propertylist_delete(phu);
81 cpl_propertylist_delete(ehu);
87 retval =
casu_imcombine(infits,NULL,10,2,1,0,100.0,
"EXPTIME",&outim,
88 &outimv,&rejmask,&rejplus,&drs,&status);
89 cpl_test_eq(status,CASU_OK);
90 cpl_test_eq(status,retval);
91 cpl_test_nonnull(outim);
92 cpl_test_null(outimv);
93 cpl_test_nonnull(rejmask);
94 cpl_test_nonnull(rejplus);
95 cpl_test_nonnull(drs);
96 cpl_test_rel(cpl_image_get_mean(outim),150.0,0.001);
97 cpl_image_delete(outim);
100 cpl_propertylist_delete(drs);
105 retval =
casu_imcombine(infits,NULL,10,2,2,0,100.0,
"EXPTIME",&outim,
106 &outimv,&rejmask,&rejplus,&drs,&status);
107 cpl_test_eq(status,CASU_OK);
108 cpl_test_eq(status,retval);
109 cpl_test_nonnull(outim);
110 cpl_test_null(outimv);
111 cpl_test_nonnull(rejmask);
112 cpl_test_nonnull(rejplus);
113 cpl_test_nonnull(drs);
114 cpl_test_rel(cpl_image_get_mean(outim),150.0,0.001);
115 cpl_image_delete(outim);
118 cpl_propertylist_delete(drs);
123 retval =
casu_imcombine(infits,NULL,10,2,3,0,100.0,
"EXPTIME",&outim,
124 &outimv,&rejmask,&rejplus,&drs,&status);
125 cpl_test_eq(status,CASU_OK);
126 cpl_test_eq(status,retval);
127 cpl_test_nonnull(outim);
128 cpl_test_null(outimv);
129 cpl_test_nonnull(rejmask);
130 cpl_test_nonnull(rejplus);
131 cpl_test_nonnull(drs);
132 cpl_test_rel(cpl_image_get_mean(outim),100.0,0.001);
133 cpl_image_delete(outim);
136 cpl_propertylist_delete(drs);
140 for (i = 0; i < 10; i++)
142 return(cpl_test_end(0));
casu_fits * casu_fits_wrap(cpl_image *im, casu_fits *model, cpl_propertylist *phu, cpl_propertylist *ehu)
void casu_fits_delete(casu_fits *p)
int casu_imcombine(casu_fits **fset, casu_fits **fsetv, int nfits, int combtype, int scaletype, int xrej, float thresh, const char *expkey, cpl_image **outimage, cpl_image **outvimage, unsigned char **rejmask, unsigned char **rejplus, cpl_propertylist **drs, int *status)
Stack images into a mean or median image with rejection.