Observatory.cc

Go to the documentation of this file.

// -*- lsst-c++ -*-

/* 
 * LSST Data Management System
 * Copyright 2008, 2009, 2010 LSST Corporation.
 * 
 * This product includes software developed by the
 * LSST Project (http://www.lsst.org/).
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * 
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * 
 * You should have received a copy of the LSST License Statement and 
 * the GNU General Public License along with this program.  If not, 
 * see <http://www.lsstcorp.org/LegalNotices/>.
 */
 
#include <sstream>
#include <cmath>

#include "lsst/pex/exceptions.h"
#include "boost/format.hpp"
#include "boost/algorithm/string.hpp"
#include "boost/tuple/tuple.hpp"

#include "lsst/afw/geom/Angle.h"
#include "lsst/afw/coord/Coord.h"
#include "lsst/afw/coord/Observatory.h"

namespace coord        = lsst::afw::coord;
namespace ex           = lsst::pex::exceptions;
namespace afwGeom      = lsst::afw::geom;



coord::Observatory::Observatory(
                                afwGeom::Angle const longitude, 
                                afwGeom::Angle const latitude,  
                                double const elevation  
                               ) :
    _latitude(latitude),
    _longitude(longitude),
    _elevation(elevation) {
}



/*
 * @brief Constructor for the observatory with lat/long as strings
 *
 * @note RA is assumed to be in DMS, not HMS!
 *
 */
coord::Observatory::Observatory(
                                std::string const longitude, 
                                std::string const latitude,  
                                double const elevation       
                               ) : 
    _latitude(dmsStringToAngle(latitude)),
    _longitude(dmsStringToAngle(longitude)),
    _elevation(elevation) {
}




afwGeom::Angle coord::Observatory::getLongitude() const {
    return _longitude;
}

afwGeom::Angle coord::Observatory::getLatitude() const {
    return _latitude;
}


void coord::Observatory::setLatitude(
                 afwGeom::Angle const latitude 
                )   {
    _latitude = latitude;
}

void coord::Observatory::setLongitude(
                                      afwGeom::Angle const longitude 
                                     ) {
    _longitude = longitude;
}


void coord::Observatory::setElevation(
                                      double const elevation 
                                     ) {
    _elevation = elevation;
}



std::string coord::Observatory::getLongitudeStr() const {
    return angleToDmsString(_longitude);
}
std::string coord::Observatory::getLatitudeStr() const {
    return angleToDmsString(_latitude);
}

std::ostream & coord::operator<<(std::ostream &os,             
                                 coord::Observatory const& obs 
                                )
{
    return os << (boost::format("%gW, %gN  %g")
                  % obs.getLatitude().asDegrees()
                  % obs.getLongitude().asDegrees()
                  % obs.getElevation()).str();
}

/************************************************************************************************************/