mapproj.h

Go to the documentation of this file.

/**
 * \file mi32/mapproj.h
 * \brief Definitions for map projection transformations
 *
 * \if NODOC
 * $Id: mapproj.h_v 1.141 2004/10/15 15:44:23 mju Exp $
 *
 * $Log: mapproj.h_v $
 * Revision 1.141  2004/10/15 15:44:23  mju
 * Remove unused latlonScaleToMeters.
 *
 * Revision 1.140  2004/08/23 17:33:36  mju
 * Define deprecated stuff when compiling DLL.
 *
 * Revision 1.139  2004/08/23 16:34:00  mju
 * Remove legacy projsys, ogcwkt, etc stuff.
 *
 * Revision 1.138  2004/07/20 20:58:33  mju
 * Accidentally removed OGC fns when still in use.
 *
 * Revision 1.137  2004/07/20 17:57:26  mju
 * Remove some obsolete methods.
 *
 * Revision 1.136  2004/07/19 17:12:29  mju
 * Remove epsg fns.
 *
 * Revision 1.135  2004/07/13 13:08:07  mju
 * Remove another obsolete prototype.
 *
 * Revision 1.134  2004/07/13 13:05:49  mju
 * Remove computeObjProjTransAffine.
 *
 * Revision 1.133  2004/04/28 16:35:17  mju
 * Add projsys_Colombia.
 *
 * Revision 1.132  2004/01/28 22:12:25  scowan
 * Moved and exported map info georeference functions.
 *
 * Revision 1.131  2003/10/01 22:13:27  dwilliss
 * doxygen
 *
 * Revision 1.130  2003/09/25 21:29:07  dwilliss
 * Don't capitalize enum any more. Genitor needed it, doxygen doesn't like it.
 *
 * Revision 1.129  2003/09/15 13:49:56  fileserver!dwilliss
 * Doxygen
 *
 * Revision 1.128  2003/09/08 21:04:47  dwilliss
 * Added doxygen start/end tags for functions gruops
 *
 * Revision 1.126  2003/02/28 17:33:59  mju
 * Add methods to test if two MPPs are equivalent
 *
 * Revision 1.125  2003/02/27 22:08:06  mju
 * Add Goode projection.
 *
 * Revision 1.123  2003/02/26 18:07:43  mju
 * Add getProjNameFull.
 *
 * Revision 1.122  2003/02/24 17:15:20  mju
 * Add datum_tokyo_grid.
 *
 * Revision 1.121  2003/02/20 18:25:53  mju
 * Add DATUM_TRANS_Grid.
 *
 * Revision 1.120  2003/02/04 18:42:17  mju
 * Add dealulpiscului datum.
 *
 * Revision 1.119  2003/02/03 20:17:01  mju
 * Add projection_ObliqueStereographic.
 * Add projsys_Italy.
 *
 * Revision 1.117  2003/01/21 14:40:37  mju
 * Add PROJSYS_Romania.
 *
 * Revision 1.116  2002/11/13 15:54:19  mju
 * When set parms for equidistant conic set both std parallels in single case.
 *
 * Revision 1.115  2002/11/05 17:46:37  mju
 * Add SetProjectionSwissConformal.
 *
 * Revision 1.114  2002/09/09 14:11:06  mju
 * Change comments re rvcobjname.
 *
 * Revision 1.113  2002/08/28 16:09:04  mju
 * In SetSystemLatLon don't change datum.
 *
 * Revision 1.112  2002/08/12 13:22:46  mju
 * Remove retired stuff.
 *
 * Revision 1.111  2002/07/15 14:57:20  mju
 * Add datum/projsys for NZGD2000/NZTM.
 *
 * Revision 1.110  2002/05/17 15:38:15  dwilliss
 * Don't leave dangling commas at the end of enums.
 *
 * Revision 1.109  2002/03/08 15:48:57  mju
 * Use differet memcpy defn for where 'amemcpy' used.
 *
 * Revision 1.108  2002/03/08 15:16:46  mju
 * Disallow memset/memcpy on MAPPROJPARM.
 *
 * Revision 1.107  2002/02/19 17:44:18  mju
 * Add PROJECTION enum.
 *
 * Revision 1.106  2002/02/18 15:31:25  mju
 * Use normal numbering for SJTSK datum.
 *
 * Revision 1.105  2002/02/18 15:28:39  mju
 * Add DATUM_SJTSK95.
 *
 * Revision 1.104  2002/02/06 18:47:49  mju
 * Add PROJ_RoyalDutch.
 *
 * Revision 1.103  2002/02/06 15:34:20  mju
 * Add PROJSYS_Dutch.
 *
 * Revision 1.102  2002/01/24 22:50:59  mju
 * Fix typo.
 *
 * Revision 1.101  2002/01/24 22:50:12  mju
 * Add PROJSYS_MichiganGeoRef.
 *
 * Revision 1.100  2001/12/20 21:43:22  mju
 * Enum docs.
 *
 * Revision 1.99  2001/12/19 14:16:52  mju
 * Attempt to get Genitor to provide proper link to nested class definitions.
 *
 * Revision 1.97  2001/11/16 14:35:39  mju
 * Add PROJSYS_Swiss.
 *
 * Revision 1.96  2001/11/15 23:28:37  mju
 * Add DATUM_Tokyo_Japan.
 *
 * Revision 1.95  2001/11/15 23:08:47  mju
 * Add PROJSYS_Japan19.
 *
 * Revision 1.94  2001/11/01 18:48:59  mju
 * Remove 'NUMPROJS' defn as obsolete.
 *
 * Revision 1.93  2001/10/30 19:14:23  mju
 * Add DATUM enums used by MapInfo.
 *
 * Revision 1.92  2001/10/19 17:56:53  mju
 * Add PROJSYS_GaussKruger3deg & Austria.
 * Add DATUM_MGI.
 *
 * Revision 1.91  2001/08/27 21:02:53  mju
 * Add PROJSYS_GhanaNational.
 *
 * Revision 1.90  2001/08/14 21:18:30  mju
 * Add SetSystemGaussKruger.
 *
 * Revision 1.89  2001/07/10 17:52:32  scowan
 * Changed function names to fit naming convension.
 *
 * Revision 1.88  2001/07/10 17:40:00  scowan
 * Added more exported functions.
 *
 * Revision 1.85  2001/06/21 17:37:26  mju
 * Add PROJSYS_AMG/MGA.
 *
 * Revision 1.84  2001/06/21 15:26:58  mju
 * Add more DATUM enums.
 *
 * Revision 1.83  2001/06/20 13:18:00  mju
 * Add ProjGetERMapper().
 *
 * Revision 1.82  2001/06/19 20:36:34  mju
 * Add new enums for datums used in ER-Mapper.
 *
 * Revision 1.80  2000/11/27 15:37:10  mju
 * Rename 'proj name' as 'oldprojname'.
 *
 * Revision 1.79  2000/11/22 17:02:05  mju
 * "Remove" setprojparmname().
 *
 * Revision 1.74  2000/09/22 15:04:29  mju
 * Add DATUM_Schwarzeck.
 *
 * Revision 1.73  2000/09/20 14:50:07  mju
 * Add op==/op!= for MAPPROJPARM.
 *
 * Revision 1.72  2000/07/18 21:36:58  mju
 * Add mapprojparm::SetSystemUPS.
 *
 * Revision 1.71  2000/07/06 17:06:17  mju
 * Define enum operators for GETPROJPARMFLAGS.
 *
 * Revision 1.70  2000/06/28 20:22:30  mju
 * Add new datum enums used in er-mapper import.
 *
 * Revision 1.63  2000/05/25 13:07:28  mju
 * Include stdlib.h if not already so hopefully will have abs() defined on the Mac.
 *
 * Revision 1.62  2000/05/24 22:56:17  dwilliss
 * use ::abs
 *
 * Revision 1.61  2000/05/12 15:18:54  mju
 * Include string.h for memset() defn.
 *
 * Revision 1.60  2000/05/11 22:12:19  mju
 * Add MAPPROJPARM ctor.
 *
 * Revision 1.55  2000/05/10 16:06:15  mju
 * Make DATUM and PROJSYS enums.
 * Add lots of inline methods to MAPPROJPARM, including methods for setting projection-specific parameters.
 *
 * Revision 1.53  2000/03/01 20:32:34  dwilliss
 * Named some parms in prototypes to make include file more self documenting.
 * (Useful for those of use who use VisualStudio to edit instead of CodeWrite)
 *
 * Revision 1.52  2000/02/29 17:01:54  mju
 * Add ELLIPSOID constants and TPARM constants.
 *
 * Revision 1.50  1999/07/21  16:48:49  mju
 * Add PROJSYS_Finnish & DATUM_KKJ.
 *
 * Revision 1.49  1999/07/21  16:32:38  mju
 * Add additional PROJSYS defns.
 *
 * Revision 1.48  1999/07/16  21:57:58  mju
 * Add DATUM_TRANS_ defns for Similarity and BursaWolfe.
 *
 * Revision 1.47  1999/07/15  21:56:57  mju
 * Add NZMG.
 *
 * Revision 1.46  1999/07/14  17:34:46  mju
 * Change PointZone to use UNICODE.
 *
 * Revision 1.45  1999/07/09  21:31:39  mju
 * Add auto-zone features.
 * Add MGRS system.
 * \endif
**/

