MISR-Toolkit/html/_misr_toolkit_test_8c_source.html

 /*===========================================================================
 =                                                                           =
 =                                  biz                                      =
 =                                                                           =
 =============================================================================

                          Jet Propulsion Laboratory
                                    MISR
                                MISR Toolkit

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

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

 #include "MisrToolkit.h"
 #include "MisrError.h"
 #include <stdio.h>
 #include <stdlib.h>
 #include <math.h>


 int main( ) {

   MTKt_status status;       /* For routine calls */
   MTKt_status status_code;  /* For error handler */
   char *err_msg[] = MTK_ERR_DESC; /* Error message descriptions */

   double lat_dd = 32.2, lon_dd = -114.5;
   double lat, lon;
   double lat_extent = 200, lon_extent = 100;
   double somx, somy;
   int path, res;
   float line, sample;
   int block;
   float blk_line, blk_sample;
   int l, s;
   int latdeg, londeg, latmin, lonmin;
   double latsec, lonsec;
   MTKt_Region region = MTKT_REGION_INIT;
   MTKt_TimeMetaData timemeta = MTKT_TIME_METADATA_INIT;
   char datetime[MTKd_DATETIME_LEN];
   char filename[200];
   char gridname[200];
   char fieldname[200];
   MTKt_DataBuffer brf_databuf = MTKT_DATABUFFER_INIT;
   MTKt_MapInfo brf_mapinfo = MTKT_MAPINFO_INIT;
   MTKt_DataBuffer hdrf_databuf = MTKT_DATABUFFER_INIT;
   MTKt_MapInfo hdrf_mapinfo = MTKT_MAPINFO_INIT;


   /* ---------------------------------------------------------- */
   /* Set a region given center latitude/longitude and an extent */
   /* ---------------------------------------------------------- */

   status = MtkSetRegionByLatLonExtent(lat_dd, lon_dd,
                       lat_extent, lon_extent, "km", &region);
   MTK_ERR_COND_JUMP(status);

   printf("region center lat/lon (dd) = (%f, %f)\n",
      region.geo.ctr.lat, region.geo.ctr.lon);
   printf("         region extent (m) = (%f, %f)\n",
           region.hextent.xlat * 2.0, region.hextent.ylon * 2.0);

   /* --------------------------------------------------------- */
   /* Read data in the region from filename/gridname/fieldname, */
   /* do some coordinate conversions, query data buffer and     */
   /* get pixel time                                            */
   /* --------------------------------------------------------- */

   //strcpy(filename, getenv("MTKHOME"));
   strcpy(filename, "");
   strcat(filename, "C:\\Mtk_testdata\\in\\");
   strcat(filename, "MISR_AM1_GRP_ELLIPSOID_GM_P037_O029058_AA_F03_0024.hdf");
   strcpy(gridname, "RedBand");
   strcpy(fieldname, "Red Brf");

   printf("Filename is: %s\n", filename);
   status = MtkReadData(filename, gridname, fieldname, region,
                &brf_databuf, &brf_mapinfo);
   MTK_ERR_COND_JUMP(status);

   status = MtkLatLonToLS(brf_mapinfo, lat_dd, lon_dd, &line, &sample);
   MTK_ERR_COND_JUMP(status);

   l = (int)round(line);
   s = (int)round(sample);

   status = MtkFileToPath(filename, &path);
   MTK_ERR_COND_JUMP(status);

   status = MtkFileGridToResolution(filename, gridname, &res);
   MTK_ERR_COND_JUMP(status);

   status = MtkLatLonToBls(path, res, lat_dd, lon_dd,
               &block, &blk_line, &blk_sample);
   MTK_ERR_COND_JUMP(status);

   status = MtkLSToLatLon(brf_mapinfo, line, sample, &lat, &lon);
   MTK_ERR_COND_JUMP(status);

   status = MtkDdToDegMinSec(lat, &latdeg, &latmin, &latsec);
   MTK_ERR_COND_JUMP(status);

   status = MtkDdToDegMinSec(lon, &londeg, &lonmin, &lonsec);
   MTK_ERR_COND_JUMP(status);

   status = MtkLSToSomXY(brf_mapinfo, line, sample, &somx, &somy);
   MTK_ERR_COND_JUMP(status);

   status = MtkTimeMetaRead(filename, &timemeta);
   MTK_ERR_COND_JUMP(status);

   status = MtkPixelTime(timemeta, somx, somy, datetime);
   MTK_ERR_COND_JUMP(status);

   printf("\n              Reading file = %s\n",filename);
   printf("              Reading grid = %s\n",gridname);
   printf("             Reading field = %s\n",fieldname);
   printf("   number of lines/samples = (%d, %d)\n",
      brf_mapinfo.nline, brf_mapinfo.nsample);
   printf("        input lat/lon (dd) = (%f, %f)\n", lat_dd, lon_dd);
   printf("               line/sample = (%f, %f)\n", line, sample);
   printf("             Brf[%d][%d] = %f\n", l, s, brf_databuf.data.f[l][s]),
   printf("         block/line/sample = (%d, %f, %f)\n",
      block, blk_line, blk_sample);
   printf("              lat/lon (dd) = (%f, %f)\n", lat, lon);
   printf("             lat/lon (dms) = (%d:%d:%f, %d:%d:%f)\n",
      latdeg, latmin, latsec, londeg, lonmin, lonsec);
   printf("                   SOM x/y = (%f, %f)\n", somx, somy);
   printf("      pixel time[%d][%d] = %s\n", l, s, datetime);

   MtkDataBufferFree(&brf_databuf);

   /* --------------------------------------------------------- */
   /* Read data in the region from filename/gridname/fieldname  */
   /* with extra dimensions, do some coordinate conversions and */
   /* query data buffer                                         */
   /* --------------------------------------------------------- */

   //strcpy(filename, getenv("MTKHOME"));
   strcpy(filename,"");
   strcat(filename, "C:\\Mtk_testdata\\in\\");
   strcat(filename, "MISR_AM1_AS_LAND_P037_O029058_F06_0017.hdf");
   strcpy(gridname, "SubregParamsLnd");
   strcpy(fieldname, "LandHDRF[3][4]");

   status = MtkReadData(filename, gridname, fieldname, region,
                &hdrf_databuf, &hdrf_mapinfo);
   MTK_ERR_COND_JUMP(status);

   status = MtkLatLonToLS(hdrf_mapinfo, lat_dd, lon_dd, &line, &sample);
   MTK_ERR_COND_JUMP(status);

   l = (int)round(line);
   s = (int)round(sample);

   status = MtkFileToPath(filename, &path);
   MTK_ERR_COND_JUMP(status);

   status = MtkFileGridToResolution(filename, gridname, &res);
   MTK_ERR_COND_JUMP(status);

   status = MtkLatLonToBls(path, res, lat_dd, lon_dd,
               &block, &blk_line, &blk_sample);
   MTK_ERR_COND_JUMP(status);

   status = MtkLSToLatLon(hdrf_mapinfo, line, sample, &lat, &lon);
   MTK_ERR_COND_JUMP(status);

   status = MtkDdToDegMinSec(lat, &latdeg, &latmin, &latsec);
   MTK_ERR_COND_JUMP(status);

   status = MtkDdToDegMinSec(lon, &londeg, &lonmin, &lonsec);
   MTK_ERR_COND_JUMP(status);

   status = MtkLatLonToSomXY(path, lat_dd, lon_dd, &somx, &somy);
   MTK_ERR_COND_JUMP(status);

   printf("\n              Reading file = %s\n",filename);
   printf("              Reading grid = %s\n",gridname);
   printf("             Reading field = %s\n",fieldname);
   printf("   number of lines/samples = (%d, %d)\n",
      hdrf_mapinfo.nline, hdrf_mapinfo.nsample);
   printf("        input lat/lon (dd) = (%f, %f)\n", lat_dd, lon_dd);
   printf("               line/sample = (%f, %f)\n", line, sample);
   printf("              HDRF[%d][%d] = %f\n", l, s, hdrf_databuf.data.f[l][s]),
   printf("         block/line/sample = (%d, %f, %f)\n",
      block, blk_line, blk_sample);
   printf("              lat/lon (dd) = (%f, %f)\n", lat, lon);
   printf("             lat/lon (dms) = (%d:%d:%f, %d:%d:%f)\n",
      latdeg, latmin, latsec, londeg, lonmin, lonsec);
   printf("                   SOM x/y = (%f, %f)\n", somx, somy);

   MtkDataBufferFree(&hdrf_databuf);

   return MTK_SUCCESS;
  ERROR_HANDLE:
   printf("Error: %s\n",err_msg[status_code]);
   MtkDataBufferFree(&hdrf_databuf);
   return status_code;
 }
