/********************************************************************** * $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 * * SPDX-License-Identifier: MIT **********************************************************************/ #include "geoconcept_syscoord.h" #include "cpl_string.h" /* -------------------------------------------------------------------- */ /* 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] = {0, 0, 0, 0, 0, 0, 0}; 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); OSRSetAxisMappingStrategy(poSR, OAMS_TRADITIONAL_GIS_ORDER); 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 */