mercator.h
上传用户:oybseng
上传日期:2015-04-27
资源大小:7831k
文件大小:9k
源码类别:

GDI/图象编程

开发平台:

Visual C++

  1. #ifndef MERCATOR_H
  2.   #define MERCATOR_H
  3. /***************************************************************************/
  4. /* RSC IDENTIFIER: MERCATOR
  5.  *
  6.  * ABSTRACT
  7.  *
  8.  *    This component provides conversions between Geodetic coordinates 
  9.  *    (latitude and longitude in radians) and Mercator projection coordinates
  10.  *    (easting and northing in meters).
  11.  *
  12.  * ERROR HANDLING
  13.  *
  14.  *    This component checks parameters for valid values.  If an invalid value
  15.  *    is found, the error code is combined with the current error code using 
  16.  *    the bitwise or.  This combining allows multiple error codes to be
  17.  *    returned. The possible error codes are:
  18.  *
  19.  *          MERC_NO_ERROR           : No errors occurred in function
  20.  *          MERC_LAT_ERROR          : Latitude outside of valid range
  21.  *                                      (-89.5 to 89.5 degrees)
  22.  *          MERC_LON_ERROR          : Longitude outside of valid range
  23.  *                                      (-180 to 360 degrees)
  24.  *          MERC_EASTING_ERROR      : Easting outside of valid range
  25.  *                                      (False_Easting +/- ~20,500,000 m,
  26.  *                                       depending on ellipsoid parameters
  27.  *                                       and Origin_Latitude)
  28.  *          MERC_NORTHING_ERROR     : Northing outside of valid range
  29.  *                                      (False_Northing +/- ~23,500,000 m,
  30.  *                                       depending on ellipsoid parameters
  31.  *                                       and Origin_Latitude)
  32.  *          MERC_ORIGIN_LAT_ERROR   : Origin latitude outside of valid range
  33.  *                                      (-89.5 to 89.5 degrees)
  34.  *          MERC_CENT_MER_ERROR     : Central meridian outside of valid range
  35.  *                                      (-180 to 360 degrees)
  36.  *          MERC_A_ERROR            : Semi-major axis less than or equal to zero
  37.  *          MERC_INV_F_ERROR        : Inverse flattening outside of valid range
  38.  *                     (250 to 350)
  39.  *
  40.  * REUSE NOTES
  41.  *
  42.  *    MERCATOR is intended for reuse by any application that performs a 
  43.  *    Mercator projection or its inverse.
  44.  *    
  45.  * REFERENCES
  46.  *
  47.  *    Further information on MERCATOR can be found in the Reuse Manual.
  48.  *
  49.  *    MERCATOR originated from :  U.S. Army Topographic Engineering Center
  50.  *                                Geospatial Information Division
  51.  *                                7701 Telegraph Road
  52.  *                                Alexandria, VA  22310-3864
  53.  *
  54.  * LICENSES
  55.  *
  56.  *    None apply to this component.
  57.  *
  58.  * RESTRICTIONS
  59.  *
  60.  *    MERCATOR has no restrictions.
  61.  *
  62.  * ENVIRONMENT
  63.  *
  64.  *    MERCATOR was tested and certified in the following environments:
  65.  *
  66.  *    1. Solaris 2.5 with GCC, version 2.8.1
  67.  *    2. Windows 95 with MS Visual C++, version 6
  68.  *
  69.  * MODIFICATIONS
  70.  *
  71.  *    Date              Description
  72.  *    ----              -----------
  73.  *    10-02-97          Original Code
  74.  *
  75.  */
  76. /***************************************************************************/
  77. /*
  78.  *                              DEFINES
  79.  */
  80.   #define MERC_NO_ERROR           0x0000
  81.   #define MERC_LAT_ERROR          0x0001
  82.   #define MERC_LON_ERROR          0x0002
  83.   #define MERC_EASTING_ERROR      0x0004
  84.   #define MERC_NORTHING_ERROR     0x0008
  85.   #define MERC_ORIGIN_LAT_ERROR   0x0010
  86.   #define MERC_CENT_MER_ERROR     0x0020
  87.   #define MERC_A_ERROR            0x0040
  88.   #define MERC_INV_F_ERROR        0x0080
  89. /***************************************************************************/
  90. /*
  91.  *                              FUNCTION PROTOTYPES
  92.  *                                for MERCATOR.C
  93.  */
  94. /* ensure proper linkage to c++ programs */
  95.   #ifdef __cplusplus
  96. extern "C" {
  97.   #endif
  98.   long Set_Mercator_Parameters(double a,      
  99.                                double f,
  100.                                double Origin_Latitude,
  101.                                double Central_Meridian,
  102.                                double False_Easting,
  103.                                double False_Northing,
  104.                                double *Scale_Factor);
  105. /*
  106.  * The function Set_Mercator_Parameters receives the ellipsoid parameters and
  107.  * Mercator projcetion parameters as inputs, and sets the corresponding state 
  108.  * variables.  It calculates and returns the scale factor.  If any errors
  109.  * occur, the error code(s) are returned by the function, otherwise 
  110.  * MERC_NO_ERROR is returned.
  111.  *
  112.  *    a                 : Semi-major axis of ellipsoid, in meters   (input)
  113.  *    f                 : Flattening of ellipsoid         (input)
  114.  *    Origin_Latitude   : Latitude in radians at which the          (input)
  115.  *                          point scale factor is 1.0
  116.  *    Central_Meridian  : Longitude in radians at the center of     (input)
  117.  *                          the projection
  118.  *    False_Easting     : A coordinate value in meters assigned to the
  119.  *                          central meridian of the projection.     (input)
  120.  *    False_Northing    : A coordinate value in meters assigned to the
  121.  *                          origin latitude of the projection       (input)
  122.  *    Scale_Factor      : Multiplier which reduces distances in the 
  123.  *                          projection to the actual distance on the
  124.  *                          ellipsoid                               (output)
  125.  */
  126.   void Get_Mercator_Parameters(double *a,
  127.                                double *f,
  128.                                double *Origin_Latitude,
  129.                                double *Central_Meridian,
  130.                                double *False_Easting,
  131.                                double *False_Northing,
  132.                                double *Scale_Factor);
  133. /*
  134.  * The function Get_Mercator_Parameters returns the current ellipsoid
  135.  * parameters, Mercator projection parameters, and scale factor.
  136.  *
  137.  *    a                 : Semi-major axis of ellipsoid, in meters   (output)
  138.  *    f                 : Flattening of ellipsoid         (output)
  139.  *    Origin_Latitude   : Latitude in radians at which the          (output)
  140.  *                          point scale factor is 1.0
  141.  *    Central_Meridian  : Longitude in radians at the center of     (output)
  142.  *                          the projection
  143.  *    False_Easting     : A coordinate value in meters assigned to the
  144.  *                          central meridian of the projection.     (output)
  145.  *    False_Northing    : A coordinate value in meters assigned to the
  146.  *                          origin latitude of the projection       (output)
  147.  *    Scale_Factor      : Multiplier which reduces distances in the 
  148.  *                          projection to the actual distance on the
  149.  *                          ellipsoid                               (output)
  150.  */
  151.   long Convert_Geodetic_To_Mercator (double Latitude,
  152.                                      double Longitude,
  153.                                      double *Easting,
  154.                                      double *Northing); 
  155. /*
  156.  * The function Convert_Geodetic_To_Mercator converts geodetic (latitude and
  157.  * longitude) coordinates to Mercator projection easting, and northing
  158.  * coordinates, according to the current ellipsoid and Mercator projection
  159.  * parameters.  If any errors occur, the error code(s) are returned by the
  160.  * function, otherwise MERC_NO_ERROR is returned.
  161.  *
  162.  *    Latitude          : Latitude (phi) in radians           (input)
  163.  *    Longitude         : Longitude (lambda) in radians       (input)
  164.  *    Easting           : Easting (X) in meters               (output)
  165.  *    Northing          : Northing (Y) in meters              (output)
  166.  */
  167.   long Convert_Mercator_To_Geodetic(double Easting,
  168.                                     double Northing,
  169.                                     double *Latitude,
  170.                                     double *Longitude);
  171. /*
  172.  * The function Convert_Mercator_To_Geodetic converts Mercator projection
  173.  * easting and northing coordinates to geodetic (latitude and longitude)
  174.  * coordinates, according to the current ellipsoid and Mercator projection
  175.  * coordinates.  If any errors occur, the error code(s) are returned by the
  176.  * function, otherwise MERC_NO_ERROR is returned.
  177.  *
  178.  *    Easting           : Easting (X) in meters                  (input)
  179.  *    Northing          : Northing (Y) in meters                 (input)
  180.  *    Latitude          : Latitude (phi) in radians              (output)
  181.  *    Longitude         : Longitude (lambda) in radians          (output)
  182.  */
  183.   #ifdef __cplusplus
  184. }
  185.   #endif
  186. #endif /* MERCATOR_H */