MtkLatLonToBls
MTKt_status MtkLatLonToBls(int path, int resolution_meters, double lat_dd, double lon_dd, int *block, float *line, float *sample)
Convert decimal degrees latitude and longitude to block, line, sample.
Definition: MtkLatLonToBls.c:34

MisrToolkit.h

MTKT_TIME_METADATA_INIT
#define MTKT_TIME_METADATA_INIT
Definition: MisrFileQuery.h:91

filename
char * filename
Definition: cdjpeg.h:133

MTKt_DataBuffer::data
MTKt_DataBufferType data
Definition: MisrUtil.h:104

MTKT_MAPINFO_INIT
#define MTKT_MAPINFO_INIT
Definition: MisrMapQuery.h:79

MtkLSToLatLon
MTKt_status MtkLSToLatLon(MTKt_MapInfo mapinfo, float line, float sample, double *lat_dd, double *lon_dd)
Convert line, sample to decimal degrees latitude and longitude.
Definition: MtkLSToLatLon.c:33

MTKt_Region::geo
MTKt_GeoCenter geo
Definition: MisrSetRegion.h:42

MTKT_REGION_INIT
#define MTKT_REGION_INIT
Definition: MisrSetRegion.h:46

MtkLSToSomXY
MTKt_status MtkLSToSomXY(MTKt_MapInfo mapinfo, float line, float sample, double *som_x, double *som_y)
Convert line, sample to SOM X, SOM Y.
Definition: MtkLSToSomXY.c:33

