MtkLatLonToPathList.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                           MtkLatLonToPathList                             =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

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

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

#include "MisrOrbitPath.h"
#include "MisrCoordQuery.h"
#include "MisrError.h"
#include "MisrProjParam.h"
#include <stdlib.h>

MTKt_status MtkLatLonToPathList(
  double lat_dd, 
  double lon_dd, 
  int *pathcnt,  
  int **pathlist  )
{
  MTKt_status status_code;      /* Return status of this function */
  MTKt_status status;           /* Return status */
  int block;            /* Block */
  float line;           /* Line */
  float sample;         /* Sample */
  int path;         /* Loop index */
  int pathlist_tmp[NPATH];  /* Temp pathlist */
  int *list = NULL;     /* Temp list */
  int cnt = 0;          /* Path count */

  if (pathcnt == NULL || pathlist == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  for (path = 1; path <= NPATH; path++) {
    status = MtkLatLonToBls(path, MAXRESOLUTION, lat_dd, lon_dd, 
           &block, &line, &sample);
    
    if (status == MTK_SUCCESS)
      pathlist_tmp[cnt++] = path;
  }

  if (cnt == 0)
    MTK_ERR_CODE_JUMP(MTK_NOT_FOUND);

  list = (int *)malloc(cnt * sizeof(int));
  if (list == NULL)
    MTK_ERR_CODE_JUMP(MTK_MALLOC_FAILED);

  for (path = 0; path < cnt; path++) {
    list[path] = pathlist_tmp[path];
  }
  *pathcnt = cnt;
  *pathlist = list;

  return MTK_SUCCESS;

ERROR_HANDLE:
  return status_code;
}