GEO_write_parameters.c

Go to the documentation of this file.

#include <stddef.h>
#include "PGS_MODIS_35251.h"
#include "GEO_product.h"
#include "GEO_output.h"
 
int GEO_write_parameters(
    MODFILE         * const geo_file,
    GEO_param_struct const  * parameter
    )
{
 
  /* Local variable declarations */
 
  #define BANDS (long)(MAX_BAND_NUMBER+1)
    struct {
      char *name;
      long rank;
      long dimensions[2];
      char data_type[DATATYPELENMAX];
      size_t data;
      char *units;
    } SDS_object[] = {
    {FOCAL_LENGTH,     1L, {BANDS}, R64,
      offsetof(GEO_param_struct, geometry_params.focal_length), "meters"},
    {BAND_POSITION,    1L, {BANDS}, R64,
      offsetof(GEO_param_struct, geometry_params.band_position),
      "frame sample time"},
    {DETECTOR_SPACE,   1L, {BANDS}, R64,
      offsetof(GEO_param_struct, geometry_params.det_space), "meters"},
    {DETECTOR_OFFSETS, 2L, {BANDS, 2L}, R64,
      offsetof(GEO_param_struct, geometry_params.det_position), "meters"},
    {T_OFFSET,         1L, {BANDS}, R64,
      offsetof(GEO_param_struct, geometry_params.F_offset), 
      "frame sample time"},
    {NUM_SAMPLES,      1L, {BANDS},     UI16,
      offsetof(GEO_param_struct, geometry_params.N_samp), NULL}
    };
    long start[] = {0L, 0L, 0L};
    int obj;
    int retval = SUCCESS;
    uint32 dim_index[] = {0,1};
    char msgbuf[PGS_SMF_MAX_MSGBUF_SIZE];
    char filefunc[] = __FILE__ ", GEO_write_parameters";
 
   
    if(parameter==NULL || geo_file==NULL)
    {
    sprintf(msgbuf, "parameter=%p, geo_file=%p",
        (void*)parameter, (void*)geo_file);
    modsmf(MODIS_E_BAD_INPUT_ARG, msgbuf, filefunc);
 
    return FAIL;
    }
 
    /* Begin output of parameters */
    for(obj=0; obj<(int)(sizeof(SDS_object)/sizeof(SDS_object[0]))
      && retval==SUCCESS; obj++)
    {
      if(createMODISarray(geo_file, SDS_object[obj].name, PARM_GRP,
      SDS_object[obj].data_type, SDS_object[obj].rank,
      SDS_object[obj].dimensions) !=MAPIOK)
      {
    retval = FAIL;
    sprintf(msgbuf, "createMODISarray(\"%s\", \"" PARM_GRP
        "\", \"%s\", %ld)", SDS_object[obj].name, SDS_object[obj].data_type,
        SDS_object[obj].rank);
        
    modsmf(MODIS_E_GEO, msgbuf, filefunc);
      }
      else if( putMODISdimname(geo_file, SDS_object[obj].name, PARM_GRP, 0,
      NUMBANDS) != MAPIOK)
      {
      retval = FAIL;
      sprintf(msgbuf, "putMODISdimname(\"%s\", \"%s\", \"" PARM_GRP
          "\", 0, \"" NUMBANDS "\")",
          geo_file->filename, SDS_object[obj].name);
      modsmf(MODIS_E_GEO, msgbuf, filefunc);
      }
      else if( putMODISarray(geo_file, SDS_object[obj].name, PARM_GRP, start,
      SDS_object[obj].dimensions, (char *)parameter+SDS_object[obj].data)
      != MAPIOK)
      {
      retval = FAIL;
      sprintf(msgbuf, "putMODISarray(\"%s\", \"%s\", \"" PARM_GRP "\")",
          geo_file->filename, SDS_object[obj].name);
      modsmf(MODIS_E_GEO, msgbuf, filefunc);
      }
      else if(SDS_object[obj].units!=NULL)
      {
      if(putMODISarinfo(geo_file, SDS_object[obj].name, PARM_GRP, MUNITS,
          TXT, (long)strlen(SDS_object[obj].units), SDS_object[obj].units)
          != MAPIOK)
      {
          retval = FAIL;
          sprintf(msgbuf, "putMODISarinfo(\"%s\", \"%s\", \"" PARM_GRP
          "\", \"" MUNITS "\", \"" TXT "\", %ld, \"%s\")",
          geo_file->filename, SDS_object[obj].name,
          (long)strlen(SDS_object[obj].units), SDS_object[obj].units);
          modsmf(MODIS_E_GEO, msgbuf, filefunc);
      }
      }
    } /* for each SDS_object    */
 
    if(putMODISfileinfo(geo_file, BAND_NUMBER, I16, 1,
    &parameter->geometry_params.band_number) != MAPIOK)
    {
    retval = FAIL;
    sprintf(msgbuf,
        "putMODISfileinfo(\"%s\", \"" BAND_NUMBER "\", \"" I16 "\", 1, %d)",
        geo_file->filename, (int)parameter->geometry_params.band_number);
    modsmf(MODIS_E_GEO, msgbuf, filefunc);
    }
 
    if(putMODISarinfo(geo_file, DETECTOR_OFFSETS, PARM_GRP, SCANDIM, UI32,
    1L, &dim_index[0]) != MAPIOK)
    {
    retval = FAIL;
    sprintf(msgbuf, "putMODISarinfo(\"%s\", \"" PARM_GRP "\", \"" SCANDIM
        "\", \"" UI32 "\", 1, %lu)", geo_file->filename,
        (unsigned long)dim_index[0]);
    modsmf(MODIS_E_GEO, msgbuf, filefunc);
    }
 
    if(putMODISarinfo(geo_file, DETECTOR_OFFSETS, PARM_GRP, TRACKDIM, UI32,
    1L, &dim_index[1]) != MAPIOK)
    {
    retval = FAIL;
    sprintf(msgbuf, "putMODISarinfo(\"%s\", \"" PARM_GRP "\", \"" TRACKDIM
        "\", \"" UI32 "\", 1, %lu)",
        geo_file->filename, (unsigned long)dim_index[1]);
    modsmf(MODIS_E_GEO, msgbuf, filefunc);
    }
 
    return retval;
}