MtkLatLonToLS
MTKt_status MtkLatLonToLS(MTKt_MapInfo mapinfo, double lat_dd, double lon_dd, float *line, float *sample)
Convert decimal degrees latitude and longitude to line, sample.
Definition: MtkLatLonToLS.c:33

MTKt_GeoCoord::lon
double lon
Definition: MisrMapQuery.h:27

MTKt_MapInfo
Map Information.
Definition: MisrMapQuery.h:65

MTKt_DataBuffer
2-dimensional Data Buffer
Definition: MisrUtil.h:98

MTKt_DataBufferType::f
MTKt_float ** f
Definition: MisrUtil.h:93

MTKt_GeoCenter::ctr
MTKt_GeoCoord ctr
Definition: MisrSetRegion.h:27

MTKt_GeoCoord::lat
double lat
Definition: MisrMapQuery.h:26

MTKT_DATABUFFER_INIT
#define MTKT_DATABUFFER_INIT
Definition: MisrUtil.h:109

MTK_SUCCESS
Definition: MisrError.h:12

MtkDataBufferFree
MTKt_status MtkDataBufferFree(MTKt_DataBuffer *databuf)
Free data buffer.
Definition: MtkDataBufferFree.c:33

MTKt_Extent::xlat
double xlat
Definition: MisrSetRegion.h:34

MTK_ERR_DESC
#define MTK_ERR_DESC
Definition: MisrError.h:91

main
int main()
Definition: MisrToolkitTest.c:25

MTKt_Region::hextent
MTKt_Extent hextent
Definition: MisrSetRegion.h:43

MtkReadData
MTKt_status MtkReadData(const char *filename, const char *gridname, const char *fieldname, MTKt_Region region, MTKt_DataBuffer *databuf, MTKt_MapInfo *mapinfo)
Reads any grid/field from any MISR product file and performs unpacking or unscaling. It also reads any MISR conventional product file.
Definition: MtkReadData.c:62

MtkSetRegionByLatLonExtent
MTKt_status MtkSetRegionByLatLonExtent(double ctr_lat_dd, double ctr_lon_dd, double lat_extent, double lon_extent, const char *extent_units, MTKt_Region *region)
Select region by latitude, longitude in decimal degrees, and extent in specified units of degrees...
Definition: MtkSetRegionByLatLonExtent.c:77

MisrError.h

MtkTimeMetaRead
MTKt_status MtkTimeMetaRead(const char *filename, MTKt_TimeMetaData *time_metadata)
Read time metadata from L1B2 Ellipsoid product file.
Definition: MtkTimeMetaRead.c:39

MTKt_TimeMetaData
Time Metadata.
Definition: MisrFileQuery.h:77

MTKt_MapInfo::nsample
int nsample
Definition: MisrMapQuery.h:72

MtkPixelTime
MTKt_status MtkPixelTime(MTKt_TimeMetaData time_metadata, double som_x, double som_y, char pixel_time[MTKd_DATETIME_LEN])
Given SOM Coordinates compute pixel time.
Definition: MtkPixelTime.c:41

MtkDdToDegMinSec
MTKt_status MtkDdToDegMinSec(double dd, int *deg, int *min, double *sec)
Convert decimal degrees to unpacked degrees, minutes, seconds.
Definition: MtkDdToDegMinSec.c:32

MTK_ERR_COND_JUMP
#define MTK_ERR_COND_JUMP(code)
Definition: MisrError.h:188

MTKt_status
MTKt_status
Definition: MisrError.h:11

MTKt_MapInfo::nline
int nline
Definition: MisrMapQuery.h:71

MTKt_Region
Region of interest.
Definition: MisrSetRegion.h:41

MtkFileToPath
MTKt_status MtkFileToPath(const char *filename, int *path)
Read path number from file.
Definition: MtkFileToPath.c:35

MTKd_DATETIME_LEN
#define MTKd_DATETIME_LEN
Definition: MisrUtil.h:149

MTKt_Extent::ylon
double ylon
Definition: MisrSetRegion.h:35

MtkLatLonToSomXY
MTKt_status MtkLatLonToSomXY(int path, double lat_dd, double lon_dd, double *som_x, double *som_y)
Convert decimal degrees latitude and longitude to SOM X, SOM Y.
Definition: MtkLatLonToSomXY.c:38

MtkFileGridToResolution
MTKt_status MtkFileGridToResolution(const char *filename, const char *gridname, int *resolution)
Get resolution of a particular grid.
Definition: MtkFileGridToResolution.c:35