MtkBlsToLatLonAry.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                            MtkBlsToLatLonAry                              =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

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

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

#include "MisrCoordQuery.h"
#include "MisrError.h"
#include "misrproj.h"

MTKt_status MtkBlsToLatLonAry(
  int path,              
  int resolution_meters, 
  int nelement,          
  const int *block,      
  const float *line,     
  const float *sample,   
  double *lat_dd,        
  double *lon_dd          )
{
  MTKt_status status_code;      /* Return status of this function */
  MTKt_status status;       /* Return status */
  double som_x;         /* SOM X */
  double som_y;         /* SOM Y */
  int i;            /* Loop index */

  if (block == NULL || line == NULL || sample == NULL ||
      lat_dd == NULL || lon_dd == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  if (nelement < 0)
    MTK_ERR_CODE_JUMP(MTK_BAD_ARGUMENT);

  for (i = 0; i < nelement; i++) {
    status = MtkBlsToSomXY(path, resolution_meters,
               block[i], line[i], sample[i], &som_x, &som_y);
    MTK_ERR_COND_JUMP(status);

    status = MtkSomXYToLatLon(path, som_x, som_y, &lat_dd[i], &lon_dd[i]);
    MTK_ERR_COND_JUMP(status);
  }

  return MTK_SUCCESS;

ERROR_HANDLE:
  return status_code;
}