#ifndef  INC_MI32_MAPPROJ_H
#define  INC_MI32_MAPPROJ_H
#define  INC_MI32_MAPPROJ_H

#ifndef  INC_MI32_STDDEFNS_H
#include <mi32/stddefns.h>
#endif

#ifndef  INC_MI32_ERRCODES_H
#include <mi32/errcodes.h>
#endif

#ifndef  INC_MI32_SIMPLEAR_H
#include <mi32/simplear.h>
#endif

#ifndef  INC_COORD_H
#include <mi32/point.h>
#endif

#ifndef  INC_MI32_MISTRING_H
   #include <mi32/mistring.h>
#endif

#ifndef  INC_STRING_H
#include <string.h>     // Includes defn for memset()
#define  INC_STRING_H
#endif

#ifndef  INC_STDLIB_H
#include <stdlib.h>
#define  INC_STDLIB_H
#endif

#ifdef GEOMDLL
   #define GEOMLIBEXPORT MI_DLLEXPORT
#else
   #define GEOMLIBEXPORT MI_DLLIMPORT
#endif

//! Forward declarations

struct CTRLPOINT3;
struct MAPPROJPARM;
class FILEPATH;

enum PROJECTION {
   PROJECTION_Undefined =                       -3,   //!< Unknown or undefined projection
   PROJECTION_None =                            -2,   //!< No projection
   PROJECTION_Arbitrary =                       -1,   //!< Arbitrary local coordinates
   PROJECTION_LatLon =                           0,   //!< Latitude/longitude
   PROJECTION_Geographic =                       0,   //!< Latitude/longitude
   PROJECTION_UTM =                              1,   //!< Universal Transverse Mercator
   PROJECTION_StatePlane =                       2,   //!< State plane
   PROJECTION_AlbersEqualAreaConic =             3,   //!< Albers Equal-Area Conic
   PROJECTION_AEAC =                             3,   //!< Albers Equal-Area Conic
   PROJECTION_LambertConformalConic =            4,   //!< Lambert Conformal Conic
   PROJECTION_LCC =                              4,   //!< Lambert Conformal Conic
   PROJECTION_Mercator =                         5,   //!< Mercator
   PROJECTION_PolarStereographic =               6,   //!< Polar Stereographic
   PROJECTION_Polyconic =                        7,   //!< Polyconic
   PROJECTION_EquidistantConic =                 8,   //!< Equidistant Conic
   PROJECTION_TransverseMercator =               9,   //!< Transverse Mercator
   PROJECTION_Stereographic =                   10,   //!< Stereographic
   PROJECTION_LambertAzimuthalEqualArea =       11,   //!< Lambert Azimuthal Equal-Area
   PROJECTION_LAEA =                            11,   //!< Lambert Azimuthal Equal-Area
   PROJECTION_AzimuthalEquidistant =            12,   //!< Azimuthal Equidistant
   PROJECTION_Gnomonic =                        13,   //!< Gnomonic
   PROJECTION_Orthographic =                    14,   //!< Orthographic
   PROJECTION_GeneralVertNearSidePerspective =  15,   //!< General Vertical Near-Side Perspective
   PROJECTION_GVNSP =                           15,   //!< General Vertical Near-Side Perspective
   PROJECTION_Sinusoidal =                      16,   //!< Sinusoidal
   PROJECTION_Equirectangular =                 17,   //!< Equirectangular
   PROJECTION_MillerCylindrical =               18,   //!< Miller Cylindrical
   PROJECTION_VanDerGrintenI =                  19,   //!< Van Der Grinten I
   PROJECTION_ObliqueMercator =                 20,   //!< Oblique Mercator (Hotine)
   PROJECTION_Robinson =                        21,   //!< Robinson
   PROJECTION_BipolarObliqueConformalConic =    22,   //!< Bipolar Oblique Conformal Conic
   PROJECTION_BOCC =                            22,   //!< Bipolar Oblique Conformal Conic
   PROJECTION_LandsatSpaceObliqueMercator =     23,   //!< Space Oblique Mercator for Landsat
   PROJECTION_LandsatSOM =                      23,   //!< Space Oblique Mercator for Landsat
   PROJECTION_Cassini =                         24,   //!< Cassini
   PROJECTION_NewZealandMapGrid =               25,   //!< New Zealand Map Grid
   PROJECTION_NZMG =                            25,   //!< New Zealand Map Grid
   PROJECTION_SwissConformalCylindrical =       26,   //!< Swiss Conformal Cylindrical
   PROJECTION_RoyalDutch =                      27,   //!< Dutch National Grid (RD)
   PROJECTION_ObliqueStereographic =            28,   //!< Oblique Stereographic
   PROJECTION_InterruptedGoodeHomolosine =      29,   //!< Interrupted Goode Homolosine
   PROJECTION_Goode =                           29,   //!< Interrupted Goode Homolosine
   };


