MtkSetRegionByPathSomUlcLrc.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                       MtkSetRegionByPathSomUlcLrc                         =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

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

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

#include "MisrSetRegion.h"
#include "MisrCoordQuery.h"
#include "MisrError.h"
#include "MisrUtil.h"
#include "MisrOrbitPath.h"

MTKt_status MtkSetRegionByPathSomUlcLrc(
  int path,         
  double ulc_som_x, 
  double ulc_som_y, 
  double lrc_som_x, 
  double lrc_som_y, 
  MTKt_Region *region  )
{
  MTKt_status status;       /* Return status */
  MTKt_status status_code;      /* Return status of this function */
  MTKt_Region rgn;      /* Region structure */
  double ctr_som_x;     /* Center of region in SOM X */
  double ctr_som_y;     /* Center of region in SOM Y */

  if (region == NULL)
    MTK_ERR_CODE_JUMP(MTK_NULLPTR);

  /* Check path bounds */
  if (path < 1 || path > 233)
    MTK_ERR_CODE_JUMP(MTK_OUTBOUNDS);

  /* Check SOM X bounds */
  if (lrc_som_x < ulc_som_x)
    MTK_ERR_CODE_JUMP(MTK_OUTBOUNDS);

  /* Check SOM Y bounds */
  if (lrc_som_y < ulc_som_y)
    MTK_ERR_CODE_JUMP(MTK_OUTBOUNDS);

  /* Determine center of region in lat/lon coordinates */

  ctr_som_x = ulc_som_x + ((lrc_som_x - ulc_som_x) / 2.0);
  ctr_som_y = ulc_som_y + ((lrc_som_y - ulc_som_y) / 2.0);

  status = MtkSomXYToLatLon(path, ctr_som_x, ctr_som_y,
                &rgn.geo.ctr.lat, &rgn.geo.ctr.lon);
  MTK_ERR_COND_JUMP(status);

  /* Half of the over all extent to measure from center */
  rgn.hextent.xlat = (lrc_som_x - ulc_som_x + MAXRESOLUTION) / 2.0;
  rgn.hextent.ylon = (lrc_som_y - ulc_som_y + MAXRESOLUTION) / 2.0;

  *region = rgn;

  return MTK_SUCCESS;

ERROR_HANDLE:
  return status_code;
}