MtkSomXYToBls.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                               MtkSomXYToBls                               =
=                                                                           =
=============================================================================

                         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 MtkSomXYToBls(
  int path,              
  int resolution_meters, 
  double som_x,          
  double som_y,          
  int *block,            
  float *line,           
  float *sample           )
{
  MTKt_status status_code;      /* Return status of this function */
  MTKt_status status;       /* Return status */
  MTKt_MisrProjParam pp;    /* Projection parameters */

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

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

  status = misr_init(pp.nblock, pp.nline, pp.nsample,
             pp.reloffset, pp.ulc, pp.lrc);
  if (status != MTK_SUCCESS)
    MTK_ERR_CODE_JUMP(MTK_MISR_PROJ_INIT_FAILED);

  status = misrfor(som_x, som_y, block, line, sample);
  if (status != MTK_SUCCESS) 
    MTK_ERR_CODE_JUMP(MTK_MISR_FORWARD_PROJ_FAILED);

  return MTK_SUCCESS;

ERROR_HANDLE:
  return status_code;
}