#define  MAXPROJ  29

#define  PROJ_Undefined    PROJECTION_Undefined
#define  PROJ_None         PROJECTION_None
#define  PROJ_Arbitrary    PROJECTION_Arbitrary
#define  PROJ_LatLon       PROJECTION_LatLon
#define  PROJ_UTM          PROJECTION_UTM
#define  PROJ_StatePlane   PROJECTION_StatePlane
#define  PROJ_AlbersConic  PROJECTION_AlbersEqualAreaConic
#define  PROJ_LambertConf  PROJECTION_LambertConformalConic
#define  PROJ_Mercator     PROJECTION_Mercator
#define  PROJ_PolarStereo  PROJECTION_PolarStereographic
#define  PROJ_Polyconic    PROJECTION_Polyconic
#define  PROJ_EquidistCon  PROJECTION_EquidistantConic
#define  PROJ_TransMercat  PROJECTION_TransverseMercator
#define  PROJ_StereoGraph  PROJECTION_Stereographic
#define  PROJ_LambertAzim  PROJECTION_LambertAzimuthalEqualArea
#define  PROJ_AzimEquidis  PROJECTION_AzimuthalEquidistant
#define  PROJ_Gnomonic     PROJECTION_Gnomonic
#define  PROJ_Orthograph   PROJECTION_Orthographic
#define  PROJ_GenVertPers  PROJECTION_GVNSP
#define  PROJ_Sinusoidal   PROJECTION_Sinusoidal
#define  PROJ_Equirectang  PROJECTION_Equirectangular
#define  PROJ_MillerCylin  PROJECTION_MillerCylindrical
#define  PROJ_VanDerGrint  PROJECTION_VanDerGrintenI
#define  PROJ_ObliqueMerc  PROJECTION_ObliqueMercator
#define  PROJ_Robinson     PROJECTION_Robinson
#define  PROJ_BipolarConic PROJECTION_BOCC
#define  PROJ_LandsatSOM   PROJECTION_LandsatSOM
#define  PROJ_Cassini      PROJECTION_Cassini
#define  PROJ_NZMG         PROJECTION_NZMG
#define  PROJ_SwissConCyl  PROJECTION_SwissConformalCylindrical
#define  PROJ_RoyalDutch   PROJECTION_RoyalDutch

#define  TPARM_SemiMajorAxis  0
#define  TPARM_SemiMinorAxis  1
#define  TPARM_FalseEasting   6
#define  TPARM_FalseNorthing  7

#define  DATUM_TRANS_BuiltIn     0     //!< Built-in transformation
#define  DATUM_TRANS_Molodensky  1     //!< Molodensky 3-parameter transformation
#define  DATUM_TRANS_Similarity  2     //!< 7-parameter similarity transformation
#define  DATUM_TRANS_BursaWolfe  3     //!< 7-parameter Bursa-Wolfe transformation
#define  DATUM_TRANS_Grid        4     //!< Grid offsets

//! Geodetic datum IDs.
enum DATUM {

   DATUM_Unspecified =           0,          //!< Unspecified or unknown
   DATUM_WGS84 =                 1,          //!< World Geodetic System 1984
   DATUM_WGS72 =                 2,          //!< World Geodetic System 1972
   DATUM_NAD27 =                 1927,       //!< North American Datum 1927
   DATUM_NAD83 =                 1983,       //!< North American Datum 1983

