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

using namespace sli;


static const int amp_n = 4; // アンプの数


static void debias(fits_image &hdu)
{
    for (int amp = 1;  amp <= amp_n;  amp++) {
        const int os_min_x = hdu.headerf("S_OSMN%d1", amp).lvalue(),
                  os_max_x = hdu.headerf("S_OSMX%d1", amp).lvalue(),
                  ef_min_x = hdu.headerf("S_EFMN%d1", amp).lvalue(),
                  ef_max_x = hdu.headerf("S_EFMX%d1", amp).lvalue();
        sli__eprintf(
            "amp = %d    os_min_x = % 5d    os_max_x = % 5d    ef_min_x = % 5d    ef_max_x = % 5d\n",
            amp, os_min_x, os_max_x, ef_min_x, ef_max_x);
    }
}


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

    fitscc fits;
    fits.read_stream(argv[1]);
    debias(fits.image(0L));
    fits.write_stream(argv[2]);

    return 0;
}