#include <stdio.h>
#include "pslib.h"

int main(void)
{
    psFits *fits = psFitsOpen("f0230_c01_606365_12k_i.fits", "r");
    psRegion region = {0, 0, 0, 0};
    psImage *image = psFitsReadImage(NULL, fits, region, 0);
    psFitsClose(fits);
    int numNaN = psImageClipNaN(image, 0.0);
    printf("Clipped %d NaN pixels.\n", numNaN);

    psStats *stats = psStatsAlloc(PS_STAT_ROBUST_MEDIAN | PS_STAT_ROBUST_STDEV);
    stats->clipSigma = 2;
    stats->clipIter = 5;
    stats = psImageStats(stats, image, NULL, 0);

    printf("%lf %lf\n", stats->robustMedian, stats->robustStdev);
    psFree(image);
    psFree(stats);
}