Coordinates.h

/*
 
MIT License
 
Copyright (c) 2017 FMI Open Development / Markus Peura, first.last@fmi.fi
 
Permission is hereby granted, free of charge, to any person obtaining a copy
of this software and associated documentation files (the "Software"), to deal
in the Software without restriction, including without limitation the rights
to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
copies of the Software, and to permit persons to whom the Software is
furnished to do so, subject to the following conditions:
 
The above copyright notice and this permission notice shall be included in all
copies or substantial portions of the Software.
 
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
SOFTWARE.
 
*/
/*
Part of Rack development has been done in the BALTRAD projects part-financed
by the European Union (European Regional Development Fund and European
Neighbourhood Partnership Instrument, Baltic Sea Region Programme 2007-2013)
*/
#ifndef RADAR__COORDINATES_ 
#define RADAR__COORDINATES_ "radar__coordinates 0.2, May 16 2011 Markus.Peura@fmi.fi"
 
#include <drain/image/CoordinatePolicy.h>
#include <math.h>
 
#include <iostream>
#include <ostream>
#include <sstream>
 
// #include <drain/util/Proj6.h>
// See also RadarProj4 @ Geometry?
#include <drain/util/Rectangle.h>
#include "Constants.h"
 
 
namespace rack {
 
extern //static
const drain::image::CoordinatePolicy polarLeftCoords;
 
extern //static
const drain::image::CoordinatePolicy limitCoords;
 
/*
class RadarProj : public drain::Proj4 {
 
public:
 
    inline
    void setSiteLocationDeg(double lon, double lat){
        std::stringstream s;
        s << "+proj=aeqd" << " +lon_0=" << lon << " +lat_0=" << lat << " +ellps=WGS84 +type=crs";
        setProjectionSrc(s.str());
    }
 
    void setSiteLocation(double lon, double lat){
        setSiteLocationDeg(lon*drain::RAD2DEG, lat*drain::RAD2DEG);
    }
 
    inline
    void setLatLonProjection(){
        setProjectionDst("+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs");
    }
 
    void determineBoundingBoxM(double range, double & xLL, double & yLL, double & xUR, double & yUR) const;
 
    inline
    void determineBoundingBoxM(double range, drain::Rectangle<double> & bbox) const {
        determineBoundingBoxM(range, bbox.lowerLeft.x, bbox.lowerLeft.y, bbox.upperRight.x, bbox.upperRight.y);
    }
 
 
 
    //void determineBoundingBoxD(double range, double & xLL, double & yLL, double & xUR, double & yUR) const;
 
 
    inline
    void XXdetermineBoundingBoxD(double range, drain::Rectangle<double> & bbox) const {
        determineBoundingBoxD(range, bbox.lowerLeft.x, bbox.lowerLeft.y, bbox.upperRight.x, bbox.upperRight.y);
    }
 
 
};
 
*/
 
//  DEPRECATED
  class Coordinates {
  public:
 
 
      Coordinates();
 
      virtual
     ~Coordinates(){};
 
    //    coordinator();
 
    void setOrigin(const double &theta,const double &phi); 
 
 
    // inline
    void setOriginDeg(const double &lat,const double &lon);
      //      origin(lat*M_PI/180.0 , lon*M_PI/180.0);
      //    }; 
  //const double &lat,const double &lon);
 
    void setBinPosition(const double &alpha, const float &range);
   
    // \param alpha is azimuth in radians, \param range in metres.
    //void setBinPosition(double alpha, float range);
   
   //  void getBoundingBox(float range,double &latMin,double &lonMin,double &latMax,double &lonMax);
 
    void info(std::ostream &ostr = std::cout);
 
    // site positional
    // double cos_theta;
 
    // (Site normal unit vector not needed as such)
    double p01, p02, p03;
 
    double p1, p2, p3;
 
    double e11, e12, e13;
 
    double e21, e22, e23;
 
    mutable double thetaBin;
 
    mutable double phiBin;
 
    inline
    double binLatitudeDeg(){ return thetaBin/M_PI*180.0;};
                //phi_bin/M_PI*180.0;};
                //
 
    inline
    double binLongitudeDeg(){ return phiBin/M_PI*180.0;};
    //    theta_bin/M_PI*180.0;};;
    //
 
    //drain::Proj6 proj;
 
    /*
    inline
    virtual std::string getProjectionString(){
        std::stringstream sstr;
        sstr << "+proj=longlat +R=" << EARTH_RADIUS << std::string(" +no_defs");
        return sstr.str();
    };
    */
  };
 
} // ::rack
 
#endif
 
// Rack