MtkSomXYToLatLonAry.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                           MtkSomXYToLatLonAry                             =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

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

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

#include "MisrCoordQuery.h"
#include "MisrUnitConv.h"
#include "MisrError.h"
#include "proj.h"
#include "gctp_prototypes.h"
#include <stddef.h>

MTKt_status MtkSomXYToLatLonAry(
  int path,            
  int nelement,        
  const double *som_x, 
  const double *som_y, 
  double *lat_dd,      
  double *lon_dd        )
{
  MTKt_status status_code;      /* Return status of this function */
  MTKt_status status;       /* Return status */
  MTKt_MisrProjParam pp;    /* Projection parameters */
  int iflg;         /* GCTP status flag */
  double lat_r;         /* Latitude in radians */
  double lon_r;         /* Longitude in radians */
  int i;            /* Loop index */
  int (*inv_trans[MAXPROJ+1])(double, double, double*, double*); /* Array of function ptrs (Not used) */

  if (som_x == NULL || som_y == NULL || lat_dd == NULL || lon_dd == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  if (nelement < 0)
    MTK_ERR_CODE_JUMP(MTK_BAD_ARGUMENT);

  status = MtkPathToProjParam(path, 0, &pp);
  MTK_ERR_COND_JUMP(status);

  inv_init(pp.projcode, pp.zonecode, pp.projparam, pp.spherecode,
       NULL, NULL, &iflg, inv_trans);
  if (iflg)
    MTK_ERR_CODE_JUMP(MTK_GCTP_INVERSE_INIT_FAILED);

  for (i = 0; i < nelement; i++) {
    iflg = sominv(som_x[i], som_y[i], &lon_r, &lat_r);
    if (iflg) 
      MTK_ERR_CODE_JUMP(MTK_GCTP_INVERSE_PROJ_FAILED);

    status = MtkRadToDd(lat_r, &lat_dd[i]);
    MTK_ERR_COND_JUMP(status);

    status = MtkRadToDd(lon_r, &lon_dd[i]);
    MTK_ERR_COND_JUMP(status);
  }

  return MTK_SUCCESS;

ERROR_HANDLE:
  return status_code;
}