#include <sli/fitscc.h>
#include <sli/mdarray_statistics.h>
#include <stdlib.h>


using namespace sli;


static void subtract_sky(fits_image &hdu, double *sky_stddev) {
}


int main(int argc, char *argv[])
{
    if (argc != 3) {
        sli__eprintf("usage: %s INPUT OUTPUT\n", argv[0]);
        exit(1);
    }

    double sky_stddev = 0.0; // 警告抑止のため0.0を代入

    fitscc fits;
    fits.read_stream(argv[1]);
    subtract_sky(fits.image(0L), &sky_stddev);
    sli__eprintf("skylevel = 0.0 +/- %f\n", sky_stddev);
    fits.write_stream(argv[2]);
    return 0;
}