   DATUM_Adindan =               101,
   DATUM_Adindan_Sudan =         107,
   DATUM_Afgooye =               108,
   DATUM_AGD66 =                 128,
   DATUM_AGD84 =                 289,
   DATUM_AinElAbd_Bahrain =      109,
   DATUM_Amersfoort =            290,
   DATUM_Anna1Astro =            111,
   DATUM_Arc1950 =               113,
   DATUM_Arc1960 =               122,
   DATUM_AscensionIsland1968 =   123,
   DATUM_AstroBeaconE =          124,
   DATUM_AstroDOS_71_4 =         125,
   DATUM_AstroTernIsland1961 =   126,
   DATUM_AstroStation1952 =      127,
   DATUM_Bahrain =               109,
   DATUM_Beijing =               324,
   DATUM_BellevueIGN =           131,
   DATUM_Bermuda1957 =           132,
   DATUM_Bogota =                134,
   DATUM_Camacupa =              291,
   DATUM_CampoInchauspe =        137,
   DATUM_CantonAstro1966 =       138,
   DATUM_Cape =                  139,
   DATUM_CapeBursaWolfe =        329,
   DATUM_CapeBursaMod =          334,
   DATUM_CapeCanaveral =         140,
   DATUM_Carthage =              141,
   DATUM_CH1903 =                331,
   DATUM_ChathamIsland1971 =     142,
   DATUM_Chua =                  143,
   DATUM_CorregoAlegre =         144,
   DATUM_DealulPiscului1970 =    417,
   DATUM_Djakarta =              146,
   DATUM_DOS1968 =               147,
   DATUM_EasterIsland1967 =      148,
   DATUM_ED50 =                  149,
   DATUM_ED50_Egypt =            153,
   DATUM_ED79 =                  162,
   DATUM_EgyptOld =              223,
   DATUM_Fahud =                 229,
   DATUM_GandajikaBase1970 =     164,
   DATUM_GDA94 =                 288,
   DATUM_GraciosaBaseSW1948 =    166,
   DATUM_Guam1963 =              167,
   DATUM_GUX1Astro =             169,
   DATUM_HitoXVIII1963 =         249,
   DATUM_Hjorsey1955 =           171,
   DATUM_HongKong1963 =          172,
   DATUM_HuTzuShan =             173,
   DATUM_Ind54 =                 154,
   DATUM_Ind60 =                 335,
   DATUM_Ind75 =                 177,
   DATUM_Indian_Bangladesh =     174,
   DATUM_Indian_India_Nepal =    175,
   DATUM_Ireland1965 =           178,
   DATUM_ISTS073Astro1969 =      180,
   DATUM_JohnstonIsland1961 =    181,
   DATUM_Kalianpur_India =       306,
   DATUM_Kandawala =             182,
   DATUM_KerguelenIsland1949 =   183,
   DATUM_Kertau =                184,
   DATUM_KKJ =                   326,
   DATUM_LC5Astro1961 =          186,
   DATUM_Leigon =                187,
   DATUM_Liberia1954 =           188,
   DATUM_Luzon =                 189,
   DATUM_Luzon_Mindanao =        190,
   DATUM_Mahe =                  191,
   DATUM_Manoca =                309,
   DATUM_Massawa =               192,
   DATUM_Merchich =              193,
   DATUM_MidwayAstro1961 =       194,
   DATUM_MGI =                   339,
   DATUM_Minna_Cameroon =        195,
   DATUM_Minna_Nigeria =         196,
   DATUM_Mporo =                 198,
   DATUM_NAD27_Alaska =          210,
   DATUM_NAD27_Bahamas =         211,
   DATUM_NAD27_Canada =          206,
   DATUM_NAD27_CanalZone =       218,
   DATUM_NAD27_Caribbean =       204,
   DATUM_NAD27_CentralAmerica =  205,
   DATUM_NAD27_Cuba =            219,
   DATUM_NAD27_Greenland =       220,
   DATUM_NAD27_Mexico =          221,
   DATUM_NAD27_SanSalvador =     212,
   DATUM_Nahrwan =               199,
   DATUM_Nahrwan_SaudiArabia =   201,
   DATUM_Nahrwan_UAE =           202,
   DATUM_NaparimaBWI =           203,
   DATUM_NewZealand1949 =        347,
   DATUM_NZGD2000 =              348,
   DATUM_NTF =                   312,
   DATUM_NWL9D =                 337,
   DATUM_ObservatorioMet1939 =   222,
   DATUM_OSGB36 =                294,
   DATUM_OldHawaii =             224,        //!< Old Hawaiian - Mean for Hawaii, Kauuai, Maui, Oahu
   DATUM_PicoDeLasNieves =       235,
   DATUM_PitcairnAstro1967 =     236,
   DATUM_PointeNoire =           238,
   DATUM_PortoSanto1936 =        239,
   DATUM_Potsdam =               325,
   DATUM_PRS92 =                 336,
   DATUM_PSAD56 =                240,
   DATUM_PuertoRico =            250,        //!< Puerto Rico datum
   DATUM_Pulkovo_Romania =       319,
   DATUM_Pulkovo_Russia =        287,
   DATUM_Qatar =                 251,
   DATUM_Qornoq =                252,
   DATUM_Reunion =               253,
   DATUM_Rome1940 =              389,
   DATUM_Rome1940_Sardinia =     387,
   DATUM_Rome1940_Sicily =       388,
   DATUM_SAD69 =                 261,
   DATUM_SantoDOS1965 =          255,
   DATUM_SaoBraz =               256,
   DATUM_SapperHill1943 =        257,
   DATUM_SelvagemGrande =        159,
   DATUM_Schwarzeck =            258,
   DATUM_SJTSK95 =               343,
   DATUM_SouthAsia =             274,
   DATUM_Tananarive =            275,
   DATUM_Timbalai =              276,
   DATUM_Tokyo_Mean =            277,
   DATUM_Tokyo_Japan =           278,
   DATUM_Tokyo_Grid =            474,
   DATUM_TristanAstro1968 =      281,
   DATUM_VitiLevu1916 =          282,
   DATUM_WakeEniwetok1960 =      283,
   DATUM_Yacare =                285,
   DATUM_Zanderij =              286,

   // The following datums have multiple regression transformations
   DATUM_MRE_AGD66 =    1001,          //!< Australian Geodetic Datum 1966
   DATUM_MRE_AGD84 =    1002,          //!< Australian Geodetic Datum 1984
   DATUM_MRE_CAI =      1003,          //!< Campo Inchauspe Datum (Argentina)
   DATUM_MRE_COA =      1004,          //!< Corrego Alegro Datum (Brazil)
   DATUM_MRE_ED50 =     1005,          //!< European Datum 1950 (Western Europe)
   DATUM_MRE_NAD27_C =  1006,          //!< North American Datum 1927 (Canada)
   DATUM_MRE_NAD27_U =  1007,          //!< North American Datum 1927 (Contental U.S.)
   DATUM_MRE_SAD69 =    1008           //!< South American Datum 1969

   };

#define  DATUM_MRE_First   1001
#define  DATUM_MRE_Last    1008

//! Commonly-used ellipsoids
#define  ELLIPSOID_GRS_1980         1
#define  ELLIPSOID_Clarke_1866      2
#define  ELLIPSOID_WGS_1972         7
#define  ELLIPSOID_WGS_1984         51

