/****************************************************************************** * $Id: ogr_srs_esri.cpp,v 1.47 2006/04/24 23:21:07 fwarmerdam Exp $ * * Project: OpenGIS Simple Features Reference Implementation * Purpose: OGRSpatialReference translation to/from ESRI .prj definitions. * Author: Frank Warmerdam, warmerdam@pobox.com * ****************************************************************************** * Copyright (c) 2000, Frank Warmerdam * * 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_esri.cpp,v $ * Revision 1.47 2006/04/24 23:21:07 fwarmerdam * fixup some common ellipsoid names (from pgao) * * Revision 1.46 2006/03/31 17:44:20 fwarmerdam * header updates * * Revision 1.45 2006/01/16 15:19:03 fwarmerdam * Fixed datum.csv error message. * * Revision 1.44 2005/12/20 21:49:57 fwarmerdam * trim log * * Revision 1.43 2005/12/20 21:48:38 fwarmerdam * Fixed last fix which seriously hosted morphToESRI(). * * Revision 1.42 2005/12/20 02:09:49 fwarmerdam * ESRI requires linear units. Ensure morphToESRI() produces them. * * Revision 1.41 2005/10/27 23:12:16 fwarmerdam * morphToESRI() now maps a PROJCS of Unknown to the PROJECTION value * * Revision 1.40 2005/08/03 16:11:18 fwarmerdam * Added additional morphToESRI() rules to: * o Change NAD27/NAD83/WGS84 GEOGCSes to GCS_ names used by ESRI. * o Change the word Unnamed to Unknown in names. * o Set exact ESRI names for UTM PROJCSes (for WGS84, NAD27 and NAD83). * * Revision 1.39 2005/05/04 14:29:20 fwarmerdam * convert Standard_Parallel_1 in ESRI Mercator to Latitude_Of_Origin. * * Revision 1.38 2005/02/18 21:52:44 fwarmerdam * Fixed support for .prj files with blank lines between the parameters. * * Revision 1.37 2005/01/13 16:32:27 fwarmerdam * added support for fipszone for stateplane * * Revision 1.36 2005/01/13 15:18:05 fwarmerdam * use SetLinearUnitsAndUpdateParameters() * * Revision 1.35 2004/09/23 16:20:40 fwarmerdam * added clode to cleanup datum mapping table: bug 613 * * Revision 1.34 2004/09/10 21:03:55 fwarmerdam * Lots of changes to map Hotine_Oblique_Mercator_Azimuth_Center in ESRI format * to Hotine_Oblique_Mercator with a rectified_grid_angle of 90 (and back). * This is a special case for the swiss oblique mercator. * See http://bugzilla.remotesensing.org/show_bug.cgi?id=423 * * Revision 1.33 2004/07/29 19:12:05 warmerda * support multiline ESRI .prj files in WKT format * * Revision 1.32 2004/05/04 13:11:49 warmerda * Added support for KRASOVSKY spheroid * * Revision 1.31 2004/04/24 15:45:04 warmerda * Added GRS80 spheroid support. * * Revision 1.30 2004/04/19 19:32:25 warmerda * added INTERNATIONAL1909 * * Revision 1.29 2004/02/25 21:14:35 warmerda * added morph of spheroid names to ESRI * * Revision 1.28 2003/11/03 21:37:34 warmerda * fixed southern hemisphere UTM support * * Revision 1.27 2003/08/14 14:36:24 warmerda * some ESRI .prj files use Central_Parallel instead of latitude_of_origin * * Revision 1.26 2003/08/03 01:30:34 warmerda * Added mapping for EUR datum as per: * http://bugzilla.remotesensing.org/show_bug.cgi?id=371 * * Revision 1.25 2003/06/23 14:49:17 warmerda * added InitDatumMappingTable, and use of gdal_datum.csv file * * Revision 1.24 2003/06/17 14:39:00 warmerda * Translate Equirectangular to/from Equidistant_Cylindrical * * Revision 1.23 2003/05/21 02:59:41 warmerda * morphToEsri() now sets specific units constants for Foot_US and Degree. * morphToEsri() converts albers parameters to match expected ESRI values. * * Revision 1.22 2003/05/08 21:52:55 warmerda * fixed C morphToESRI func, added ESRI unit remapping * * Revision 1.21 2003/02/25 04:53:16 warmerda * Added support for the LAMBERT projection. Added support for defining a * GEOGCS from the SPHEROID if there is no known DATUM. Fixed bug with * GREATBRITIAN_GRID. * * Revision 1.20 2003/02/14 22:15:04 warmerda * expand tabs * * Revision 1.19 2003/01/24 20:15:34 warmerda * added polar stereographic support */ #include "ogr_spatialref.h" #include "ogr_p.h" #include "cpl_csv.h" CPL_CVSID("$Id: ogr_srs_esri.cpp,v 1.47 2006/04/24 23:21:07 fwarmerdam Exp $"); static char *apszProjMapping[] = { "Albers", SRS_PT_ALBERS_CONIC_EQUAL_AREA, "Cassini", SRS_PT_CASSINI_SOLDNER, "Hotine_Oblique_Mercator_Azimuth_Natural_Origin", SRS_PT_HOTINE_OBLIQUE_MERCATOR, "Hotine_Oblique_Mercator_Azimuth_Center", SRS_PT_HOTINE_OBLIQUE_MERCATOR, "Lambert_Conformal_Conic", SRS_PT_LAMBERT_CONFORMAL_CONIC_2SP, "Lambert_Conformal_Conic", SRS_PT_LAMBERT_CONFORMAL_CONIC_1SP, "Van_der_Grinten_I", SRS_PT_VANDERGRINTEN, SRS_PT_TRANSVERSE_MERCATOR, SRS_PT_TRANSVERSE_MERCATOR, "Gauss_Kruger", SRS_PT_TRANSVERSE_MERCATOR, "Mercator", SRS_PT_MERCATOR_1SP, "Equidistant_Cylindrical", SRS_PT_EQUIRECTANGULAR, NULL, NULL }; static char *apszAlbersMapping[] = { SRS_PP_CENTRAL_MERIDIAN, SRS_PP_LONGITUDE_OF_CENTER, SRS_PP_LATITUDE_OF_ORIGIN, SRS_PP_LATITUDE_OF_CENTER, "Central_Parallel", SRS_PP_LATITUDE_OF_CENTER, NULL, NULL }; static char *apszMercatorMapping[] = { SRS_PP_STANDARD_PARALLEL_1, SRS_PP_LATITUDE_OF_ORIGIN, NULL, NULL }; static char **papszDatumMapping = NULL; static char *apszDefaultDatumMapping[] = { "6267", "North_American_1927", SRS_DN_NAD27, "6269", "North_American_1983", SRS_DN_NAD83, NULL, NULL, NULL }; static char *apszUnitMapping[] = { "Meter", "meter", "Meter", "metre", "Foot", "foot", "Foot", "feet", "Foot_US", SRS_UL_US_FOOT, "Degree", "degree", "Degree", "degrees", "Degree", SRS_UA_DEGREE, "Radian", SRS_UA_RADIAN, NULL, NULL }; /* -------------------------------------------------------------------- */ /* Table relating USGS and ESRI state plane zones. */ /* -------------------------------------------------------------------- */ static int anUsgsEsriZones[] = { 101, 3101, 102, 3126, 201, 3151, 202, 3176, 203, 3201, 301, 3226, 302, 3251, 401, 3276, 402, 3301, 403, 3326, 404, 3351, 405, 3376, 406, 3401, 407, 3426, 501, 3451, 502, 3476, 503, 3501, 600, 3526, 700, 3551, 901, 3601, 902, 3626, 903, 3576, 1001, 3651, 1002, 3676, 1101, 3701, 1102, 3726, 1103, 3751, 1201, 3776, 1202, 3801, 1301, 3826, 1302, 3851, 1401, 3876, 1402, 3901, 1501, 3926, 1502, 3951, 1601, 3976, 1602, 4001, 1701, 4026, 1702, 4051, 1703, 6426, 1801, 4076, 1802, 4101, 1900, 4126, 2001, 4151, 2002, 4176, 2101, 4201, 2102, 4226, 2103, 4251, 2111, 6351, 2112, 6376, 2113, 6401, 2201, 4276, 2202, 4301, 2203, 4326, 2301, 4351, 2302, 4376, 2401, 4401, 2402, 4426, 2403, 4451, 2500, 0, 2501, 4476, 2502, 4501, 2503, 4526, 2600, 0, 2601, 4551, 2602, 4576, 2701, 4601, 2702, 4626, 2703, 4651, 2800, 4676, 2900, 4701, 3001, 4726, 3002, 4751, 3003, 4776, 3101, 4801, 3102, 4826, 3103, 4851, 3104, 4876, 3200, 4901, 3301, 4926, 3302, 4951, 3401, 4976, 3402, 5001, 3501, 5026, 3502, 5051, 3601, 5076, 3602, 5101, 3701, 5126, 3702, 5151, 3800, 5176, 3900, 0, 3901, 5201, 3902, 5226, 4001, 5251, 4002, 5276, 4100, 5301, 4201, 5326, 4202, 5351, 4203, 5376, 4204, 5401, 4205, 5426, 4301, 5451, 4302, 5476, 4303, 5501, 4400, 5526, 4501, 5551, 4502, 5576, 4601, 5601, 4602, 5626, 4701, 5651, 4702, 5676, 4801, 5701, 4802, 5726, 4803, 5751, 4901, 5776, 4902, 5801, 4903, 5826, 4904, 5851, 5001, 6101, 5002, 6126, 5003, 6151, 5004, 6176, 5005, 6201, 5006, 6226, 5007, 6251, 5008, 6276, 5009, 6301, 5010, 6326, 5101, 5876, 5102, 5901, 5103, 5926, 5104, 5951, 5105, 5976, 5201, 6001, 5200, 6026, 5200, 6076, 5201, 6051, 5202, 6051, 5300, 0, 5400, 0 }; void OGREPSGDatumNameMassage( char ** ppszDatum ); /************************************************************************/ /* RemapSpheroidName() */ /* */ /* Convert Spheroid name to ESRI style name */ /************************************************************************/ static const char* RemapSpheroidName(const char* pszName) { if (strcmp(pszName, "WGS 84") == 0) return "WGS 1984"; if (strcmp(pszName, "WGS 72") == 0) return "WGS 1972"; return pszName; } /************************************************************************/ /* ESRIToUSGSZone() */ /* */ /* Convert ESRI style state plane zones to USGS style state */ /* plane zones. */ /************************************************************************/ static int ESRIToUSGSZone( int nESRIZone ) { int nPairs = sizeof(anUsgsEsriZones) / (2*sizeof(int)); int i; for( i = 0; i < nPairs; i++ ) { if( anUsgsEsriZones[i*2+1] == nESRIZone ) return anUsgsEsriZones[i*2]; } return 0; } /************************************************************************/ /* MorphNameToESRI() */ /* */ /* Make name ESRI compatible. Convert spaces and special */ /* characters to underscores and then strip down. */ /************************************************************************/ static void MorphNameToESRI( char ** ppszName ) { int i, j; char *pszName = *ppszName; /* -------------------------------------------------------------------- */ /* Translate non-alphanumeric values to underscores. */ /* -------------------------------------------------------------------- */ for( i = 0; pszName[i] != '\0'; i++ ) { if( !(pszName[i] >= 'A' && pszName[i] <= 'Z') && !(pszName[i] >= 'a' && pszName[i] <= 'z') && !(pszName[i] >= '0' && pszName[i] <= '9') ) { pszName[i] = '_'; } } /* -------------------------------------------------------------------- */ /* Remove repeated and trailing underscores. */ /* -------------------------------------------------------------------- */ for( i = 1, j = 0; pszName[i] != '\0'; i++ ) { if( pszName[j] == '_' && pszName[i] == '_' ) continue; pszName[++j] = pszName[i]; } if( pszName[j] == '_' ) pszName[j] = '\0'; else pszName[j+1] = '\0'; } /************************************************************************/ /* CleanESRIDatumMappingTable() */ /************************************************************************/ CPL_C_START void CleanupESRIDatumMappingTable() { if( papszDatumMapping == NULL ) return; if( papszDatumMapping != apszDefaultDatumMapping ) { CSLDestroy( papszDatumMapping ); papszDatumMapping = NULL; } } CPL_C_END /************************************************************************/ /* InitDatumMappingTable() */ /************************************************************************/ static void InitDatumMappingTable() { if( papszDatumMapping != NULL ) return; /* -------------------------------------------------------------------- */ /* Try to open the datum.csv file. */ /* -------------------------------------------------------------------- */ const char *pszFilename = CSVFilename("gdal_datum.csv"); FILE * fp = VSIFOpen( pszFilename, "rb" ); /* -------------------------------------------------------------------- */ /* Use simple default set if we can't find the file. */ /* -------------------------------------------------------------------- */ if( fp == NULL ) { papszDatumMapping = apszDefaultDatumMapping; return; } /* -------------------------------------------------------------------- */ /* Figure out what fields we are interested in. */ /* -------------------------------------------------------------------- */ char **papszFieldNames = CSVReadParseLine( fp ); int nDatumCodeField = CSLFindString( papszFieldNames, "DATUM_CODE" ); int nEPSGNameField = CSLFindString( papszFieldNames, "DATUM_NAME" ); int nESRINameField = CSLFindString( papszFieldNames, "ESRI_DATUM_NAME" ); CSLDestroy( papszFieldNames ); if( nDatumCodeField == -1 || nEPSGNameField == -1 || nESRINameField == -1 ) { CPLError( CE_Failure, CPLE_AppDefined, "Failed to find required field in gdal_datum.csv in InitDatumMappingTable(), using default table setup." ); papszDatumMapping = apszDefaultDatumMapping; return; } /* -------------------------------------------------------------------- */ /* Read each line, adding a detail line for each. */ /* -------------------------------------------------------------------- */ int nMappingCount = 0; const int nMaxDatumMappings = 1000; char **papszFields; papszDatumMapping = (char **)CPLCalloc(sizeof(char*),nMaxDatumMappings*3); for( papszFields = CSVReadParseLine( fp ); papszFields != NULL; papszFields = CSVReadParseLine( fp ) ) { int nFieldCount = CSLCount(papszFields); CPLAssert( nMappingCount+1 < nMaxDatumMappings ); if( MAX(nEPSGNameField,MAX(nDatumCodeField,nESRINameField)) < nFieldCount && nMaxDatumMappings > nMappingCount+1 ) { papszDatumMapping[nMappingCount*3+0] = CPLStrdup( papszFields[nDatumCodeField] ); papszDatumMapping[nMappingCount*3+1] = CPLStrdup( papszFields[nESRINameField] ); papszDatumMapping[nMappingCount*3+2] = CPLStrdup( papszFields[nEPSGNameField] ); OGREPSGDatumNameMassage( &(papszDatumMapping[nMappingCount*3+2]) ); nMappingCount++; } CSLDestroy( papszFields ); } VSIFClose( fp ); papszDatumMapping[nMappingCount*3+0] = NULL; papszDatumMapping[nMappingCount*3+1] = NULL; papszDatumMapping[nMappingCount*3+2] = NULL; } /************************************************************************/ /* OSRImportFromESRI() */ /************************************************************************/ OGRErr OSRImportFromESRI( OGRSpatialReferenceH hSRS, char **papszPrj ) { return ((OGRSpatialReference *) hSRS)->importFromESRI( papszPrj ); } /************************************************************************/ /* 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 importFromESRI(). */ /************************************************************************/ static double OSR_GDV( char **papszNV, const char * pszField, double dfDefaultValue ) { int iLine; if( papszNV == NULL || papszNV[0] == NULL ) return dfDefaultValue; if( EQUALN(pszField,"PARAM_",6) ) { int nOffset; for( iLine = 0; papszNV[iLine] != NULL && !EQUALN(papszNV[iLine],"Paramet",7); iLine++ ) {} for( nOffset=atoi(pszField+6); papszNV[iLine] != NULL && nOffset > 0; iLine++ ) { if( strlen(papszNV[iLine]) > 0 ) nOffset--; } while( papszNV[iLine] != NULL && strlen(papszNV[iLine]) == 0 ) iLine++; if( papszNV[iLine] != NULL ) { char **papszTokens, *pszLine = papszNV[iLine]; double dfValue; int i; // Trim comments. for( i=0; pszLine[i] != '\0'; i++ ) { if( pszLine[i] == '/' && pszLine[i+1] == '*' ) pszLine[i] = '\0'; } papszTokens = CSLTokenizeString(papszNV[iLine]); if( CSLCount(papszTokens) == 3 ) { dfValue = ABS(atof(papszTokens[0])) + atof(papszTokens[1]) / 60.0 + atof(papszTokens[2]) / 3600.0; if( atof(papszTokens[0]) < 0.0 ) dfValue *= -1; } else if( CSLCount(papszTokens) > 0 ) dfValue = atof(papszTokens[0]); else dfValue = dfDefaultValue; CSLDestroy( papszTokens ); return dfValue; } else return dfDefaultValue; } else { for( iLine = 0; papszNV[iLine] != NULL && !EQUALN(papszNV[iLine],pszField,strlen(pszField)); iLine++ ) {} if( papszNV[iLine] == NULL ) return dfDefaultValue; else return atof( papszNV[iLine] + strlen(pszField) ); } } /************************************************************************/ /* OSR_GDS() */ /************************************************************************/ static const char*OSR_GDS( char **papszNV, const char * pszField, const char *pszDefaultValue ) { int iLine; if( papszNV == NULL || papszNV[0] == NULL ) return pszDefaultValue; for( iLine = 0; papszNV[iLine] != NULL && !EQUALN(papszNV[iLine],pszField,strlen(pszField)); iLine++ ) {} if( papszNV[iLine] == NULL ) return pszDefaultValue; else { static char szResult[80]; char **papszTokens; papszTokens = CSLTokenizeString(papszNV[iLine]); if( CSLCount(papszTokens) > 1 ) strncpy( szResult, papszTokens[1], sizeof(szResult)); else strncpy( szResult, pszDefaultValue, sizeof(szResult)); CSLDestroy( papszTokens ); return szResult; } } /************************************************************************/ /* importFromESRI() */ /************************************************************************/ /** * Import coordinate system from ESRI .prj format(s). * * This function will read the text loaded from an ESRI .prj file, and * translate it into an OGRSpatialReference definition. This should support * many (but by no means all) old style (Arc/Info 7.x) .prj files, as well * as the newer pseudo-OGC WKT .prj files. Note that new style .prj files * are in OGC WKT format, but require some manipulation to correct datum * names, and units on some projection parameters. This is addressed within * importFromESRI() by an automatical call to morphFromESRI(). * * Currently only GEOGRAPHIC, UTM, STATEPLANE, GREATBRITIAN_GRID, ALBERS, * EQUIDISTANT_CONIC, and TRANSVERSE (mercator) projections are supported * from old style files. * * At this time there is no equivelent exportToESRI() method. Writing old * style .prj files is not supported by OGRSpatialReference. However the * morphToESRI() and exportToWkt() methods can be used to generate output * suitable to write to new style (Arc 8) .prj files. * * This function is the equilvelent of the C function OSRImportFromESRI(). * * @param papszPrj NULL terminated list of strings containing the definition. * * @return OGRERR_NONE on success or an error code in case of failure. */ OGRErr OGRSpatialReference::importFromESRI( char **papszPrj ) { if( papszPrj == NULL || papszPrj[0] == NULL ) return OGRERR_CORRUPT_DATA; /* -------------------------------------------------------------------- */ /* ArcGIS and related products now use a varient of Well Known */ /* Text. Try to recognise this and ingest it. WKT is usually */ /* all on one line, but we will accept multi-line formats and */ /* concatenate. */ /* -------------------------------------------------------------------- */ if( EQUALN(papszPrj[0],"GEOGCS",6) || EQUALN(papszPrj[0],"PROJCS",6) || EQUALN(papszPrj[0],"LOCAL_CS",8) ) { char *pszWKT, *pszWKT2; OGRErr eErr; int i; pszWKT = CPLStrdup(papszPrj[0]); for( i = 1; papszPrj[i] != NULL; i++ ) { pszWKT = (char *) CPLRealloc(pszWKT,strlen(pszWKT)+strlen(papszPrj[i])+1); strcat( pszWKT, papszPrj[i] ); } pszWKT2 = pszWKT; eErr = importFromWkt( &pszWKT2 ); CPLFree( pszWKT ); if( eErr == OGRERR_NONE ) eErr = morphFromESRI(); return eErr; } /* -------------------------------------------------------------------- */ /* Operate on the basis of the projection name. */ /* -------------------------------------------------------------------- */ const char *pszProj = OSR_GDS( papszPrj, "Projection", NULL ); if( pszProj == NULL ) { CPLDebug( "OGR_ESRI", "Can't find Projection\n" ); return OGRERR_CORRUPT_DATA; } else if( EQUAL(pszProj,"GEOGRAPHIC") ) { } else if( EQUAL(pszProj,"utm") ) { if( (int) OSR_GDV( papszPrj, "zone", 0.0 ) != 0 ) { double dfYShift = OSR_GDV( papszPrj, "Yshift", 0.0 ); SetUTM( (int) OSR_GDV( papszPrj, "zone", 0.0 ), dfYShift == 0.0 ); } else { double dfCentralMeridian, dfRefLat; int nZone; dfCentralMeridian = OSR_GDV( papszPrj, "PARAM_1", 0.0 ); dfRefLat = OSR_GDV( papszPrj, "PARAM_2", 0.0 ); nZone = (int) ((dfCentralMeridian+183) / 6.0 + 0.0000001); SetUTM( nZone, dfRefLat >= 0.0 ); } } else if( EQUAL(pszProj,"STATEPLANE") ) { int nZone = (int) OSR_GDV( papszPrj, "zone", 0.0 ); if( nZone != 0 ) nZone = ESRIToUSGSZone( nZone ); else nZone = (int) OSR_GDV( papszPrj, "fipszone", 0.0 ); if( nZone != 0 ) { if( EQUAL(OSR_GDS( papszPrj, "Datum", "NAD83"),"NAD27") ) SetStatePlane( nZone, FALSE ); else SetStatePlane( nZone, TRUE ); } } else if( EQUAL(pszProj,"GREATBRITIAN_GRID") || EQUAL(pszProj,"GREATBRITAIN_GRID") ) { const char *pszWkt = "PROJCS[\"OSGB 1936 / British National Grid\",GEOGCS[\"OSGB 1936\",DATUM[\"OSGB_1936\",SPHEROID[\"Airy 1830\",6377563.396,299.3249646]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",49],PARAMETER[\"central_meridian\",-2],PARAMETER[\"scale_factor\",0.999601272],PARAMETER[\"false_easting\",400000],PARAMETER[\"false_northing\",-100000],UNIT[\"metre\",1]]"; importFromWkt( (char **) &pszWkt ); } else if( EQUAL(pszProj,"ALBERS") ) { SetACEA( OSR_GDV( papszPrj, "PARAM_1", 0.0 ), OSR_GDV( papszPrj, "PARAM_2", 0.0 ), OSR_GDV( papszPrj, "PARAM_4", 0.0 ), OSR_GDV( papszPrj, "PARAM_3", 0.0 ), OSR_GDV( papszPrj, "PARAM_5", 0.0 ), OSR_GDV( papszPrj, "PARAM_6", 0.0 ) ); } else if( EQUAL(pszProj,"LAMBERT") ) { SetLCC( OSR_GDV( papszPrj, "PARAM_1", 0.0 ), OSR_GDV( papszPrj, "PARAM_2", 0.0 ), OSR_GDV( papszPrj, "PARAM_4", 0.0 ), OSR_GDV( papszPrj, "PARAM_3", 0.0 ), OSR_GDV( papszPrj, "PARAM_5", 0.0 ), OSR_GDV( papszPrj, "PARAM_6", 0.0 ) ); } else if( EQUAL(pszProj,"EQUIDISTANT_CONIC") ) { int nStdPCount = (int) OSR_GDV( papszPrj, "PARAM_1", 0.0 ); if( nStdPCount == 1 ) { SetEC( OSR_GDV( papszPrj, "PARAM_2", 0.0 ), OSR_GDV( papszPrj, "PARAM_2", 0.0 ), OSR_GDV( papszPrj, "PARAM_4", 0.0 ), OSR_GDV( papszPrj, "PARAM_3", 0.0 ), OSR_GDV( papszPrj, "PARAM_5", 0.0 ), OSR_GDV( papszPrj, "PARAM_6", 0.0 ) ); } else { SetEC( OSR_GDV( papszPrj, "PARAM_2", 0.0 ), OSR_GDV( papszPrj, "PARAM_3", 0.0 ), OSR_GDV( papszPrj, "PARAM_5", 0.0 ), OSR_GDV( papszPrj, "PARAM_4", 0.0 ), OSR_GDV( papszPrj, "PARAM_5", 0.0 ), OSR_GDV( papszPrj, "PARAM_7", 0.0 ) ); } } else if( EQUAL(pszProj,"TRANSVERSE") ) { SetTM( OSR_GDV( papszPrj, "PARAM_2", 0.0 ), OSR_GDV( papszPrj, "PARAM_3", 0.0 ), OSR_GDV( papszPrj, "PARAM_1", 0.0 ), OSR_GDV( papszPrj, "PARAM_4", 0.0 ), OSR_GDV( papszPrj, "PARAM_5", 0.0 ) ); } else if( EQUAL(pszProj,"POLAR") ) { SetPS( OSR_GDV( papszPrj, "PARAM_2", 0.0 ), OSR_GDV( papszPrj, "PARAM_1", 0.0 ), 1.0, OSR_GDV( papszPrj, "PARAM_3", 0.0 ), OSR_GDV( papszPrj, "PARAM_4", 0.0 ) ); } else { CPLDebug( "OGR_ESRI", "Unsupported projection: %s", pszProj ); SetLocalCS( pszProj ); } /* -------------------------------------------------------------------- */ /* Try to translate the datum/spheroid. */ /* -------------------------------------------------------------------- */ if( !IsLocal() && GetAttrNode( "GEOGCS" ) == NULL ) { const char *pszDatum; pszDatum = OSR_GDS( papszPrj, "Datum", ""); if( EQUAL(pszDatum,"NAD27") || EQUAL(pszDatum,"NAD83") || EQUAL(pszDatum,"WGS84") || EQUAL(pszDatum,"WGS72") ) { SetWellKnownGeogCS( pszDatum ); } else if( EQUAL( pszDatum, "EUR" ) ) { SetWellKnownGeogCS( "EPSG:4230" ); } else { const char *pszSpheroid; pszSpheroid = OSR_GDS( papszPrj, "Spheroid", ""); if( EQUAL(pszSpheroid,"INT1909") || EQUAL(pszSpheroid,"INTERNATIONAL1909") ) { OGRSpatialReference oGCS; oGCS.importFromEPSG( 4022 ); CopyGeogCSFrom( &oGCS ); } else if( EQUAL(pszSpheroid,"AIRY") ) { OGRSpatialReference oGCS; oGCS.importFromEPSG( 4001 ); CopyGeogCSFrom( &oGCS ); } else if( EQUAL(pszSpheroid,"CLARKE1866") ) { OGRSpatialReference oGCS; oGCS.importFromEPSG( 4008 ); CopyGeogCSFrom( &oGCS ); } else if( EQUAL(pszSpheroid,"GRS80") ) { OGRSpatialReference oGCS; oGCS.importFromEPSG( 4019 ); CopyGeogCSFrom( &oGCS ); } else if( EQUAL(pszSpheroid,"KRASOVSKY") || EQUAL(pszSpheroid,"KRASSOVSKY") ) { OGRSpatialReference oGCS; oGCS.importFromEPSG( 4024 ); CopyGeogCSFrom( &oGCS ); } else { // If we don't know, default to WGS84 so there is something there. SetWellKnownGeogCS( "WGS84" ); } } } /* -------------------------------------------------------------------- */ /* Linear units translation */ /* -------------------------------------------------------------------- */ if( IsLocal() || IsProjected() ) { const char *pszValue; pszValue = OSR_GDS( papszPrj, "Units", NULL ); if( pszValue == NULL ) SetLinearUnitsAndUpdateParameters( SRS_UL_METER, 1.0 ); else if( EQUAL(pszValue,"FEET") ) SetLinearUnitsAndUpdateParameters( SRS_UL_FOOT, atof(SRS_UL_FOOT_CONV) ); else SetLinearUnitsAndUpdateParameters( pszValue, 1.0 ); } return OGRERR_NONE; } /************************************************************************/ /* morphToESRI() */ /************************************************************************/ /** * Convert in place from ESRI WKT format. * * The value notes of this coordinate system as modified in various manners * to adhere more closely to the WKT standard. This mostly involves * translating a variety of ESRI names for projections, arguments and * datums to "standard" names, as defined by Adam Gawne-Cain's reference * translation of EPSG to WKT for the CT specification. * * This does the same as the C function OSRMorphToESRI(). * * @return OGRERR_NONE unless something goes badly wrong. */ OGRErr OGRSpatialReference::morphToESRI() { OGRErr eErr; /* -------------------------------------------------------------------- */ /* Fixup ordering, missing linear units, etc. */ /* -------------------------------------------------------------------- */ eErr = Fixup(); if( eErr != OGRERR_NONE ) return eErr; /* -------------------------------------------------------------------- */ /* Strip all CT parameters (AXIS, AUTHORITY, TOWGS84, etc). */ /* -------------------------------------------------------------------- */ eErr = StripCTParms(); if( eErr != OGRERR_NONE ) return eErr; if( GetRoot() == NULL ) return OGRERR_NONE; /* -------------------------------------------------------------------- */ /* There is a special case for Hotine Oblique Mercator to split */ /* out the case with an angle to rectified grid. Bug 423 */ /* -------------------------------------------------------------------- */ const char *pszProjection = GetAttrValue("PROJECTION"); if( pszProjection != NULL && EQUAL(pszProjection,SRS_PT_HOTINE_OBLIQUE_MERCATOR) && fabs(GetProjParm(SRS_PP_AZIMUTH, 0.0 )-90) < 0.0001 && fabs(GetProjParm(SRS_PP_RECTIFIED_GRID_ANGLE, 0.0 )-90) < 0.0001 ) { SetNode( "PROJCS|PROJECTION", "Hotine_Oblique_Mercator_Azimuth_Center" ); /* ideally we should strip out of the rectified_grid_angle */ } /* -------------------------------------------------------------------- */ /* Translate PROJECTION keywords that are misnamed. */ /* -------------------------------------------------------------------- */ GetRoot()->applyRemapper( "PROJECTION", apszProjMapping+1, apszProjMapping, 2 ); /* -------------------------------------------------------------------- */ /* Translate DATUM keywords that are misnamed. */ /* -------------------------------------------------------------------- */ InitDatumMappingTable(); GetRoot()->applyRemapper( "DATUM", papszDatumMapping+2, papszDatumMapping+1, 3 ); /* -------------------------------------------------------------------- */ /* Very specific handling for some well known geographic */ /* coordinate systems. */ /* -------------------------------------------------------------------- */ OGR_SRSNode *poGeogCS = GetAttrNode( "GEOGCS" ); if( poGeogCS != NULL ) { const char *pszGeogCSName = poGeogCS->GetChild(0)->GetValue(); const char *pszAuthName = GetAuthorityName("GEOGCS"); const char *pszUTMPrefix = NULL; int nGCSCode = -1; if( pszAuthName != NULL && EQUAL(pszAuthName,"EPSG") ) nGCSCode = atoi(GetAuthorityCode("GEOGCS")); if( nGCSCode == 4326 || EQUAL(pszGeogCSName,"WGS84") || EQUAL(pszGeogCSName,"WGS 84") ) { poGeogCS->GetChild(0)->SetValue( "GCS_WGS_1984" ); pszUTMPrefix = "WGS_1984"; } else if( nGCSCode == 4267 || EQUAL(pszGeogCSName,"NAD27") || EQUAL(pszGeogCSName,"NAD 27") ) { poGeogCS->GetChild(0)->SetValue( "GCS_North_American_1927" ); pszUTMPrefix = "NAD_1927"; } else if( nGCSCode == 4269 || EQUAL(pszGeogCSName,"NAD83") || EQUAL(pszGeogCSName,"NAD 83") ) { poGeogCS->GetChild(0)->SetValue( "GCS_North_American_1983" ); pszUTMPrefix = "NAD_1983"; } /* -------------------------------------------------------------------- */ /* Force Unnamed to Unknown for most common locations. */ /* -------------------------------------------------------------------- */ static char *apszUnknownMapping[] = { "Unknown", "Unnamed", NULL, NULL }; GetRoot()->applyRemapper( "PROJCS", apszUnknownMapping+1, apszUnknownMapping+0, 2 ); GetRoot()->applyRemapper( "GEOGCS", apszUnknownMapping+1, apszUnknownMapping+0, 2 ); GetRoot()->applyRemapper( "DATUM", apszUnknownMapping+1, apszUnknownMapping+0, 2 ); GetRoot()->applyRemapper( "SPHEROID", apszUnknownMapping+1, apszUnknownMapping+0, 2 ); GetRoot()->applyRemapper( "PRIMEM", apszUnknownMapping+1, apszUnknownMapping+0, 2 ); /* -------------------------------------------------------------------- */ /* If the PROJCS name is unset, use the PROJECTION name in */ /* place of unknown, or unnamed. At the request of Peng Gao. */ /* -------------------------------------------------------------------- */ OGR_SRSNode *poNode; poNode = GetAttrNode( "PROJCS" ); if( poNode != NULL ) poNode = poNode->GetChild( 0 ); if( poNode != NULL && ( EQUAL(poNode->GetValue(),"unnamed") || EQUAL(poNode->GetValue(),"unknown") || EQUAL(poNode->GetValue(),"") ) ) { if( GetAttrValue( "PROJECTION", 0 ) != NULL ) poNode->SetValue( GetAttrValue( "PROJECTION", 0 ) ); } /* -------------------------------------------------------------------- */ /* Prepare very specific PROJCS names for UTM coordinate */ /* systems. */ /* -------------------------------------------------------------------- */ int bNorth, nZone; nZone = GetUTMZone( &bNorth ); if( nZone > 0 && pszUTMPrefix != NULL ) { char szUTMName[128]; if( bNorth ) sprintf( szUTMName, "%s_UTM_Zone_%dN", pszUTMPrefix, nZone ); else sprintf( szUTMName, "%s_UTM_Zone_%dS", pszUTMPrefix, nZone ); OGR_SRSNode *poProjCS = GetAttrNode( "PROJCS" ); if( poProjCS != NULL ) poProjCS->GetChild(0)->SetValue( szUTMName ); } } /* -------------------------------------------------------------------- */ /* Translate UNIT keywords that are misnamed, or even the wrong */ /* case. */ /* -------------------------------------------------------------------- */ GetRoot()->applyRemapper( "UNIT", apszUnitMapping+1, apszUnitMapping, 2 ); /* -------------------------------------------------------------------- */ /* reset constants for decimal degrees to the exact string ESRI */ /* expects when encountered to ensure a matchup. */ /* -------------------------------------------------------------------- */ OGR_SRSNode *poUnit = GetAttrNode( "GEOGCS|UNIT" ); if( poUnit != NULL && poUnit->GetChildCount() >= 2 && ABS(GetAngularUnits()-0.0174532925199433) < 0.00000000001 ) { poUnit->GetChild(0)->SetValue("Degree"); poUnit->GetChild(1)->SetValue("0.017453292519943295"); } /* -------------------------------------------------------------------- */ /* Make sure we reproduce US Feet exactly too. */ /* -------------------------------------------------------------------- */ poUnit = GetAttrNode( "PROJCS|UNIT" ); if( poUnit != NULL && poUnit->GetChildCount() >= 2 && ABS(GetLinearUnits()- 0.30480060960121924) < 0.000000000000001) { poUnit->GetChild(0)->SetValue("Foot_US"); poUnit->GetChild(1)->SetValue("0.30480060960121924"); } /* -------------------------------------------------------------------- */ /* Remap parameters used for Albers and Mercator. */ /* -------------------------------------------------------------------- */ pszProjection = GetAttrValue("PROJECTION"); if( pszProjection != NULL && EQUAL(pszProjection,"Albers") ) GetRoot()->applyRemapper( "PARAMETER", apszAlbersMapping + 1, apszAlbersMapping + 0, 2 ); if( pszProjection != NULL && EQUAL(pszProjection,"Mercator") ) GetRoot()->applyRemapper( "PARAMETER", apszAlbersMapping + 1, apszAlbersMapping + 0, 2 ); /* -------------------------------------------------------------------- */ /* Convert SPHEROID name to use underscores instead of spaces. */ /* -------------------------------------------------------------------- */ OGR_SRSNode *poSpheroid; poSpheroid = GetAttrNode( "SPHEROID" ); if( poSpheroid != NULL ) poSpheroid = poSpheroid->GetChild(0); if( poSpheroid != NULL ) { char *pszNewValue = CPLStrdup(RemapSpheroidName(poSpheroid->GetValue())); MorphNameToESRI( &pszNewValue ); poSpheroid->SetValue( pszNewValue ); CPLFree( pszNewValue ); } /* -------------------------------------------------------------------- */ /* Try to insert a D_ in front of the datum name. */ /* -------------------------------------------------------------------- */ OGR_SRSNode *poDatum; poDatum = GetAttrNode( "DATUM" ); if( poDatum != NULL ) poDatum = poDatum->GetChild(0); if( poDatum != NULL ) { if( !EQUALN(poDatum->GetValue(),"D_",2) ) { char *pszNewValue; pszNewValue = (char *) CPLMalloc(strlen(poDatum->GetValue())+3); strcpy( pszNewValue, "D_" ); strcat( pszNewValue, poDatum->GetValue() ); poDatum->SetValue( pszNewValue ); CPLFree( pszNewValue ); } } return OGRERR_NONE; } /************************************************************************/ /* OSRMorphToESRI() */ /************************************************************************/ OGRErr OSRMorphToESRI( OGRSpatialReferenceH hSRS ) { return ((OGRSpatialReference *) hSRS)->morphToESRI(); } /************************************************************************/ /* morphFromESRI() */ /* */ /* modify this definition from the ESRI definition of WKT to */ /* the "Standard" definition. */ /************************************************************************/ /** * Convert in place to ESRI WKT format. * * The value nodes of this coordinate system as modified in various manners * more closely map onto the ESRI concept of WKT format. This includes * renaming a variety of projections and arguments, and stripping out * nodes note recognised by ESRI (like AUTHORITY and AXIS). * * This does the same as the C function OSRMorphFromESRI(). * * @return OGRERR_NONE unless something goes badly wrong. */ OGRErr OGRSpatialReference::morphFromESRI() { OGRErr eErr = OGRERR_NONE; if( GetRoot() == NULL ) return OGRERR_NONE; /* -------------------------------------------------------------------- */ /* Translate DATUM keywords that are oddly named. */ /* -------------------------------------------------------------------- */ InitDatumMappingTable(); GetRoot()->applyRemapper( "DATUM", papszDatumMapping+1, papszDatumMapping+2, 3 ); /* -------------------------------------------------------------------- */ /* Try to remove any D_ in front of the datum name. */ /* -------------------------------------------------------------------- */ OGR_SRSNode *poDatum; poDatum = GetAttrNode( "DATUM" ); if( poDatum != NULL ) poDatum = poDatum->GetChild(0); if( poDatum != NULL ) { if( EQUALN(poDatum->GetValue(),"D_",2) ) { char *pszNewValue = CPLStrdup( poDatum->GetValue() + 2 ); poDatum->SetValue( pszNewValue ); CPLFree( pszNewValue ); } } /* -------------------------------------------------------------------- */ /* Split Lambert_Conformal_Conic into 1SP or 2SP form. */ /* */ /* See bugzilla.remotesensing.org/show_bug.cgi?id=187 */ /* -------------------------------------------------------------------- */ const char *pszProjection = GetAttrValue("PROJECTION"); if( pszProjection != NULL && EQUAL(pszProjection,"Lambert_Conformal_Conic") ) { if( GetProjParm( "Scale_Factor", 2.0 ) == 2.0 ) SetNode( "PROJCS|PROJECTION", SRS_PT_LAMBERT_CONFORMAL_CONIC_2SP ); else SetNode( "PROJCS|PROJECTION", SRS_PT_LAMBERT_CONFORMAL_CONIC_1SP ); } /* -------------------------------------------------------------------- */ /* If we are remapping Hotine_Oblique_Mercator_Azimuth_Center */ /* add a rectified_grid_angle parameter - to match the azimuth */ /* I guess. */ /* -------------------------------------------------------------------- */ if( pszProjection != NULL && EQUAL(pszProjection,"Hotine_Oblique_Mercator_Azimuth_Center") ) { SetProjParm( SRS_PP_RECTIFIED_GRID_ANGLE , GetProjParm( SRS_PP_AZIMUTH, 0.0 ) ); FixupOrdering(); } /* -------------------------------------------------------------------- */ /* Remap albers and Mercator parameters. */ /* -------------------------------------------------------------------- */ if( pszProjection != NULL && EQUAL(pszProjection,"Albers") ) GetRoot()->applyRemapper( "PARAMETER", apszAlbersMapping + 0, apszAlbersMapping + 1, 2 ); if( pszProjection != NULL && EQUAL(pszProjection,"Mercator") ) GetRoot()->applyRemapper( "PARAMETER", apszMercatorMapping + 0, apszMercatorMapping + 1, 2 ); /* -------------------------------------------------------------------- */ /* Translate PROJECTION keywords that are misnamed. */ /* -------------------------------------------------------------------- */ GetRoot()->applyRemapper( "PROJECTION", apszProjMapping, apszProjMapping+1, 2 ); /* -------------------------------------------------------------------- */ /* Translate DATUM keywords that are misnamed. */ /* -------------------------------------------------------------------- */ InitDatumMappingTable(); GetRoot()->applyRemapper( "DATUM", papszDatumMapping+1, papszDatumMapping+2, 3 ); return eErr; } /************************************************************************/ /* OSRMorphFromESRI() */ /************************************************************************/ OGRErr OSRMorphFromESRI( OGRSpatialReferenceH hSRS ) { return ((OGRSpatialReference *) hSRS)->morphFromESRI(); }