/********************************************************************** * $Id: geoconcept_syscoord.c$ * * Name: geoconcept_syscoord.c * Project: OpenGIS Simple Features Reference Implementation * Purpose: Implements translation between Geoconcept SysCoord * and OGRSpatialRef format * Language: C * ********************************************************************** * Copyright (c) 2007, Geoconcept and IGN * Copyright (c) 2008-2010, Even Rouault * * Permission is hereby granted, free of charge, to any person obtaining a * copy of this software and associated documentation files (the "Software"), * to deal in the Software without restriction, including without limitation * the rights to use, copy, modify, merge, publish, distribute, sublicense, * and/or sell copies of the Software, and to permit persons to whom the * Software is furnished to do so, subject to the following conditions: * * The above copyright notice and this permission notice shall be included * in all copies or substantial portions of the Software. * * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER * DEALINGS IN THE SOFTWARE. **********************************************************************/ #include "geoconcept_syscoord.h" #include "cpl_string.h" CPL_CVSID("$Id: geoconcept_syscoord.c 7eb5fc0b4dd0a6df3b94dccf36730770d3e11427 2018-04-05 17:48:43 +0200 Even Rouault $") /* -------------------------------------------------------------------- */ /* GCSRS globals */ /* -------------------------------------------------------------------- */ /* * The following information came from GEO CONCEPT PROJECTION files, * aka GCP files. * A lot of information has been added to these GCP. There are mostly * noticed as FIXME in the source. */ static const GCSysCoord gk_asSysCoordList[]= /* * pszSysCoordName, pszUnit, dfPM, dfLambda0, dfPhi0, dfk0, dfX0, dfY0, dfPhi1, dfPhi2, nDatumID, nProjID, coordSystemID, timeZoneValue * * #12, #14, #15, #17 : parameters listed below are "generic" ... * * Geoconcept uses cos(lat_ts) as scale factor, but cos(lat_ts)==cos(-lat_ts) : I then set dfPhi1 with lat_ts */ { {"Lambert 2 extended", NULL, 2.337229166667, 0.000000000, 46.80000000,0.99987742000, 600000.000, 2200000.000, 0.0, 0.0, 13, 2, 1,-1}, {"Lambert 1", NULL, 2.337229166667, 0.000000000, 49.50000000,0.99987734000, 600000.000, 200000.000, 0.0, 0.0, 13, 2, 2,-1}, {"Lambert 2", NULL, 2.337229166667, 0.000000000, 46.80000000,0.99987742000, 600000.000, 200000.000, 0.0, 0.0, 13, 2, 3,-1}, {"Lambert 3", NULL, 2.337229166667, 0.000000000, 44.10000000,0.99987750000, 600000.000, 200000.000, 0.0, 0.0, 13, 2, 4,-1}, {"Lambert 4", NULL, 2.337229166667, 0.000000000, 42.16500000,0.99994471000, 234.358, 185861.369, 0.0, 0.0, 13, 2, 5,-1}, {"Bonne NTF", NULL, 2.337222222222, 0.000000000, 48.86000000,1.00000000000, 0.000, 0.000, 0.0, 0.0, 1, 3, 11,-1}, {"UTM Nord - ED50", NULL, 0.000000000000, 0.000000000, 0.00000000,0.99960000000, 500000.000, 0.000, 0.0, 0.0, 14, 1, 12, 0}, {"Plate carr""\xe9""e", NULL, 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 11, 4, 13,-1}, {"MGRS (Military UTM)", NULL, 0.000000000000, 0.000000000, 0.00000000,0.99960000000, 0.000, 0.000, 0.0, 0.0, 4, 11, 14,-1}, {"UTM Sud - WGS84", NULL, 0.000000000000, 0.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0, 4, 1, 15, 0}, {"National GB projection", NULL, 0.000000000000, -2.000000000, 49.00000000,0.99960127170, 400000.000, -100000.000, 0.0, 0.0, 12, 12, 16,-1}, {"UTM Nord - WGS84", NULL, 0.000000000000, 0.000000000, 0.00000000,0.99960000000, 500000.000, 0.000, 0.0, 0.0, 4, 1, 17, 0}, {"UTM Nord - WGS84", NULL, 0.000000000000, 0.000000000, 0.00000000,0.99960000000, 500000.000, 0.000, 0.0, 0.0,9990, 1, 17, 0}, {"Lambert 2 ""\xe9""tendu - sans grille", NULL, 2.337229166667, 0.000000000, 46.80000000,0.99987742000, 600000.000, 2200000.000, 0.0, 0.0, 1, 2, 91,-1}, {"Lambert 1 - sans grille", NULL, 2.337229166667, 0.000000000, 49.50000000,0.99987734000, 600000.000, 200000.000, 0.0, 0.0, 1, 2, 92,-1}, {"Lambert 2 - sans grille", NULL, 2.337229166667, 0.000000000, 46.80000000,0.99987742000, 600000.000, 200000.000, 0.0, 0.0, 1, 2, 93,-1}, {"Lambert 3 - sans grille", NULL, 2.337229166667, 0.000000000, 44.10000000,0.99987750000, 600000.000, 200000.000, 0.0, 0.0, 1, 2, 94,-1}, {"Lambert 4 - sans grille", NULL, 2.337229166667, 0.000000000, 42.16500000,0.99994471000, 234.358, 185861.369, 0.0, 0.0, 1, 2, 95,-1}, {"(Long/Lat) NTF", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 1, 0, 100,-1}, {"(Long/Lat) WGS84", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 4, 0, 101,-1}, {"(Long/Lat) ED50", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 14, 0, 102,-1}, {"(Long/Lat) Australian 1984", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 7, 0, 103,-1}, {"(Long/Lat) Airy", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 12, 0, 104,-1}, {"(Long/Lat) NTF Paris (gr)", "gr", 2.337229166667, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 1, 0, 105,-1}, {"(Long/Lat) WGS 72", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 3, 0, 107,-1}, {"Geoportail MILLER", NULL, 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 4, 24, 222,-1}, {"IGN-RRAFGUADU20", NULL, 0.000000000000, -63.000000000, 0.00000000,0.99960000000, 500000.000, 0.000, 0.0, 0.0,9984, 1, 501,-1},/* FIXME does not exist in IGNF, use IGN-UTM20W84GUAD instead */ {"IGN-RRAFMARTU20", NULL, 0.000000000000, -63.000000000, 0.00000000,0.99960000000, 500000.000, 0.000, 0.0, 0.0,9984, 1, 502,-1},/* FIXME does not exist in IGNF, use IGN-UTM20W84MART instead, never reached cause identical to 501:-1 */ {"IGN-RGM04UTM38S", NULL, 0.000000000000, 45.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 503,-1},/* FIXME 5030 datum changed into 9984 */ {"IGN-RGR92UTM40S", NULL, 0.000000000000, 57.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 504,-1}, {"IGN-UTM22RGFG95", NULL, 0.000000000000, -51.000000000, 0.00000000,0.99960000000, 500000.000, 0.000, 0.0, 0.0,9984, 1, 505,-1}, {"IGN-UTM01SWG84", NULL, 0.000000000000,-177.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 506,-1},/* never reached cause identical to 15:1 */ {"IGN-RGSPM06U21", NULL, 0.000000000000, -57.000000000, 0.00000000,0.99960000000, 500000.000, 0.000, 0.0, 0.0,9984, 1, 507,-1}, {"IGN-RGPFUTM5S", NULL, 0.000000000000,-153.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 508,-1}, {"IGN-RGPFUTM6S", NULL, 0.000000000000,-147.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 509,-1}, {"IGN-RGPFUTM7S", NULL, 0.000000000000,-141.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 510,-1}, {"IGN-CROZ63UTM39S", NULL, 0.000000000000, 51.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9983, 1, 511,-1}, {"IGN-WGS84UTM1S", NULL, 0.000000000000,-177.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0, 4, 1, 512,-1}, {"IGN-RGNCUTM57S", NULL, 0.000000000000, 159.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 513,-1}, {"IGN-RGNCUTM58S", NULL, 0.000000000000, 165.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 514,-1}, {"IGN-RGNCUTM59S", NULL, 0.000000000000, 171.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9984, 1, 515,-1}, {"IGN-KERG62UTM42S", NULL, 0.000000000000, 69.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0,9988, 1, 516,-1}, {"IGN-REUN47GAUSSL", NULL, 0.000000000000, 55.533333333,-21.11666667,1.00000000000, 160000.000, 50000.000, 0.0, 0.0, 2, 19, 520,-1}, {"Lambert 1 Carto", NULL, 2.337229166667, 0.000000000, 49.50000000,0.99987734000, 600000.000, 1200000.000, 0.0, 0.0, 13, 2, 1002,-1}, {"Lambert 2 Carto", NULL, 2.337229166667, 0.000000000, 46.80000000,0.99987742000, 600000.000, 2200000.000, 0.0, 0.0, 13, 2, 1003,-1},/* never reached cause identical to 1:-1 */ {"Lambert 3 Carto", NULL, 2.337229166667, 0.000000000, 44.10000000,0.99987750000, 600000.000, 3200000.000, 0.0, 0.0, 13, 2, 1004,-1}, {"Lambert 4 Carto", NULL, 2.337229166667, 0.000000000, 42.16500000,0.99994471000, 234.358, 4185861.369, 0.0, 0.0, 13, 2, 1005,-1}, {"Lambert 93", NULL, 0.000000000000, 3.000000000, 46.50000000,0.00000000000, 700000.000, 6600000.000, 44.0, 49.0,9984, 18, 1006,-1}, {"IGN-RGNCLAM", NULL, 0.000000000000, 166.000000000,-21.30000000,0.00000000000, 400000.000, 300000.000,-20.4,-22.2,9984, 18, 1007,-1},/* Added in GCP */ {"Lambert 1 Carto - sans grille", NULL, 2.337229166667, 0.000000000, 49.50000000,0.99987734000, 600000.000, 1200000.000, 0.0, 0.0, 1, 2, 1092,-1}, {"Lambert 2 Carto - sans grille", NULL, 2.337229166667, 0.000000000, 46.80000000,0.99987742000, 600000.000, 2200000.000, 0.0, 0.0, 1, 2, 1093,-1}, {"Lambert 3 Carto - sans grille", NULL, 2.337229166667, 0.000000000, 44.10000000,0.99987750000, 600000.000, 3200000.000, 0.0, 0.0, 1, 2, 1094,-1}, {"Lambert 4 Carto - sans grille", NULL, 2.337229166667, 0.000000000, 42.16500000,0.99994471000, 234.358, 185861.369, 0.0, 0.0, 1, 2, 1095,-1}, {"Suisse", NULL, 0.000000000000, 7.439583333, 46.95240556,1.00000000000, 600000.000, 200000.000, 0.0, 0.0, 2, 25, 1556,-1}, {"Geoportail France", NULL, 0.000000000000, 0.000000000, 0.00000000,0.68835457569, 0.000, 0.000, 46.5, 0.0,9984, 26, 2012,-1}, {"Geoportail Antilles", NULL, 0.000000000000, 0.000000000, 0.00000000,0.96592582629, 0.000, 0.000, 15.0, 0.0,9984, 26, 2016,-1}, {"Geoportail Guyane", NULL, 0.000000000000, 0.000000000, 0.00000000,0.99756405026, 0.000, 0.000, 4.0, 0.0,9984, 26, 2017,-1}, {"Geoportail Reunion", NULL, 0.000000000000, 0.000000000, 0.00000000,0.93358042649, 0.000, 0.000,-21.0, 0.0,9984, 26, 2018,-1}, {"Geoportail Mayotte", NULL, 0.000000000000, 0.000000000, 0.00000000,0.97814760073, 0.000, 0.000,-12.0, 0.0,9984, 26, 2019,-1}, {"Geoportail ST Pierre et Miquelon",NULL, 0.000000000000, 0.000000000, 0.00000000,0.68199836006, 0.000, 0.000, 47.0, 0.0,9984, 26, 2020,-1}, {"Geoportail Nouvelle Caledonie", NULL, 0.000000000000, 0.000000000, 0.00000000,0.92718385456, 0.000, 0.000,-22.0, 0.0,9984, 26, 2021,-1}, {"Geoportail Wallis", NULL, 0.000000000000, 0.000000000, 0.00000000,0.97029572627, 0.000, 0.000,-14.0, 0.0,9984, 26, 2022,-1}, {"Geoportail Polynesie", NULL, 0.000000000000, 0.000000000, 0.00000000,0.96592582628, 0.000, 0.000,-15.0, 0.0,9984, 26, 2023,-1}, {"Mercator sur sph""\xe8""re WGS84", NULL, 0.000000000000, 0.000000000, 0.00000000,1.00000000000, 0.000, 0.000, 0.0, 0.0,2015, 21, 2027,-1}, {"(Long/Lat) RGF 93", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 13, 0, 2028,-1}, {"(Long/Lat) ITRS-89", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9984, 0, 2028,-1}, {"Geoportail Crozet", NULL, 0.000000000000, 0.000000000, 0.00000000,0.69465837046, 0.000, 0.000,-46.0, 0.0,9984, 26, 2040,-1},/* FIXME : wrong scale factor was 0.69088241108 */ {"Geoportail Kerguelen", NULL, 0.000000000000, 0.000000000, 0.00000000,0.64944804833, 0.000, 0.000,-49.5, 0.0,9984, 26, 2042,-1},/* FIXME : wrong scale factor was 0.67815966987 */ {"Lambert CC 42", NULL, 0.000000000000, 3.000000000, 42.00000000,0.00000000000,1700000.000, 1200000.000, 41.2, 42.8,9984, 18, 2501,-1}, {"Lambert CC 43", NULL, 0.000000000000, 3.000000000, 43.00000000,0.00000000000,1700000.000, 2200000.000, 42.2, 43.8,9984, 18, 2502,-1}, {"Lambert CC 44", NULL, 0.000000000000, 3.000000000, 44.00000000,0.00000000000,1700000.000, 3200000.000, 43.2, 44.8,9984, 18, 2503,-1}, {"Lambert CC 45", NULL, 0.000000000000, 3.000000000, 45.00000000,0.00000000000,1700000.000, 4200000.000, 44.2, 45.8,9984, 18, 2504,-1}, {"Lambert CC 46", NULL, 0.000000000000, 3.000000000, 46.00000000,0.00000000000,1700000.000, 5200000.000, 45.2, 46.8,9984, 18, 2505,-1}, {"Lambert CC 47", NULL, 0.000000000000, 3.000000000, 47.00000000,0.00000000000,1700000.000, 6200000.000, 46.2, 47.8,9984, 18, 2506,-1}, {"Lambert CC 48", NULL, 0.000000000000, 3.000000000, 48.00000000,0.00000000000,1700000.000, 7200000.000, 47.2, 48.8,9984, 18, 2507,-1}, {"Lambert CC 49", NULL, 0.000000000000, 3.000000000, 49.00000000,0.00000000000,1700000.000, 8200000.000, 48.2, 49.8,9984, 18, 2508,-1}, {"Lambert CC 50", NULL, 0.000000000000, 3.000000000, 50.00000000,0.00000000000,1700000.000, 9200000.000, 49.2, 50.8,9984, 18, 2509,-1}, {"(Long/Lat) IGN-RGM04GEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9984, 0,10001,-1}, {"(Long/Lat) IGN-RGFG95GEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9984, 0,10002,-1},/* never reached, identical to 10001:-1 */ {"(Long/Lat) IGN-WGS84RRAFGEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9984, 0,10003,-1},/* never reached, identical to 10001:-1 */ {"(Long/Lat) IGN-RGR92GEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9984, 0,10004,-1},/* never reached, identical to 10001:-1 */ {"(Long/Lat) IGN-WGS84G", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 4, 0,10005,-1}, {"(Long/Lat) CROZ63GEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, 4, 0,10006,-1}, {"(Long/Lat) RGSPM06GEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9984, 0,10007,-1},/* never reached, identical to 10001:-1 */ {"(Long/Lat) RGPFGEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9984, 0,10008,-1},/* never reached, identical to 10001:-1 */ {"(Long/Lat) RGNCGEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9984, 0,10009,-1},/* never reached, identical to 10001:-1 */ {"(Long/Lat) KER62GEO", "d", 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0,9988, 0,10010,-1}, {"UTM Sud - ED50", NULL, 0.000000000000, 0.000000000, 0.00000000,0.99960000000, 500000.000,10000000.000, 0.0, 0.0, 14, 1,99912, 0},/* FIXME allow retrieving 12:0 - See _findSysCoord_GCSRS() */ {NULL, NULL, 0.000000000000, 0.000000000, 0.00000000,0.00000000000, 0.000, 0.000, 0.0, 0.0, -1, -1, -1,-1} }; static const GCProjectionInfo gk_asProjList[]= /* * pszProjName, nSphere, nProjID */ { {"Geographic shift", 0, 0}, {"UTM", 0, 1}, {"Lambert Conform Conic", 0, 2}, {"Bonne", 0, 3}, {"Plate carr""\xe9""e", 0, 4}, {"MGRS (Military UTM)", 0, 11}, {"Transversal Mercator", 0, 12}, {"Lambert secant", 0, 18}, {"Gauss Laborde", 1, 19}, {"Polyconic", 0, 20}, {"Direct Mercator", 0, 21}, {"Stereographic oblic", 1, 22}, {"Miller", 0, 24}, {"Mercator oblic", 1, 25}, {"Equi rectangular", 1, 26}, {NULL, 0, -1} }; static const GCDatumInfo gk_asDatumList[]= /* * pszDatumName, dfShiftX, dfShiftY, dfShiftZ, dfRotX, dfRotY, dfRotZ, dfScaleFactor, dfFA, dfFlattening, nEllipsoidID, nDatumID */ /* * Wrong dx, dy, dz : * IGN-RGM04GEO, was -217, -216, 67 * IGN-RGFG95GEO, was -2, -2, 2 * IGN-RGSPM06GEO, was -125.593, 143.763, -194.558 * * #1 and #13 are identical * #8, #11, #2015 are spherical views of #4 * #5030, #5031 and #5032 are identical * FIXME : #5030, #5031, #5032 are ITRS89 compliant, so "compatible" with #4, better use #9999 as ellipsoid * FIXME : #9999 to #9986 added */ { {"NTF (Clarke 1880)", -168.0000, -60.0000, 320.0000, 0.00000, 0.00000, 0.00000, 0.0, -112.200,-54.7388e-6, 3, 1}, {"ED50 France (International 1909)", -84.0000, -97.0000,-117.0000, 0.00000, 0.00000, 0.00000, 0.0, -251.000,-14.1927e-6, 5, 2}, {"WGS 72", 0.0000, 12.0000, 6.0000, 0.00000, 0.00000, 0.00000, 0.0, 2.000, 0.0312e-6, 6, 3}, {"WGS_1984", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, 9999, 4}, {"ED 79", -83.0000, -95.0000,-116.0000, 0.00000, 0.00000, 0.00000, 0.0, -251.000,-14.1927e-6, 5, 5}, {"Australian Geodetic 1966", -133.0000, -48.0000, 148.0000, 0.00000, 0.00000, 0.00000, 0.0, -23.000, -0.0081e-6, 7, 6}, {"Australian Geodetic 1984", -134.0000, -48.0000, 149.0000, 0.00000, 0.00000, 0.00000, 0.0, -23.000, -0.0081e-6, 7, 7}, {"Sphere", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, 1, 8}, {"Sphere DCW", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, 1, 11}, {"Airy", 375.0000,-111.0000, 431.0000, 0.00000, 0.00000, 0.00000, 0.0, 573.604, 11.96002325e-6, 8, 12}, {"NTF-Grille", -168.0000, -60.0000, 320.0000, 0.00000, 0.00000, 0.00000, 0.0, -112.200,-54.7388e-6, 3, 13}, {"ED50 (International 1909)", -87.0000, -98.0000,-121.0000, 0.00000, 0.00000, 0.00000, 0.0, -251.000,-14.1927e-6, 5, 14}, {"WGS 84 sur sphere", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, 1,2015}, {"IGN-RGM04GEO", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, 4,5030}, {"IGN-RGFG95GEO", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, 4,5031}, {"IGN-RGSPM06GEO", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, 4,5032}, {"IGN-WALL78", 253.0000,-133.0000,-127.0000, 0.00000, 0.00000, 0.00000, 0.0, -251.000,-14.1927e-6, 5,9999},/* FIXME */ {"IGN-TAHA", 72.4380, 345.9180, 79.4860,-1.60450,-0.88230, -0.55650, 1.3746e-6, -251.000,-14.1927e-6, 5,9998},/* FIXME */ {"IGN-MOOREA87", 215.9820, 149.5930, 176.2290, 3.26240, 1.69200, 1.15710, 10.47730e-6,-251.000,-14.1927e-6, 5,9997},/* FIXME */ {"IGN-TAHI51", 162.0000, 117.0000, 154.0000, 0.00000, 0.00000, 0.00000, 0.0, -251.000,-14.1927e-6, 5,9996},/* FIXME */ {"IGN-NUKU72", 165.7320, 216.7200, 180.5050,-0.64340,-0.45120, -0.07910, 7.42040e-6,-251.000,-14.1927e-6, 5,9995},/* FIXME */ {"IGN-IGN63", 410.7210, 55.0490, 80.7460,-2.57790,-2.35140, -0.66640, 17.33110e-6,-251.000,-14.1927e-6, 5,9994},/* FIXME */ {"IGN-MART38", 126.9260, 547.9390, 130.4090,-2.78670, 5.16124, -0.85844, 13.82265e-6,-251.000,-14.1927e-6, 5,9993},/* FIXME */ {"IGN-GUAD48", -472.2900, -5.6300,-304.1200, 0.43620,-0.83740, 0.25630, 1.89840e-6,-251.000,-14.1927e-6, 5,9992},/* FIXME */ {"IGN-GUADFM49", 136.5960, 248.1480,-429.7890, 0.00000, 0.00000, 0.00000, 0.0, -251.000,-14.1927e-6, 5,9991},/* FIXME */ {"IGN-STPM50", -95.5930, 573.7630, 173.4420,-0.96020, 1.25100, -1.39180, 42.62650e-6, -69.400,-37.2957e-6, 2,9990},/* FIXME */ {"IGN-CSG67", -193.0660, 236.9930, 105.4470, 0.48140,-0.80740, 0.12760, 1.56490e-6,-251.000,-14.1927e-6, 5,9989},/* FIXME */ {"IGN-KERG62", 145.0000,-187.0000, 103.0000, 0.00000, 0.00000, 0.00000, 0.0, -251.000,-14.1927e-6, 5,9988},/* FIXME */ {"IGN-REUN47", 789.5240,-626.4860, -89.9040, 0.60060,76.79460,-10.57880,-32.32410e-6,-251.000,-14.1927e-6, 5,9987},/* FIXME */ {"IGN-MAYO50", -599.9280,-275.5520,-195.6650, 0.08350, 0.47150, -0.06020,-49.28140e-6,-251.000,-14.1927e-6, 5,9986},/* FIXME */ {"IGN-TAHI79", 221.5250, 152.9480, 176.7680, 2.38470, 1.38960, 0.87700, 11.47410e-6,-251.000,-14.1927e-6, 5,9985},/* FIXME */ {"ITRS-89", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, 4,9984}, {"IGN-CROZ63", 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, -251.000,-14.1927e-6, 5,9983},/* FIXME added cause the Bursa-Wolf parameters are not known */ {NULL, 0.0000, 0.0000, 0.0000, 0.00000, 0.00000, 0.00000, 0.0, 0.000, 0.0, -1, -1} }; static const GCSpheroidInfo gk_asSpheroidList[]= /* * pszSpheroidName, dfA, dfE, nEllipsoidID * * cause Geoconcept assimilates WGS84 and GRS80, WGS84 is added to the list */ { {"Sphere", 6378137.0000, 0.00000000000000, 1}, {"Clarke 1866", 6378206.4000, 0.08227185423947, 2},/* Wrong, semi-major was 6378249.4000 */ {"Clarke 1880", 6378249.2000, 0.08248325676300, 3},/* Wrong, eccentricity was 0.082483256945 */ {"GRS 80", 6378137.0000, 0.08181919104300, 4},/* Wrong, eccentricity was 0.081819191060 */ {"International 1909", 6378388.0000, 0.08199188997900, 5}, {"WGS 72", 6378135.0000, 0.08181881201777, 6}, {"Australian National", 6378160.0000, 0.08182017998700, 7}, {"Airy", 6377563.3960, 0.08167337387420, 8}, {"WGS 84", 6378137.0000, 0.08181919084262,9999}, {NULL, 0, 0, -1} }; /* -------------------------------------------------------------------- */ /* GCSRS API Prototypes */ /* -------------------------------------------------------------------- */ /* -------------------------------------------------------------------- */ static int GCSRSAPI_CALL _areCompatibleSpheroids_GCSRS ( int id1, int id2 ) { if( id1==id2 ) return TRUE; switch( id1 ) { case 4 : case 9999 : switch( id2 ) { case 4 : case 9999 : return TRUE; default : break; } break; default : break; } return FALSE; }/* _areCompatibleSpheroids_GCSRS */ /* -------------------------------------------------------------------- */ static int GCSRSAPI_CALL _areCompatibleDatums_GCSRS ( int id1, int id2 ) { if( id1==id2 ) return TRUE; switch( id1 ) { case 1 : /* NTF */ case 13 : switch( id2 ) { case 1 : case 13 : return TRUE; default : break; } break; case 2 : /* ED50 */ case 14 : case 9983 : case 9985 : case 9986 : case 9987 : case 9989 : case 9991 : case 9992 : case 9993 : case 9994 : case 9995 : case 9997 : case 9998 : case 9999 : switch( id2 ) { case 2 : case 14 : case 9983 : case 9985 : case 9986 : case 9987 : case 9989 : case 9991 : case 9992 : case 9993 : case 9994 : case 9995 : case 9997 : case 9998 : case 9999 : return TRUE; default : break; } break; case 4 : /* WGS84 - ITRS89 */ case 8 : case 11 : case 2015 : case 5030 : case 5031 : case 5032 : case 9984 : switch( id2 ) { case 4 : case 8 : case 11 : case 2015 : case 5030 : case 5031 : case 5032 : case 9984 : return TRUE; default : break; } break; default : break; } return FALSE; }/* _areCompatibleDatums_GCSRS */ #define CPLDebugSpheroid_GCSRS(e) \ CPLDebug( "GEOCONCEPT", "SemiMajor:%.4f;Eccentricity:%.10f;",\ GetInfoSpheroidSemiMajor_GCSRS(e),\ GetInfoSpheroidExcentricity_GCSRS(e)\ ); /* -------------------------------------------------------------------- */ static const GCSpheroidInfo GCSRSAPI_CALL1(*) _findSpheroid_GCSRS ( double a, double rf ) { int iSpheroid, iResol= 0, nResol= 2; const GCSpheroidInfo* ell; double e, p[]= {1e-10, 1e-8}; /* f = 1 - sqrt(1 - e^2) */ e= (rf == 0.0) ? 0.0 : 1.0/rf; e= sqrt(e*(2.0-e)); ell_relax: for( iSpheroid= 0, ell= &(gk_asSpheroidList[0]); GetInfoSpheroidID_GCSRS(ell)!=-1; iSpheroid++, ell= &(gk_asSpheroidList[iSpheroid]) ) { if( fabs(GetInfoSpheroidSemiMajor_GCSRS(ell) - a) > 1e-4 ) continue; if( fabs(GetInfoSpheroidExcentricity_GCSRS(ell) - e) > p[iResol] ) continue; break; } if( GetInfoSpheroidID_GCSRS(ell)==-1 && iResol!=nResol-1 ) { iResol++; goto ell_relax; } return ell; }/* _findSpheroid_GCSRS */ #define CPLDebugDatum_GCSRS(d) \ CPLDebug( "GEOCONCEPT", "ID:%d;ShiftX:%.4f;ShiftY:%.4f;ShiftZ:%.4f;DiffA:%.4f;DiffFlattening:%.7f;",\ GetInfoDatumID_GCSRS((d)),\ GetInfoDatumShiftX_GCSRS((d)),\ GetInfoDatumShiftY_GCSRS((d)),\ GetInfoDatumShiftZ_GCSRS((d)),\ GetInfoDatumDiffA_GCSRS((d)),\ GetInfoDatumDiffFlattening_GCSRS((d))\ ); /* -------------------------------------------------------------------- */ static const GCDatumInfo GCSRSAPI_CALL1(*) _findDatum_GCSRS ( double dx, double dy, double dz, double a, double f ) { int iDatum, bRelax= FALSE; const GCDatumInfo* datum; datum_relax: for( iDatum= 0, datum= &(gk_asDatumList[0]); GetInfoDatumID_GCSRS(datum)!=-1; iDatum++, datum= &(gk_asDatumList[iDatum]) ) { if( !bRelax ) { if( fabs(GetInfoDatumShiftX_GCSRS(datum) - dx) > 1e-4 ) continue; if( fabs(GetInfoDatumShiftY_GCSRS(datum) - dy) > 1e-4 ) continue; if( fabs(GetInfoDatumShiftZ_GCSRS(datum) - dz) > 1e-4 ) continue; } if( fabs(GetInfoDatumDiffA_GCSRS(datum) - (6378137.0000-a)) > 1e-4 ) continue; if( fabs(GetInfoDatumDiffFlattening_GCSRS(datum) - (0.003352779565406696648-f)) > 1e-7 ) continue; break; } if( GetInfoDatumID_GCSRS(datum)==-1 && !bRelax ) { /* * FIXME : when both nadgrids and towgs84 are defined, bursa-wolf parameters are lost ! * if the projection and the ellipsoid are known, one can retrieve the datum * Try relaxed search ... */ bRelax= TRUE; goto datum_relax; } return datum; }/* _findDatum_GCSRS */ /* -------------------------------------------------------------------- */ static const GCProjectionInfo GCSRSAPI_CALL1(*) _findProjection_GCSRS ( const char* p, double lat_ts ) { int iProj; const GCProjectionInfo* proj; for( iProj= 0, proj= &(gk_asProjList[0]); GetInfoProjID_GCSRS(proj)!=-1; iProj++, proj= &(gk_asProjList[iProj]) ) { if( iProj==0 && p==NULL) break; if( iProj==1 && ( EQUAL(p,SRS_PT_TRANSVERSE_MERCATOR) || EQUAL(p,SRS_PT_TRANSVERSE_MERCATOR_SOUTH_ORIENTED) ) ) break; if( iProj==2 && EQUAL(p,SRS_PT_LAMBERT_CONFORMAL_CONIC_1SP) ) break; if( iProj==3 && EQUAL(p,SRS_PT_BONNE) ) break; if( iProj==4 && EQUAL(p,SRS_PT_EQUIRECTANGULAR) && lat_ts==0.0 ) break; /* FIXME : iProj==6 ? */ if( iProj==7 && ( EQUAL(p,SRS_PT_LAMBERT_CONFORMAL_CONIC_2SP) || EQUAL(p,SRS_PT_LAMBERT_CONFORMAL_CONIC_2SP_BELGIUM) ) ) break; if( iProj==8 && EQUAL(p,SRS_PT_GAUSSSCHREIBERTMERCATOR) ) break; if( iProj==9 && EQUAL(p,SRS_PT_POLYCONIC) ) break; /* FIXME if( iProj==10 && ( EQUAL(p,SRS_PT_MERCATOR_1SP) || EQUAL(p,SRS_PT_MERCATOR_2SP) ) ) break; */ if( iProj==11 && ( EQUAL(p,SRS_PT_OBLIQUE_STEREOGRAPHIC) || EQUAL(p,SRS_PT_POLAR_STEREOGRAPHIC) ) ) break; if( iProj==12 && EQUAL(p,SRS_PT_MILLER_CYLINDRICAL) ) break; /* FIXME if( iProj==13 && ( EQUAL(p,SRS_PT_HOTINE_OBLIQUE_MERCATOR) || EQUAL(p,SRS_PT_HOTINE_OBLIQUE_MERCATOR_TWO_POINT_NATURAL_ORIGIN) || EQUAL(p,SRS_PT_LABORDE_OBLIQUE_MERCATOR) ) ) break; */ if( iProj==14 && EQUAL(p,SRS_PT_EQUIRECTANGULAR) && lat_ts!=0.0 ) break; } return proj; }/* _findProjection_GCSRS */ #define CPLDebugSysCoord_GCSRS(m,s) \ CPLDebug( "GEOCONCEPT", "[%s]ID=%d;Zone=%d;DatumID=%d;ProjID=%d;PrimeMeridian=%.10f;CentralMeridian=%.10f;LatitudeOfOrigin=%.10f;StandardParallel1=%.10f;StandardParallel2=%.10f;ScaleFactor=%.10f;FalseEasting=%.10f;FalseNorthing=%.10f;",\ (m)? (m):"",\ GetSysCoordSystemID_GCSRS((s)),\ GetSysCoordTimeZone_GCSRS((s)),\ GetSysCoordDatumID_GCSRS((s)),\ GetSysCoordProjID_GCSRS((s)),\ GetSysCoordPrimeMeridian_GCSRS((s)),\ GetSysCoordCentralMeridian_GCSRS((s)),\ GetSysCoordLatitudeOfOrigin_GCSRS((s)),\ GetSysCoordStandardParallel1_GCSRS((s)),\ GetSysCoordStandardParallel2_GCSRS((s)),\ GetSysCoordScaleFactor_GCSRS((s)),\ GetSysCoordFalseEasting_GCSRS((s)),\ GetSysCoordFalseNorthing_GCSRS((s))\ ); /* -------------------------------------------------------------------- */ static GCSysCoord GCSRSAPI_CALL1(*) _findSysCoord_GCSRS ( GCSysCoord* theSysCoord ) { int iSysCoord, bestSysCoord= -1; const GCSysCoord* gcsc; if( !theSysCoord) return NULL; SetSysCoordSystemID_GCSRS(theSysCoord, -1); SetSysCoordTimeZone_GCSRS(theSysCoord, -1); CPLDebugSysCoord_GCSRS(NULL,theSysCoord); for( iSysCoord= 0, gcsc= &(gk_asSysCoordList[0]); GetSysCoordSystemID_GCSRS(gcsc)!=-1; iSysCoord++, gcsc= &(gk_asSysCoordList[iSysCoord]) ) { if( !_areCompatibleDatums_GCSRS(GetSysCoordDatumID_GCSRS(gcsc), GetSysCoordDatumID_GCSRS(theSysCoord)) ) continue; if( GetSysCoordProjID_GCSRS(gcsc) != GetSysCoordProjID_GCSRS(theSysCoord) ) continue; if( fabs(GetSysCoordPrimeMeridian_GCSRS(gcsc) - GetSysCoordPrimeMeridian_GCSRS(theSysCoord) ) > 1e-8 ) continue; if( fabs(GetSysCoordCentralMeridian_GCSRS(gcsc) - GetSysCoordCentralMeridian_GCSRS(theSysCoord) ) > 1e-8 ) { /* UTM family: central meridian is the 6* zone - 183 (in degrees) */ if( GetSysCoordProjID_GCSRS(gcsc) == 1 && /* generic UTM definition */ GetSysCoordCentralMeridian_GCSRS(gcsc)==0.0 ) { /* go on */ } else { continue; } } if( fabs(GetSysCoordLatitudeOfOrigin_GCSRS(gcsc) - GetSysCoordLatitudeOfOrigin_GCSRS(theSysCoord) ) > 1e-8 ) continue; if( fabs(GetSysCoordStandardParallel1_GCSRS(gcsc) - GetSysCoordStandardParallel1_GCSRS(theSysCoord) ) > 1e-8 ) continue; if( fabs(GetSysCoordStandardParallel2_GCSRS(gcsc) - GetSysCoordStandardParallel2_GCSRS(theSysCoord) ) > 1e-8 ) continue; if( fabs(GetSysCoordScaleFactor_GCSRS(gcsc) - GetSysCoordScaleFactor_GCSRS(theSysCoord) ) > 1e-8 ) continue; if( fabs(GetSysCoordFalseEasting_GCSRS(gcsc) - GetSysCoordFalseEasting_GCSRS(theSysCoord) ) > 1e-4 ) continue; if( fabs(GetSysCoordFalseNorthing_GCSRS(gcsc) - GetSysCoordFalseNorthing_GCSRS(theSysCoord) ) > 1e-4 ) continue; /* Found a candidate : */ if( bestSysCoord==-1) { bestSysCoord= iSysCoord; } else { switch( GetSysCoordProjID_GCSRS(gcsc) ) { case 0:/* long/lat */ if( GetSysCoordDatumID_GCSRS(gcsc)==GetSysCoordDatumID_GCSRS(theSysCoord) && GetSysCoordDatumID_GCSRS(&(gk_asSysCoordList[bestSysCoord]))!=GetSysCoordDatumID_GCSRS(theSysCoord)) /* exact match */ { bestSysCoord= iSysCoord; } break; case 1: /* UTM family: central meridian is the 6* zone - 183 (in degrees) */ if( GetSysCoordCentralMeridian_GCSRS(gcsc)!=0.0 && GetSysCoordDatumID_GCSRS(gcsc)==GetSysCoordDatumID_GCSRS(theSysCoord) && GetSysCoordDatumID_GCSRS(&(gk_asSysCoordList[bestSysCoord]))!=GetSysCoordDatumID_GCSRS(theSysCoord)) /* exact match */ { bestSysCoord= iSysCoord; } break; default : break; } } } /* Seems to be the right Geoconcept system: */ if( bestSysCoord>=0 ) { gcsc= &(gk_asSysCoordList[bestSysCoord]); switch( GetSysCoordSystemID_GCSRS(gcsc) ) { case 99912 : /* hack */ SetSysCoordSystemID_GCSRS(theSysCoord, 12); break; default : SetSysCoordSystemID_GCSRS(theSysCoord, GetSysCoordSystemID_GCSRS(gcsc)); break; } SetSysCoordTimeZone_GCSRS(theSysCoord, GetSysCoordTimeZone_GCSRS(gcsc)); if( GetSysCoordName_GCSRS(gcsc) ) SetSysCoordName_GCSRS(theSysCoord, GetSysCoordName_GCSRS(gcsc)); if( GetSysCoordUnit_GCSRS(gcsc) ) SetSysCoordUnit_GCSRS(theSysCoord, GetSysCoordUnit_GCSRS(gcsc)); } return theSysCoord; }/* _findSysCoord_GCSRS */ /* -------------------------------------------------------------------- */ static void GCSRSAPI_CALL _InitSysCoord_GCSRS ( GCSysCoord* theSysCoord ) { SetSysCoordSystemID_GCSRS(theSysCoord, -1); SetSysCoordTimeZone_GCSRS(theSysCoord, -1); SetSysCoordName_GCSRS(theSysCoord, NULL); SetSysCoordUnit_GCSRS(theSysCoord, NULL); SetSysCoordCentralMeridian_GCSRS(theSysCoord, 0.0); SetSysCoordLatitudeOfOrigin_GCSRS(theSysCoord, 0.0); SetSysCoordStandardParallel1_GCSRS(theSysCoord, 0.0); SetSysCoordStandardParallel2_GCSRS(theSysCoord, 0.0); SetSysCoordScaleFactor_GCSRS(theSysCoord, 0.0); SetSysCoordFalseEasting_GCSRS(theSysCoord, 0.0); SetSysCoordFalseNorthing_GCSRS(theSysCoord, 0.0); SetSysCoordDatumID_GCSRS(theSysCoord, -1); SetSysCoordProjID_GCSRS(theSysCoord, -1); SetSysCoordPrimeMeridian_GCSRS(theSysCoord, 0); }/* _InitSysCoord_GCSRS */ /* -------------------------------------------------------------------- */ GCSysCoord GCSRSAPI_CALL1(*) CreateSysCoord_GCSRS ( int srsid, int nTimezone ) { int iSysCoord; GCSysCoord* theSysCoord; const GCSysCoord* gcsc; if( !(theSysCoord= VSI_MALLOC_VERBOSE(sizeof(GCSysCoord))) ) { return NULL; } _InitSysCoord_GCSRS(theSysCoord); if( srsid>=0) { for( iSysCoord= 0, gcsc= &(gk_asSysCoordList[0]); GetSysCoordSystemID_GCSRS(gcsc)!=-1; iSysCoord++, gcsc= &(gk_asSysCoordList[iSysCoord]) ) { if( srsid==GetSysCoordSystemID_GCSRS(gcsc) ) { SetSysCoordSystemID_GCSRS(theSysCoord, srsid); SetSysCoordTimeZone_GCSRS(theSysCoord, nTimezone); if( GetSysCoordName_GCSRS(gcsc) ) SetSysCoordName_GCSRS(theSysCoord, GetSysCoordName_GCSRS(gcsc)); if( GetSysCoordUnit_GCSRS(gcsc) ) SetSysCoordUnit_GCSRS(theSysCoord, GetSysCoordUnit_GCSRS(gcsc)); SetSysCoordCentralMeridian_GCSRS(theSysCoord, GetSysCoordCentralMeridian_GCSRS(gcsc)); SetSysCoordLatitudeOfOrigin_GCSRS(theSysCoord, GetSysCoordLatitudeOfOrigin_GCSRS(gcsc)); SetSysCoordStandardParallel1_GCSRS(theSysCoord, GetSysCoordStandardParallel1_GCSRS(gcsc)); SetSysCoordStandardParallel2_GCSRS(theSysCoord, GetSysCoordStandardParallel2_GCSRS(gcsc)); SetSysCoordScaleFactor_GCSRS(theSysCoord, GetSysCoordScaleFactor_GCSRS(gcsc)); SetSysCoordFalseEasting_GCSRS(theSysCoord, GetSysCoordFalseEasting_GCSRS(gcsc)); SetSysCoordFalseNorthing_GCSRS(theSysCoord, GetSysCoordFalseNorthing_GCSRS(gcsc)); SetSysCoordDatumID_GCSRS(theSysCoord, GetSysCoordDatumID_GCSRS(gcsc)); SetSysCoordProjID_GCSRS(theSysCoord, GetSysCoordProjID_GCSRS(gcsc)); break; } } } return theSysCoord; }/* CreateSysCoord_GCSRS */ /* -------------------------------------------------------------------- */ static void GCSRSAPI_CALL _ReInitSysCoord_GCSRS ( GCSysCoord* theSysCoord ) { _InitSysCoord_GCSRS(theSysCoord); }/* _ReInitSysCoord_GCSRS */ /* -------------------------------------------------------------------- */ void GCSRSAPI_CALL DestroySysCoord_GCSRS ( GCSysCoord** theSysCoord ) { _ReInitSysCoord_GCSRS(*theSysCoord); CPLFree(*theSysCoord); *theSysCoord= NULL; }/* DestroySysCoord_GCSRS */ /* -------------------------------------------------------------------- */ GCSysCoord GCSRSAPI_CALL1(*) OGRSpatialReference2SysCoord_GCSRS ( OGRSpatialReferenceH poSR ) { char* pszProj4= NULL; const GCSpheroidInfo* ell= NULL; const GCDatumInfo* datum= NULL; const GCProjectionInfo* gcproj= NULL; double a, rf, f, p[7]; GCSysCoord* syscoord= NULL; if( !poSR ) return NULL; pszProj4= NULL; OSRExportToProj4(poSR, &pszProj4); if( !pszProj4 ) pszProj4= CPLStrdup(""); CPLDebug("GEOCONCEPT", "SRS : %s", pszProj4); if( !(syscoord= CreateSysCoord_GCSRS(-1,-1)) ) { goto onError; } SetSysCoordPrimeMeridian_GCSRS(syscoord, OSRGetPrimeMeridian(poSR,NULL)); a= OSRGetSemiMajor(poSR,NULL); rf= OSRGetInvFlattening(poSR,NULL); ell= _findSpheroid_GCSRS(a, rf); if( GetInfoSpheroidID_GCSRS(ell)==-1 ) { CPLDebug("GEOCONCEPT", "Unsupported ellipsoid : %.4f %.10f", a, rf); goto onError; } CPLDebug("GEOCONCEPT", "ellipsoid found : %s", GetInfoSpheroidName_GCSRS(ell)); OSRGetTOWGS84(poSR,p,7); f= 1.0 - sqrt(1.0 - GetInfoSpheroidExcentricity_GCSRS(ell)*GetInfoSpheroidExcentricity_GCSRS(ell)); datum= _findDatum_GCSRS(p[0], p[1], p[2], GetInfoSpheroidSemiMajor_GCSRS(ell), f); if( GetInfoDatumID_GCSRS(datum)==-1 ) { CPLDebug("GEOCONCEPT", "Unsupported datum : %.4f %.4f; %.4f a=%.4f rf=%.10f", p[0], p[1], p[2], a, rf); goto onError; } /* FIXME : WGS 84 and GRS 80 assimilation by Geoconcept : */ if( GetInfoSpheroidID_GCSRS(ell)==4 ) /* GRS 80 */ { datum= &(gk_asDatumList[31]); } else if( GetInfoSpheroidID_GCSRS(ell)==9999 ) /* WGS 84 */ { datum= &(gk_asDatumList[3]); } CPLDebug("GEOCONCEPT", "datum found : %s", GetInfoDatumName_GCSRS(datum)); SetSysCoordDatumID_GCSRS(syscoord, GetInfoDatumID_GCSRS(datum)); gcproj= _findProjection_GCSRS(OSRIsGeographic(poSR)? NULL : OSRGetAttrValue(poSR, "PROJECTION", 0), OSRGetProjParm(poSR,SRS_PP_PSEUDO_STD_PARALLEL_1,0.0,NULL)); if( GetInfoProjID_GCSRS(gcproj)==-1 ) { CPLDebug("GEOCONCEPT", "Unsupported projection : %s", OSRIsGeographic(poSR)? "GEOCS":OSRGetAttrValue(poSR, "PROJECTION", 0)); goto onError; } CPLDebug("GEOCONCEPT", "projection : %s", GetInfoProjName_GCSRS(gcproj)); SetSysCoordProjID_GCSRS(syscoord, GetInfoProjID_GCSRS(gcproj)); /* then overwrite them with projection specific parameters ... */ if( OSRIsProjected(poSR) ) { double v; SetSysCoordPrimeMeridian_GCSRS(syscoord, OSRGetPrimeMeridian(poSR,NULL)); SetSysCoordCentralMeridian_GCSRS(syscoord, OSRGetProjParm(poSR,SRS_PP_CENTRAL_MERIDIAN,0.0,NULL)); SetSysCoordLatitudeOfOrigin_GCSRS(syscoord, OSRGetProjParm(poSR,SRS_PP_LATITUDE_OF_ORIGIN,0.0,NULL)); SetSysCoordStandardParallel1_GCSRS(syscoord, OSRGetProjParm(poSR,SRS_PP_STANDARD_PARALLEL_1,0.0,NULL)); SetSysCoordStandardParallel2_GCSRS(syscoord, OSRGetProjParm(poSR,SRS_PP_STANDARD_PARALLEL_2,0.0,NULL)); SetSysCoordFalseEasting_GCSRS(syscoord, OSRGetProjParm(poSR,SRS_PP_FALSE_EASTING,0.0,NULL)); SetSysCoordFalseNorthing_GCSRS(syscoord, OSRGetProjParm(poSR,SRS_PP_FALSE_NORTHING,0.0,NULL)); if( (v= OSRGetProjParm(poSR,SRS_PP_SCALE_FACTOR,0.0,NULL))!= 0.0 ) { SetSysCoordScaleFactor_GCSRS(syscoord, v); } if( (v= OSRGetProjParm(poSR,SRS_PP_PSEUDO_STD_PARALLEL_1,0.0,NULL))!= 0.0 ) { /* should be SRS_PT_EQUIRECTANGULAR : */ SetSysCoordScaleFactor_GCSRS(syscoord, cos(v*M_PI/180.0)); SetSysCoordStandardParallel1_GCSRS(syscoord, v);/* allow keeping lat_ts sign */ } } /* Retrieve the syscoord : */ if( !_findSysCoord_GCSRS(syscoord) ) { CPLDebug("GEOCONCEPT", "invalid syscoord ?!"); goto onError; } if( GetSysCoordSystemID_GCSRS(syscoord)==-1 ) { CPLDebug("GEOCONCEPT", "Cannot find syscoord"); goto onError; } /* when SRS_PT_TRANSVERSE_MERCATOR, get zone : */ if( GetSysCoordTimeZone_GCSRS(syscoord)==0 ) { int pbNorth= 1; SetSysCoordTimeZone_GCSRS(syscoord, OSRGetUTMZone(poSR,&pbNorth)); } if( pszProj4 ) { CPLFree(pszProj4); } CPLDebug( "GEOCONCEPT", "SysCoord value: %d:%d", GetSysCoordSystemID_GCSRS(syscoord), GetSysCoordTimeZone_GCSRS(syscoord) ); return syscoord; onError: if( pszProj4 ) { CPLDebug( "GEOCONCEPT", "Unhandled spatial reference system '%s'.", pszProj4); CPLFree(pszProj4); } if( syscoord ) { DestroySysCoord_GCSRS(&syscoord); } return NULL; }/* OGRSpatialReference2SysCoord_GCSRS */ /* -------------------------------------------------------------------- */ OGRSpatialReferenceH GCSRSAPI_CALL SysCoord2OGRSpatialReference_GCSRS ( GCSysCoord* syscoord ) { OGRSpatialReferenceH poSR; const GCDatumInfo* datum= NULL; const GCSpheroidInfo* ell= NULL; int i; double f; poSR= OSRNewSpatialReference(NULL); if( syscoord && GetSysCoordSystemID_GCSRS(syscoord)!=-1 ) { switch( GetSysCoordProjID_GCSRS(syscoord) ) { case 0 : /* long/lat */ break; case 1: /* UTM */ case 11: /* MGRS */ case 12: /* TM */ OSRSetTM(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordScaleFactor_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 2: /* LCC 1SP */ OSRSetLCC1SP(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordScaleFactor_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 3: /* Bonne */ OSRSetBonne(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 4: /* Plate Caree */ OSRSetEquirectangular(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 18: /* LCC 2SP */ OSRSetLCC(poSR, GetSysCoordStandardParallel1_GCSRS(syscoord), GetSysCoordStandardParallel2_GCSRS(syscoord), GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 19: /* Gauss Schreiber : Reunion */ OSRSetGaussSchreiberTMercator(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordScaleFactor_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 20: /* Polyconic */ OSRSetPolyconic(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 21: /* Direct Mercator */ OSRSetMercator(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordScaleFactor_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 22: /* Stereographic oblic */ OSRSetOS(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordScaleFactor_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 24: /* Miller */ OSRSetMC(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; case 26: /* Equi rectangular */ OSRSetEquirectangular2(poSR, GetSysCoordLatitudeOfOrigin_GCSRS(syscoord), GetSysCoordCentralMeridian_GCSRS(syscoord), GetSysCoordStandardParallel1_GCSRS(syscoord), GetSysCoordFalseEasting_GCSRS(syscoord), GetSysCoordFalseNorthing_GCSRS(syscoord)); break; default : break; } if( GetSysCoordProjID_GCSRS(syscoord)>0 ) OSRSetProjCS(poSR, GetSysCoordName_GCSRS(syscoord)); for( i= 0, datum= &(gk_asDatumList[0]); GetInfoDatumID_GCSRS(datum)!=-1; i++, datum= &(gk_asDatumList[i]) ) { if( GetInfoDatumID_GCSRS(datum)==GetSysCoordDatumID_GCSRS(syscoord) ) break; } for( i= 0, ell= &(gk_asSpheroidList[0]); GetInfoSpheroidID_GCSRS(ell)!=-1; i++, ell= &(gk_asSpheroidList[i]) ) { if( _areCompatibleSpheroids_GCSRS(GetInfoSpheroidID_GCSRS(ell), GetInfoDatumSpheroidID_GCSRS(datum)) ) break; } /* FIXME : WGS 84 and GRS 80 assimilation by Geoconcept : */ if( GetInfoDatumID_GCSRS(datum)==4 ) /* WGS 84 */ { ell= &(gk_asSpheroidList[8]); } else if( GetInfoDatumID_GCSRS(datum)==9984 ) /* GRS 80 */ { ell= &(gk_asSpheroidList[3]); } f= 1.0 - sqrt(1.0 - GetInfoSpheroidExcentricity_GCSRS(ell)*GetInfoSpheroidExcentricity_GCSRS(ell)); OSRSetGeogCS(poSR, GetSysCoordProjID_GCSRS(syscoord)!=0 || !GetSysCoordName_GCSRS(syscoord)? "unnamed":GetSysCoordName_GCSRS(syscoord), GetInfoDatumID_GCSRS(datum)>=0? GetInfoDatumName_GCSRS(datum):"unknown", GetInfoSpheroidID_GCSRS(ell)>=0? GetInfoSpheroidName_GCSRS(ell):"unknown", GetInfoSpheroidID_GCSRS(ell)>=0? GetInfoSpheroidSemiMajor_GCSRS(ell):6378137.0, GetInfoSpheroidID_GCSRS(ell)>=0? (f==0? 0:1/f):298.257223563, "Greenwich", GetSysCoordPrimeMeridian_GCSRS(syscoord), SRS_UA_DEGREE, CPLAtof(SRS_UA_DEGREE_CONV)); /* As Geoconcept uses Molodensky, we've got only 3 out of 7 params for Bursa-Wolf : */ /* the 4 missing Bursa-Wolf parameters have been added to the gk_asDatumList ! */ if( GetInfoProjID_GCSRS(syscoord)>0 && GetInfoDatumID_GCSRS(datum)!=-1 ) { OSRSetTOWGS84(poSR, GetInfoDatumShiftX_GCSRS(datum), GetInfoDatumShiftY_GCSRS(datum), GetInfoDatumShiftZ_GCSRS(datum), GetInfoDatumRotationX_GCSRS(datum), GetInfoDatumRotationY_GCSRS(datum), GetInfoDatumRotationZ_GCSRS(datum), 1e6*GetInfoDatumScaleFactor_GCSRS(datum)); } } /* -------------------------------------------------------------------- */ /* Report on translation. */ /* -------------------------------------------------------------------- */ { char* pszWKT; OSRExportToWkt(poSR,&pszWKT); if( pszWKT!=NULL ) { CPLDebug( "GEOCONCEPT", "This SysCoord value: %d:%d was translated to : %s", syscoord ? GetSysCoordSystemID_GCSRS(syscoord) : -1, syscoord ? GetSysCoordTimeZone_GCSRS(syscoord) : -1, pszWKT ); CPLFree( pszWKT ); } } return poSR; }/* SysCoord2OGRSpatialReference_GCSRS */