//! Projection system IDs.
enum PROJSYS {
   PROJSYS_UserDefined =      0,
   PROJSYS_LatLon =           1,       //!< Latitude-Longitude
   PROJSYS_Geographic =       2,       //!< Same as latitude-longitude
   PROJSYS_StatePlane1927 =   3,       //!< United States State Plane 1927
   PROJSYS_StatePlane1983 =   4,       //!< United States State Plane 1983
   PROJSYS_UTM =              5,       //!< Universal Transverse Mercator
   PROJSYS_UPS =              6,       //!< Universal Polar Stereographic
   PROJSYS_GaussKruger =      7,       //!< Gauss-Kruger (6-degree zones)
   PROJSYS_MGRS =             8,       //!< Military Grid Reference System
   PROJSYS_NZMG =             9,       //!< New Zealand Map Grid
   PROJSYS_BritishNational =  10,      //!< British National Grid
   PROJSYS_German =           11,      //!< Germany
   PROJSYS_Finnish =          12,      //!< Finnish National Grid
   PROJSYS_Namibian =         13,      //!< Namibian Map System
   PROJSYS_AMG =              14,      //!< Australian Map Grid
   PROJSYS_MGA =              15,      //!< Map Grid of Australia
   PROJSYS_GhanaNational =    16,      //!< Ghana National Grid
   PROJSYS_GaussKruger3deg =  17,      //!< Gauss-Kruger (3-degree zones)
   PROJSYS_Austria =          18,      //!< Austria
   PROJSYS_Japan19 =          19,      //!< Japan-19
   PROJSYS_Swiss =            20,      //!< Swiss
   PROJSYS_MichiganGeoRef =   21,      //!< Michigan GeoRef
   PROJSYS_Dutch =            22,      //!< Dutch National Grid
   PROJSYS_NZTM =             23,      //!< New Zealand Transverse Mercator
   PROJSYS_Romania =          24,      //!< Romania
   PROJSYS_Italy =            25,      //!< Italy
   PROJSYS_Colombia =         26,      //!< Colombia
   PROJSYS_WisconsinCounty =  101,     //!< Wisconsin County Coordinate System

   PROJSYS_COUNT
   };

//----------------------------------------------------------------------------
//    Function Prototypes                                                     
//----------------------------------------------------------------------------

//!:Associate with "Map Projection Functions"
//!\addtogroup mapproj Map Projection Functions
//!@{

GEOMLIBEXPORT ERRVALUE MapInfoReadTabGeoreference (
   const FILEPATH& filepath, 
   MAPPROJPARM& projparm, 
   SIMPLE_ARRAY<CTRLPOINT3>& points
   );
   
GEOMLIBEXPORT ERRVALUE MapInfoExportGeoreference (
   const MAPPROJPARM& projparm, 
   int xyunits, 
   int& MapInfoProj, 
   int& MapInfoDatum, 
   double MapInfoProjParm[6], 
   char ProjStr[256]
   );

extern "C" {

#if !defined(NO_DEPRECATED) || defined(GEOMDLL)
//! Compute distance in meters between two points in specified projection.
//! \deprecated Use SPATREF classes.
//! @return Distance between points in meters.
//! If the coordinate system is Latitude-Longitude the distance along a geodetic
//! is returned (actual distance over curved surface).
GEOMLIBEXPORT double ProjDistToMeters (
   const MAPPROJPARM* proj,            //!< Projection parameters, will use GRS-1980 ellipsoid if NULL
   double x1,                          //!< X coordinate of first point (in meters or degrees)
   double y1,                          //!< Y coordinate of first point
   double x2,                          //!< X coordinate of second point
   double y2                           //!< Y coordinate of second point
   );
#endif

//! Set ER-Mapper projection and datum names from MAPPROJPARM.
GEOMLIBEXPORT ERRVALUE ProjGetERMapper (
   const MAPPROJPARM& projparm,        //!< Projection parameters
   char *ProjectionName,               //!< ER-Mapper projection name returned, must be at least 16 characters
   char *DatumName                     //!< ER-Mapper datum name returned, must be at least 16 characters
   );

//! Set projection parameters from ER-Mapper Projection and Datum names.
GEOMLIBEXPORT ERRVALUE ProjSetERMapper (
   MAPPROJPARM& projparm,              //!< Projection parameters to fill in
   const char *ProjectionName,         //!< ER-Mapper projection name
   const char *DatumName               //!< ER-Mapper datum name
   );

}  // Extern "C"

GEOMLIBEXPORT void ProjGetArcPrj (
   const MAPPROJPARM& ProjParm,
   MISTRINGLIST& PrjList
   );

GEOMLIBEXPORT bool ProjSetArcPrj (
   MAPPROJPARM& ProjParm,
   const MISTRINGLIST& PrjList,
   double& xshift,
   double& yshift,
   double& scale
   );

//!@}

//----------------------------------------------------------------------------
//!      Map projection parameter structure                                      
//----------------------------------------------------------------------------

#define  PROJZ_ZONE     0
#define  PROJZ_SYSTEM   1
#define  PROJZ_DATUM    2
#define  PROJZ_ALTPROJ  3

//! Parameters for a specific map projection, coordinate system and geodetic datum.
struct MAPPROJPARM {
   // This structure is used extensively within other structure and in DLLs.
   // Therefore, members may not be added or removed, nor may virtual methods be added.

   #ifndef GENERATING_DOXYGEN_OUTPUT

   INT16 mapproj;
   UNICODE oldprojname[35];      //!< No longer valid, use GetName() instead
   INT16 zone[4];
   double tparm[15];

   #endif // GENERATING_DOXYGEN_OUTPUT

   //! Default constructor.
   MAPPROJPARM (
      ) { Clear(); }

   //! Clear all parameters and reset to Arbitrary Local coordinates.
   void Clear (
      ) { memset(static_cast<void*>(this),0,sizeof(*this)); mapproj = PROJECTION_Arbitrary; }

