GLIMSView Source Code Documentation

hdfeosfrmt.h

Go to the documentation of this file.

#ifndef __HDFEOSFRMT_H
#define __HDFEOSFRMT_H

#include "imageformat.h"
#include "odlparser.h"

#include <hdf.h>
#include <mfhdf.h>
#include <HdfEosDef.h>
#include <proj_api.h>

#include <string>
#include <cstring>
#include <vector>
#include <cmath>

#include "TraceLog.h"

#ifndef M_PI
#define M_PI 3.1415926
#endif

class HdfEosFrmt : public ImageFormat {
protected:
  int32   mSWfid;

  typedef enum {
    UNSUPPORTED, 
    UTM          
  } Projection ;

  typedef struct {
    std::string                 swathname;
    std::string                 gfldname;
    std::vector<int32>          dims;
    std::vector<std::string>    dimnames;
    int32                       type;
    void*                       geomap;
  } GeoField ;

  typedef struct {
    std::string                 swathname;
    std::string                 fieldname;
    std::vector<int32>          dims;
    std::vector<std::string>    dimnames;
    int32                       type;
    GeoField*                   geofldlat;
    GeoField*                   geofldlon;
    std::vector<int32>          latoff;
    std::vector<int32>          lonoff;
    std::vector<int32>          latinc;
    std::vector<int32>          loninc;
  } DataField ;

  std::vector<DataField>  mDFList;
  std::vector<GeoField> mGFList;
  Projection              mProj;

  void      parseGeoFlds( int32 swid, std::string swname );

  std::string     getHdfGblAttr( const std::string &fname,
                                 const std::string &attrname ) const;

  std::string   getAquiDateTime( const std::string &fname ) const;

  std::string   getProjection( const std::string &fname ) const;

  std::string   getOrigID( const std::string &fname ) const;

  std::string   getInstZenith( const std::string &fname ) const;

  std::string   getCloudPct( const std::string &fname ) const;

  std::string   getSunAzim( const std::string &fname ) const;

  std::string   getSunElev( const std::string &fname ) const;

  std::string   getImageAzim( const std::string &fname ) const;


  //std::string   getImageID( const std::string &fname ) const;

  std::string   getInstID( const std::string &fname ) const;

  std::string   getInstAz( const std::string &fname ) const;


  void          getBandMetaData( ) const;

  std::vector<std::string>
  splitList( char *clist, char delim ) {
    int l, r;
    std::vector<std::string> ret;
    std::string list = clist;
    l = 0;
    while ( l < (int)list.length()-1 ) {
      r = list.find( delim, l );
      if ( r == -1 ) r = list.length();
      ret.push_back( list.substr( l, r-l ) );
      l = r+1;
    }
    return ret;
  }

  template<class DType>
  void findrange( int32 swid,
                  std::string field,
                  DType &lwr,
                  DType &upr,
                  std::vector<int32> dims );

  template<class DType>
  void maxmin( DType* data, int sz, DType &lwr, DType &upr ) {
    for ( int i=0; i < sz; i++ ) {
      if ( lwr > data[i] )
        lwr = data[i];
      if ( upr < data[i] )
        upr = data[i];
    }
  }

  template<class DType>
  void bilinterp( DType  *map,
                  std::vector<int32>  dims,
                  std::vector<int32>  offset,
                  std::vector<int32>  inc,
                  int    x,
                  int    y,
                  double &res );

  void getGeoFields( int iband );
  void getIntersect( int &track,
                     int &p1,
                     int p2,
                     int iband,
                     double searchVal,
                     bool isLon,
                     bool horizontal=true );

  void trackLineTo( int &x,
                    int &y,
                    int iband,
                    double trackVal,
                    double searchVal );

  /*
  void getULCorners( );
  void xyToLL( int iband, double &x, double &y );
  void xyToMeter( int iband, double &x, double &y );
  void meterToLL( int iband, double &x, double &y );
  void llToMeter( int iband, double &lon, double &lat );
  void llToXY( int iband, double &lon, double &lat );
  projPJ  createProj( int iband );
  */

  bool between( double a, double b, double c ) {
    if ( a < b )
      return( c > a && c < b );
    else
      return( c < a && c > b );
  }

public:
  ~HdfEosFrmt( );

  virtual bool open( std:: string fname );

  virtual void getLL( int iband, double &i, double &j, bool toWGS84 );
  virtual void getLL( int iband, double &i, double &j ) {
    if ( iband >= 0 && iband < (int)mBandSet.size() &&
         mBandSet[iband].proj != "" ) {
      xyToLL( iband, i, j );
      return;
    }
    getLL( iband, i, j, true );
  }

  virtual void getXY( int iband, double &i, double &j );
  virtual void* getBand( int iband,
                         int x1,
                         int y1,
                         int w1,
                         int h1,
                         int w2,
                         int h2,
                         ScaleType = SCALE_BI ) const;
};

#endif

---

Home | Search | Disclaimers & Privacy | Contact Us GLIMSView Maintainer: dsoltesz@usgs.gov