NASA Ocean Color

ocssw V2022

 /*******************************************************************
  
    hdfio.c - general routines for hdf file input, output
  
    Contents:
      hdfio_open - open the file
      hdfio_rd_gattr - read a global attribute
      hdfio_rd_sd - read a science dataset
      hdfio_close - close a hdf file
  
  *******************************************************************/
  
 #include <string.h>
 #include "hdfio.h"
  
 int hdfio_open(char *fname, hdfio_struc *hdfinfo)
 /*******************************************************************
  
    hdfio_open
  
    purpose: open a general hdf file
  
    Returns type: int - return status: 0 is good
  
    Parameters: (in calling order)
       Type              Name            I/O     Description
       ----              ----            ---     -----------
       char *            fname            I      name of file to open
       hdfio_struc *     hdfinfo          O      information  struct
                                                 about the opened file
  
    Modification history:
       Programmer        Date            Description of change
       ----------        ----            ---------------------
       W. Robinson, SAIC 5 Aug 2004      Original development
  
  *******************************************************************/
  {
  
     if ((hdfinfo->fid = Hopen(fname, DFACC_RDONLY, 0)) < 0) {
         printf("hdfio_open: Failed on the Hopen of \n%s\n", fname);
         return -1;
     }
  
     /*
      *  SDstart opens the hdf interface and inits the SD interface, wierd.
      *  apparently, if both SD and HDF are used, both these need to
      *  be called.
      */
     if ((hdfinfo->sdfid = SDstart(fname, DFACC_RDONLY)) < 0) {
         printf("hdfio_open: Failure at SDstart of \n%s\n", fname);
         hdfio_close(*hdfinfo);
         return -1;
     }
  
     /*
      *  ok, it's open.
      */
     return 0;
 }
  
 int hdfio_rd_gattr(hdfio_struc hdfinfo, char *name, int32 n_type,
         int32 count, void *data)
 /*******************************************************************
  
    hdfio_rd_gattr
  
    purpose: read a global attribute to a waiting array (if size correct
               and type is correct)
  
    Returns type: int - return status: 0 is good
  
    Parameters: (in calling order)
       Type              Name            I/O     Description
       ----              ----            ---     -----------
       hdfio_struc       hdfinfo          I      information  struct
                                                 about the opened file
       char *            name             I      name of glabal attr to read
       int32             n_type           I      number type expected using
                                                 HDF nomenclature
       int32             count            I      count of values
       void*             data             O      array of data from the
                                                 global attribute 
  
    Modification history:
       Programmer        Date            Description of change
       ----------        ----            ---------------------
       W. Robinson, SAIC 5 Aug 2004      Original development
  
  *******************************************************************/
  {
     int32 attr_index;
     char a_name[200];
     int32 a_n_type, a_count;
  
     /*
      *  get the index for that attribute name
      */
     if ((attr_index = SDfindattr(hdfinfo.sdfid, name))
             == -1) {
         printf("hdfio_rd_gattr: Error in SDfindattr for attribute: '%s'\n",
                 name);
         return -1;
     }
     /*
      *  find out about the attribute and be sure it is the type and size expected
      */
     if (SDattrinfo(hdfinfo.sdfid, attr_index, a_name, &a_n_type, &a_count)
             == -1) {
         printf("hdfio_rd_gattr: Error in SDattrinfo for attribute: '%s'\n",
                 name);
         return -1;
     }
  
     if (a_n_type != n_type) {
         printf(
                 "hdfio_rd_gattr: attribute: '%s', difference in number type found\n",
                 name);
         printf("          expected: %d, read: %d\n", (int) n_type, (int) a_n_type);
         return -1;
     }
  
     if (((n_type != DFNT_CHAR) && (a_count != count)) ||
             ((n_type == DFNT_CHAR) && (count < a_count))) {
         printf("hdfio_rd_gattr: attribute: '%s', difference in count found\n",
                 name);
         printf("             expected: %d, read: %d\n", (int) count, (int) a_count);
         return -1;
     }
  
     /*
      *  read in the stuff
      */
     if (SDreadattr(hdfinfo.sdfid, attr_index, data) == -1) {
         printf("hdfio_rd_gattr: attribute: '%s', Error reading data\n",
                 name);
         return -1;
     }
     return 0;
 }
  
 int hdfio_rd_sd(hdfio_struc hdfinfo, char *arr_name, void *array)
 /*******************************************************************
  
    hdfio_rd_sd
  
    purpose: read a full science dataset to a float array
  
    Returns type: int - return status: 0 is good
  
    Parameters: (in calling order)
       Type              Name            I/O     Description
       ----              ----            ---     -----------
       hdfio_struc       hdfinfo          I      hdf information  struct
                                                 about the opened file
       char *            arr_name         I      name of science dataset 
                                                 to read
       void *            array            O      array of data from the sds
                                                 already allocated outside
  
    Modification history:
       Programmer        Date            Description of change
       ----------        ----            ---------------------
       W. Robinson       9-Feb-1995      Original development
       W. Robinson, SAIC  5 Aug 2004     update this version to do more general 
                                         void array type
  
  *******************************************************************/
  {
     int j;
     int32 index, rank, sdid, numbertype, nattrs, data_dims[5];
     int32 start[5], edge[5];
     char name[H4_MAX_NC_NAME]; /*  MAX_NC_NAME is max # chars for this */
  
     /*
      *  SDnametoindex is used to zero in on the data we want,
      *  ie. get the index of the data
      */
     if ((index = SDnametoindex(hdfinfo.sdfid, arr_name)) < 0) {
         printf("hdfio_rd_sd: couldn't find %s in the dataset\n",
                 arr_name);
         return -1;
     }
  
     /*
      *  next, get the SD ID for the data
      */
     if ((sdid = SDselect(hdfinfo.sdfid, index)) < 0) {
         printf("hdfio_rd_sd: Failed in SDselect for item %s\n",
                 arr_name);
         return -1;
     }
  
     /*
      *  now, get the size information of the data item
      */
     if (SDgetinfo(sdid, name, &rank, data_dims, &numbertype,
             &nattrs) < 0) {
         printf("hdfio_rd_sd: Failed in SDgetinfo for item %s\n",
                 arr_name);
         return -1;
     }
  
     /*
      *  compute the read controls for whole array
      */
  
     for (j = 0; j < rank; j++) {
         start[j] = 0;
         edge[j] = data_dims[j];
     }
  
     /*
      *  read the data
      */
     if (SDreaddata(sdid, start, NULL, edge, array) < 0) {
         printf("hdfio_rd_sd: failure to read data for item %s\n", arr_name);
         return -1;
     }
     return 0;
 }
  
 void hdfio_close(hdfio_struc hdfinfo)
 /*******************************************************************
  
    hdfio_close
  
    purpose: close the hdf dataset 
  
    Returns type: void - nothing
  
    Parameters: (in calling order)
       Type              Name            I/O     Description
       ----              ----            ---     -----------
       hdfio_struc       hdfinfo          O      information  struct
                                                 about the opened file
  
    Modification history:
       Programmer        Date            Description of change
       ----------        ----            ---------------------
       W. Robinson, SAIC, 5 Aug 2004     Original development
  
  *******************************************************************/
  {
  
     /*
      *  just close the file up
      */
     Hclose(hdfinfo.fid);
     SDend(hdfinfo.sdfid);
 }