   //! Get datum ID.
   DATUM GetDatum (
      ) const { return (static_cast<DATUM>(zone[PROJZ_DATUM])); }

#ifndef NO_DEPRECATED
   //! Compute distance in meters between specified points.
   //! \deprecated Use SPATREF classes.
   //! @return Distance in meters.
   //! If the coordinate system is Latitude-Longitude the distance along a geodetic
   //! is returned (actual distance over curved surface).
   DEPRECATED double GetDistanceInMeters (
      double x1,
      double y1,
      double x2,
      double y2
      ) const { return (ProjDistToMeters(this,x1,y1,x2,y2)); }
#endif

#ifndef NO_DEPRECATED
   //! Determine distance in meters between specified points.
   //! \deprecated Use SPATREF classes.
   //! @return Distance in meters.
   //! If the coordinate system is Latitude-Longitude the distance along a geodetic
   //! is returned (actual distance over curved surface).
   DEPRECATED double GetDistanceInMeters (
      const DPOINT2D& point1,
      const DPOINT2D& point2
      ) const { return (ProjDistToMeters(this,point1.x,point1.y,point2.x,point2.y)); }
#endif

   //! Get projection ID.
   PROJECTION GetProjection (
      ) const { return (static_cast<PROJECTION>(mapproj)); }

   //! Get sub-projection ID.
   INT16 GetProjSub (
      ) const { return (zone[PROJZ_ALTPROJ]); }

   //! Get system ID.
   PROJSYS GetSystem (
      ) const { return (static_cast<PROJSYS>(zone[PROJZ_SYSTEM])); }

   //! Get zone ID.
   INT16 GetZone (
      ) const { return (zone[PROJZ_ZONE]); }

   //! Set datum ID.
   void SetDatum (
      DATUM datum
      ) { zone[PROJZ_DATUM] = static_cast<INT16>(datum); }

   //! Set projection to Albers Equal-Area Conic.
   ERRVALUE SetProjectionAEAC (
      double LatStdParallel1,          //!< Latitude of first standard parallel (degrees)
      double LatStdParallel2,          //!< Latitude of second standard parallel (degrees)
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double LatOrigin = 0.0,          //!< Latitude of origin (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_AEAC;
      tparm[2] = LatStdParallel1;
      tparm[3] = LatStdParallel2;
      tparm[4] = LonCentralMeridian;
      tparm[5] = LatOrigin;
      tparm[TPARM_FalseEasting] = FalseEasting;
      tparm[TPARM_FalseNorthing] = FalseNorthing;
      return (0);
      }

   //! Set projection to arbitrary.
   void SetProjectionArbitrary (
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      mapproj = PROJECTION_Arbitrary;
      return;
      }

   //! Set projection to Azimuthal Equidistant.
   ERRVALUE SetProjectionAzimuthalEquidistant (
      double LonCenter,                //!< Longitude of center of projection (degrees)
      double LatCenter,                //!< Latitude of center of projection (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj4567(PROJECTION_AzimuthalEquidistant,LonCenter,LatCenter,FalseEasting,FalseNorthing));
      }

   //! Set projection to Bipolar Oblique Conformal Conic.
   ERRVALUE SetProjectionBOCC (
      bool skewed = true               //!< true if skewed, false if not
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_BOCC;
      tparm[2] = !skewed;
      return (0);
      }

   //! Set projection to Cassini.
   ERRVALUE SetProjectionCassini (
      double LonCenter,                //!< Longitude of center of projection (degrees)
      double LatCenter,                //!< Latitude of center of projection (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj4567(PROJECTION_Cassini,LonCenter,LatCenter,FalseEasting,FalseNorthing));
      }

   //! Set projection to Equidistant Conic with single standard parallel.
   ERRVALUE SetProjectionEquidistantConic1SP (
      double LatStdParallel,           //!< Latitude of standard parallel (degrees)
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double LatOrigin = 0.0,          //!< Latitude of origin (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_EquidistantConic;
      tparm[2] = tparm[3] = LatStdParallel;
      tparm[4] = LonCentralMeridian;
      tparm[5] = LatOrigin;
      tparm[TPARM_FalseEasting] = FalseEasting;
      tparm[TPARM_FalseNorthing] = FalseNorthing;
      tparm[8] = 0.0;
      return (0);
      }

   //! Set projection to Equidistant Conic with two standard parallels.
   ERRVALUE SetProjectionEquidistantConic2SP (
      double LatStdParallel1,          //!< Latitude of first standard parallel (degrees)
      double LatStdParallel2,          //!< Latitude of second standard parallel (degrees)
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double LatOrigin = 0.0,          //!< Latitude of origin (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_EquidistantConic;
      tparm[2] = LatStdParallel1;
      tparm[3] = LatStdParallel2;
      tparm[4] = LonCentralMeridian;
      tparm[5] = LatOrigin;
      tparm[TPARM_FalseEasting] = FalseEasting;
      tparm[TPARM_FalseNorthing] = FalseNorthing;
      tparm[8] = 1.0;
      return (0);
      }

   //! Set projection to Equirectangular.
   ERRVALUE SetProjectionEquirectangular (
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double LatTrueScale,             //!< Latitude of true scale (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj4567(PROJECTION_Equirectangular,LonCentralMeridian,LatTrueScale,FalseEasting,FalseNorthing));
      }

   //! Set projection to Gnomonic.
   ERRVALUE SetProjectionGnomonic (
      double LonCenter,                //!< Longitude of center of projection (degrees)
      double LatCenter,                //!< Latitude of center of projection (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj4567(PROJECTION_Gnomonic,LonCenter,LatCenter,FalseEasting,FalseNorthing));
      }

   //! Set projection to General Vertical Near-Side Perspective.
   ERRVALUE SetProjectionGVNSP (
      double Height,                   //!< Height above sphere of reference (meters)
      double LonCenter,                //!< Longitude of center of projection (degrees)
      double LatCenter,                //!< Latitude of center of projection (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj24567(PROJECTION_GVNSP,Height,LonCenter,LatCenter,FalseEasting,FalseNorthing));
      }

   //! Set projection to Lambert Azimuthal Equal-Area.
   ERRVALUE SetProjectionLAEA (
      double LonCenter,                //!< Longitude of center of projection (degrees)
      double LatCenter,                //!< Latitude of center of projection (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj4567(PROJECTION_LAEA,LonCenter,LatCenter,FalseEasting,FalseNorthing));
      }

