cfitsio.cpp

Go to the documentation of this file.

/**
 * @file
 * @ingroup daq_ocm_libfits
 * @copyright 2021 ESO - European Southern Observatory
 *
 * @brief Contains definitions for daq/fits/cfitsio.hpp
 */
#include <daq/fits/cfitsio.hpp>
 
#include <cassert>
 
#include <fmt/format.h>
 
namespace daq::fits {
 
CfitsioError::CfitsioError(int status, char const* message)
    : std::runtime_error(message), m_status(status) {
}
CfitsioError::CfitsioError(int status, std::string const& message)
    : std::runtime_error(message), m_status(status) {
}
 
void DefaultClose(fitsfile* ptr) noexcept {
    if (ptr == nullptr) {
        return;
    }
    int status = 0;
    fits_close_file(ptr, &status);
}
 
int CfitsioError::GetStatus() const noexcept {
    return m_status;
}
 
void InitPrimaryHduNoImage(fitsfile* ptr) {
    /*
    SIMPLE  =                    T / file does conform to FITS standard
    BITPIX  =                    8 / number of bits per data pixel
    NAXIS   =                    0 / number of data axes
    */
    assert(ptr);
    int status = 0;
    int simple = true;
    char bitpix = 8;
    char naxis = 0;
    int num_hdus = 0;
    (void)fits_get_num_hdus(ptr, &num_hdus, &status);
    if (num_hdus != 0) {
        throw CfitsioError(HEADER_NOT_EMPTY, "Cannot safely initialize an non-empty FITS file");
    }
    (void)fits_create_hdu(ptr, &status);
    (void)fits_write_key(ptr,
                         TLOGICAL,
                         "SIMPLE",
                         &simple,
                         "file does conform to FITS standard",
                         &status);
    (void)fits_write_key(ptr,
                         TBYTE,
                         "BITPIX",
                         &bitpix,
                         "number of bits per data pixel",
                         &status);
    (void)fits_write_key(ptr,
                         TBYTE,
                         "NAXIS",
                         &naxis,
                         "file provides primary HDU keywords only",
                         &status);
    if (status != 0) {
        char err_text[31]; // returns maximum 30 characters
        fits_get_errstatus(status, &err_text[0]);
        err_text[30] = '\0';
        throw CfitsioError(status, fmt::format("failed to initialize primary HDU: {}",
                                               err_text));
    }
}
 
UniqueFitsFile CreateEmpty(char const* filename) {
    fitsfile* ptr = nullptr;
    int status = 0;
    fits_create_file(&ptr, const_cast<char*>(filename), &status);
    if (status != 0 || ptr == nullptr) {
        throw CfitsioError(status, fmt::format("failed to create FITS file '{}'",
                                               filename));
    }
    return UniqueFitsFile(ptr, DefaultClose);
    /*, +[](fitsfile* p) {
        int status = 0;
        fits_close_file(p, &status);
    }*/
}
 
}  // namespace daq::fits