index

an array had not been initialized Several spelling and grammar corrections were which is read from the appropriate MCF the above metadata values were hard coded A problem calculating the average background DN for SWIR bands when the moon is in the space view port was corrected The new algorithm used to calculate the average background DN for all reflective bands when the moon is in the space view port is now the same as the algorithm employed by the thermal bands For non SWIR changes in the averages are typically less than Also for non SWIR the black body DNs remain a backup in case the SV DNs are not available For SWIR the changes in computed averages were larger because the old which used the black body suffered from contamination by the micron leak As a consequence of the if SV DNs are not available for the SWIR the EV pixels will not be the granule time is used to identify the appropriate tables within the set given for one LUT the first two or last two tables respectively will be used for the interpolation If there is only one LUT in the set of it will be treated as a constant LUT The manner in which Earth View data is checked for saturation was changed Previously the raw Earth View DNs and Space View DNs were checked against the lookup table values contained in the table dn_sat The change made is to check the raw Earth and Space View DNs to be sure they are less than the maximum saturation value and to check the Space View subtracted Earth View dns against a set of values contained in the new lookup table dn_sat_ev The metadata configuration and ASSOCIATEDINSTRUMENTSHORTNAME from the MOD02HKM product The same metatdata with extensions and were removed from the MOD021KM and MOD02OBC products ASSOCIATEDSENSORSHORTNAME was set to MODIS in all products These changes are reflected in new File Specification which users may consult for exact the pow functions were eliminated in Emissive_Cal and Emissive bands replaced by more efficient code Other calculations throughout the code were also made more efficient Aside from a few round off there was no difference to the product The CPU time decreased by about for a day case and for a night case A minor bug in calculating the uncertainty index for emissive bands was corrected The frame index(0-based) was previously being used the frame number(1-based) should have been used. There were only a few minor changes to the uncertainty index(maximum of 1 digit). 3. Some inefficient arrays(Sigma_RVS_norm_sq) were eliminated and some code lines in Preprocess_L1A_Data were moved into Process_OBCEng_Emiss. There were no changes to the product. Required RAM was reduced by 20 MB. Now

int j

Definition: decode_rs.h:73

hdfio.h

MDN.setup.name

name

Definition: setup.py:14

NULL

#define NULL

Definition: decode_rs.h:63

hdfio_close

void hdfio_close(hdfio_struc hdfinfo)

Definition: hdfio.c:223