   //! Set projection to Lambert Conformal Conic.
   ERRVALUE SetProjectionLCC (
      double LatStdParallel1,          //!< Latitude of first standard parallel (degrees)
      double LatStdParallel2,          //!< Latitude of second standard parallel (degrees)
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double LatOrigin = 0.0,          //!< Latitude of origin (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0,      //!< False northing (meters)
      double ScaleFactor = 1.0         //!< Scale factor
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_LCC;
      tparm[2] = LatStdParallel1;
      tparm[3] = LatStdParallel2;
      tparm[4] = LonCentralMeridian;
      tparm[5] = LatOrigin;
      tparm[TPARM_FalseEasting] = FalseEasting;
      tparm[TPARM_FalseNorthing] = FalseNorthing;
      tparm[8] = ScaleFactor;
      return (0);
      }

   //! Set projection to Landsat Space Oblique Mercator.
   ERRVALUE SetProjectionLandsatSOM (
      int satellite,                   //!< Satellite (1-5)
      int path                         //!< Satellite path (1-251)
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_LandsatSOM;
      tparm[2] = satellite;
      tparm[3] = path;
      return (0);
      }

   //! Set projection to Mercator.
   ERRVALUE SetProjectionMercator (
      double ScaleFactor,              //!< Scale factor at central meridian
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double LatTrueScale,             //!< Latitude of true scale (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj24567(PROJECTION_Mercator,ScaleFactor,LonCentralMeridian,LatTrueScale,FalseEasting,FalseNorthing));
      }

   //! Set projection to MillerCylindrical.
   ERRVALUE SetProjectionMillerCylindrical (
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj467(PROJECTION_MillerCylindrical,LonCentralMeridian,FalseEasting,FalseNorthing));
      }

   //! Set projection to Oblique Mercator (format A).
   ERRVALUE SetProjectionObliqueMercatorA (
      double ScaleFactor,              //!< Scale factor at center
      double LatOrigin,                //!< Latitude of origin (degrees)
      double LonGeodeticPoint1,        //!< Longitude of first point defining central geodetic line of projection
      double LatGeodeticPoint1,        //!< Latitude of first point defining central geodetic line of projection
      double LonGeodeticPoint2,        //!< Longitude of second point defining central geodetic line of projection
      double LatGeodeticPoint2,        //!< Latitude of second point defining central geodetic line of projection
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_ObliqueMercator;
      tparm[2] = ScaleFactor;
      tparm[5] = LatOrigin;
      tparm[TPARM_FalseEasting] = FalseEasting;
      tparm[TPARM_FalseNorthing] = FalseNorthing;
      tparm[8] = LonGeodeticPoint1;
      tparm[9] = LatGeodeticPoint1;
      tparm[10] = LonGeodeticPoint2;
      tparm[11] = LatGeodeticPoint2;
      return (0);
      }

   //! Set projection to Oblique Mercator (format B).
   ERRVALUE SetProjectionObliqueMercatorB (
      double ScaleFactor,              //!< Scale factor at center
      double AzimuthAngle,             //!< Angle of azimuth east of north for central line of projection
      double LonAzimuthPoint,          //!< Longitude of point along central line of projection where azimuth is measured (degrees)
      double LatOrigin,                //!< Latitude of origin (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_ObliqueMercator;
      tparm[2] = ScaleFactor;
      tparm[3] = AzimuthAngle;
      tparm[4] = LonAzimuthPoint;
      tparm[5] = LatOrigin;
      tparm[TPARM_FalseEasting] = FalseEasting;
      tparm[TPARM_FalseNorthing] = FalseNorthing;
      return (0);
      }

   //! Set projection to Orthographic.
   ERRVALUE SetProjectionOrthographic (
      double LonCenter,                //!< Longitude of center of projection (degrees)
      double LatCenter,                //!< Latitude of center of projection (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj4567(PROJECTION_Orthographic,LonCenter,LatCenter,FalseEasting,FalseNorthing));
      }

   //! Set projection to Polar Stereographic.
   ERRVALUE SetProjectionPolarStereo (
      double ScaleFactor,              //!< Scale factor at center
      double LonDownward,              //!< Longitude directed straight down from pole
      double LatTrueScale,             //!< Latitude of true scale (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj24567(PROJECTION_PolarStereographic,ScaleFactor,LonDownward,LatTrueScale,FalseEasting,FalseNorthing));
      }

   //! Set projection to Polyconic.
   ERRVALUE SetProjectionPolyconic (
      double ScaleFactor,              //!< Scale factor at center of projection
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double LatOrigin = 0.0,          //!< Latitude of origin (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj24567(PROJECTION_Polyconic,ScaleFactor,LonCentralMeridian,LatOrigin,FalseEasting,FalseNorthing));
      }

   //! Set projection to Robinson.
   ERRVALUE SetProjectionRobinson (
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj467(PROJECTION_Robinson,LonCentralMeridian,FalseEasting,FalseNorthing));
      }

   //! Set projection to Sinusoidal.
   ERRVALUE SetProjectionSinusoidal (
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj467(PROJECTION_Sinusoidal,LonCentralMeridian,FalseEasting,FalseNorthing));
      }

   //! Set projection to Stereographic.
   ERRVALUE SetProjectionStereographic (
      double ScaleFactor,              //!< Scale factor at center of projection
      double LonCenter,                //!< Longitude of center of projection (degrees)
      double LatCenter,                //!< Latitude of center of projection (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj24567(PROJECTION_Stereographic,ScaleFactor,LonCenter,LatCenter,FalseEasting,FalseNorthing));
      }

   //! Set projection to Swiss Conformal.
   void SetProjectionSwissConformal (
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = PROJECTION_SwissConformalCylindrical;
      tparm[6] = FalseEasting;
      tparm[7] = FalseNorthing;
      }
      

   //! Set projection to Transverse Mercator.
   ERRVALUE SetProjectionTransverseMercator (
      double ScaleFactor,              //!< Scale factor at center of projection
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double LatOrigin = 0.0,          //!< Latitude of origin (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj24567(PROJECTION_TransverseMercator,ScaleFactor,LonCentralMeridian,LatOrigin,FalseEasting,FalseNorthing));
      }

   //! Set projection to Van der Grinten I.
   ERRVALUE SetProjectionVanDerGrintenI (
      double LonCentralMeridian,       //!< Longitude of central meridian (degrees)
      double FalseEasting = 0.0,       //!< False easting (meters)
      double FalseNorthing = 0.0       //!< False northing (meters)
      ) {
      return (SetProj467(PROJECTION_VanDerGrintenI,LonCentralMeridian,FalseEasting,FalseNorthing));
      }

   //! Set coordinate system to Gauss-Kruger.
   ERRVALUE SetSystemGaussKruger (
      INT16 gkzone,                    //!< Zone (1-60, or -1 - -60 as alternative for South zones)
      bool southzone = false           //!< South zones
      ) {
      if (gkzone < -60 || gkzone == 0 || gkzone > 60) return (EBadFuncParm);
      zone[PROJZ_SYSTEM] = PROJSYS_GaussKruger;
      zone[PROJZ_ZONE] = static_cast<INT16>(abs(gkzone));
      zone[PROJZ_ALTPROJ] = 1;
      mapproj = PROJECTION_TransverseMercator;
      tparm[2] = 1.0;
      tparm[4] = (zone[PROJZ_ZONE] <= 30) ? 6 * zone[PROJZ_ZONE] - 3 : 6 * zone[PROJZ_ZONE] - 363;
      tparm[5] = 0.0;
      tparm[TPARM_FalseEasting] = 500000.0;
      if (gkzone < 0 || southzone) {
         tparm[TPARM_FalseNorthing] = 10000000.0;
         }
      else {
         tparm[TPARM_FalseNorthing] = 0.0;
         }
      return (0);
      }

   //! Set coordinate system to Latitude-Longitude.
   void SetSystemLatLon (
      ) {
      memset(tparm,0,sizeof(tparm));
      zone[PROJZ_SYSTEM] = PROJSYS_LatLon;
      zone[PROJZ_ZONE] = 0;
      zone[PROJZ_ALTPROJ] = 0;
      mapproj = PROJECTION_LatLon;
      }

   //! Set parameters for Universal Transverse Mercator projection system.
   void SetSystemNZMG (
      ) {
      memset(static_cast<void*>(this),0,sizeof(*this));
      zone[PROJZ_SYSTEM] = PROJSYS_NZMG;
      zone[PROJZ_DATUM] = DATUM_NewZealand1949;
      mapproj = PROJECTION_NZMG;
      }

   //! Set parameters for Universal Polar Stereographic projection system.
   void SetSystemUPS (
      bool southpole = false
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UPS;
      SetProjectionPolarStereo(0.994,0.0,(southpole)?-90.0:90.0,2000000.0,2000000.0);
      return;
      }

   //! Set parameters for Universal Transverse Mercator projection system.
   ERRVALUE SetSystemUTM (
      INT16 utmzone,                   //!< Zone (1-60, or -1 - -60 as alternative for South zones)
      bool southzone = false           //!< South zones
      ) {
      if (utmzone < -60 || utmzone == 0 || utmzone > 60) return (EBadFuncParm);
      zone[PROJZ_SYSTEM] = PROJSYS_UTM;
      zone[PROJZ_ZONE] = static_cast<INT16>(abs(utmzone));
      mapproj = PROJECTION_TransverseMercator;
      tparm[2] = .9996;
      tparm[4] = 6 * zone[PROJZ_ZONE] - 183;
      tparm[5] = 0.0;
      tparm[TPARM_FalseEasting] = 500000.0;
      if (utmzone < 0 || southzone) {
         tparm[TPARM_FalseNorthing] = 10000000.0;
         }
      else {
         tparm[TPARM_FalseNorthing] = 0.0;
         }
      return (0);
      }

