NASA Ocean Color

ocssw V2022

 #ifndef _GEOLOCATE_OCI_H_
 #define _GEOLOCATE_OCI_H_
  
 #include "common.h"
  
 #include <netcdf>
  
 using namespace netCDF;
 using namespace netCDF::exceptions;
 /*
 #define SWAP_2(x) ( (((x) & 0xff) << 8) | ((unsigned short)(x) >> 8) )
  
 #define SWAP_4(x) ( ((x) << 24) | \
          (((x) << 8) & 0x00ff0000) | \
          (((x) >> 8) & 0x0000ff00) | \
          ((x) >> 24) )
 */
 typedef float quat_array[4];
 typedef float orb_array[3];
  
 double constexpr PI = 3.14159265358979323846;
 double constexpr RADEG = 180 / PI;
 double constexpr DTOR = PI / 180;
  
 float constexpr focal_length = 45.184;
  
  
 class geoFile {
   //  netCDF::NcFile *geofile;
  
   std::string fileName;
   
   int ngrps;
   int ndims;
   
   netCDF::NcDim ncDims[1000];
   
  public:
   geoFile();
   ~geoFile();
  
   netCDF::NcFile *geofile;
   netCDF::NcGroup ncGrps[10];
  
   int createFile( const char* filename, const char* cdlfile,
                   size_t sdim, int *ncid, int *gid,
                   uint32_t bbb, uint32_t rbb,
                   int16_t pcdim, int16_t psdim, size_t nswband, int32_t *rta_nadir,
                   bool radiance);
  
   int parseDims( std::string dimString, std::vector<netCDF::NcDim>& varDims);
  
   int close();
 };
  
 int get_nbands( uint32_t ntaps, int16_t jagg[16], uint32_t& nbb);
  
 int j2000_to_ecr(int32_t iyr, int32_t idy, double sec, double ecmat[3][3]);
 int j2000_to_mod(int32_t iyr, int32_t idy, double sec, double j2mod[3][3]);
 int get_nut(int32_t iyr, int32_t idy, double xnut[3][3]);
 int get_ut1(int32_t iyr, int32_t idy, double &ut1utc);
 int ephparms(double t, double &xls, double &gs, double &xlm, double &omega);
 int nutate(double t, double xls, double gs, double xlm, double omega,
            double &dpsi, double &eps, double &epsm);
 int gha2000(int32_t iyr, double day, double &gha);
 int mtoq(double rm[3][3], double q[4]);
 int qprod(double q1[4], float q2[4], double q3[4]);
 int qprod(float q1[4], float q2[4], float q3[4]);
 int orb_interp(size_t n_SC_rec, size_t sdim,
                double *torb, orb_array *p, orb_array *v,
                double *time, orb_array *posi, orb_array *veli);
 int q_interp(size_t n_SC_rec, size_t sdim, double *tq, quat_array *q,
              double *time, quat_array *qi);
 int tilt_interp(size_t n_tilts, size_t sdim, double *ttilt, float *tiltin,
                 double *time, float *tilt);
 int l_sun(size_t sdim, int32_t iyr, int32_t iday,
           double *sec, orb_array *sunr, double *rs);
 int sun2000(size_t sdim, int32_t iyr, int32_t idy,
             double *sec, orb_array *sun, double *rs);
 int qtom(float quat[4], double rm[3][3]);
 int scan_ell(float p[3], double sm[3][3], double coef[10]);
  
 int oci_geonav(const char* dem_file, float pos[3], float vel[3],
                double smat[3][3], double coef[10],
                float sunr[3], float **pview, size_t npix, double *delt,
                float *xlat, float *xlon, short *solz, short *sola,
                short *senz, short *sena, short *height,
                float& lonwest, float& loneast,
                float& latsouth, float& latnorth);
  
 int mat2rpy(float pos[3], float vel[3], double smat[3][3], float rpy[3]);
  
  
  
 inline
 int expandEnvVar( std::string *sValue) {
   if ( (*sValue).find_first_of( "$" ) == std::string::npos) return 0;
   std::string::size_type posEndIdx = (*sValue).find_first_of( "/" );
   if ( posEndIdx == std::string::npos) return 0;
   const std::string envVar = sValue->substr (1, posEndIdx - 1);
   char *envVar_str = getenv(envVar.c_str());
   if (envVar_str == 0x0) {
     printf("Environment variable: %s not defined.\n", sValue->c_str());
     exit(1);
   }
   *sValue = envVar_str + (*sValue).substr( posEndIdx);
  
   return 0;
 }
  
  
 geoFile::geoFile() {
  
 }
  
  
 geoFile::~geoFile() {
   
 }
  
  
 #endif  // _GEOLOCATE_OCI_H_

focal_length

constexpr float focal_length

Definition: geolocate_oci.h:25

data_t q

Definition: decode_rs.h:74

data_t t[NROOTS+1]

Definition: decode_rs.h:77

get_nbands

int get_nbands(uint32_t ntaps, int16_t jagg[16], uint32_t &nbb)

nutate

int nutate(double t, double xls, double gs, double xlm, double omega, double &dpsi, double &eps, double &epsm)

