MtkFileCoreMetaDataRaw.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                          MtkFileCoreMetaDataRaw                           =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

            Copyright 2006, California Institute of Technology.
                           ALL RIGHTS RESERVED.
                 U.S. Government Sponsorship acknowledged.

============================================================================*/

#include "MisrFileQuery.h"
#include "MisrError.h"
#include "MisrUtil.h"
#include <hdf.h>

MTKt_status MtkFileCoreMetaDataRaw(
  const char *filename,    
  char **coremeta           )
{
  MTKt_status status;      /* Return status */

  status = MtkFileCoreMetaDataRawNC(filename, coremeta); // try netCDF
  if (status != MTK_NETCDF_OPEN_FAILED) return status;

  return MtkFileCoreMetaDataRawHDF(filename, coremeta); // try HDF
}

MTKt_status MtkFileCoreMetaDataRawNC(
  const char *filename,    
  char **coremeta           )
{
  MTKt_status status_code; /* Return status of this function */
  MTKt_status status;      /* Return status */

  if (filename == NULL) return MTK_NULLPTR;

  /* Open file */
  int ncid = 0;
  {
    int nc_status = nc_open(filename, NC_NOWRITE, &ncid);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_OPEN_FAILED);
  }

  /* Read grid attribute */
  status = MtkFileCoreMetaDataRawNcid(ncid, coremeta);
  MTK_ERR_COND_JUMP(status);

  /* Close file */
  {
    int nc_status = nc_close(ncid);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_CLOSE_FAILED);
  }
  ncid = 0;

  return MTK_SUCCESS;

 ERROR_HANDLE:
  if (ncid != 0) nc_close(ncid);
  return status_code;
}

MTKt_status MtkFileCoreMetaDataRawHDF(
  const char *filename,    
  char **coremeta           )
{
  MTKt_status status_code; /* Return status of this function */
  MTKt_status status;      /* Return status */
  intn hdf_status;
  int32 sds_id = FAIL;

  if (filename == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  /* Open HDF File */
  hdf_status = sds_id = SDstart(filename, DFACC_READ);
  if (hdf_status == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDSTART_FAILED);

  /* Read coremetadata. */
  status = MtkFileCoreMetaDataRawFid(sds_id, coremeta);
  MTK_ERR_COND_JUMP(status);

  /* Close HDF File */
  hdf_status = SDend(sds_id);
  if (hdf_status == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDEND_FAILED);
  sds_id = FAIL;

  return MTK_SUCCESS;

ERROR_HANDLE:
  if (sds_id != FAIL)
    SDend(sds_id);

  return status_code;
}

MTKt_status MtkFileCoreMetaDataRawFid(
  int32 sds_id,            
  char **coremeta           )
{
  MTKt_status status_code; /* Return status of this function */
  intn hdf_status;
  int32 attr_index;
  char attr_name[80];
  int32 data_type;
  int32 count;
  char *attr_buf = NULL;

  if (coremeta == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  /* Find attribute index */
  hdf_status = attr_index = SDfindattr(sds_id, "coremetadata");
  if (hdf_status == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDFINDATTR_FAILED);

  /* Get attribute information */
  hdf_status = SDattrinfo(sds_id, attr_index, attr_name, &data_type, &count);
  if (hdf_status == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDATTRINFO_FAILED);

  /* Allocate Memory */
  attr_buf = (char*)malloc((count + 1) * sizeof(char));
  if (attr_buf == NULL)
    MTK_ERR_CODE_JUMP(MTK_MALLOC_FAILED);

  /* Read attribute */
  hdf_status = SDreadattr(sds_id, attr_index, (VOIDP)attr_buf);
  if (hdf_status == FAIL)
    MTK_ERR_CODE_JUMP(MTK_HDF_SDREADATTR_FAILED);

  attr_buf[count] = '\0';

  *coremeta = attr_buf;

  return MTK_SUCCESS;

ERROR_HANDLE:
  if (attr_buf != NULL)
    free(attr_buf);

  return status_code;
}

MTKt_status MtkFileCoreMetaDataRawNcid(
  int ncid,               
  char **coremeta           )
{
  MTKt_status status_code; /* Return status of this function */
  char *attr_buf = NULL;

  if (coremeta == NULL) {
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);
  }

  const char *group_name = "HDFEOS INFORMATION";
  const char *var_name = "coremetadata";

  {
    int group_id;
    int nc_status = nc_inq_ncid(ncid, group_name, &group_id);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);

    MTKt_ncvarid var;
    int status = MtkNCVarId(group_id, var_name, &var);
    MTK_ERR_COND_JUMP(status);

    int ndims;
    nc_status = nc_inq_varndims(var.gid, var.varid, &ndims);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);

    if (ndims != 1) {
      MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);
    }

    int dimid;
    nc_status = nc_inq_vardimid(var.gid, var.varid, &dimid);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);

    size_t len;
    nc_status = nc_inq_dimlen(group_id, dimid, &len);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);

    /* Allocate Memory */
    attr_buf = (char*)calloc((len + 1), sizeof(char));
    if (attr_buf == NULL)
      MTK_ERR_CODE_JUMP(MTK_MALLOC_FAILED);
    
    nc_status = nc_get_var_text(var.gid, var.varid, attr_buf);
    if (nc_status != NC_NOERR) MTK_ERR_CODE_JUMP(MTK_NETCDF_READ_FAILED);
  }

  *coremeta = attr_buf;

  return MTK_SUCCESS;

ERROR_HANDLE:
  if (attr_buf != NULL)
    free(attr_buf);

  return status_code;
}