private:
   #ifndef GENERATING_DOXYGEN_OUTPUT
   ERRVALUE SetProj24567 (
      INT16 proj,
      double parm2,
      double parm4,
      double parm5,
      double parm6,
      double parm7
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = proj;
      tparm[2] = parm2;
      tparm[4] = parm4;
      tparm[5] = parm5;
      tparm[6] = parm6;
      tparm[7] = parm7;
      return (0);
      }

   ERRVALUE SetProj4567 (
      INT16 proj,
      double parm4,
      double parm5,
      double parm6,
      double parm7
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = proj;
      tparm[4] = parm4;
      tparm[5] = parm5;
      tparm[6] = parm6;
      tparm[7] = parm7;
      return (0);
      }

   ERRVALUE SetProj467 (
      INT16 proj,
      double parm4,
      double parm6,
      double parm7
      ) {
      zone[PROJZ_SYSTEM] = PROJSYS_UserDefined;
      zone[PROJZ_ZONE] = 0;
      mapproj = proj;
      tparm[4] = parm4;
      tparm[6] = parm6;
      tparm[7] = parm7;
      return (0);
      }

   CHECKSIZE(80+15*8);                 //! Make compiler validate structure size
   #endif // GENERATING_DOXYGEN_OUTPUT
   };

//! Comparison for equality.
inline bool operator== (
   const MAPPROJPARM& lhs,
   const MAPPROJPARM& rhs
   ) {
   return (lhs.mapproj == rhs.mapproj && !memcmp(lhs.zone,rhs.zone,3*sizeof(lhs.zone[0])) && !memcmp(lhs.tparm,rhs.tparm,sizeof(lhs.tparm)));
   }

//! Comparison for inequality.
inline bool operator!= (
   const MAPPROJPARM& lhs,
   const MAPPROJPARM& rhs
   ) {
   return (!operator==(lhs,rhs));
   }

#ifndef GENERATING_DOXYGEN_OUTPUT
//! Prevent unsafe usage.
#ifdef memcpy
//! Assume defined as 'amemcpy' so need definition to match.
extern void * (*amemcpy)(MAPPROJPARM*, MAPPROJPARM*, int);
extern void * (*amemcpy)(MAPPROJPARM*, const MAPPROJPARM*, int);
#else
void memcpy(MAPPROJPARM*, MAPPROJPARM*, int);
void memcpy(MAPPROJPARM*, const MAPPROJPARM*, int);
#endif
void memset(MAPPROJPARM*, int, int);
#endif // GENERATING_DOXYGEN_OUTPUT


/*----------------------------------------------------------------------------*/

#ifdef GEOMLIBEXPORT
#undef GEOMLIBEXPORT
#endif

#endif   // INC_MI32_MAPPROJ_H

---