conversion.hh

Go to the documentation of this file.

/*
 * Copyright (c) 1999-2003 Bert Kampes
 * Copyright (c) 1999-2003 Delft University of Technology, The Netherlands
 *
 * This file is part of Doris, the Delft o-o radar interferometric software.
 *
 * Doris 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.
 *
 * Doris 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, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 * Publications that contain results produced by the Doris software should
 * contain an acknowledgment. (For example: The interferometric processing
 * was performed using the freely available Doris software package developed
 * by the Delft Institute for Earth-Oriented Space Research (DEOS), Delft
 * University of Technology, or include a reference to: Bert Kampes and
 * Stefania Usai. \"Doris: The Delft Object-oriented Radar Interferometric
 * software.\" In: proceedings 2nd ITC ORS symposium, August 1999. (cdrom)).
 *
 */
/****************************************************************
 * $Source: /users/kampes/DEVELOP/DORIS/doris/src/RCS/conversion.hh,v $
 * $Revision: 3.10 $
 * $Date: 2005/03/07 10:38:40 $
 * $Author: kampes $
 *
 * Some conversion routines
 ****************************************************************/


#ifndef CONVERSION_H                    // guard
#define CONVERSION_H

#include "constants.hh"                 // typedefs
#include "readinput.hh"                 // input structs
#include "orbitbk.hh"                   // my orbit class
#include "slcimage.hh"                  // my slc image class



// ====== Inline ======
inline real8 deg2rad(real8 x)   {return x * PI / 180.0;}
inline real4 deg2rad(real4 x)   {return x * PI / 180.0;}
inline real8 rad2deg(real8 x)   {return x * 180.0 / PI;}
inline real4 rad2deg(real4 x)   {return x * 180.0 / PI;}


// these routines are now members of class slcimage.h
// Bert Kampes, 05-Mar-2005
// ______ Convert line number to azimuth time (1 is first line) ______
//inline real8 line2ta(real8 line, real8 ta1, real8 prf)
//  {return ta1+(line-1.0)/prf;}
//
// ______ Convert pixel number to range time (1 is first pixel) ______
//inline real8 pix2tr(real8 pixel, real8 tr1, real8 rangesamplingratex2)
//  {return tr1+(pixel-1.0)/rangesamplingratex2;}
//
// ______ Convert pixel number to range (1 is first pixel) ______
//inline real8 pix2range(real8 pixel, real8 tr1, real8 rangesamplingratex2)
//  {return SOL*pix2tr(pixel,tr1,rangesamplingratex2);}
//
// ______ Convert azimuth time to line number (1 is first line) ______
//inline real8 ta2line(real8 azitime, real8 ta1, real8 prf)
//  {return 1.0+prf*(azitime-ta1);}
//
// ______ Convert range time to pixel number (1 is first pixel) ______
//inline real8 tr2pix(real8 rangetime, real8 tr1, real8 rangesamplingratex2)
//  {return 1.0+rangesamplingratex2*(rangetime-tr1);}
//
// ______ Convert range pixel to fDC (1 is first pixel) ______
// Bert Kampes, 03-Mar-2005
//inline real8 pix2fdc(real8 pixel, const slcimage &image)
//  {const real8 tau=(pixel-1.0)/(image.rsr2x/2.0);// two-way time
//   return image.f_DC_a0+image.f_DC_a1*tau+image.f_DC_a2*tau*tau;}

// ______ polar to xyz ??? ______
void pol2xyz(
        cn    &xyz,
        real8  lat,
        real8  lon,
        real8  hei=0.);

// ______ xyz to polar ??? NOT OK ______
void xyz2pol(
        const cn &xyz, 
        real8    &lat,
        real8    &lon,
        real8    &hei);

// ______ Convert xyz cartesian coordinates to ellipsoid coordinates ______
void xyz2ell(
        const input_ell &ell,
        const cn        &xyz,
        real8           &phi,
        real8           &lambda,
        real8           &height);

// ______ Convert xyz cartesian to ellipsoid coordinates phi,lambda (no h) ______
void xyz2ell(
        const input_ell &ell,
        const cn        &xyz,
        real8           &phi,
        real8           &lambda);

// ______ Convert ellipsoidal coordinates to xyz cartesian ______
void ell2xyz(
        const input_ell &ell,
        cn &xyz,
        real8 phi,
        real8 lambda,
        real8 height);

// ______ Dump tiepoint INFO ______
void tiepoint(
        const input_gen     &generalinput,
        const slcimage      &master,
        const slcimage      &slave,
              orbit         &masterorbit,
              orbit         &slaveorbit,
        const input_ell     &ellips);

// ====== Inline functions ======
// ______ For resampling output ci2, this might be slow, make function
//        i16 myround(real4) to avoid overhead real() ______
inline compli16 cr4toci2(complr4 x)
  { return compli16((x.real() > 0) ? short(x.real()+.5) : short(x.real()-.5),
                    (x.imag() > 0) ? short(x.imag()+.5) : short(x.imag()-.5));}

#endif // CONVERSION_H

---