slcimage.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/slcimage.hh,v $   *
 * $Revision: 3.12 $                                            *
 * $Date: 2005/04/11 13:47:45 $                                 *
 * $Author: kampes $                                            *
 *                                                              *
 * The slcimage class contains the definition of the data and   *
 * functions for slc images.                                    *
 * Data mainly public because this used to be a struct and I    *
 * did not want to change the other code.                       *
 * functions for updating the data like prf, reading data into  *
 * a matrix, etc.                                               *
 #%// BK 14-Aug-2000                                            *
 ****************************************************************/


#ifndef SLCIMAGE_H
#define SLCIMAGE_H


#include "constants.hh"                 // typedefs


// ====== Define template functions (no member no friend) ======
// ______ (matrix class is declared way below) ______
template <class Type> class matrix;




// ====== Struct slcimage: information on master/slave ======
class slcimage                          // info per image
  {
  private:
    // copy constructor; copy matrices etc as public.
    slcimage(const slcimage& img)
      {};// prevent copy constructor usage by privatizing

  public:
    char        file[EIGHTY];           // current filename
    char        utc1[25];               // string for getorb
    int16       sensor;                 // flag for ers/asar/rsat/jers [BK 21.08.2002]
    int16       sar_processor;          // VMP or ATLANTIS or TUDELFT [BK 04.2005]
    int16       formatflag;             // flag for current file format
    cn          approxcentreoriginal;   //  wgs84? sphere x,y,z
    window      originalwindow;         // linelo etc. for normalization (master only).
    window      currentwindow;          // linelo etc. of current on file.
    real8       prf;                    // pulse repetition frequency (Hz)
    real8       abw;                    // azimuth band width (Hz)
    real8       rsr2x;                  // 2 times range sampling rate (Hz)
    real8       rbw;                    // range band width (Hz)
    real8       t_azi1;                 // sec. of day of first azimuth line
    real8       t_range1;               // one way time (s) to first range pixel
    real8       wavelength;             // meters
    // real8 hamming_azi;               // weighting function
    // real8 hamming_range;             // weighting function
    // ______ Xtrack f_DC coefficients from leader, not used for now, ______
    // ______ better use matrix for polyval?
    real8       f_DC_a0;                // constant term Hz
    real8       f_DC_a1;                // linear term Hz/s
    real8       f_DC_a2;                // quadratic term Hz/s/s

    // ______ offset = X(l,p) - X(L,P) ______
    // ______ Where l,p are in the local slave coordinate system and ______
    // ______ where L,P are in the local master coordinate system ______
    // ______ These variables are stored in the slaveinfo variable only ______
    int32       coarseoffsetL;          // offset in line (azimuth) direction
    int32       coarseoffsetP;          // offset in pixel (range) direction

  
    // ______ Public functions of class ______
    // ___Constructor/Destructor___
    slcimage();// set defaults
    ~slcimage()
      {
      TRACE_FUNCTION("~slcimage() (BK 06-Mar-2005)");
      ;// nothing to destruct ?
      }// dealloc

    // ---- Helper ----
    inline void showdata() const
      {DEBUG << "\tcurrent file: \t" << file
           << "\n\tformatflag:   \t" << formatflag
           << "\n\tcurrent win:  \t"
           << currentwindow.linelo << " " << currentwindow.linehi << " "
           << currentwindow.pixlo  << " " << currentwindow.pixhi; 
       DEBUG.print();
      }

    // ______ Add a bias 1 way time range ______
    inline void add_rg_t_error(const real8 dt) {t_range1+=dt;}
    inline void add_az_t_error(const real8 dt) {t_azi1+=dt;}

    // ______ Read readfiles section from resultfile ______
    void fillslcimage(const char *file);
    void updateslcimage(const char *file, const char *iden);    
    //____RaffaeleNutricato START MODIFICATION SECTION 1
    void updateslcimageML(const char *file, const char *iden);  
    //____RaffaeleNutricato END MODIFICATION SECTION 1
    // ______ Read matrix from file (all FORMATS) ______
    //void readdata(matrix<complr4> &Result,window win); 
    //void readdata(matrix<complr4> &Result,window win) const; 
    matrix<complr4> readdata(window win) const; 

    // ______ Convert line number to azimuth time (1 is first line) ______
    inline real8 line2ta(const real8 line) const
      {return t_azi1+(line-1.0)/prf;}

    // ______ Convert pixel number to range time (1 is first pixel) ______
    inline real8 pix2tr(const real8 pixel) const
      {return t_range1+(pixel-1.0)/rsr2x;}

    // ______ Convert pixel number to range (1 is first pixel) ______
    inline real8 pix2range(const real8 pixel) const
      {return SOL*pix2tr(pixel);}

    // ______ Convert azimuth time to line number (1 is first line) ______
    inline real8 ta2line(const real8 azitime) const
      {return 1.0+prf*(azitime-t_azi1);}

    // ______ Convert range time to pixel number (1 is first pixel) ______
    inline real8 tr2pix(const real8 rangetime) const
      {return 1.0+rsr2x*(rangetime-t_range1);}

    // ______ Convert range pixel to fDC (1 is first pixel) ______
    // Bert Kampes, 03-Mar-2005
    inline real8 pix2fdc(const real8 pixel) const
      {const real8 tau=(pixel-1.0)/(rsr2x/2.0);// two-way time
       return f_DC_a0+f_DC_a1*tau+f_DC_a2*sqr(tau);}


  }; // END class slcimage


#endif // SLCIMAGE_H

---