a_unwarp.cxx

Go to the documentation of this file.

/*
Copyright 2000-12 Pierre SMARS (smars@yuntech.edu.tw)
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 2 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 GNU General Public License
along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "a_unwarp.h"
//---------------------------------------------------------------------------
void a_unwarp::unwarp(double u, double v, double & un, double& vn)
{
        double du = u-u0_;
        double dv = v-v0_;
        double r2 = du*du+dv*dv;
        un = u + K1_*du*r2 + K2_*du*r2*r2 + P1_*(2*du*du+r2) + 2*P2_*du*dv;
        vn = v + K1_*dv*r2 + K2_*dv*r2*r2 + P2_*(2*dv*dv+r2) + 2*P1_*du*dv;
}
//---------------------------------------------------------------------------
std::istream& operator>> (std::istream& i, a_unwarp& l)
{
        i >> l.u0_ >> l.v0_;
        i >> l.K1_ >> l.K2_;
        i >> l.P1_ >> l.P2_;
        return i;
}
//---------------------------------------------------------------------------
std::ostream& operator<< (std::ostream& o, const a_unwarp& l)
{
        o << l.u0_ << " " << l.v0_ << std::endl;
        o << l.K1_ << " " << l.K2_ << std::endl;
        o << l.P1_ << " " << l.P2_ << std::endl;
        return o;
}