Definition: geolocate_hawkeye.cpp:604

day

int32_t day

Definition: atrem_corl1v2.h:308

get_ut1

int get_ut1(int32_t iyr, int32_t idy, double &ut1utc)

Definition: geolocate_hawkeye.cpp:635

ephparms

int ephparms(double t, double &xls, double &gs, double &xlm, double &omega)

Definition: geolocate_hawkeye.cpp:575

mtoq

int mtoq(double rm[3][3], double q[4])

Definition: geolocate_hawkeye.cpp:786

createFile

int createFile(const char *filename, const char *cdlfile, size_t sdim, int *ncid, int *gid)

Definition: hawkeyeUtil.cpp:39

color_dtdb.v

Definition: color_dtdb.py:430

qprod

int qprod(double q1[4], float q2[4], double q3[4])

Definition: geolocate_hawkeye.cpp:823

pos

float32 * pos

Definition: l1_czcs_hdf.c:35

orb_array

float orb_array[3]

Definition: geolocate_oci.h:19

constexpr double PI

Definition: geolocate_oci.h:21

dpsi

double dpsi

Definition: sun2000.c:6

eps

double eps

Definition: gha2000.c:3

string

@ string

Definition: GEO_read_param_file.c:64

RADEG

constexpr double RADEG

Definition: geolocate_oci.h:22

hico.proc_hico.solz

solz

Definition: proc_hico.py:240

hico.proc_hico.senz

senz

Definition: proc_hico.py:239

sun2000

int sun2000(size_t sdim, int32_t iyr, int32_t idy, double *sec, orb_array *sun, double *rs)

Definition: geolocate_oci.cpp:1362

quat_array

float quat_array[4]

Definition: geolocate_oci.h:18

expandEnvVar

int expandEnvVar(std::string *sValue)

Definition: geolocate_oci.h:96

j2000_to_ecr

int j2000_to_ecr(int32_t iyr, int32_t idy, double sec, double ecmat[3][3])

Definition: geolocate_hawkeye.cpp:455

geoFile::geofile

netCDF::NcFile * geofile

Definition: geolocate_oci.h:42

omega

data_t omega[NROOTS+1]

Definition: decode_rs.h:77

hico.proc_hico.sena

sena

Definition: proc_hico.py:241

j2000_to_mod

int j2000_to_mod(int32_t iyr, int32_t idy, double sec, double j2mod[3][3])

Definition: geolocate_hawkeye.cpp:512

oci_geonav

int oci_geonav(const char *dem_file, float pos[3], float vel[3], double smat[3][3], double coef[10], float sunr[3], float **pview, size_t npix, double *delt, float *xlat, float *xlon, short *solz, short *sola, short *senz, short *sena, short *height, float &lonwest, float &loneast, float &latsouth, float &latnorth)

Definition: geolocate_oci.cpp:1486

tilt_interp

int tilt_interp(size_t n_tilts, size_t sdim, double *ttilt, float *tiltin, double *time, float *tilt)

Definition: geolocate_oci.cpp:1306

scan_ell

int scan_ell(float p[3], double sm[3][3], double coef[10])

Definition: geolocate_hawkeye.cpp:1078

filename

char filename[FILENAME_MAX]

Definition: atrem_corl1.h:122

parseDims

int parseDims(int ncid, int ndims, string dimString, int *numDims, int *dimid, int *varDims)

Definition: hawkeyeUtil.cpp:339

xlon

float xlon[LAC_PIXEL_NUM]

Definition: l1a_seawifs.c:93

orb_interp

int orb_interp(size_t n_SC_rec, size_t sdim, double *torb, orb_array *p, orb_array *v, double *time, orb_array *posi, orb_array *veli)

Definition: geolocate_hawkeye.cpp:845

hico.proc_hico.sola

sola

Definition: proc_hico.py:242

DTOR

constexpr double DTOR

Definition: geolocate_oci.h:23

time

this program makes no use of any feature of the SDP Toolkit that could generate such a then geolocation is calculated at that and then aggregated up to Resolved feature request Bug by adding three new int8 SDSs for each high resolution offsets between the high resolution geolocation and a bi linear interpolation extrapolation of the positions This can be used to reconstruct the high resolution geolocation Resolved Bug by delaying cumulation of gflags until after validation of derived products Resolved Bug by setting Latitude and Longitude to the correct fill resolving to support Near Real Time because they may be unnecessary if use of entrained ephemeris and attitude data is turned resolving bug report Corrected to filter out Aqua attitude records with missing status helping resolve bug MOD_PR03 will still correctly write scan and pixel data that does not depend upon the start time

Definition: HISTORY.txt:248

l_sun

int l_sun(size_t sdim, int32_t iyr, int32_t iday, double *sec, orb_array *sunr, double *rs)

Definition: geolocate_oci.cpp:1337

get_nut

int get_nut(int32_t iyr, int32_t idy, double xnut[3][3])

Definition: geolocate_hawkeye.cpp:547

common.h

mat2rpy

int mat2rpy(float pos[3], float vel[3], double smat[3][3], float rpy[3])

Definition: geolocate_oci.cpp:1701

geoFile