MtkReadData.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                               MtkReadData                                 =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

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

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

#include "MisrReadData.h"
#include "MisrFileQuery.h"
#include "MisrError.h"
#include <stdlib.h>
#include <string.h>
#include <hdf.h>
#include <HdfEosDef.h>
#include <netcdf.h>

MTKt_status MtkReadData(
  const char *filename,     
  const char *gridname,     
  const char *fieldname,    
  MTKt_Region region,       
  MTKt_DataBuffer *databuf, 
  MTKt_MapInfo *mapinfo      )
{
  MTKt_status status;       /* Return status */

  status = MtkReadDataNC(filename, gridname, fieldname, region, databuf, mapinfo);  // try netCDF
  if (status != MTK_NETCDF_OPEN_FAILED) return status;

  return MtkReadDataHDF(filename, gridname, fieldname, region, databuf, mapinfo);  // try HDF
}

MTKt_status MtkReadDataHDF(
  const char *filename,     
  const char *gridname,     
  const char *fieldname,    
  MTKt_Region region,       
  MTKt_DataBuffer *databuf, 
  MTKt_MapInfo *mapinfo      )
{
  MTKt_status status_code;  /* Return status code for error macros */
  MTKt_status status;       /* Return status */
  int32 fid = FAIL;     /* HDF-EOS File id */
  intn hdfstatus;       /* HDF return status */

  if (filename == NULL) return MTK_NULLPTR;

  /* Open file. */
  fid = GDopen((char*)filename, DFACC_READ);
  if (fid == FAIL) MTK_ERR_CODE_JUMP(MTK_HDFEOS_GDOPEN_FAILED);

  /* Read data. */
  status = MtkReadDataFid(fid, gridname, fieldname, region, databuf, mapinfo);
  MTK_ERR_COND_JUMP(status);

  /* Close file. */
  hdfstatus = GDclose(fid);
  if (hdfstatus == FAIL) MTK_ERR_CODE_JUMP(MTK_HDFEOS_GDCLOSE_FAILED);
  fid = FAIL;

  return MTK_SUCCESS;

 ERROR_HANDLE:
  if (fid != FAIL) GDclose(fid);
  return status_code;
}

MTKt_status MtkReadDataNC(
  const char *filename,     
  const char *gridname,     
  const char *fieldname,    
  MTKt_Region region,       
  MTKt_DataBuffer *databuf, 
  MTKt_MapInfo *mapinfo      )
{
  MTKt_status status_code;  /* Return status code for error macros */
  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 data. */
  status = MtkReadDataNcid(ncid, gridname, fieldname, region, databuf, mapinfo);
  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;
  /* Close file. */

  return MTK_SUCCESS;

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

MTKt_status MtkReadDataFid(
  int32 fid,                
  const char *gridname,     
  const char *fieldname,    
  MTKt_Region region,       
  MTKt_DataBuffer *databuf, 
  MTKt_MapInfo *mapinfo      )
{
  MTKt_status status_code;  /* Return status code for error macros */
  MTKt_status status;       /* Return status */
  MTKt_FileType filetype;   /* File type */

  /* Determine file type */
  status = MtkFileTypeFid(fid, &filetype);
  MTK_ERR_COND_JUMP(status);

  switch (filetype) {

  case MTK_AGP:
  case MTK_GP_GMP:
  case MTK_GRP_RCCM:
  case MTK_PP:
    status = MtkReadRawFid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_GRP_ELLIPSOID_GM:
  case MTK_GRP_ELLIPSOID_LM:
  case MTK_GRP_TERRAIN_GM:
  case MTK_GRP_TERRAIN_LM:
    status = MtkReadL1B2Fid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_AS_AEROSOL:
    status = MtkReadRawFid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_AS_LAND:
    status = MtkReadL2LandFid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_TC_ALBEDO:
    status = MtkReadRawFid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_TC_CLASSIFIERS:
    status = MtkReadRawFid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_TC_STEREO:
    status = MtkReadRawFid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_TC_CLOUD:
    status = MtkReadL2TCCloudFid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_CMV_NRT:
    status = MtkReadRawFid(fid, gridname, fieldname, region,
               databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_HR_BRF:
  case MTK_HR_RPV:
  case MTK_HR_TIP:
  case MTK_CONVENTIONAL:
    status = MtkReadConvFid(fid, gridname, fieldname, region,
                databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_UNKNOWN:
  default:
    return MTK_FILETYPE_NOT_SUPPORTED;
  }

  return MTK_SUCCESS;

 ERROR_HANDLE:

  return status_code;
}

MTKt_status MtkReadDataNcid(
  int ncid,                
  const char *gridname,     
  const char *fieldname,    
  MTKt_Region region,       
  MTKt_DataBuffer *databuf, 
  MTKt_MapInfo *mapinfo      )
{
  MTKt_status status_code;  /* Return status code for error macros */
  MTKt_status status;       /* Return status */
  MTKt_FileType filetype;   /* File type */

  /* Determine file type */
  status = MtkFileTypeNcid(ncid, &filetype);
  MTK_ERR_COND_JUMP(status);

  switch (filetype) {

  case MTK_AS_AEROSOL:
    status = MtkReadRawNcid(ncid, gridname, fieldname, region, databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_AS_LAND:
    status = MtkReadL2LandNcid(ncid, gridname, fieldname, region, databuf, mapinfo);
    MTK_ERR_COND_JUMP(status);
    break;

  case MTK_UNKNOWN:
  default:
    return MTK_FILETYPE_NOT_SUPPORTED;
  }

  return MTK_SUCCESS;

 ERROR_HANDLE:

  return status_code;
}