RectangularRegion.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 "lsst/pex/exceptions.h"

#include "lsst/ap/CircularRegion.h"
#include "lsst/ap/RectangularRegion.h"
#include "lsst/ap/SpatialUtil.h"


namespace ex = lsst::pex::exceptions;

lsst::ap::RectangularRegion::RectangularRegion(
    double const minRa,
    double const maxRa,
    double const minDec,
    double const maxDec
) :
    _minRa(minRa),
    _maxRa(maxRa),
    _minDec(minDec),
    _maxDec(maxDec)
{
    if (minRa < 0.0 || minRa >= 360.0 || maxRa < 0.0 || maxRa >= 360.0) {
        throw LSST_EXCEPT(ex::RangeError,
                          "right ascension must be in range [0, 360) degrees");
    }
    if (minDec < -90.0 || minDec > 90.0 || maxDec < -90.0 || maxDec > 90.0) {
        throw LSST_EXCEPT(ex::RangeError,
                          "declination must be in range [-90, 90] degrees");
    }
    if (maxDec < minDec) {
        throw LSST_EXCEPT(ex::InvalidParameterError,
                          "minimum declination greater than maximum declination");
    }
}


lsst::ap::RectangularRegion::RectangularRegion(
    double const centerRa,
    double const centerDec,
    double const radius
) {
    fromCircle(centerRa, centerDec, radius);
}


lsst::ap::RectangularRegion::RectangularRegion(CircularRegion const & region) {
    fromCircle(region.getCenterRa(), region.getCenterDec(), region.getRadius());
}


void lsst::ap::RectangularRegion::fromCircle(
    double const ra,
    double const dec,
    double const radius
) {
    if (ra < 0.0 || ra >= 360.0) {
        throw LSST_EXCEPT(ex::RangeError,
                          "right ascension must be in range [0, 360) degrees");
    }
    if (dec < -90.0 || dec > 90.0) {
        throw LSST_EXCEPT(ex::RangeError,
                          "declination must be in range  [-90, 90] degrees");
    }
    if (radius < 0.0 || radius > 90.0) {
        throw LSST_EXCEPT(ex::RangeError,
                          "circle radius must be in range  [0, 90] degrees");
    }
    double alpha = maxAlpha(radius, dec);
    _minRa = ra - alpha;
    if (_minRa < 0.0) {
        _minRa += 360.0;
    }
    _maxRa = ra + alpha;
    if (_maxRa >= 360.0) {
        _maxRa -= 360.0;
    }
    _minDec = dec - radius;
    if (_minDec < -90.0) {
        _minDec = -90.0;
    }
    _maxDec = dec + radius;
    if (_maxDec > 90.0) {
        _maxDec = 90.0;
    }
}