lsst.afw.display Namespace Reference

Namespaces
	ds9
	utils

Functions
template<typename ImageT >
void	writeBasicFits (int fd, ImageT const &data, image::Wcs const *Wcs, char const *title)
template<typename ImageT >
void	writeBasicFits (std::string const &filename, ImageT const &data, image::Wcs const *Wcs, char const *title)

Function Documentation

template<typename ImageT >

void lsst::afw::display::writeBasicFits	(	int	fd,
		ImageT const &	data,
		image::Wcs const *	Wcs,
		char const *	title
	)

Definition at line 361 of file simpleFits.cc.

                     {
    /*
     * Allocate cards for FITS headers
     */
    std::list<Card> cards;
    /*
     * What sort if image is it?
     */
    int bitpix = lsst::afw::fits::getBitPix<typename ImageT::Pixel>();
    if (bitpix == 20) {                 // cfitsio for "Unsigned short"
        cards.push_back(Card("BZERO",  32768.0, ""));
        cards.push_back(Card("BSCALE", 1.0,     ""));
        bitpix = 16;
    } else if (bitpix == 0) {
        throw LSST_EXCEPT(lsst::pex::exceptions::RuntimeError, "Unsupported image type");
    }
    /*
     * Generate WcsA, pixel coordinates, allowing for X0 and Y0
     */
    std::string wcsName = "A";
    cards.push_back(Card(str(boost::format("CRVAL1%s") % wcsName),
                         data.getX0(), "(output) Column pixel of Reference Pixel"));
    cards.push_back(Card(str(boost::format("CRVAL2%s") % wcsName),
                         data.getY0(), "(output) Row pixel of Reference Pixel"));
    cards.push_back(Card(str(boost::format("CRPIX1%s") % wcsName), 1.0,
                         "Column Pixel Coordinate of Reference"));
    cards.push_back(Card(str(boost::format("CRPIX2%s") % wcsName), 1.0, "Row Pixel Coordinate of Reference"));
    cards.push_back(Card(str(boost::format("CTYPE1%s") % wcsName), "LINEAR", "Type of projection"));
    cards.push_back(Card(str(boost::format("CTYPE1%s") % wcsName), "LINEAR", "Type of projection"));
    cards.push_back(Card(str(boost::format("CUNIT1%s") % wcsName), "PIXEL", "Column unit"));
    cards.push_back(Card(str(boost::format("CUNIT2%s") % wcsName), "PIXEL", "Row unit"));
    /*
     * Now WcsB, so that pixel (0,0) is correctly labelled (but ignoring XY0)
     */
    wcsName = "B";
    cards.push_back(Card(str(boost::format("CRVAL1%s") % wcsName), 0,
                         "(output) Column pixel of Reference Pixel"));
    cards.push_back(Card(str(boost::format("CRVAL2%s") % wcsName), 0,
                         "(output) Row pixel of Reference Pixel"));
    cards.push_back(Card(str(boost::format("CRPIX1%s") % wcsName), 1.0,
                         "Column Pixel Coordinate of Reference"));
    cards.push_back(Card(str(boost::format("CRPIX2%s") % wcsName), 1.0, "Row Pixel Coordinate of Reference"));
    cards.push_back(Card(str(boost::format("CTYPE1%s") % wcsName), "LINEAR", "Type of projection"));
    cards.push_back(Card(str(boost::format("CTYPE1%s") % wcsName), "LINEAR", "Type of projection"));
    cards.push_back(Card(str(boost::format("CUNIT1%s") % wcsName), "PIXEL", "Column unit"));
    cards.push_back(Card(str(boost::format("CUNIT2%s") % wcsName), "PIXEL", "Row unit"));

    if (title) {
        cards.push_back(Card("OBJECT", title, "Image being displayed"));
    }
    /*
     * Was there something else?
     */
    if (Wcs != NULL) {
        typedef std::vector<std::string> NameList;

        image::Wcs::Ptr newWcs = Wcs->clone(); //Create a copy
        newWcs->shiftReferencePixel(-data.getX0(), -data.getY0());
        
        lsst::daf::base::PropertySet::Ptr metadata = newWcs->getFitsMetadata();

        NameList paramNames = metadata->paramNames();
        
        for (NameList::const_iterator i = paramNames.begin(), end = paramNames.end(); i != end; ++i) {
            if (*i == "SIMPLE" ||
                *i == "BITPIX" ||
                *i == "NAXIS" ||
                *i == "NAXIS1" ||
                *i == "NAXIS2" ||
                *i == "XTENSION" ||
                *i == "PCOUNT" ||
                *i == "GCOUNT"
               ) {
                continue;
            }
            std::type_info const &type = metadata->typeOf(*i);
            if (type == typeid(bool)) {
                cards.push_back(Card(*i, metadata->get<bool>(*i)));
            } else if (type == typeid(int)) {
                cards.push_back(Card(*i, metadata->get<int>(*i)));
            } else if (type == typeid(float)) {
                cards.push_back(Card(*i, metadata->get<float>(*i)));
            } else if (type == typeid(double)) {
                cards.push_back(Card(*i, metadata->get<double>(*i)));
            } else {
                cards.push_back(Card(*i, metadata->get<std::string>(*i)));
            }
        }
    }
    /*
     * Basic FITS stuff
     */
    const int naxis = 2;                // == NAXIS
    int naxes[naxis];                   /* values of NAXIS1 etc */
    naxes[0] = data.getWidth();
    naxes[1] = data.getHeight();
    
    write_fits_hdr(fd, bitpix, naxis, naxes, cards, 1);
    for (int y = 0; y != data.getHeight(); ++y) {
        if (write_fits_data(fd, bitpix, (char *)(data.row_begin(y)), (char *)(data.row_end(y))) < 0){
            throw LSST_EXCEPT(lsst::pex::exceptions::RuntimeError,
                              (boost::format("Error writing data for row %d") % y).str());
        }
    }

    pad_to_fits_record(fd, data.getWidth()*data.getHeight(), bitpix);
}   

template<typename ImageT >

void lsst::afw::display::writeBasicFits	(	std::string const &	filename,
		ImageT const &	data,
		image::Wcs const *	Wcs,
		char const *	title
	)

Definition at line 476 of file simpleFits.cc.

                     {
    int fd;
    if ((filename.c_str())[0] == '|') {         // a command
        const char *cmd = filename.c_str() + 1;
        while (isspace(*cmd)) {
            cmd++;
        }

        fd = fileno(popen(cmd, "w"));
    } else {
        fd = creat(filename.c_str(), 777);
    }

    if (fd < 0) {
        throw LSST_EXCEPT(lsst::pex::exceptions::RuntimeError,
                          (boost::format("Cannot open \"%s\"") % filename).str());
    }

    try {
        writeBasicFits(fd, data, Wcs, title);
    } catch(lsst::pex::exceptions::Exception &) {
        (void)close(fd);
        throw;
    }

    (void)close(fd);
}