/****************************************************************************** * $Id: ogr_srs_proj4.cpp,v 1.63 2006/04/21 03:33:23 fwarmerdam Exp $ * * Project: OpenGIS Simple Features Reference Implementation * Purpose: OGRSpatialReference interface to PROJ.4. * Author: Frank Warmerdam * ****************************************************************************** * Copyright (c) 1999, Les Technologies SoftMap Inc. * * 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. ****************************************************************************** * * $Log: ogr_srs_proj4.cpp,v $ * Revision 1.63 2006/04/21 03:33:23 fwarmerdam * Per bug 811 we should not actually be adjusting the central meridian or * other longitudes. They are relative to the geogcs prime meridian in * WKT and PROJ.4 format. Ugg. * * Revision 1.62 2006/04/21 02:44:35 fwarmerdam * Fixed bogota meridian sign (west, not east). * * Revision 1.61 2006/03/21 15:12:22 fwarmerdam * Fixed support for +R parameter. * * Revision 1.60 2005/12/13 15:57:28 dron * Export scale coefficient when exporting to stereographic projection. * * Revision 1.59 2005/12/01 04:59:46 fwarmerdam * added two point equidistant support * * Revision 1.58 2005/11/10 04:12:40 fwarmerdam * The last change related to: * http://bugzilla.remotesensing.org/show_bug.cgi?id=980 * * Revision 1.57 2005/11/10 04:12:14 fwarmerdam * Changed the translation for Oblique Stereographic from proj=stere to * proj=sterea. * * Revision 1.56 2005/09/29 20:22:26 dmorissette * Improved support for MapInfo modified TM projections in exportToProj4() * (Anthony D - MITAB bug 1155) * * Revision 1.55 2005/09/05 20:42:52 fwarmerdam * ensure that non-geographic, non-projected srses return empty string * * Revision 1.54 2005/04/06 00:02:05 fwarmerdam * various osr and oct functions now stdcall * * Revision 1.53 2005/02/28 15:01:56 fwarmerdam * Applied patch to add +towgs84 parameter if one is well known for a given * EPSG GEOGCS, and no datum shifting info is available in the source defn. * See http://bugzilla.remotesensing.org/show_bug.cgi?id=784 * * Revision 1.52 2005/02/09 20:25:39 fwarmerdam * Changed default height for GEOS projection. * * Revision 1.51 2005/01/05 21:02:33 fwarmerdam * added Goode Homolosine * * Revision 1.50 2004/11/11 18:28:45 fwarmerdam * added Bonne projection support * * Revision 1.49 2004/09/10 21:04:51 fwarmerdam * Various fixes for swiss oblique mercator (somerc) support. * http://bugzilla.remotesensing.org/show_bug.cgi?id=423 * * Revision 1.48 2004/09/10 20:29:16 fwarmerdam * Re bug 423: fixed up SRS_PT_SWISS_OBLIQUE_CYLINDRICAL to use somerc. * * Revision 1.47 2004/05/12 14:27:47 warmerda * added translation for nzgd49 and potsdam using epsg codes * * Revision 1.46 2004/05/04 18:49:43 warmerda * Fixed handling of prime meridian. In PROJ.4 strings the central meridian * or other longitude parameters are relative to the prime meridian. Apply * appropriate transformations going and coming from WKT where that is not * the case. * * Revision 1.45 2004/02/05 17:07:59 dron * Support for HOM projection, specified by two points on centerline. * * Revision 1.44 2003/06/27 19:02:04 warmerda * added OSRProj4Tokenize(), fixes plus related tokenizing bug in bug 355 * * Revision 1.43 2003/06/27 17:54:12 warmerda * allow k_0 in place of k parsing proj4, bug 355 * * Revision 1.42 2003/06/21 23:25:03 warmerda * added +towgs84 support in importFromProj4() * * Revision 1.41 2003/05/28 18:17:38 warmerda * treat meter or m as meter * * Revision 1.40 2003/05/20 18:09:36 warmerda * fixed so that importFromProj4() will transform linear parameter units * * Revision 1.39 2003/02/13 19:27:28 warmerda * always append no_defs to avoid getting hosted by defaults file * * Revision 1.38 2002/12/15 23:42:59 warmerda * added initial support for normalizing proj params * * Revision 1.37 2002/12/14 19:50:21 warmerda * implement support for NZMG (New Zealand Map Grid) * * Revision 1.36 2002/12/10 04:05:04 warmerda * fixed some translation problems with prime meridians * * Revision 1.35 2002/12/09 23:03:45 warmerda * Fixed translation of paris prime meridian. * * Revision 1.34 2002/12/09 18:55:07 warmerda * moved DMS stuff to gdal/port * * Revision 1.33 2002/12/09 16:12:14 warmerda * added +pm= support * * Revision 1.32 2002/12/03 15:38:10 warmerda * dont write linear units to geograpic srs but avoid freaking out * * Revision 1.31 2002/12/01 20:19:21 warmerda * use %.16g to preserve precision as per bug 184 * * Revision 1.30 2002/11/29 21:23:38 warmerda * implement LCC1SP PROJ.4 to WKT translation * * Revision 1.29 2002/11/29 20:56:37 warmerda * added LCC1SP conversion to PROJ.4, missing reverse * * Revision 1.28 2002/11/25 03:28:16 warmerda * added/improved documentation * * Revision 1.27 2002/07/09 14:49:07 warmerda * Improved translation from/to polar stereographic as per * http://bugzilla.remotesensing.org/show_bug.cgi?id=172 * * Revision 1.26 2002/06/11 18:02:03 warmerda * add PROJ.4 normalization and EPSG support * * Revision 1.25 2002/04/18 14:22:45 warmerda * made OGRSpatialReference and co 'const correct' * * Revision 1.24 2002/03/08 20:17:52 warmerda * assume WGS84 if ellipse info missing * * Revision 1.23 2002/01/18 04:48:48 warmerda * clean tabs, and newlines from input to importProj4 * * Revision 1.22 2001/12/12 20:18:06 warmerda * added support for units in proj.4 to WKT * * Revision 1.21 2001/11/07 22:08:51 warmerda * dont put + in to_meter value * * Revision 1.20 2001/11/02 22:21:15 warmerda * fixed memory leak * * Revision 1.19 2001/10/11 19:23:30 warmerda * fixed datum names * * Revision 1.18 2001/07/18 05:03:05 warmerda * added CPL_CVSID * * Revision 1.17 2001/03/02 04:37:43 danmo * Return empty string for LOCAL_CS in exportToProj4(). * * Revision 1.16 2001/01/22 14:00:28 warmerda * added untested support for Swiss Oblique Cylindrical * * Revision 1.15 2001/01/19 21:10:46 warmerda * replaced tabs * * Revision 1.14 2000/12/05 23:10:34 warmerda * Added cassini support * * Revision 1.13 2000/12/05 17:46:16 warmerda * Use +R_A for VanDerGrinten and miller * * Revision 1.12 2000/09/25 20:22:21 warmerda * ensure pszPROJ4Units is initialized * * Revision 1.11 2000/09/25 15:46:56 warmerda * added sphere * * Revision 1.10 2000/08/28 20:13:23 warmerda * added importFromProj4 * * Revision 1.9 2000/07/11 01:02:06 warmerda * added ExportToProj4() * * Revision 1.8 2000/07/09 20:49:21 warmerda * added +datum support * * Revision 1.7 2000/03/06 02:23:54 warmerda * don't use +datum syntax * * Revision 1.6 1999/12/22 15:39:43 warmerda * fix to differentiate WGS variants * * Revision 1.5 1999/12/13 16:29:59 warmerda * Added improved units, and ellipse support. * * Revision 1.4 1999/12/08 16:34:05 warmerda * added five or six more projections * * Revision 1.1 1999/07/29 17:29:15 warmerda * New * */ #include "ogr_spatialref.h" #include "ogr_p.h" #include "cpl_conv.h" extern int EPSGGetWGS84Transform( int nGeogCS, double *padfTransform ); CPL_CVSID("$Id: ogr_srs_proj4.cpp,v 1.63 2006/04/21 03:33:23 fwarmerdam Exp $"); /* -------------------------------------------------------------------- */ /* The following list comes from osrs/proj/src/pj_ellps.c */ /* ... please update from time to time. */ /* -------------------------------------------------------------------- */ static const char *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)", 0, 0, 0, 0, }; /************************************************************************/ /* 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. */ /************************************************************************/ char **OSRProj4Tokenize( const char *pszFull ) { char *pszStart = NULL; char *pszFullWrk; char **papszTokens = NULL; int i; if( pszFull == NULL ) return NULL; pszFullWrk = CPLStrdup( pszFull ); for( i=0; pszFullWrk[i] != '\0'; i++ ) { switch( pszFullWrk[i] ) { case '+': if( i == 0 || pszFullWrk[i-1] == '\0' ) { if( pszStart != NULL ) papszTokens = CSLAddString( papszTokens, pszStart ); pszStart = pszFullWrk + i + 1; } break; case ' ': case '\t': case '\n': pszFullWrk[i] = '\0'; break; default: break; } } if( pszStart != NULL && strlen(pszStart) > 0 ) papszTokens = CSLAddString( papszTokens, pszStart ); CPLFree( pszFullWrk ); return papszTokens; } /************************************************************************/ /* OSRImportFromProj4() */ /************************************************************************/ OGRErr OSRImportFromProj4( OGRSpatialReferenceH hSRS, const char *pszProj4 ) { return ((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; pszValue = CSLFetchNameValue( papszNV, pszField ); // special hack to use k_0 if available. if( pszValue == NULL && EQUAL(pszField,"k") ) pszValue = CSLFetchNameValue( papszNV, "k_0" ); if( pszValue == NULL ) return dfDefaultValue; else return CPLDMSToDec(pszValue); } /************************************************************************/ /* importFromProj4() */ /************************************************************************/ /** * 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 equivelents, 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" * * This method is the equivelent 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 ) { char **papszNV = NULL; char **papszTokens; int i; char *pszCleanCopy; /* -------------------------------------------------------------------- */ /* Strip any newlines or other "funny" stuff that might occur */ /* if this string just came from reading a file. */ /* -------------------------------------------------------------------- */ pszCleanCopy = CPLStrdup( pszProj4 ); for( i = 0; pszCleanCopy[i] != '\0'; i++ ) { if( pszCleanCopy[i] == 10 || pszCleanCopy[i] == 13 || pszCleanCopy[i] == 9 ) pszCleanCopy[i] = ' '; } /* -------------------------------------------------------------------- */ /* Try to normalize the definition. This should expand +init= */ /* clauses and so forth. */ /* -------------------------------------------------------------------- */ char *pszNormalized; pszNormalized = OCTProj4Normalize( pszCleanCopy ); CPLFree( pszCleanCopy ); /* -------------------------------------------------------------------- */ /* If we have an EPSG based init string, try to process it */ /* directly to get the fullest possible EPSG definition. */ /* -------------------------------------------------------------------- */ if( strstr(pszNormalized,"init=epsg:") != NULL ) { OGRErr eErr; const char *pszNumber = strstr(pszNormalized,"init=epsg:") + 10; eErr = importFromEPSG( atoi(pszNumber) ); if( eErr == OGRERR_NONE ) { CPLFree( pszNormalized ); return eErr; } } /* -------------------------------------------------------------------- */ /* Parse the PROJ.4 string into a cpl_string.h style name/value */ /* list. */ /* -------------------------------------------------------------------- */ papszTokens = OSRProj4Tokenize( pszNormalized ); CPLFree( pszNormalized ); for( i = 0; papszTokens != NULL && papszTokens[i] != NULL; i++ ) { char *pszEqual = strstr(papszTokens[i],"="); if( pszEqual == NULL ) 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 dfFromGreenwich = 0.0; int nPMCode = -1; if( pszPM != NULL ) { if( EQUAL(pszPM,"lisbon") ) { dfFromGreenwich = CPLDMSToDec( "9d07'54.862\"W" ); nPMCode = 8902; } else if( EQUAL(pszPM,"paris") ) { dfFromGreenwich = CPLDMSToDec( "2d20'14.025\"E" ); nPMCode = 8903; } else if( EQUAL(pszPM,"bogota") ) { dfFromGreenwich = CPLDMSToDec( "74d04'51.3\"W" ); nPMCode = 8904; } else if( EQUAL(pszPM,"madrid") ) { dfFromGreenwich = CPLDMSToDec( "3d41'16.48\"W" ); nPMCode = 8905; } else if( EQUAL(pszPM,"rome") ) { dfFromGreenwich = CPLDMSToDec( "12d27'8.4\"E" ); nPMCode = 8906; } else if( EQUAL(pszPM,"bern") ) { dfFromGreenwich = CPLDMSToDec( "7d26'22.5\"E" ); nPMCode = 8907; } else if( EQUAL(pszPM,"jakarta") ) { dfFromGreenwich = CPLDMSToDec( "106d48'27.79\"E" ); nPMCode = 8908; } else if( EQUAL(pszPM,"ferro") ) { dfFromGreenwich = CPLDMSToDec( "17d40'W" ); nPMCode = 8909; } else if( EQUAL(pszPM,"brussels") ) { dfFromGreenwich = CPLDMSToDec( "4d22'4.71\"E" ); nPMCode = 8910; } else if( EQUAL(pszPM,"stockholm") ) { dfFromGreenwich = CPLDMSToDec( "18d3'29.8\"E" ); nPMCode = 8911; } else if( EQUAL(pszPM,"athens") ) { dfFromGreenwich = CPLDMSToDec( "23d42'58.815\"E" ); nPMCode = 8912; } else if( EQUAL(pszPM,"oslo") ) { dfFromGreenwich = CPLDMSToDec( "10d43'22.5\"E" ); nPMCode = 8913; } else { dfFromGreenwich = CPLDMSToDec( pszPM ); pszPM = "unnamed"; } } else pszPM = "Greenwich"; /* -------------------------------------------------------------------- */ /* Operate on the basis of the projection name. */ /* -------------------------------------------------------------------- */ const char *pszProj = CSLFetchNameValue(papszNV,"proj"); if( pszProj == NULL ) { 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,"bonne") ) { SetBonne( OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ) + dfFromGreenwich, 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 ) + dfFromGreenwich, 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 ) + dfFromGreenwich, 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 ) + dfFromGreenwich, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj,"tmerc") ) { SetTM( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ) + dfFromGreenwich, 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") ) { SetUTM( (int) OSR_GDV( papszNV, "zone", 0.0 ), (int) OSR_GDV( papszNV, "south", 1.0 ) ); } else if( EQUAL(pszProj,"merc") ) { SetMercator( OSR_GDV( papszNV, "lat_ts", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ) + dfFromGreenwich, 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") && 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 ) + dfFromGreenwich, 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") && 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 ) + dfFromGreenwich, OSR_GDV( papszNV, "k", 1.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUALN(pszProj,"stere",5) /* mostly sterea */ && CSLFetchNameValue(papszNV,"k") != NULL ) { SetOS( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 ) + dfFromGreenwich, 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 ) + dfFromGreenwich, 1.0, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj,"eqc") ) { SetEquirectangular( OSR_GDV( papszNV, "lat_ts", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 )+dfFromGreenwich, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj,"gnom") ) { SetGnomonic( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj,"eck4") ) { SetEckertIV( OSR_GDV( papszNV, "lon_0", 0.0 )+dfFromGreenwich, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj,"eck6") ) { SetEckertVI( OSR_GDV( papszNV, "lon_0", 0.0 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, 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 )+dfFromGreenwich, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj,"geos") ) { SetGEOS( OSR_GDV( papszNV, "lon_0", 0.0 )+dfFromGreenwich, OSR_GDV( papszNV, "h", 35785831.0 ), OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else if( EQUAL(pszProj,"lcc") ) { if( OSR_GDV(papszNV, "lat_0", 0.0 ) == OSR_GDV(papszNV, "lat_1", 0.0 ) ) { /* 1SP form */ SetLCC1SP( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 )+dfFromGreenwich, 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 )+dfFromGreenwich, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } } else if( EQUAL(pszProj,"omerc") ) { SetHOM( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lonc", 0.0 )+dfFromGreenwich, OSR_GDV( papszNV, "alpha", 0.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,"somerc") ) { SetHOM( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 )+dfFromGreenwich, 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,"tpeqd") ) { SetTPED( OSR_GDV( papszNV, "lat_0", 0.0 ), OSR_GDV( papszNV, "lon_0", 0.0 )+dfFromGreenwich, OSR_GDV( papszNV, "lat_1", 0.0 ), OSR_GDV( papszNV, "lon_1", 0.0 )+dfFromGreenwich, OSR_GDV( papszNV, "x_0", 0.0 ), OSR_GDV( papszNV, "y_0", 0.0 ) ); } else { CPLDebug( "OGR_PROJ4", "Unsupported projection: %s", pszProj ); CSLDestroy( papszNV ); return OGRERR_CORRUPT_DATA; } /* -------------------------------------------------------------------- */ /* Try to translate the datum. */ /* -------------------------------------------------------------------- */ const char *pszValue; int bFullyDefined = FALSE; pszValue = CSLFetchNameValue(papszNV, "datum"); if( pszValue == NULL ) { /* do nothing */ } else if( (EQUAL(pszValue,"NAD27") || EQUAL(pszValue,"NAD83") || EQUAL(pszValue,"WGS84") || EQUAL(pszValue,"WGS72")) && dfFromGreenwich == 0.0 ) { SetWellKnownGeogCS( pszValue ); bFullyDefined = TRUE; } else if( EQUAL(pszValue,"potsdam") ) { OGRSpatialReference oGCS; oGCS.importFromEPSG( 4314 ); CopyGeogCSFrom( &oGCS ); bFullyDefined = TRUE; } else if( EQUAL(pszValue,"nzgd49") ) { OGRSpatialReference oGCS; oGCS.importFromEPSG( 4272 ); CopyGeogCSFrom( &oGCS ); bFullyDefined = TRUE; } else { /* we don't recognise the datum, and ignore it */ } /* -------------------------------------------------------------------- */ /* Set the ellipsoid information. */ /* -------------------------------------------------------------------- */ double dfSemiMajor, dfInvFlattening, dfSemiMinor; pszValue = CSLFetchNameValue(papszNV, "ellps"); if( pszValue != NULL && !bFullyDefined ) { for( i = 0; ogr_pj_ellps[i] != NULL; i += 4 ) { if( !EQUAL(ogr_pj_ellps[i],pszValue) ) continue; CPLAssert( EQUALN(ogr_pj_ellps[i+1],"a=",2) ); dfSemiMajor = atof(ogr_pj_ellps[i+1]+2); if( EQUALN(ogr_pj_ellps[i+2],"rf=",3) ) dfInvFlattening = atof(ogr_pj_ellps[i+2]+3); else { CPLAssert( EQUALN(ogr_pj_ellps[i+2],"b=",2) ); dfSemiMinor = atof(ogr_pj_ellps[i+2]+2); if( ABS(dfSemiMajor/dfSemiMinor) - 1.0 < 0.0000000000001 ) dfInvFlattening = 0.0; else dfInvFlattening = -1.0 / (dfSemiMinor/dfSemiMajor - 1.0); } SetGeogCS( ogr_pj_ellps[i+3], "unknown", ogr_pj_ellps[i], dfSemiMajor, dfInvFlattening, pszPM, 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 ) { CPLDebug( "OGR_PROJ4", "Can't find ellipse definition in:\n%s", pszProj4 ); CSLDestroy( papszNV ); return OGRERR_UNSUPPORTED_SRS; } if( dfInvFlattening == -1.0 ) { if( ABS(dfSemiMajor/dfSemiMinor) - 1.0 < 0.0000000000001 ) dfInvFlattening = 0.0; else dfInvFlattening = -1.0 / (dfSemiMinor/dfSemiMajor - 1.0); } SetGeogCS( "unnamed ellipse", "unknown", "unnamed", dfSemiMajor, dfInvFlattening, pszPM, dfFromGreenwich ); bFullyDefined = TRUE; } /* -------------------------------------------------------------------- */ /* Handle TOWGS84 conversion. */ /* -------------------------------------------------------------------- */ pszValue = CSLFetchNameValue(papszNV, "towgs84"); if(pszValue!=NULL) { char **papszToWGS84 = CSLTokenizeStringComplex( pszValue, ",", FALSE, TRUE ); if( CSLCount(papszToWGS84) >= 7 ) SetTOWGS84( atof(papszToWGS84[0]), atof(papszToWGS84[1]), atof(papszToWGS84[2]), atof(papszToWGS84[3]), atof(papszToWGS84[4]), atof(papszToWGS84[5]), atof(papszToWGS84[6]) ); else if( CSLCount(papszToWGS84) >= 3 ) SetTOWGS84( atof(papszToWGS84[0]), atof(papszToWGS84[1]), atof(papszToWGS84[2]) ); else CPLError( CE_Warning, CPLE_AppDefined, "Seemingly corrupt +towgs84 option (%s), ignoring.", pszValue ); CSLDestroy(papszToWGS84); } /* -------------------------------------------------------------------- */ /* Linear units translation */ /* -------------------------------------------------------------------- */ if( IsProjected() || IsLocal() ) { pszValue = CSLFetchNameValue(papszNV, "to_meter"); if( pszValue != NULL && atof(pszValue) > 0.0 ) { SetLinearUnits( "unknown", atof(pszValue) ); } else if( (pszValue = CSLFetchNameValue(papszNV, "units")) != NULL ) { if( EQUAL(pszValue,"meter" ) || EQUAL(pszValue,"m") ) SetLinearUnits( SRS_UL_METER, 1.0 ); else if( EQUAL(pszValue,"us-ft" ) ) SetLinearUnits( SRS_UL_US_FOOT, atof(SRS_UL_US_FOOT_CONV) ); else if( EQUAL(pszValue,"ft" ) ) SetLinearUnits( SRS_UL_FOOT, atof(SRS_UL_FOOT_CONV) ); else if( EQUAL(pszValue,"yd" ) ) SetLinearUnits( pszValue, 0.9144 ); else if( EQUAL(pszValue,"us-yd" ) ) SetLinearUnits( pszValue, 0.914401828803658 ); else // This case is untranslatable. Should add all proj.4 unts SetLinearUnits( pszValue, 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" ); int i; for( 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)); } } CSLDestroy( papszNV ); return OGRERR_NONE; } /************************************************************************/ /* OSRExportToProj4() */ /************************************************************************/ OGRErr CPL_STDCALL OSRExportToProj4( OGRSpatialReferenceH hSRS, char ** ppszReturn ) { *ppszReturn = NULL; return ((OGRSpatialReference *) hSRS)->exportToProj4( ppszReturn ); } /************************************************************************/ /* exportToProj4() */ /************************************************************************/ /** * 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. * * This method is the equivelent 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 { char szProj4[512]; const char *pszProjection = GetAttrValue("PROJECTION"); szProj4[0] = '\0'; /* -------------------------------------------------------------------- */ /* Get the prime meridian info. */ /* -------------------------------------------------------------------- */ const OGR_SRSNode *poPRIMEM = GetAttrNode( "PRIMEM" ); double dfFromGreenwich = 0.0; if( poPRIMEM != NULL && poPRIMEM->GetChildCount() >= 2 && atof(poPRIMEM->GetChild(1)->GetValue()) != 0.0 ) { dfFromGreenwich = atof(poPRIMEM->GetChild(1)->GetValue()); } /* ==================================================================== */ /* Handle the projection definition. */ /* ==================================================================== */ if( pszProjection == NULL && IsGeographic() ) { sprintf( szProj4+strlen(szProj4), "+proj=longlat " ); } else if( pszProjection == NULL && !IsGeographic() ) { // LOCAL_CS, or incompletely initialized coordinate systems. *ppszProj4 = CPLStrdup(""); return OGRERR_NONE; } else if( EQUAL(pszProjection,SRS_PT_CYLINDRICAL_EQUAL_AREA) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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; int nZone = GetUTMZone( &bNorth ); if( nZone != 0 ) { if( bNorth ) sprintf( szProj4+strlen(szProj4), "+proj=utm +zone=%d ", nZone ); else sprintf( szProj4+strlen(szProj4),"+proj=utm +zone=%d +south ", nZone ); } else sprintf( szProj4+strlen(szProj4), "+proj=tmerc +lat_0=%.16g +lon_0=%.16g +k=%f +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_MERCATOR_1SP) ) { sprintf( szProj4+strlen(szProj4), "+proj=merc +lat_ts=%.16g +lon_0=%.16g +k=%f +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_OBLIQUE_STEREOGRAPHIC) ) { sprintf( szProj4+strlen(szProj4), "+proj=sterea +lat_0=%.16g +lon_0=%.16g +k=%f +x_0=%.16g +y_0=%.16g ", // "+proj=stere +lat_0=%.16g +lon_0=%.16g +k=%f +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) ) { sprintf( szProj4+strlen(szProj4), "+proj=stere +lat_0=%.16g +lon_0=%.16g +k=%f +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 ) sprintf( 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 sprintf( 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) ) { sprintf( szProj4+strlen(szProj4), "+proj=eqc +lat_ts=%.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_GNOMONIC) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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_IV) ) { sprintf( 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_VI) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( szProj4+strlen(szProj4), "+proj=robin +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_VANDERGRINTEN) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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_GEOSTATIONARY_SATELLITE) ) { sprintf( 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) ) { sprintf( 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) ) { sprintf( 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) ) { /* not clear how ProjParm[3] - angle from rectified to skewed grid - should be applied ... see the +not_rot flag for PROJ.4. Just ignoring for now. */ /* 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 ) { sprintf( 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 { sprintf( 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) ); } } else if( EQUAL(pszProjection, SRS_PT_HOTINE_OBLIQUE_MERCATOR_TWO_POINT_NATURAL_ORIGIN) ) { sprintf( 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_LATITUDE_OF_POINT_1,0.0), GetNormProjParm(SRS_PP_LONGITUDE_OF_POINT_1,0.0), GetNormProjParm(SRS_PP_LATITUDE_OF_POINT_2,0.0), GetNormProjParm(SRS_PP_LONGITUDE_OF_POINT_2,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) ) { sprintf( 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,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_TWO_POINT_EQUIDISTANT) ) { sprintf( szProj4+strlen(szProj4), "+proj=tpeqd +lat_0=%.16g +lon_0=%.16g " "+lat_1=%.16g +lon_1=%.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) ); } /* Note: This never really gets used currently. See bug 423 */ else if( EQUAL(pszProjection,SRS_PT_SWISS_OBLIQUE_CYLINDRICAL) ) { sprintf( 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) ); } /* -------------------------------------------------------------------- */ /* Handle earth model. For now we just always emit the user */ /* defined ellipsoid parameters. */ /* -------------------------------------------------------------------- */ double dfSemiMajor = GetSemiMajor(); double dfInvFlattening = GetInvFlattening(); const char *pszPROJ4Ellipse = NULL; const char *pszDatum = GetAttrValue("DATUM"); if( ABS(dfSemiMajor-6378249.145) < 0.01 && ABS(dfInvFlattening-293.465) < 0.0001 ) { pszPROJ4Ellipse = "clrk80"; /* Clark 1880 */ } else if( ABS(dfSemiMajor-6378245.0) < 0.01 && ABS(dfInvFlattening-298.3) < 0.0001 ) { pszPROJ4Ellipse = "krass"; /* Krassovsky */ } else if( ABS(dfSemiMajor-6378388.0) < 0.01 && ABS(dfInvFlattening-297.0) < 0.0001 ) { pszPROJ4Ellipse = "intl"; /* International 1924 */ } else if( ABS(dfSemiMajor-6378160.0) < 0.01 && ABS(dfInvFlattening-298.25) < 0.0001 ) { pszPROJ4Ellipse = "aust_SA"; /* Australian */ } else if( ABS(dfSemiMajor-6377397.155) < 0.01 && ABS(dfInvFlattening-299.1528128) < 0.0001 ) { pszPROJ4Ellipse = "bessel"; /* Bessel 1841 */ } else if( ABS(dfSemiMajor-6377483.865) < 0.01 && ABS(dfInvFlattening-299.1528128) < 0.0001 ) { pszPROJ4Ellipse = "bess_nam"; /* Bessel 1841 (Namibia / Schwarzeck)*/ } else if( ABS(dfSemiMajor-6378160.0) < 0.01 && ABS(dfInvFlattening-298.247167427) < 0.0001 ) { pszPROJ4Ellipse = "GRS67"; /* GRS 1967 */ } else if( ABS(dfSemiMajor-6378137) < 0.01 && ABS(dfInvFlattening-298.257222101) < 0.000001 ) { pszPROJ4Ellipse = "GRS80"; /* GRS 1980 */ } else if( ABS(dfSemiMajor-6378206.4) < 0.01 && ABS(dfInvFlattening-294.9786982) < 0.0001 ) { pszPROJ4Ellipse = "clrk66"; /* Clarke 1866 */ } else if( ABS(dfSemiMajor-6378206.4) < 0.01 && ABS(dfInvFlattening-294.9786982) < 0.0001 ) { pszPROJ4Ellipse = "mod_airy"; /* Modified Airy */ } else if( ABS(dfSemiMajor-6377563.396) < 0.01 && ABS(dfInvFlattening-299.3249646) < 0.0001 ) { pszPROJ4Ellipse = "airy"; /* Modified Airy */ } else if( ABS(dfSemiMajor-6378200) < 0.01 && ABS(dfInvFlattening-298.3) < 0.0001 ) { pszPROJ4Ellipse = "helmert"; /* Helmert 1906 */ } else if( ABS(dfSemiMajor-6378155) < 0.01 && ABS(dfInvFlattening-298.3) < 0.0001 ) { pszPROJ4Ellipse = "fschr60m"; /* Modified Fischer 1960 */ } else if( ABS(dfSemiMajor-6377298.556) < 0.01 && ABS(dfInvFlattening-300.8017) < 0.0001 ) { pszPROJ4Ellipse = "evrstSS"; /* Everest (Sabah & Sarawak) */ } else if( ABS(dfSemiMajor-6378165.0) < 0.01 && ABS(dfInvFlattening-298.3) < 0.0001 ) { pszPROJ4Ellipse = "WGS60"; } else if( ABS(dfSemiMajor-6378145.0) < 0.01 && ABS(dfInvFlattening-298.25) < 0.0001 ) { pszPROJ4Ellipse = "WGS66"; } else if( ABS(dfSemiMajor-6378135.0) < 0.01 && ABS(dfInvFlattening-298.26) < 0.0001 ) { pszPROJ4Ellipse = "WGS72"; } else if( ABS(dfSemiMajor-6378137.0) < 0.01 && ABS(dfInvFlattening-298.257223563) < 0.000001 ) { pszPROJ4Ellipse = "WGS84"; } else if( EQUAL(pszDatum,"North_American_Datum_1927") ) { // pszPROJ4Ellipse = "clrk66:+datum=nad27"; /* NAD 27 */ pszPROJ4Ellipse = "clrk66"; } else if( EQUAL(pszDatum,"North_American_Datum_1983") ) { // pszPROJ4Ellipse = "GRS80:+datum=nad83"; /* NAD 83 */ pszPROJ4Ellipse = "GRS80"; } if( pszPROJ4Ellipse == NULL ) sprintf( szProj4+strlen(szProj4), "+a=%.16g +b=%.16g ", GetSemiMajor(), GetSemiMinor() ); else sprintf( szProj4+strlen(szProj4), "+ellps=%s ", pszPROJ4Ellipse ); /* -------------------------------------------------------------------- */ /* Translate the datum. */ /* -------------------------------------------------------------------- */ const char *pszPROJ4Datum = NULL; const OGR_SRSNode *poTOWGS84 = GetAttrNode( "TOWGS84" ); char szTOWGS84[256]; int nEPSGDatum = -1; const char *pszAuthority; int nEPSGGeogCS = -1; const char *pszGeogCSAuthority; pszAuthority = GetAuthorityName( "DATUM" ); if( pszAuthority != NULL && EQUAL(pszAuthority,"EPSG") ) nEPSGDatum = atoi(GetAuthorityCode( "DATUM" )); pszGeogCSAuthority = GetAuthorityName( "GEOGCS" ); if( pszGeogCSAuthority != NULL && EQUAL(pszGeogCSAuthority,"EPSG") ) nEPSGGeogCS = atoi(GetAuthorityCode( "GEOGCS" )); if( pszDatum == NULL ) /* nothing */; else if( EQUAL(pszDatum,SRS_DN_NAD27) || nEPSGDatum == 6267 ) pszPROJ4Datum = "+datum=NAD27"; else if( EQUAL(pszDatum,SRS_DN_NAD83) || nEPSGDatum == 6269 ) pszPROJ4Datum = "+datum=NAD83"; else if( EQUAL(pszDatum,SRS_DN_WGS84) || nEPSGDatum == 6326 ) pszPROJ4Datum = "+datum=WGS84"; else if( nEPSGDatum == 6314 ) pszPROJ4Datum = "+datum=potsdam"; else if( nEPSGDatum == 6272 ) pszPROJ4Datum = "+datum=nzgd49"; else if( poTOWGS84 != NULL ) { if( poTOWGS84->GetChildCount() > 2 && (poTOWGS84->GetChildCount() < 6 || EQUAL(poTOWGS84->GetChild(3)->GetValue(),"") && EQUAL(poTOWGS84->GetChild(4)->GetValue(),"") && EQUAL(poTOWGS84->GetChild(5)->GetValue(),"") && EQUAL(poTOWGS84->GetChild(6)->GetValue(),"")) ) { sprintf( szTOWGS84, "+towgs84=%s,%s,%s", poTOWGS84->GetChild(0)->GetValue(), poTOWGS84->GetChild(1)->GetValue(), poTOWGS84->GetChild(2)->GetValue() ); pszPROJ4Datum = szTOWGS84; } else if( poTOWGS84->GetChildCount() > 6 ) { sprintf( szTOWGS84, "+towgs84=%s,%s,%s,%s,%s,%s,%s", poTOWGS84->GetChild(0)->GetValue(), poTOWGS84->GetChild(1)->GetValue(), poTOWGS84->GetChild(2)->GetValue(), poTOWGS84->GetChild(3)->GetValue(), poTOWGS84->GetChild(4)->GetValue(), poTOWGS84->GetChild(5)->GetValue(), poTOWGS84->GetChild(6)->GetValue() ); pszPROJ4Datum = szTOWGS84; } } else if( nEPSGGeogCS != -1 ) { double padfTransform[7]; if( EPSGGetWGS84Transform( nEPSGGeogCS, padfTransform ) ) { sprintf( szTOWGS84, "+towgs84=%f,%f,%f,%f,%f,%f,%f", padfTransform[0], padfTransform[1], padfTransform[2], padfTransform[3], padfTransform[4], padfTransform[5], padfTransform[6] ); pszPROJ4Datum = szTOWGS84; } } if( pszPROJ4Datum != NULL ) { strcat( szProj4, pszPROJ4Datum ); strcat( szProj4, " " ); } /* -------------------------------------------------------------------- */ /* Is there prime meridian info to apply? */ /* -------------------------------------------------------------------- */ if( poPRIMEM != NULL && poPRIMEM->GetChildCount() >= 2 && atof(poPRIMEM->GetChild(1)->GetValue()) != 0.0 ) { const char *pszAuthority = GetAuthorityName( "PRIMEM" ); char szPMValue[128]; int nCode = -1; if( pszAuthority != NULL && EQUAL(pszAuthority,"EPSG") ) nCode = atoi(GetAuthorityCode( "PRIMEM" )); switch( nCode ) { case 8902: strcpy( szPMValue, "lisbon" ); break; case 8903: strcpy( szPMValue, "paris" ); break; case 8904: strcpy( szPMValue, "bogota" ); break; case 8905: strcpy( szPMValue, "madrid" ); break; case 8906: strcpy( szPMValue, "rome" ); break; case 8907: strcpy( szPMValue, "bern" ); break; case 8908: strcpy( szPMValue, "jakarta" ); break; case 8909: strcpy( szPMValue, "ferro" ); break; case 8910: strcpy( szPMValue, "brussels" ); break; case 8911: strcpy( szPMValue, "stockholm" ); break; case 8912: strcpy( szPMValue, "athens" ); break; case 8913: strcpy( szPMValue, "oslo" ); break; default: sprintf( szPMValue, "%.16g", dfFromGreenwich ); } sprintf( szProj4+strlen(szProj4), "+pm=%s ", szPMValue ); } /* -------------------------------------------------------------------- */ /* Handle linear units. */ /* -------------------------------------------------------------------- */ const char *pszPROJ4Units=NULL; char *pszLinearUnits = NULL; double dfLinearConv; dfLinearConv = GetLinearUnits( &pszLinearUnits ); if( strstr(szProj4,"longlat") != NULL ) pszPROJ4Units = NULL; else if( dfLinearConv == 1.0 ) pszPROJ4Units = "m"; else if( dfLinearConv == 1000.0 ) pszPROJ4Units = "km"; else if( dfLinearConv == 0.0254 ) pszPROJ4Units = "in"; else if( EQUAL(pszLinearUnits,SRS_UL_FOOT) ) pszPROJ4Units = "ft"; else if( EQUAL(pszLinearUnits,"IYARD") || dfLinearConv == 0.9144 ) pszPROJ4Units = "yd"; else if( dfLinearConv == 0.001 ) pszPROJ4Units = "mm"; else if( dfLinearConv == 0.01 ) pszPROJ4Units = "cm"; else if( EQUAL(pszLinearUnits,SRS_UL_US_FOOT) ) pszPROJ4Units = "us-ft"; else if( EQUAL(pszLinearUnits,SRS_UL_NAUTICAL_MILE) ) pszPROJ4Units = "kmi"; else if( EQUAL(pszLinearUnits,"Mile") || EQUAL(pszLinearUnits,"IMILE") ) pszPROJ4Units = "mi"; else { sprintf( szProj4+strlen(szProj4), "+to_meter=%.16g ", dfLinearConv ); } if( pszPROJ4Units != NULL ) sprintf( szProj4+strlen(szProj4), "+units=%s ", pszPROJ4Units ); /* -------------------------------------------------------------------- */ /* Add the no_defs flag to ensure that no values from */ /* proj_def.dat are implicitly used with our definitions. */ /* -------------------------------------------------------------------- */ sprintf( szProj4+strlen(szProj4), "+no_defs " ); *ppszProj4 = CPLStrdup( szProj4 ); return OGRERR_NONE; }