MisrMapQuery.h

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                              MisrMapQuery                                 =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

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

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

#ifndef MISRMAPQUERY_H
#define MISRMAPQUERY_H

#include "MisrError.h"
#include "MisrProjParam.h"
#include "MisrUtil.h"

typedef struct {
  double lat;           
  double lon;           
} MTKt_GeoCoord;

#define MTKT_GEOCOORD_INIT { 0.0, 0.0 }

typedef struct {
  double x;         
  double y;         
} MTKt_SomCoord;

#define MTKT_SOMCOORD_INIT { 0.0, 0.0 }

typedef struct {
  MTKt_GeoCoord ulc;        
  MTKt_GeoCoord urc;        
  MTKt_GeoCoord ctr;        
  MTKt_GeoCoord lrc;        
  MTKt_GeoCoord llc;        
} MTKt_GeoRegion;

#define MTKT_GEOREGION_INIT { MTKT_GEOCOORD_INIT, MTKT_GEOCOORD_INIT, \
                              MTKT_GEOCOORD_INIT, MTKT_GEOCOORD_INIT, \
                              MTKT_GEOCOORD_INIT }

typedef struct {
  int path;         
  MTKt_SomCoord ulc;        
  MTKt_SomCoord ctr;        
  MTKt_SomCoord lrc;        
} MTKt_SomRegion;

#define MTKT_SOMREGION_INIT { 0, MTKT_SOMCOORD_INIT, MTKT_SOMCOORD_INIT, \
                              MTKT_SOMCOORD_INIT }

typedef struct {
  int path;         
  int start_block;      
  int end_block;        
  int resolution;       
  int resfactor;        
  int nline;            
  int nsample;          
  MTKt_boolean pixelcenter; 
  MTKt_SomRegion som;       
  MTKt_GeoRegion geo;       
  MTKt_MisrProjParam pp;    
} MTKt_MapInfo;

#define MTKT_MAPINFO_INIT { 0, 0, 0, 0, 0, 0, 0, MTK_TRUE, \
                            MTKT_SOMREGION_INIT, MTKT_GEOREGION_INIT, \
                            MTKT_MISRPROJPARAM_INIT }

typedef enum { 
  MTKe_ORIGIN_UL = 0,    
  MTKe_ORIGIN_UR = 1,    
  MTKe_ORIGIN_LL = 2,    
  MTKe_ORIGIN_LR = 3     
} MTKt_OriginCode;

typedef enum {
  MTKe_PIX_REG_CENTER = 0,  
  MTKe_PIX_REG_CORNER = 1   
} MTKt_PixRegCode;

typedef struct {
  double min_x;         
  double min_y;         
  double max_x;         
  double max_y;         
  int size_line;        
  int size_sample;      
  double resolution_x;          
  double resolution_y;          
  double tline[4];      
  double tsample[4];        
  MTKt_OriginCode origin_code;  
  MTKt_PixRegCode pix_reg_code; 
} MTKt_GenericMapInfo;


#define MTKT_GENERICMAPINFO_INIT { 0, 0, 0, 0, 0, 0, 0, 0, {0,0,0}, {0,0,0}, MTKe_ORIGIN_UL, MTKe_PIX_REG_CENTER }

typedef struct {
  int proj_code;
  int sphere_code;
  int zone_code;
  double proj_param[15];
} MTKt_GCTPProjInfo;

#define MTKT_GCTPPROJINFO_INIT { 0, 0, 0, {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0} }

MTKt_status MtkLSToSomXY( MTKt_MapInfo mapinfo,
              float line,
              float sample,
              double *som_x,
              double *som_y );

MTKt_status MtkLSToSomXYAry( MTKt_MapInfo mapinfo,
                 int nelement,
                 const float *line,
                 const float *sample,
                 double *som_x,
                 double *som_y );

MTKt_status MtkSomXYToLS( MTKt_MapInfo mapinfo,
              double som_x,
              double som_y,
              float *line,
              float *sample );

MTKt_status MtkSomXYToLSAry( MTKt_MapInfo mapinfo,
                 int nelement,
                 const double *som_x,
                 const double *som_y,
                 float *line,
                 float *sample );

MTKt_status MtkLatLonToLS( MTKt_MapInfo mapinfo,
               double lat_dd,
               double lon_dd,
               float *line,
               float *sample );

MTKt_status MtkLatLonToLSAry( MTKt_MapInfo mapinfo,
                  int nelement,
                  const double *lat_dd,
                  const double *lon_dd,
                  float *line,
                  float *sample );

MTKt_status MtkLSToLatLon( MTKt_MapInfo mapinfo,
               float line,
               float sample,
               double *lat_dd,
               double *lon_dd );

MTKt_status MtkLSToLatLonAry( MTKt_MapInfo mapinfo,
                  int nelement,
                  const float *line,
                  const float *sample,
                  double *lat_dd,
                  double *lon_dd );

MTKt_status MtkCreateLatLon( MTKt_MapInfo mapinfo,
                 MTKt_DataBuffer *latbuf,
                 MTKt_DataBuffer *lonbuf );

MTKt_status MtkGenericMapInfo(
  double Min_x, 
  double Min_y,  
  double Resolution_x, 
  double Resolution_y, 
  int Number_pixel_x,  
  int Number_pixel_y,  
  MTKt_OriginCode Origin_code, 
  MTKt_PixRegCode Pix_reg_code, 
  MTKt_GenericMapInfo *Map_info );

MTKt_status MtkGCTPProjInfo(
  int Proj_code,  
  int Sphere_code,
  int Zone_code,  
  double Proj_param[15],
  MTKt_GCTPProjInfo *Proj_info );

MTKt_status MtkGCTPCreateLatLon(
  const MTKt_GenericMapInfo *Map_info, 
  const MTKt_GCTPProjInfo *Proj_info,
  MTKt_DataBuffer *Latitude, 
  MTKt_DataBuffer *Longitude );

MTKt_status MtkGenericMapInfoRead(
  const char *Filename,  
  MTKt_GenericMapInfo *Map_info 
);

MTKt_status MtkGCTPProjInfoRead(
  const char *Filename,  
  MTKt_GCTPProjInfo *Proj_info 
);

MTKt_status MtkChangeMapResolution(
  const MTKt_MapInfo *Map_info_in,  
  int Resolution,            
  MTKt_MapInfo *Map_info_out  
);

#endif /* MISRMAPQUERY_H */