EVOLUTION-MANAGER

Edit File: ogr_srs_proj4.cpp

/****************************************************************************** * * Project: OpenGIS Simple Features Reference Implementation * Purpose: OGRSpatialReference interface to PROJ.4. * Author: Frank Warmerdam <warmerdam@pobox.com> * ****************************************************************************** * Copyright (c) 1999, Les Technologies SoftMap Inc. * Copyright (c) 2008-2013, Even Rouault <even dot rouault at mines-paris dot org> * Copyright (c) 2014, Kyle Shannon <kyle at pobox dot com> * * Permission is hereby granted, free of charge, to any person obtaining a * copy of this software and associated documentation files (the "Software"), * to deal in the Software without restriction, including without limitation * the rights to use, copy, modify, merge, publish, distribute, sublicense, * and/or sell copies of the Software, and to permit persons to whom the * Software is furnished to do so, subject to the following conditions: * * The above copyright notice and this permission notice shall be included * in all copies or substantial portions of the Software. * * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER * DEALINGS IN THE SOFTWARE. ****************************************************************************/ #include "cpl_port.h" #include "ogr_srs_api.h" #include <cmath> #include <cstddef> #include <cstdlib> #include <cstring> #include <memory> #include <string> #include <vector> #include "cpl_conv.h" #include "cpl_error.h" #include "cpl_string.h" #include "ogr_core.h" #include "ogr_p.h" #include "ogr_spatialref.h" extern int EPSGGetWGS84Transform( int nGeogCS, std::vector<CPLString>& asTransform ); CPL_CVSID("$Id: ogr_srs_proj4.cpp 825a6a02f5f7031274c5e717aee82b24ff02d9b8 2018-08-18 18:54:43 +0200 Even Rouault $") /* -------------------------------------------------------------------- */ /* The following list comes from osrs/proj/src/pj_ellps.c. */ /* Please update from time to time. */ /* -------------------------------------------------------------------- */ static const char * const ogr_pj_ellps[] = { "MERIT", "a=6378137.0", "rf=298.257", "MERIT 1983", "SGS85", "a=6378136.0", "rf=298.257", "Soviet Geodetic System 85", "GRS80", "a=6378137.0", "rf=298.257222101", "GRS 1980(IUGG, 1980)", "IAU76", "a=6378140.0", "rf=298.257", "IAU 1976", "airy", "a=6377563.396", "b=6356256.910", "Airy 1830", "APL4.9", "a=6378137.0.", "rf=298.25", "Appl. Physics. 1965", "NWL9D", "a=6378145.0.", "rf=298.25", "Naval Weapons Lab., 1965", "mod_airy", "a=6377340.189", "b=6356034.446", "Modified Airy", "andrae", "a=6377104.43", "rf=300.0", "Andrae 1876 (Den., Iclnd.)", "aust_SA", "a=6378160.0", "rf=298.25", "Australian Natl & S. Amer. 1969", "GRS67", "a=6378160.0", "rf=298.2471674270", "GRS 67(IUGG 1967)", "bessel", "a=6377397.155", "rf=299.1528128", "Bessel 1841", "bess_nam", "a=6377483.865", "rf=299.1528128", "Bessel 1841 (Namibia)", "clrk66", "a=6378206.4", "b=6356583.8", "Clarke 1866", "clrk80", "a=6378249.145", "rf=293.4663", "Clarke 1880 mod.", "CPM", "a=6375738.7", "rf=334.29", "Comm. des Poids et Mesures 1799", "delmbr", "a=6376428.", "rf=311.5", "Delambre 1810 (Belgium)", "engelis", "a=6378136.05", "rf=298.2566", "Engelis 1985", "evrst30", "a=6377276.345", "rf=300.8017", "Everest 1830", "evrst48", "a=6377304.063", "rf=300.8017", "Everest 1948", "evrst56", "a=6377301.243", "rf=300.8017", "Everest 1956", "evrst69", "a=6377295.664", "rf=300.8017", "Everest 1969", "evrstSS", "a=6377298.556", "rf=300.8017", "Everest (Sabah & Sarawak)", "fschr60", "a=6378166.", "rf=298.3", "Fischer (Mercury Datum) 1960", "fschr60m", "a=6378155.", "rf=298.3", "Modified Fischer 1960", "fschr68", "a=6378150.", "rf=298.3", "Fischer 1968", "helmert", "a=6378200.", "rf=298.3", "Helmert 1906", "hough", "a=6378270.0", "rf=297.", "Hough", "intl", "a=6378388.0", "rf=297.", "International 1909 (Hayford)", "krass", "a=6378245.0", "rf=298.3", "Krassovsky, 1942", "kaula", "a=6378163.", "rf=298.24", "Kaula 1961", "lerch", "a=6378139.", "rf=298.257", "Lerch 1979", "mprts", "a=6397300.", "rf=191.", "Maupertius 1738", "new_intl", "a=6378157.5", "b=6356772.2", "New International 1967", "plessis", "a=6376523.", "b=6355863.", "Plessis 1817 (France)", "SEasia", "a=6378155.0", "b=6356773.3205", "Southeast Asia", "walbeck", "a=6376896.0", "b=6355834.8467", "Walbeck", "WGS60", "a=6378165.0", "rf=298.3", "WGS 60", "WGS66", "a=6378145.0", "rf=298.25", "WGS 66", "WGS72", "a=6378135.0", "rf=298.26", "WGS 72", "WGS84", "a=6378137.0", "rf=298.257223563", "WGS 84", "sphere", "a=6370997.0", "b=6370997.0", "Normal Sphere (r=6370997)", nullptr, nullptr, nullptr, nullptr, }; typedef struct { const char* pszPJ; const char* pszOGR; int nEPSG; int nGCS; } OGRProj4Datum; // Derived from proj/src/pj_datum.c. // WGS84, NAD27 and NAD83 are directly hard-coded in the code. static const OGRProj4Datum ogr_pj_datums[] = { { "GGRS87", "Greek_Geodetic_Reference_System_1987", 4121, 6121}, { "potsdam", "Deutsches_Hauptdreiecksnetz", 4314, 6314}, { "carthage", "Carthage", 4223, 6223}, { "hermannskogel", "Militar_Geographische_Institut", 4312, 6312}, { "ire65", "TM65", 4299, 6299}, { "nzgd49", "New_Zealand_Geodetic_Datum_1949", 4272, 6272}, { "OSGB36", "OSGB_1936", 4277, 6277} }; typedef struct { const char* pszProj4PMName; const char* pszWKTPMName; const char* pszFromGreenwich; int nPMCode; } OGRProj4PM; // Derived from pj_datums.c. static const OGRProj4PM ogr_pj_pms [] = { { "greenwich", "Greenwich", "0dE", 8901 }, { "lisbon", "Lisbon", "9d07'54.862\"W", 8902 }, { "paris", "Paris", "2d20'14.025\"E", 8903 }, { "bogota", "Bogota", "74d04'51.3\"W", 8904 }, { "madrid", "Madrid", "3d41'16.58\"W", 8905 }, { "rome", "Rome", "12d27'8.4\"E", 8906 }, { "bern", "Bern", "7d26'22.5\"E", 8907 }, { "jakarta", "Jakarta", "106d48'27.79\"E", 8908 }, { "ferro", "Ferro", "17d40'W", 8909 }, { "brussels", "Brussels", "4d22'4.71\"E", 8910 }, { "stockholm", "Stockholm", "18d3'29.8\"E", 8911 }, { "athens", "Athens", "23d42'58.815\"E", 8912 }, { "oslo", "Oslo", "10d43'22.5\"E", 8913 } }; static const char* OGRGetProj4Datum( const char* pszDatum, int nEPSGDatum ) { for( size_t i = 0; i < CPL_ARRAYSIZE(ogr_pj_datums); i++ ) { if( nEPSGDatum == ogr_pj_datums[i].nGCS || EQUAL(pszDatum, ogr_pj_datums[i].pszOGR) ) { return ogr_pj_datums[i].pszPJ; } } return nullptr; } static const OGRProj4PM* OGRGetProj4PMFromProj4Name( const char* pszProj4PMName ) { for( size_t i = 0; i < CPL_ARRAYSIZE(ogr_pj_pms); i++ ) { if( EQUAL(pszProj4PMName, ogr_pj_pms[i].pszProj4PMName) ) { return &ogr_pj_pms[i]; } } return nullptr; } static const OGRProj4PM* OGRGetProj4PMFromCode( int nPMCode ) { for( size_t i = 0; i < CPL_ARRAYSIZE(ogr_pj_pms); i++ ) { if( nPMCode == ogr_pj_pms[i].nPMCode ) { return &ogr_pj_pms[i]; } } return nullptr; } static const OGRProj4PM* OGRGetProj4PMFromVal( double dfVal ) { for( unsigned int i = 0; i < CPL_ARRAYSIZE(ogr_pj_pms); i++ ) { // TODO(schwehr): Use an almost equal call. if( fabs(dfVal - CPLDMSToDec(ogr_pj_pms[i].pszFromGreenwich)) < 1e-10 ) { return &ogr_pj_pms[i]; } } return nullptr; } typedef struct { const char* pszWKTName; const char* pszValueInMeter; const char* pszProjName; } LinearUnitsStruct; static const LinearUnitsStruct asLinearUnits [] = { { SRS_UL_METER, "1.0", "m" }, { SRS_UL_METER, "1.0", "meter" }, // alias of former { SRS_UL_METER, "1.0", "metre" }, // alias of former { "metre", "1.0", "m" }, // alias of former // Leave as 'kilometre' instead of SRS_UL_KILOMETER due to historical usage. { "kilometre", SRS_UL_KILOMETER_CONV, "km" }, { SRS_UL_KILOMETER, SRS_UL_KILOMETER_CONV, "km" }, // alias of former { SRS_UL_DECIMETER, SRS_UL_DECIMETER_CONV, "dm" }, { SRS_UL_CENTIMETER, SRS_UL_CENTIMETER_CONV, "cm" }, { SRS_UL_MILLIMETER, SRS_UL_MILLIMETER_CONV, "mm" }, // Leave as 'Foot (International)' or SRS_UL_FOOT instead of SRS_UL_INTL_FOOT // due to historical usage. { SRS_UL_FOOT, SRS_UL_FOOT_CONV, "ft" }, { SRS_UL_INTL_FOOT, SRS_UL_INTL_FOOT_CONV, "ft" }, // alias of former { SRS_UL_US_FOOT, SRS_UL_US_FOOT_CONV, "us-ft" }, { SRS_UL_INDIAN_FOOT, SRS_UL_INDIAN_FOOT_CONV, "ind-ft" }, { SRS_UL_INTL_NAUT_MILE, SRS_UL_INTL_NAUT_MILE_CONV, "kmi" }, { SRS_UL_NAUTICAL_MILE, SRS_UL_NAUTICAL_MILE_CONV, "kmi" }, // alias of former { SRS_UL_INTL_STAT_MILE, SRS_UL_INTL_STAT_MILE_CONV, "mi" }, { "Mile", SRS_UL_INTL_STAT_MILE_CONV, "mi" }, // alias of former { "IMILE", SRS_UL_INTL_STAT_MILE_CONV, "mi" }, // alias of former { SRS_UL_US_STAT_MILE, SRS_UL_US_STAT_MILE_CONV, "us-mi"}, { SRS_UL_INTL_LINK, SRS_UL_INTL_LINK_CONV, "link" }, { SRS_UL_LINK, SRS_UL_LINK_CONV, "link" }, // alias of former { SRS_UL_INTL_YARD, SRS_UL_INTL_YARD_CONV, "yd" }, { "IYARD", SRS_UL_INTL_YARD_CONV, "yd" }, // alias of former { SRS_UL_US_YARD, SRS_UL_US_YARD_CONV, "us-yd" }, { SRS_UL_INDIAN_YARD, SRS_UL_INDIAN_YARD_CONV, "ind-yd" }, { SRS_UL_INTL_INCH, SRS_UL_INTL_INCH_CONV, "in" }, { SRS_UL_US_INCH, SRS_UL_US_INCH_CONV, "us-in" }, { SRS_UL_INTL_FATHOM, SRS_UL_INTL_FATHOM_CONV, "fath" }, { SRS_UL_INTL_CHAIN, SRS_UL_INTL_CHAIN_CONV, "ch" }, { SRS_UL_US_CHAIN, SRS_UL_US_CHAIN_CONV, "us-ch" }, { SRS_UL_INDIAN_CHAIN, SRS_UL_INDIAN_CHAIN_CONV, "ind-ch" }, // { SRS_UL_ROD, SRS_UL_ROD_CONV, "????" } }; /************************************************************************/ /* GetLinearFromLinearConvOrName() */ /************************************************************************/ static const LinearUnitsStruct *GetLinearFromLinearConvOrName( double dfLinearConv, const char *pszLinearUnits ) { for( size_t i = 0; i < CPL_ARRAYSIZE(asLinearUnits); i++ ) { if( (pszLinearUnits != nullptr && EQUAL(pszLinearUnits, asLinearUnits[i].pszWKTName)) || fabs(dfLinearConv - CPLAtof(asLinearUnits[i].pszValueInMeter)) < 0.00000001 ) { return &(asLinearUnits[i]); } } return nullptr; } /************************************************************************/ /* GetLinearFromProjName() */ /************************************************************************/ static const LinearUnitsStruct* GetLinearFromProjName( const char* pszProjName ) { for( size_t i = 0; i < CPL_ARRAYSIZE(asLinearUnits); i++ ) { if( EQUAL(pszProjName, asLinearUnits[i].pszProjName) ) { return &(asLinearUnits[i]); } } return nullptr; } /************************************************************************/ /* OSRProj4Tokenize() */ /* */ /* Custom tokenizing function for PROJ.4 strings. The main */ /* reason we can't just use CSLTokenizeString is to handle */ /* strings with a + sign in the exponents of parameter values. */ /************************************************************************/ static char **OSRProj4Tokenize( const char *pszFull ) { if( pszFull == nullptr ) return nullptr; char *pszFullWrk = CPLStrdup( pszFull ); char *pszStart = nullptr; char **papszTokens = nullptr; for( int i=0; pszFullWrk[i] != '\0'; i++ ) { switch( pszFullWrk[i] ) { case '+': if( i == 0 || pszFullWrk[i-1] == '\0' ) { if( pszStart != nullptr ) { if( strstr(pszStart, "=") != nullptr ) papszTokens = CSLAddString( papszTokens, pszStart ); else { CPLString osAsBoolean = pszStart; osAsBoolean += "=yes"; papszTokens = CSLAddString( papszTokens, osAsBoolean ); } } pszStart = pszFullWrk + i + 1; } break; case ' ': case '\t': case '\n': pszFullWrk[i] = '\0'; break; default: break; } } if( pszStart != nullptr && strlen(pszStart) > 0 ) papszTokens = CSLAddString( papszTokens, pszStart ); CPLFree( pszFullWrk ); return papszTokens; } /************************************************************************/ /* OSRImportFromProj4() */ /************************************************************************/ /** * \brief Import PROJ.4 coordinate string. * * This function is the same as OGRSpatialReference::importFromProj4(). */ OGRErr OSRImportFromProj4( OGRSpatialReferenceH hSRS, const char *pszProj4 ) { VALIDATE_POINTER1( hSRS, "OSRImportFromProj4", OGRERR_FAILURE ); return reinterpret_cast<OGRSpatialReference *>(hSRS)-> importFromProj4( pszProj4 ); } /************************************************************************/ /* OSR_GDV() */ /* */ /* Fetch a particular parameter out of the parameter list, or */ /* the indicated default if it isn't available. This is a */ /* helper function for importFromProj4(). */ /************************************************************************/ static double OSR_GDV( char **papszNV, const char * pszField, double dfDefaultValue ) { const char *pszValue = CSLFetchNameValue( papszNV, pszField ); // Special hack to use k_0 if available. if( pszValue == nullptr && EQUAL(pszField, "k") ) pszValue = CSLFetchNameValue( papszNV, "k_0" ); if( pszValue == nullptr ) return dfDefaultValue; return CPLDMSToDec(pszValue); } /************************************************************************/ /* importFromProj4() */ /************************************************************************/ /** * \brief Import PROJ.4 coordinate string. * * The OGRSpatialReference is initialized from the passed PROJ.4 style * coordinate system string. In addition to many +proj formulations which * have OGC equivalents, it is also possible to import "+init=epsg:n" style * definitions. These are passed to importFromEPSG(). Other init strings * (such as the state plane zones) are not currently supported. * * Example: * pszProj4 = "+proj=utm +zone=11 +datum=WGS84" * * Some parameters, such as grids, recognized by PROJ.4 may not be well * understood and translated into the OGRSpatialReference model. It is possible * to add the +wktext parameter which is a special keyword that OGR recognized * as meaning "embed the entire PROJ.4 string in the WKT and use it literally * when converting back to PROJ.4 format". * * For example: * "+proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000 +y_0=6023150 +ellps=intl * +units=m +nadgrids=nzgd2kgrid0005.gsb +wktext" * * will be translated as : * \code * PROJCS["unnamed", * GEOGCS["International 1909 (Hayford)", * DATUM["unknown", * SPHEROID["intl",6378388,297]], * PRIMEM["Greenwich",0], * UNIT["degree",0.0174532925199433]], * PROJECTION["New_Zealand_Map_Grid"], * PARAMETER["latitude_of_origin",-41], * PARAMETER["central_meridian",173], * PARAMETER["false_easting",2510000], * PARAMETER["false_northing",6023150], * UNIT["Meter",1], * EXTENSION["PROJ4","+proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000 * +y_0=6023150 +ellps=intl +units=m +nadgrids=nzgd2kgrid0005.gsb +wktext"]] * \endcode * * Special processing for 'etmerc' (GDAL &gt;= 1.10 ): if +proj=etmerc is found * in the passed string, the SRS built will use the WKT representation for a * standard Transverse Mercator, but will aso include a PROJ4 EXTENSION node to * preserve the etmerc projection method. * * For example: * "+proj=etmerc +lat_0=0 +lon_0=9 +k=0.9996 +units=m +x_0=500000 +datum=WGS84" * * will be translated as : * \code * PROJCS["unnamed", * GEOGCS["WGS 84", * DATUM["WGS_1984", * SPHEROID["WGS 84",6378137,298.257223563, * AUTHORITY["EPSG","7030"]], * TOWGS84[0,0,0,0,0,0,0], * AUTHORITY["EPSG","6326"]], * PRIMEM["Greenwich",0, * AUTHORITY["EPSG","8901"]], * UNIT["degree",0.0174532925199433, * AUTHORITY["EPSG","9108"]], * AUTHORITY["EPSG","4326"]], * PROJECTION["Transverse_Mercator"], * PARAMETER["latitude_of_origin",0], * PARAMETER["central_meridian",9], * PARAMETER["scale_factor",0.9996], * PARAMETER["false_easting",500000], * PARAMETER["false_northing",0], * UNIT["Meter",1], * EXTENSION["PROJ4","+proj=etmerc +lat_0=0 +lon_0=9 +k=0.9996 +units=m +x_0=500000 +datum=WGS84 +nodefs"]] * \endcode * * This method is the equivalent of the C function OSRImportFromProj4(). * * @param pszProj4 the PROJ.4 style string. * * @return OGRERR_NONE on success or OGRERR_CORRUPT_DATA on failure. */ OGRErr OGRSpatialReference::importFromProj4( const char * pszProj4 ) { /* -------------------------------------------------------------------- */ /* Clear any existing definition. */ /* -------------------------------------------------------------------- */ Clear(); /* -------------------------------------------------------------------- */ /* Strip any newlines or other "funny" stuff that might occur */ /* if this string just came from reading a file. */ /* -------------------------------------------------------------------- */ char *pszCleanCopy = CPLStrdup( pszProj4 ); for( int i = 0; pszCleanCopy[i] != '\0'; i++ ) { if( pszCleanCopy[i] == 10 || pszCleanCopy[i] == 13 || pszCleanCopy[i] == 9 ) pszCleanCopy[i] = ' '; } const char* pszInitEpsgCleanCopy = strstr(pszCleanCopy, "init=epsg:"); bool bSetAuthorityCode = true; // If there's an override, then drop the authority code. if( pszInitEpsgCleanCopy != nullptr && strchr(pszInitEpsgCleanCopy, '+') != nullptr ) { bSetAuthorityCode = false; } /* -------------------------------------------------------------------- */ /* Try to normalize the definition. This should expand +init= */ /* clauses and so forth. */ /* -------------------------------------------------------------------- */ char *pszNormalized = OCTProj4Normalize( pszCleanCopy ); // Workaround proj.4 bug (#239) by manually re-adding no_off/no_uoff. if( strstr(pszCleanCopy, "+no_off") != nullptr && strstr(pszNormalized, "+no_off") == nullptr ) { char* pszTmp = CPLStrdup(CPLSPrintf("%s +no_off", pszNormalized)); CPLFree(pszNormalized); pszNormalized = pszTmp; } else if( strstr(pszCleanCopy, "+no_uoff") != nullptr && strstr(pszNormalized, "+no_uoff") == nullptr ) { char* pszTmp = CPLStrdup(CPLSPrintf("%s +no_uoff", pszNormalized)); CPLFree(pszNormalized); pszNormalized = pszTmp; } CPLFree( pszCleanCopy ); /* -------------------------------------------------------------------- */ /* If we have an EPSG based init string, and no existing +proj */ /* portion then try to normalize into into a PROJ.4 string. */ /* This can happen if the proj.4 epsg dictionary is missing. */ /* -------------------------------------------------------------------- */ const char* pszInitEpsg = strstr(pszNormalized, "init=epsg:"); if( pszInitEpsg != nullptr && strstr(pszNormalized, "proj=") == nullptr ) { const char *pszNumber = pszInitEpsg + strlen("init=epsg:"); const OGRErr eErr = importFromEPSG( atoi(pszNumber) ); if( eErr != OGRERR_NONE || strchr(pszNumber, '+') == nullptr ) { CPLFree( pszNormalized ); return eErr; } const int nIdx = GetRoot()->FindChild("AUTHORITY"); if( nIdx >= 0 ) { GetRoot()->DestroyChild( nIdx ); } } /* -------------------------------------------------------------------- */ /* Parse the PROJ.4 string into a cpl_string.h style name/value */ /* list. */ /* -------------------------------------------------------------------- */ char **papszTokens = OSRProj4Tokenize( pszNormalized ); CPLFree( pszNormalized ); char **papszNV = nullptr; for( int i = 0; papszTokens != nullptr && papszTokens[i] != nullptr; i++ ) { char *pszEqual = strstr(papszTokens[i], "="); if( pszEqual == nullptr ) papszNV = CSLAddNameValue(papszNV, papszTokens[i], "" ); else { pszEqual[0] = '\0'; papszNV = CSLAddNameValue( papszNV, papszTokens[i], pszEqual+1 ); } } CSLDestroy( papszTokens ); /* -------------------------------------------------------------------- */ /* Extract the prime meridian, if there is one set. */ /* -------------------------------------------------------------------- */ const char *pszPM = CSLFetchNameValue( papszNV, "pm" ); double l_dfFromGreenwich = 0.0; if( pszPM != nullptr ) { const OGRProj4PM* psProj4PM = OGRGetProj4PMFromProj4Name(pszPM); if( psProj4PM ) { l_dfFromGreenwich = CPLDMSToDec(psProj4PM->pszFromGreenwich); pszPM = psProj4PM->pszWKTPMName; } else { l_dfFromGreenwich = CPLDMSToDec( pszPM ); pszPM = "unnamed"; } } else { pszPM = "Greenwich"; } /* -------------------------------------------------------------------- */ /* Operate on the basis of the projection name. */ /* -------------------------------------------------------------------- */ const char *pszProj = CSLFetchNameValue(papszNV, "proj"); bool bAddProj4Extension = false; if( pszProj == nullptr ) { CPLDebug( "OGR_PROJ4", "Can't find +proj= in:\n%s", pszProj4 ); CSLDestroy( papszNV ); return OGRERR_CORRUPT_DATA; } else if( EQUAL(pszProj, "longlat") || EQUAL(pszProj, "latlong") ) { } else if( EQUAL(pszProj, "geocent") ) { SetGeocCS( "Geocentric" ); } else if( EQUAL(pszProj, "bonne") ) { SetBonne( OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "cass") ) { SetCS( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "nzmg") ) { SetNZMG( OSR_GDV( papszNV, "lat_0", -41.0 ), OSR_GDV( papszNV, "lon_0", 173.0 ), OSR_GDV( papszNV, "x_0", 2510000.0 ), OSR_GDV( papszNV, "y_0", 6023150.0 ) ); } else if( EQUAL(pszProj, "cea") ) { SetCEA( OSR_GDV( papszNV, "lat_ts", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "tmerc") ) { const char *pszAxis = CSLFetchNameValue( papszNV, "axis" ); if( pszAxis == nullptr || !EQUAL(pszAxis, "wsu") ) SetTM( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); else SetTMSO( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } // For etmerc, we translate it into standard TM for the WKT // point of view, but make sure that the original proj.4 // definition is preserved for accurate reprojection. else if( EQUAL(pszProj, "etmerc") && CSLFetchNameValue( papszNV, "axis" ) == nullptr ) { bAddProj4Extension = true; SetTM( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "utm") ) { const double dfZone = OSR_GDV(papszNV, "zone", 0.0); const double dfSouth = OSR_GDV(papszNV, "south", 1.0); if( dfZone > 60.0 || dfZone < 0.0 || CPLIsNan(dfZone) || dfSouth < 0.0 || dfSouth > 1.0 || CPLIsNan(dfSouth) ) { CSLDestroy(papszNV); return OGRERR_CORRUPT_DATA; } SetUTM(static_cast<int>(dfZone), static_cast<int>(dfSouth)); } else if( EQUAL(pszProj, "merc") // 2SP form. && OSR_GDV(papszNV, "lat_ts", 1000.0) < 999.0 ) { SetMercator2SP( OSR_GDV( papszNV, "lat_ts", 0.0 ), 0.0, OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "merc") ) // 1SP form. { SetMercator( 0.0, OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "stere") && std::abs(OSR_GDV( papszNV, "lat_0", 0.0 ) - 90) < 0.001 ) { SetPS( OSR_GDV( papszNV, "lat_ts", 90.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "stere") && std::abs(OSR_GDV( papszNV, "lat_0", 0.0 ) + 90) < 0.001 ) { SetPS( OSR_GDV( papszNV, "lat_ts", -90.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "sterea") ) { SetOS( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "stere") ) { SetStereographic( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "eqc") ) { if( OSR_GDV( papszNV, "lat_ts", 0.0 ) != 0.0 ) SetEquirectangular2( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "lat_ts", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); else SetEquirectangular( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "gstmerc") ) { SetGaussSchreiberTMercator( OSR_GDV( papszNV, "lat_0", -21.116666667 ), OSR_GDV( papszNV, "lon_0", 55.53333333309), OSR_GDV( papszNV, "k_0", 1.0 ), OSR_GDV( papszNV, "x_0", 160000.000 ), OSR_GDV( papszNV, "y_0", 50000.000 ) ); } else if( EQUAL(pszProj, "gnom") ) { SetGnomonic( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "ortho") ) { SetOrthographic( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "laea") ) { SetLAEA( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "aeqd") ) { SetAE( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "eqdc") ) { SetEC( OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lat_2", 0.0 ), OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "mill") ) { SetMC( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "moll") ) { SetMollweide( OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "eck1") || EQUAL(pszProj, "eck2") || EQUAL(pszProj, "eck3") || EQUAL(pszProj, "eck4") || EQUAL(pszProj, "eck5") || EQUAL(pszProj, "eck6") ) { SetEckert( pszProj[3] - '0', // TODO(schwehr): Symbolic const for 3. OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "poly") ) { SetPolyconic( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "aea") ) { SetACEA( OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lat_2", 0.0 ), OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "robin") ) { SetRobinson( OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "vandg") ) { SetVDG( OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "sinu") ) { SetSinusoidal( OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "gall") ) { SetGS( OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "goode") ) { SetGH( OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "igh") ) { SetIGH(); } else if( EQUAL(pszProj, "geos") ) { SetGEOS( OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "h", 35785831.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); if( EQUAL(CSLFetchNameValueDef( papszNV, "sweep", "y" ), "x") ) bAddProj4Extension = true; } else if( EQUAL(pszProj, "lcc") ) { if( OSR_GDV(papszNV, "lat_0", 0.0 ) == OSR_GDV(papszNV, "lat_1", 0.0 ) && CSLFetchNameValue( papszNV, "lat_2" ) == nullptr ) { // 1SP form. SetLCC1SP( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "k_0", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else { // 2SP form. SetLCC( OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lat_2", 0.0 ), OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } } else if( EQUAL(pszProj, "omerc") ) { if( CSLFetchNameValue(papszNV, "no_uoff") != nullptr || CSLFetchNameValue(papszNV, "no_off") != nullptr ) { // From PJ_omerc, when alpha is defined but not gamma // the default gamma value is alpha. // if( alp || gam ) // { // if( alp ) // { // gamma0 = asin(sin(alpha_c) / D); // if( !gam ) // gamma = alpha_c; SetHOM( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lonc", 0.0 ), OSR_GDV( papszNV, "alpha", 0.0 ), OSR_GDV( papszNV, "gamma", OSR_GDV(papszNV, "alpha", 0.0 )), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( CSLFetchNameValue(papszNV, "lat_1") && CSLFetchNameValue(papszNV, "lon_1") && CSLFetchNameValue(papszNV, "lat_2") && CSLFetchNameValue(papszNV, "lon_2") ) { SetHOM2PNO( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lon_1", 0.0 ), OSR_GDV( papszNV, "lat_2", 0.0 ), OSR_GDV( papszNV, "lon_2", 0.0 ), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else { SetHOMAC( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lonc", 0.0 ), OSR_GDV( papszNV, "alpha", 0.0 ), OSR_GDV( papszNV, "gamma", OSR_GDV(papszNV, "alpha", 0.0 )), OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } } else if( EQUAL(pszProj, "somerc") ) { SetHOMAC( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), 90.0, 90.0, OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "krovak") ) { SetKrovak( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), 30.28813972222222, // Alpha. ignored by PROJ 78.5, // Pseudo_standard_parallel_1. ignored by PROJ OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "imw_p") ) { SetIWMPolyconic( OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lat_2", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "wag1") ) { SetWagner( 1, 0.0, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "wag2") ) { SetWagner( 2, 0.0, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "wag3") ) { SetWagner( 3, OSR_GDV( papszNV, "lat_ts", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "wag4") ) { SetWagner( 4, 0.0, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "wag5") ) { SetWagner( 5, 0.0, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "wag6") ) { SetWagner( 6, 0.0, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "wag7") ) { SetWagner( 7, 0.0, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj, "qsc") ) { SetQSC( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ) ); } else if( EQUAL(pszProj, "sch") ) { SetSCH( OSR_GDV( papszNV, "plat_0", 0.0 ), OSR_GDV( papszNV, "plon_0", 0.0 ), OSR_GDV( papszNV, "phdg_0", 0.0 ), OSR_GDV( papszNV, "h_0", 0.0) ); } else if( EQUAL(pszProj, "tpeqd") ) { SetTPED( OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lon_1", 0.0 ), OSR_GDV( papszNV, "lat_2", 0.0 ), OSR_GDV( papszNV, "lon_2", 0.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( strstr(pszProj4, "wktext") != nullptr ) { // Fake out a projected coordinate system for otherwise // unrecognized projections for which we are already planning // to embed the actual PROJ.4 string via extension node. SetProjection( "custom_proj4" ); } else { CPLDebug( "OGR_PROJ4", "Unsupported projection: %s", pszProj ); CSLDestroy( papszNV ); return OGRERR_CORRUPT_DATA; } /* -------------------------------------------------------------------- */ /* Try to translate the datum. */ /* -------------------------------------------------------------------- */ bool bFullyDefined = false; const char *pszValue = CSLFetchNameValue(papszNV, "datum"); if( pszValue == nullptr ) { // Do nothing. } else if( (EQUAL(pszValue, "NAD27") || EQUAL(pszValue, "NAD83") || EQUAL(pszValue, "WGS84") || EQUAL(pszValue, "WGS72")) && l_dfFromGreenwich == 0.0 ) { SetWellKnownGeogCS( pszValue ); bFullyDefined = true; } else { for( size_t i = 0; i < CPL_ARRAYSIZE(ogr_pj_datums); i++ ) { if( EQUAL(pszValue, ogr_pj_datums[i].pszPJ) ) { OGRSpatialReference oGCS; if( oGCS.importFromEPSG( ogr_pj_datums[i].nEPSG ) == OGRERR_NONE && CopyGeogCSFrom( &oGCS ) == OGRERR_NONE ) { bFullyDefined = true; } break; } } // If we don't recognise the datum, we ignore it. } /* -------------------------------------------------------------------- */ /* Set the ellipsoid information. */ /* -------------------------------------------------------------------- */ double dfSemiMajor = 0.0; double dfInvFlattening = 0.0; double dfSemiMinor = 0.0; pszValue = CSLFetchNameValue(papszNV, "ellps"); if( pszValue != nullptr && !bFullyDefined ) { for( int i = 0; ogr_pj_ellps[i] != nullptr; i += 4 ) { if( !EQUAL(ogr_pj_ellps[i], pszValue) ) continue; CPLAssert( STARTS_WITH_CI(ogr_pj_ellps[i+1], "a=") ); dfSemiMajor = CPLAtof(ogr_pj_ellps[i+1]+2); if( STARTS_WITH_CI(ogr_pj_ellps[i+2], "rf=") ) dfInvFlattening = CPLAtof(ogr_pj_ellps[i+2]+3); else { CPLAssert( STARTS_WITH_CI(ogr_pj_ellps[i+2], "b=") ); dfSemiMinor = CPLAtof(ogr_pj_ellps[i+2]+2); dfInvFlattening = OSRCalcInvFlattening(dfSemiMajor, dfSemiMinor); } SetGeogCS( ogr_pj_ellps[i+3], "unknown", ogr_pj_ellps[i], dfSemiMajor, dfInvFlattening, pszPM, l_dfFromGreenwich ); bFullyDefined = true; break; } } if( !bFullyDefined ) { dfSemiMajor = OSR_GDV( papszNV, "a", 0.0 ); if( dfSemiMajor == 0.0 ) { dfSemiMajor = OSR_GDV( papszNV, "R", 0.0 ); if( dfSemiMajor != 0.0 ) { dfSemiMinor = -1.0; dfInvFlattening = 0.0; } else { CPLDebug("OGR_PROJ4", "Can't find ellipse definition, default to WGS84:\n%s", pszProj4); dfSemiMajor = SRS_WGS84_SEMIMAJOR; dfSemiMinor = -1.0; dfInvFlattening = SRS_WGS84_INVFLATTENING; } } else { dfSemiMinor = OSR_GDV( papszNV, "b", -1.0 ); dfInvFlattening = OSR_GDV( papszNV, "rf", -1.0 ); if( dfSemiMinor == -1.0 && dfInvFlattening == -1.0 ) { const double dfFlattening = OSR_GDV( papszNV, "f", -1.0 ); if( dfFlattening == 0.0 ) dfSemiMinor = dfSemiMajor; else if( dfFlattening != -1.0 ) dfInvFlattening = 1.0 / dfFlattening; } } if( dfSemiMinor == -1.0 && dfInvFlattening == -1.0 ) { CPLDebug( "OGR_PROJ4", "Can't find ellipse definition in:\n%s", pszProj4 ); CSLDestroy( papszNV ); return OGRERR_UNSUPPORTED_SRS; } if( dfInvFlattening == -1.0 ) { dfInvFlattening = OSRCalcInvFlattening(dfSemiMajor, dfSemiMinor); } SetGeogCS( "unnamed ellipse", "unknown", "unnamed", dfSemiMajor, dfInvFlattening, pszPM, l_dfFromGreenwich ); // bFullyDefined = true; } /* -------------------------------------------------------------------- */ /* Handle TOWGS84 conversion. */ /* -------------------------------------------------------------------- */ pszValue = CSLFetchNameValue(papszNV, "towgs84"); // Make sure that we do not include a useless TOWGS84 node if the datum // is WGS84. Depending on the proj.4 version +datum=WGS84 might be // expanded with a +towgs84=0,0,0 or not. if( pszValue != nullptr && !( EQUAL(CSLFetchNameValueDef(papszNV, "datum", ""), "WGS84") && EQUAL(pszValue, "0,0,0")) ) { char **papszToWGS84 = CSLTokenizeStringComplex( pszValue, ",", FALSE, TRUE ); if( CSLCount(papszToWGS84) >= 7 ) SetTOWGS84( CPLAtof(papszToWGS84[0]), CPLAtof(papszToWGS84[1]), CPLAtof(papszToWGS84[2]), CPLAtof(papszToWGS84[3]), CPLAtof(papszToWGS84[4]), CPLAtof(papszToWGS84[5]), CPLAtof(papszToWGS84[6]) ); else if( CSLCount(papszToWGS84) >= 3 ) SetTOWGS84( CPLAtof(papszToWGS84[0]), CPLAtof(papszToWGS84[1]), CPLAtof(papszToWGS84[2]) ); else CPLError( CE_Warning, CPLE_AppDefined, "Seemingly corrupt +towgs84 option (%s), ignoring.", pszValue ); CSLDestroy(papszToWGS84); } /* -------------------------------------------------------------------- */ /* Handle nadgrids via an extension node. */ /* -------------------------------------------------------------------- */ pszValue = CSLFetchNameValue(papszNV, "nadgrids"); if( pszValue != nullptr ) { SetExtension( "DATUM", "PROJ4_GRIDS", pszValue ); FixupOrdering(); } /* -------------------------------------------------------------------- */ /* Linear units translation */ /* -------------------------------------------------------------------- */ if( IsProjected() || IsLocal() || IsGeocentric() ) { pszValue = CSLFetchNameValue(papszNV, "to_meter"); if( pszValue != nullptr && CPLAtofM(pszValue) > 0.0 ) { const double dfValue = CPLAtofM(pszValue); const LinearUnitsStruct* psLinearUnits = GetLinearFromLinearConvOrName( dfValue, pszValue ); if( psLinearUnits != nullptr ) { SetLinearUnits( psLinearUnits->pszWKTName, CPLAtof(psLinearUnits->pszValueInMeter) ); } else SetLinearUnits( "unknown", CPLAtofM(pszValue) ); } // All units reported by cs2cs -lu are supported, fall back to meter. else if( (pszValue = CSLFetchNameValue(papszNV, "units")) != nullptr ) { const LinearUnitsStruct* psLinearUnits = GetLinearFromProjName( pszValue ); if( psLinearUnits != nullptr ) { SetLinearUnits( psLinearUnits->pszWKTName, CPLAtof(psLinearUnits->pszValueInMeter) ); } else { // This case is untranslatable. Should add all proj.4 unts. SetLinearUnits( !EQUAL(pszValue, "AUTHORITY") ? pszValue : "unknown", 1.0 ); } } } /* -------------------------------------------------------------------- */ /* Adjust linear parameters into PROJCS units if the linear */ /* units are not meters. */ /* -------------------------------------------------------------------- */ if( GetLinearUnits() != 1.0 && IsProjected() ) { OGR_SRSNode *poPROJCS = GetAttrNode( "PROJCS" ); for( int i = 0; i < poPROJCS->GetChildCount(); i++ ) { OGR_SRSNode *poParm = poPROJCS->GetChild(i); if( !EQUAL(poParm->GetValue(), "PARAMETER") || poParm->GetChildCount() != 2 ) continue; const char *pszParmName = poParm->GetChild(0)->GetValue(); if( IsLinearParameter(pszParmName) ) SetNormProjParm(pszParmName, GetProjParm(pszParmName)); } } /* -------------------------------------------------------------------- */ /* Handle geoidgrids via an extension node and COMPD_CS. */ /* -------------------------------------------------------------------- */ pszValue = CSLFetchNameValue(papszNV, "geoidgrids"); OGR_SRSNode *poVERT_CS = nullptr; if( pszValue != nullptr ) { OGR_SRSNode *poHorizSRS = GetRoot()->Clone(); Clear(); CPLString osName = poHorizSRS->GetChild(0)->GetValue(); osName += " + "; osName += "Unnamed Vertical Datum"; SetNode( "COMPD_CS", osName ); GetRoot()->AddChild( poHorizSRS ); poVERT_CS = new OGR_SRSNode( "VERT_CS" ); GetRoot()->AddChild( poVERT_CS ); poVERT_CS->AddChild( new OGR_SRSNode( "Unnamed" ) ); CPLString osTarget = GetRoot()->GetValue(); osTarget += "|VERT_CS|VERT_DATUM"; SetNode( osTarget, "Unnamed" ); poVERT_CS->GetChild(1)->AddChild( new OGR_SRSNode( "2005" ) ); SetExtension( osTarget, "PROJ4_GRIDS", pszValue ); } /* -------------------------------------------------------------------- */ /* Handle vertical units. */ /* -------------------------------------------------------------------- */ if( poVERT_CS != nullptr ) { const char *pszUnitName = nullptr; const char *pszUnitConv = nullptr; pszValue = CSLFetchNameValue(papszNV, "vto_meter"); if( pszValue != nullptr && CPLAtofM(pszValue) > 0.0 ) { const double dfValue = CPLAtofM(pszValue); const LinearUnitsStruct* psLinearUnits = GetLinearFromLinearConvOrName( dfValue, pszValue ); if( psLinearUnits != nullptr ) { pszUnitName = psLinearUnits->pszWKTName; pszUnitConv = psLinearUnits->pszValueInMeter; } else { pszUnitName = "unknown"; pszUnitConv = pszValue; } } else if( (pszValue = CSLFetchNameValue(papszNV, "vunits")) != nullptr ) { const LinearUnitsStruct* psLinearUnits = GetLinearFromProjName( pszValue ); if( psLinearUnits != nullptr ) { pszUnitName = psLinearUnits->pszWKTName; pszUnitConv = psLinearUnits->pszValueInMeter; } else { pszUnitName = "unknown"; pszUnitConv = pszValue; } } if( pszUnitName != nullptr ) { OGR_SRSNode *poUnits = new OGR_SRSNode( "UNIT" ); poUnits->AddChild( new OGR_SRSNode( pszUnitName ) ); poUnits->AddChild( new OGR_SRSNode( pszUnitConv ) ); poVERT_CS->AddChild( poUnits ); } } // Add AXIS to VERT_CS node. if( poVERT_CS != nullptr ) { OGR_SRSNode *poAxis = new OGR_SRSNode( "AXIS" ); poAxis->AddChild( new OGR_SRSNode( "Up" ) ); poAxis->AddChild( new OGR_SRSNode( "UP" ) ); poVERT_CS->AddChild( poAxis ); } /* -------------------------------------------------------------------- */ /* do we want to insert a PROJ.4 EXTENSION item? */ /* -------------------------------------------------------------------- */ if( strstr(pszProj4, "wktext") != nullptr || bAddProj4Extension ) SetExtension( GetRoot()->GetValue(), "PROJ4", pszProj4 ); /* -------------------------------------------------------------------- */ /* Preserve authority (for example IGNF) */ /* -------------------------------------------------------------------- */ const char *pszINIT = CSLFetchNameValue(papszNV, "init"); const char *pszColumn = nullptr; if( bSetAuthorityCode && pszINIT != nullptr && (pszColumn = strchr(pszINIT, ':')) != nullptr && GetRoot()->FindChild( "AUTHORITY" ) < 0 ) { CPLString osAuthority; osAuthority.assign(pszINIT, pszColumn - pszINIT); osAuthority.toupper(); OGR_SRSNode* poAuthNode = new OGR_SRSNode( "AUTHORITY" ); poAuthNode->AddChild( new OGR_SRSNode( osAuthority ) ); poAuthNode->AddChild( new OGR_SRSNode( pszColumn + 1 ) ); GetRoot()->AddChild( poAuthNode ); } CSLDestroy( papszNV ); return OGRERR_NONE; } /************************************************************************/ /* OSRExportToProj4() */ /************************************************************************/ /** * \brief Export coordinate system in PROJ.4 format. * * This function is the same as OGRSpatialReference::exportToProj4(). */ OGRErr CPL_STDCALL OSRExportToProj4( OGRSpatialReferenceH hSRS, char ** ppszReturn ) { VALIDATE_POINTER1( hSRS, "OSRExportToProj4", OGRERR_FAILURE ); *ppszReturn = nullptr; return OGRSpatialReference::FromHandle(hSRS)->exportToProj4( ppszReturn ); } /************************************************************************/ /* exportToProj4() */ /************************************************************************/ // TODO(schwehr): Un-macrofy this. #define SAFE_PROJ4_STRCAT(szNewStr) do { \ if( CPLStrlcat(szProj4, szNewStr, sizeof(szProj4)) >= sizeof(szProj4) ) { \ CPLError(CE_Failure, CPLE_AppDefined, \ "String overflow when formatting proj.4 string"); \ *ppszProj4 = CPLStrdup(""); \ return OGRERR_FAILURE; \ } } while( false ); /** * \brief Export coordinate system in PROJ.4 format. * * Converts the loaded coordinate reference system into PROJ.4 format * to the extent possible. The string returned in ppszProj4 should be * deallocated by the caller with CPLFree() when no longer needed. * * LOCAL_CS coordinate systems are not translatable. An empty string * will be returned along with OGRERR_NONE. * * Special processing for Transverse Mercator with GDAL &gt;= 1.10 and PROJ * &gt;= 4.8 : * If the OSR_USE_ETMERC configuration option is set to YES, the PROJ.4 * definition built from the SRS will use the 'etmerc' projection method, * rather than the default 'tmerc'. This will give better accuracy (at the * expense of computational speed) when reprojection occurs near the edges * of the validity area for the projection. * Starting with GDAL &gt;= 2.2, setting OSR_USE_ETMERC to NO will expand to the * 'tmerc' projection method (useful with PROJ &gt;= 4.9.3, where utm uses * etmerc) * * This method is the equivalent of the C function OSRExportToProj4(). * * @param ppszProj4 pointer to which dynamically allocated PROJ.4 definition * will be assigned. * * @return OGRERR_NONE on success or an error code on failure. */ OGRErr OGRSpatialReference::exportToProj4( char ** ppszProj4 ) const { if( GetRoot() == nullptr ) { *ppszProj4 = CPLStrdup(""); CPLError( CE_Failure, CPLE_NotSupported, "No translation for an empty SRS to PROJ.4 format is known."); return OGRERR_UNSUPPORTED_SRS; } /* -------------------------------------------------------------------- */ /* Do we have a PROJ.4 override definition? */ /* -------------------------------------------------------------------- */ const char *pszPredefProj4 = GetExtension( GetRoot()->GetValue(), "PROJ4", nullptr ); if( pszPredefProj4 != nullptr ) { *ppszProj4 = CPLStrdup( pszPredefProj4 ); return OGRERR_NONE; } /* -------------------------------------------------------------------- */ /* Get the prime meridian info. */ /* -------------------------------------------------------------------- */ const OGR_SRSNode *poPRIMEM = GetAttrNode( "PRIMEM" ); double l_dfFromGreenwich = 0.0; if( poPRIMEM != nullptr && poPRIMEM->GetChildCount() >= 2 && CPLAtof(poPRIMEM->GetChild(1)->GetValue()) != 0.0 ) { l_dfFromGreenwich = CPLAtof(poPRIMEM->GetChild(1)->GetValue()); } /* ==================================================================== */ /* Handle the projection definition. */ /* ==================================================================== */ const char *pszProjection = GetAttrValue("PROJECTION"); char szProj4[512] = {}; // TODO(schwehr): Cleanup CPLsnprintf calls to do less pointer arithmetic. if( pszProjection == nullptr && IsGeographic() ) { CPLsnprintf(szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=longlat "); } else if( IsGeocentric() ) { CPLsnprintf(szProj4 + strlen(szProj4) , sizeof(szProj4) - strlen(szProj4), "+proj=geocent "); } else if( pszProjection == nullptr && !IsGeographic() ) { // LOCAL_CS, or incompletely initialized coordinate systems. *ppszProj4 = CPLStrdup(""); return OGRERR_NONE; } else if( EQUAL(pszProjection, SRS_PT_CYLINDRICAL_EQUAL_AREA) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=cea +lon_0=%.16g +lat_ts=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_BONNE) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=bonne +lon_0=%.16g +lat_1=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_CASSINI_SOLDNER) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=cass +lat_0=%.16g +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_NEW_ZEALAND_MAP_GRID) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=nzmg +lat_0=%.16g +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_TRANSVERSE_MERCATOR) || EQUAL(pszProjection, SRS_PT_TRANSVERSE_MERCATOR_MI_21) || EQUAL(pszProjection, SRS_PT_TRANSVERSE_MERCATOR_MI_22) || EQUAL(pszProjection, SRS_PT_TRANSVERSE_MERCATOR_MI_23) || EQUAL(pszProjection, SRS_PT_TRANSVERSE_MERCATOR_MI_24) || EQUAL(pszProjection, SRS_PT_TRANSVERSE_MERCATOR_MI_25) ) { int bNorth = FALSE; const int nZone = GetUTMZone( &bNorth ); const char* pszUseETMERC = CPLGetConfigOption("OSR_USE_ETMERC", nullptr); if( pszUseETMERC && CPLTestBool(pszUseETMERC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=etmerc +lat_0=%.16g +lon_0=%.16g +k=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( pszUseETMERC && !CPLTestBool(pszUseETMERC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=tmerc +lat_0=%.16g +lon_0=%.16g +k=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( nZone != 0 ) { if( bNorth ) CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=utm +zone=%d ", nZone ); else CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=utm +zone=%d +south ", nZone ); } else CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=tmerc +lat_0=%.16g +lon_0=%.16g +k=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_TRANSVERSE_MERCATOR_SOUTH_ORIENTED) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=tmerc +lat_0=%.16g +lon_0=%.16g +k=%.16g " "+x_0=%.16g +y_0=%.16g +axis=wsu ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_MERCATOR_1SP) ) { if( GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0) == 0.0 ) CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=merc +lon_0=%.16g +k=%.16g +x_0=%.16g " "+y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); else if( GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0) == 1.0 ) CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=merc +lon_0=%.16g +lat_ts=%.16g +x_0=%.16g " "+y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); else { CPLError(CE_Failure, CPLE_NotSupported, "Mercator_1SP with scale != 1.0 and " "latitude of origin != 0, not supported by PROJ.4."); *ppszProj4 = CPLStrdup(""); return OGRERR_UNSUPPORTED_SRS; } } else if( EQUAL(pszProjection, SRS_PT_MERCATOR_2SP) ) { if( GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0) == 0.0 ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=merc +lon_0=%.16g +lat_ts=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else { CPLError(CE_Failure, CPLE_NotSupported, "Mercator_2SP with " "latitude of origin != 0, not supported by PROJ.4."); *ppszProj4 = CPLStrdup(""); return OGRERR_UNSUPPORTED_SRS; } } else if( EQUAL(pszProjection, SRS_PT_MERCATOR_AUXILIARY_SPHERE) ) { // This is EPSG:3875 Pseudo Mercator. No point in trying to parse the // rest of the parameters, since we know pretty much everything at this // stage. CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=merc +a=%.16g +b=%.16g +lat_ts=%.16g " "+lon_0=%.16g +x_0=%.16g +y_0=%.16g +k=%.16g +units=m " "+nadgrids=@null +wktext +no_defs", GetSemiMajor(), GetSemiMajor(), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0) ); *ppszProj4 = CPLStrdup( szProj4 ); return OGRERR_NONE; } else if( EQUAL(pszProjection, SRS_PT_OBLIQUE_STEREOGRAPHIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=sterea +lat_0=%.16g +lon_0=%.16g +k=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_STEREOGRAPHIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=stere +lat_0=%.16g +lon_0=%.16g +k=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_POLAR_STEREOGRAPHIC) ) { if( GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0) >= 0.0 ) CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=stere +lat_0=90 +lat_ts=%.16g +lon_0=%.16g " "+k=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 90.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); else CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=stere +lat_0=-90 +lat_ts=%.16g +lon_0=%.16g " "+k=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, -90.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_EQUIRECTANGULAR) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=eqc +lat_ts=%.16g +lat_0=%.16g +lon_0=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_GAUSSSCHREIBERTMERCATOR) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=gstmerc +lat_0=%.16g +lon_0=%.16g " "+k_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, -21.116666667), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 55.53333333309), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 160000.000), GetNormProjParm(SRS_PP_FALSE_NORTHING, 50000.000) ); } else if( EQUAL(pszProjection, SRS_PT_GNOMONIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=gnom +lat_0=%.16g +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ORTHOGRAPHIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=ortho +lat_0=%.16g +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_LAMBERT_AZIMUTHAL_EQUAL_AREA) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=laea +lat_0=%.16g +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_AZIMUTHAL_EQUIDISTANT) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=aeqd +lat_0=%.16g +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_EQUIDISTANT_CONIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=eqdc +lat_0=%.16g +lon_0=%.16g +lat_1=%.16g +lat_2=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_CENTER, 0.0), GetNormProjParm(SRS_PP_LONGITUDE_OF_CENTER, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_2, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_MILLER_CYLINDRICAL) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=mill +lat_0=%.16g +lon_0=%.16g +x_0=%.16g +y_0=%.16g +R_A ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_MOLLWEIDE) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=moll +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ECKERT_I) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=eck1 +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ECKERT_II) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=eck2 +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ECKERT_III) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=eck3 +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ECKERT_IV) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=eck4 +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ECKERT_V) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=eck5 +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ECKERT_VI) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=eck6 +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_POLYCONIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=poly +lat_0=%.16g +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ALBERS_CONIC_EQUAL_AREA) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=aea +lat_1=%.16g +lat_2=%.16g +lat_0=%.16g +lon_0=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_2, 0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_ROBINSON) ) { // Workaround a bug in proj.4 : // https://github.com/OSGeo/proj.4/commit/ // bc7453d1a75aab05bdff2c51ed78c908e3efa3cd const double dfLon0 = GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0); const double dfX0 = GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0); const double dfY0 = GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0); if( CPLIsNan(dfLon0) || CPLIsNan(dfX0) || CPLIsNan(dfY0) ) { return OGRERR_FAILURE; } CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=robin +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", dfLon0, dfX0, dfY0 ); } else if( EQUAL(pszProjection, SRS_PT_VANDERGRINTEN) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=vandg +lon_0=%.16g +x_0=%.16g +y_0=%.16g +R_A ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_SINUSOIDAL) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=sinu +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LONGITUDE_OF_CENTER, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_GALL_STEREOGRAPHIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=gall +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_GOODE_HOMOLOSINE) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=goode +lon_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_IGH) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=igh " ); } else if( EQUAL(pszProjection, SRS_PT_GEOSTATIONARY_SATELLITE) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=geos +lon_0=%.16g +h=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SATELLITE_HEIGHT, 35785831.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_LAMBERT_CONFORMAL_CONIC_2SP) || EQUAL(pszProjection, SRS_PT_LAMBERT_CONFORMAL_CONIC_2SP_BELGIUM) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=lcc +lat_1=%.16g +lat_2=%.16g +lat_0=%.16g +lon_0=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_2, 0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_LAMBERT_CONFORMAL_CONIC_1SP) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=lcc +lat_1=%.16g +lat_0=%.16g +lon_0=%.16g " "+k_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_HOTINE_OBLIQUE_MERCATOR) ) { // Special case for swiss oblique mercator: see bug 423. if( fabs(GetNormProjParm(SRS_PP_AZIMUTH, 0.0) - 90.0) < 0.0001 && fabs(GetNormProjParm(SRS_PP_RECTIFIED_GRID_ANGLE, 0.0) - 90.0) < 0.0001 ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=somerc +lat_0=%.16g +lon_0=%.16g " "+k_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=omerc +lat_0=%.16g +lonc=%.16g +alpha=%.16g " "+k=%.16g +x_0=%.16g +y_0=%.16g +no_uoff ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_AZIMUTH, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); // RSO variant: http://trac.osgeo.org/proj/ticket/62 // Note that gamma is only supported by PROJ 4.8.0 and later. if( GetNormProjParm(SRS_PP_RECTIFIED_GRID_ANGLE, 1000.0) != 1000.0 ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+gamma=%.16g ", GetNormProjParm(SRS_PP_RECTIFIED_GRID_ANGLE, 1000.0)); } } } else if( EQUAL(pszProjection, SRS_PT_HOTINE_OBLIQUE_MERCATOR_AZIMUTH_CENTER) ) { // Special case for swiss oblique mercator: see bug 423. if( fabs(GetNormProjParm(SRS_PP_AZIMUTH, 0.0) - 90.0) < 0.0001 && fabs(GetNormProjParm(SRS_PP_RECTIFIED_GRID_ANGLE, 0.0) - 90.0) < 0.0001 ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=somerc +lat_0=%.16g +lon_0=%.16g " "+k_0=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=omerc +lat_0=%.16g +lonc=%.16g +alpha=%.16g " "+k=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_AZIMUTH, 0.0), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); // RSO variant: http://trac.osgeo.org/proj/ticket/62 // Note that gamma is only supported by PROJ 4.8.0 and later. if( GetNormProjParm(SRS_PP_RECTIFIED_GRID_ANGLE, 1000.0) != 1000.0 ) { CPLsnprintf(szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+gamma=%.16g ", GetNormProjParm(SRS_PP_RECTIFIED_GRID_ANGLE, 1000.0)); } } } else if( EQUAL(pszProjection, SRS_PT_HOTINE_OBLIQUE_MERCATOR_TWO_POINT_NATURAL_ORIGIN) ) { // Not really clear which of Point_1/1st_Point convention is the // "normalized" one, so accept both. CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=omerc +lat_0=%.16g " "+lon_1=%.16g +lat_1=%.16g +lon_2=%.16g +lat_2=%.16g " "+k=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_LONGITUDE_OF_POINT_1, GetNormProjParm(SRS_PP_LONGITUDE_OF_1ST_POINT, 0.0)), GetNormProjParm(SRS_PP_LATITUDE_OF_POINT_1, GetNormProjParm(SRS_PP_LATITUDE_OF_1ST_POINT, 0.0)), GetNormProjParm(SRS_PP_LONGITUDE_OF_POINT_2, GetNormProjParm(SRS_PP_LONGITUDE_OF_2ND_POINT, 0.0)), GetNormProjParm(SRS_PP_LATITUDE_OF_POINT_2, GetNormProjParm(SRS_PP_LATITUDE_OF_2ND_POINT, 0.0)), GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_KROVAK) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=krovak +lat_0=%.16g +lon_0=%.16g +alpha=%.16g " "+k=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_CENTER, 0.0), GetNormProjParm(SRS_PP_LONGITUDE_OF_CENTER, 0.0), GetNormProjParm(SRS_PP_AZIMUTH, 30.28813972222222), // ignored by PROJ GetNormProjParm(SRS_PP_SCALE_FACTOR, 1.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_TWO_POINT_EQUIDISTANT) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=tpeqd +lat_1=%.16g +lon_1=%.16g " "+lat_2=%.16g +lon_2=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_1ST_POINT, 0.0), GetNormProjParm(SRS_PP_LONGITUDE_OF_1ST_POINT, 0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_2ND_POINT, 0.0), GetNormProjParm(SRS_PP_LONGITUDE_OF_2ND_POINT, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_IMW_POLYCONIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=imw_p +lat_1=%.16g +lat_2=%.16g +lon_0=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_1ST_POINT, 0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_2ND_POINT, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_WAGNER_I) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wag1 +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_WAGNER_II) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wag2 +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_WAGNER_III) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wag3 +lat_ts=%.16g +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_WAGNER_IV) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wag4 +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_WAGNER_V) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wag5 +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_WAGNER_VI) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wag6 +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_WAGNER_VII) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wag7 +x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_QSC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=qsc +lat_0=%.16g +lon_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_SCH) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=sch +plat_0=%.16g +plon_0=%.16g +phdg_0=%.16g +h_0=%.16g ", GetNormProjParm(SRS_PP_PEG_POINT_LATITUDE, 0.0), GetNormProjParm(SRS_PP_PEG_POINT_LONGITUDE, 0.0), GetNormProjParm(SRS_PP_PEG_POINT_HEADING, 0.0), GetNormProjParm(SRS_PP_PEG_POINT_HEIGHT, 0.0) ); } // Note: This never really gets used currently. See bug 423. else if( EQUAL(pszProjection, SRS_PT_SWISS_OBLIQUE_CYLINDRICAL) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=somerc +lat_0=%.16g +lon_0=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_AITOFF) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=aitoff +lat_0=%.16g +lon_0=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0) ); } else if( EQUAL(pszProjection, SRS_PT_WINKEL_I) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wink1 +lat_0=%.16g +lon_0=%.16g lat_ts=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 45.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0)); } else if( EQUAL(pszProjection, SRS_PT_WINKEL_II) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wink2 +lat_0=%.16g +lon_0=%.16g +lat_1=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 40.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0)); } else if( EQUAL(pszProjection, SRS_PT_WINKEL_TRIPEL) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=wintri +lat_0=%.16g +lon_0=%.16g +lat_1=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_STANDARD_PARALLEL_1, 40.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0)); } else if( EQUAL(pszProjection, SRS_PT_CRASTER_PARABOLIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=crast +lat_0=%.16g +lon_0=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0)); } else if( EQUAL(pszProjection, SRS_PT_LOXIMUTHAL) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=loxim +lon_0=%.16g +lat_1=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 40.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0)); } else if( EQUAL(pszProjection, SRS_PT_QUARTIC_AUTHALIC) ) { CPLsnprintf( szProj4 + strlen(szProj4), sizeof(szProj4) - strlen(szProj4), "+proj=qua_aut +lat_0=%.16g +lon_0=%.16g " "+x_0=%.16g +y_0=%.16g ", GetNormProjParm(SRS_PP_LATITUDE_OF_ORIGIN, 0.0), GetNormProjParm(SRS_PP_CENTRAL_MERIDIAN, 0.0), GetNormProjParm(SRS_PP_FALSE_EASTING, 0.0), GetNormProjParm(SRS_PP_FALSE_NORTHING, 0.0)); } else { CPLError( CE_Failure, CPLE_NotSupported, "No translation for %s to PROJ.4 format is known.", pszProjection ); *ppszProj4 = CPLStrdup(""); return OGRERR_UNSUPPORTED_SRS; } /* -------------------------------------------------------------------- */ /* Handle earth model. For now we just always emit the user */ /* defined ellipsoid parameters. */ /* -------------------------------------------------------------------- */ const double dfSemiMajor = GetSemiMajor(); const double dfInvFlattening = GetInvFlattening(); const char *pszPROJ4Ellipse = nullptr; const char *pszDatum = GetAttrValue("DATUM"); if( std::abs(dfSemiMajor-6378249.145) < 0.01 && std::abs(dfInvFlattening-293.465) < 0.0001 ) { pszPROJ4Ellipse = "clrk80"; // Clark 1880 } else if( std::abs(dfSemiMajor-6378245.0) < 0.01 && std::abs(dfInvFlattening-298.3) < 0.0001 ) { pszPROJ4Ellipse = "krass"; // Krassovsky } else if( std::abs(dfSemiMajor-6378388.0) < 0.01 && std::abs(dfInvFlattening-297.0) < 0.0001 ) { pszPROJ4Ellipse = "intl"; // International 1924 } else if( std::abs(dfSemiMajor-6378160.0) < 0.01 && std::abs(dfInvFlattening-298.25) < 0.0001 ) { pszPROJ4Ellipse = "aust_SA"; // Australian } else if( std::abs(dfSemiMajor-6377397.155) < 0.01 && std::abs(dfInvFlattening-299.1528128) < 0.0001 ) { pszPROJ4Ellipse = "bessel"; // Bessel 1841 } else if( std::abs(dfSemiMajor-6377483.865) < 0.01 && std::abs(dfInvFlattening-299.1528128) < 0.0001 ) { pszPROJ4Ellipse = "bess_nam"; // Bessel 1841 (Namibia / Schwarzeck) } else if( std::abs(dfSemiMajor-6378160.0) < 0.01 && std::abs(dfInvFlattening-298.247167427) < 0.0001 ) { pszPROJ4Ellipse = "GRS67"; // GRS 1967 } else if( std::abs(dfSemiMajor-6378137) < 0.01 && std::abs(dfInvFlattening-298.257222101) < 0.000001 ) { pszPROJ4Ellipse = "GRS80"; // GRS 1980 } else if( std::abs(dfSemiMajor-6378206.4) < 0.01 && std::abs(dfInvFlattening-294.9786982) < 0.0001 ) { pszPROJ4Ellipse = "clrk66"; // Clarke 1866 } else if( std::abs(dfSemiMajor-6377340.189) < 0.01 && std::abs(dfInvFlattening-299.3249646) < 0.0001 ) { pszPROJ4Ellipse = "mod_airy"; // Modified Airy } else if( std::abs(dfSemiMajor-6377563.396) < 0.01 && std::abs(dfInvFlattening-299.3249646) < 0.0001 ) { pszPROJ4Ellipse = "airy"; // Airy } else if( std::abs(dfSemiMajor-6378200) < 0.01 && std::abs(dfInvFlattening-298.3) < 0.0001 ) { pszPROJ4Ellipse = "helmert"; // Helmert 1906 } else if( std::abs(dfSemiMajor-6378155) < 0.01 && std::abs(dfInvFlattening-298.3) < 0.0001 ) { pszPROJ4Ellipse = "fschr60m"; // Modified Fischer 1960 } else if( std::abs(dfSemiMajor-6377298.556) < 0.01 && std::abs(dfInvFlattening-300.8017) < 0.0001 ) { pszPROJ4Ellipse = "evrstSS"; // Everest (Sabah & Sarawak) } else if( std::abs(dfSemiMajor-6378165.0) < 0.01 && std::abs(dfInvFlattening-298.3) < 0.0001 ) { pszPROJ4Ellipse = "WGS60"; } else if( std::abs(dfSemiMajor-6378145.0) < 0.01 && std::abs(dfInvFlattening-298.25) < 0.0001 ) { pszPROJ4Ellipse = "WGS66"; } else if( std::abs(dfSemiMajor-6378135.0) < 0.01 && std::abs(dfInvFlattening-298.26) < 0.0001 ) { pszPROJ4Ellipse = "WGS72"; } else if( std::abs(dfSemiMajor-6378137.0) < 0.01 && std::abs(dfInvFlattening-298.257223563) < 0.000001 ) { pszPROJ4Ellipse = "WGS84"; } else if( pszDatum != nullptr && EQUAL(pszDatum, "North_American_Datum_1927") ) { // pszPROJ4Ellipse = "clrk66:+datum=nad27"; // NAD 27 pszPROJ4Ellipse = "clrk66"; } else if( pszDatum != nullptr && EQUAL(pszDatum, "North_American_Datum_1983") ) { // pszPROJ4Ellipse = "GRS80:+datum=nad83"; // NAD 83 pszPROJ4Ellipse = "GRS80"; } char szEllipseDef[128] = {}; if( pszPROJ4Ellipse == nullptr ) CPLsnprintf( szEllipseDef, sizeof(szEllipseDef), "+a=%.16g +b=%.16g ", GetSemiMajor(), GetSemiMinor() ); else CPLsnprintf( szEllipseDef, sizeof(szEllipseDef), "+ellps=%s ", pszPROJ4Ellipse ); /* -------------------------------------------------------------------- */ /* Translate the datum. */ /* -------------------------------------------------------------------- */ const char *pszPROJ4Datum = nullptr; const OGR_SRSNode *poTOWGS84 = GetAttrNode( "TOWGS84" ); int nEPSGDatum = -1; int nEPSGGeogCS = -1; const char *pszProj4Grids = GetExtension( "DATUM", "PROJ4_GRIDS" ); const char *pszAuthority = GetAuthorityName( "DATUM" ); if( pszAuthority != nullptr && EQUAL(pszAuthority, "EPSG") ) nEPSGDatum = atoi(GetAuthorityCode( "DATUM" )); const char *pszGeogCSAuthority = GetAuthorityName( "GEOGCS" ); if( pszGeogCSAuthority != nullptr && EQUAL(pszGeogCSAuthority, "EPSG") ) nEPSGGeogCS = atoi(GetAuthorityCode( "GEOGCS" )); if( pszDatum == nullptr ) { // Nothing. } else if( EQUAL(pszDatum, SRS_DN_NAD27) || nEPSGDatum == 6267 ) { pszPROJ4Datum = "NAD27"; } else if( EQUAL(pszDatum, SRS_DN_NAD83) || nEPSGDatum == 6269 ) { pszPROJ4Datum = "NAD83"; } else if( EQUAL(pszDatum, SRS_DN_WGS84) || nEPSGDatum == 6326 ) { pszPROJ4Datum = "WGS84"; } else if( (pszPROJ4Datum = OGRGetProj4Datum(pszDatum, nEPSGDatum)) != nullptr ) { // Nothing. } if( pszProj4Grids != nullptr ) { SAFE_PROJ4_STRCAT( szEllipseDef ); szEllipseDef[0] = '\0'; SAFE_PROJ4_STRCAT( "+nadgrids=" ); SAFE_PROJ4_STRCAT( pszProj4Grids ); SAFE_PROJ4_STRCAT( " " ); pszPROJ4Datum = nullptr; } if( pszPROJ4Datum == nullptr || CPLTestBool(CPLGetConfigOption("OVERRIDE_PROJ_DATUM_WITH_TOWGS84", "YES")) ) { if( poTOWGS84 != nullptr ) { if( poTOWGS84->GetChildCount() >= 3 && (poTOWGS84->GetChildCount() < 7 || (EQUAL(poTOWGS84->GetChild(3)->GetValue(), "") && EQUAL(poTOWGS84->GetChild(4)->GetValue(), "") && EQUAL(poTOWGS84->GetChild(5)->GetValue(), "") && EQUAL(poTOWGS84->GetChild(6)->GetValue(), ""))) ) { SAFE_PROJ4_STRCAT( szEllipseDef ); szEllipseDef[0] = '\0'; SAFE_PROJ4_STRCAT( "+towgs84="); for( int iChild = 0; iChild < 3; iChild++ ) { if( iChild > 0 ) SAFE_PROJ4_STRCAT( "," ); SAFE_PROJ4_STRCAT(poTOWGS84->GetChild(iChild)->GetValue()); } SAFE_PROJ4_STRCAT( " " ); pszPROJ4Datum = nullptr; } else if( poTOWGS84->GetChildCount() >= 7) { SAFE_PROJ4_STRCAT( szEllipseDef ); szEllipseDef[0] = '\0'; SAFE_PROJ4_STRCAT( "+towgs84="); for( int iChild = 0; iChild < 7; iChild++ ) { if( iChild > 0 ) SAFE_PROJ4_STRCAT( "," ); SAFE_PROJ4_STRCAT(poTOWGS84->GetChild(iChild)->GetValue()); } SAFE_PROJ4_STRCAT( " " ); pszPROJ4Datum = nullptr; } } // If we don't know the datum, trying looking up TOWGS84 parameters // based on the EPSG GCS code. else if( nEPSGGeogCS != -1 && pszPROJ4Datum == nullptr ) { std::vector<CPLString> asBursaTransform; if( EPSGGetWGS84Transform( nEPSGGeogCS, asBursaTransform ) ) { char szTOWGS84[256] = {}; CPLsnprintf( szTOWGS84, sizeof(szTOWGS84), "+towgs84=%s,%s,%s,%s,%s,%s,%s", asBursaTransform[0].c_str(), asBursaTransform[1].c_str(), asBursaTransform[2].c_str(), asBursaTransform[3].c_str(), asBursaTransform[4].c_str(), asBursaTransform[5].c_str(), asBursaTransform[6].c_str() ); SAFE_PROJ4_STRCAT( szEllipseDef ); szEllipseDef[0] = '\0'; SAFE_PROJ4_STRCAT( szTOWGS84 ); SAFE_PROJ4_STRCAT( " " ); pszPROJ4Datum = nullptr; } } } if( pszPROJ4Datum != nullptr ) { SAFE_PROJ4_STRCAT( "+datum=" ); SAFE_PROJ4_STRCAT( pszPROJ4Datum ); SAFE_PROJ4_STRCAT( " " ); } else { // The ellipsedef may already have been appended and will now // be empty, otherwise append now. SAFE_PROJ4_STRCAT( szEllipseDef ); // szEllipseDef[0] = '\0'; } /* -------------------------------------------------------------------- */ /* Is there prime meridian info to apply? */ /* -------------------------------------------------------------------- */ if( poPRIMEM != nullptr && poPRIMEM->GetChildCount() >= 2 && CPLAtof(poPRIMEM->GetChild(1)->GetValue()) != 0.0 ) { pszAuthority = GetAuthorityName( "PRIMEM" ); const int nCode = pszAuthority != nullptr && EQUAL(pszAuthority, "EPSG") ? atoi(GetAuthorityCode( "PRIMEM" )) : -1; const OGRProj4PM* psProj4PM = nullptr; if( nCode > 0 ) psProj4PM = OGRGetProj4PMFromCode(nCode); if( psProj4PM == nullptr ) psProj4PM = OGRGetProj4PMFromVal(l_dfFromGreenwich); char szPMValue[128] = {}; if( psProj4PM != nullptr ) { strcpy( szPMValue, psProj4PM->pszProj4PMName ); } else { CPLsnprintf( szPMValue, sizeof(szPMValue), "%.16g", l_dfFromGreenwich ); } SAFE_PROJ4_STRCAT( "+pm=" ); SAFE_PROJ4_STRCAT( szPMValue ); SAFE_PROJ4_STRCAT( " " ); } /* -------------------------------------------------------------------- */ /* Handle linear units. */ /* -------------------------------------------------------------------- */ const char *pszPROJ4Units=nullptr; const char *pszLinearUnits = nullptr; double dfLinearConv = GetLinearUnits( &pszLinearUnits ); if( strstr(szProj4, "longlat") != nullptr ) { pszPROJ4Units = nullptr; } else { const LinearUnitsStruct* psLinearUnits = GetLinearFromLinearConvOrName( dfLinearConv, pszLinearUnits ); if( psLinearUnits != nullptr ) pszPROJ4Units = psLinearUnits->pszProjName; else pszPROJ4Units = nullptr; if( pszPROJ4Units == nullptr ) { char szLinearConv[128] = {}; CPLsnprintf( szLinearConv, sizeof(szLinearConv), "%.16g", dfLinearConv ); SAFE_PROJ4_STRCAT( "+to_meter=" ); SAFE_PROJ4_STRCAT( szLinearConv ); SAFE_PROJ4_STRCAT( " " ); } } if( pszPROJ4Units != nullptr ) { SAFE_PROJ4_STRCAT( "+units="); SAFE_PROJ4_STRCAT( pszPROJ4Units ); SAFE_PROJ4_STRCAT( " " ); } /* -------------------------------------------------------------------- */ /* If we have vertical datum grids, attach them to the proj.4 string. */ /* -------------------------------------------------------------------- */ const char *pszProj4Geoids = GetExtension( "VERT_DATUM", "PROJ4_GRIDS" ); if( pszProj4Geoids != nullptr ) { SAFE_PROJ4_STRCAT( "+geoidgrids=" ); SAFE_PROJ4_STRCAT( pszProj4Geoids ); SAFE_PROJ4_STRCAT( " " ); } /* -------------------------------------------------------------------- */ /* Handle vertical units, but only if we have them. */ /* -------------------------------------------------------------------- */ const OGR_SRSNode *poVERT_CS = GetRoot()->GetNode( "VERT_CS" ); const OGR_SRSNode *poVUNITS = nullptr; if( poVERT_CS != nullptr ) poVUNITS = poVERT_CS->GetNode( "UNIT" ); if( poVUNITS != nullptr && poVUNITS->GetChildCount() >= 2 ) { pszPROJ4Units = nullptr; dfLinearConv = CPLAtof( poVUNITS->GetChild(1)->GetValue() ); const LinearUnitsStruct* psLinearUnits = GetLinearFromLinearConvOrName(dfLinearConv, poVUNITS->GetChild(0)->GetValue()); if( psLinearUnits != nullptr ) pszPROJ4Units = psLinearUnits->pszProjName; else pszPROJ4Units = nullptr; if( pszPROJ4Units == nullptr ) { char szLinearConv[128] = {}; CPLsnprintf( szLinearConv, sizeof(szLinearConv), "%.16g", dfLinearConv ); SAFE_PROJ4_STRCAT( "+vto_meter=" ); SAFE_PROJ4_STRCAT( szLinearConv ); SAFE_PROJ4_STRCAT( " " ); } else { SAFE_PROJ4_STRCAT( "+vunits="); SAFE_PROJ4_STRCAT( pszPROJ4Units ); SAFE_PROJ4_STRCAT( " " ); } } /* -------------------------------------------------------------------- */ /* Add the no_defs flag to ensure that no values from */ /* proj_def.dat are implicitly used with our definitions. */ /* -------------------------------------------------------------------- */ SAFE_PROJ4_STRCAT( "+no_defs " ); *ppszProj4 = CPLStrdup( szProj4 ); return OGRERR_NONE; }