Fix wrong EPSG conversion code for UTM south

When generating a Conversion object for a UTM south conversion (for example when importing from PROJ.4 or WKT1), we currently use a wrong offset. Fixes https://gis.stackexchange.com/questions/482037/what-is-epsg17056
rouault · Jun 5, 2024 · 355a7c7 · 355a7c7
1 parent e0e1482
commit 355a7c7
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/src/iso19111/operation/conversion.cpp b/src/iso19111/operation/conversion.cpp
@@ -251,7 +251,7 @@ getUTMConversionProperty(const util::PropertyMap &properties, int zone,
         conversionName += (north ? 'N' : 'S');
 
         return createMapNameEPSGCode(conversionName,
-                                     (north ? 16000 : 17000) + zone);
+                                     (north ? 16000 : 16100) + zone);
     } else {
         return properties;
     }
@@ -1760,7 +1760,8 @@ ConversionNNPtr Conversion::createPopularVisualisationPseudoMercator(
  * sphere at centerLat.
  *
  * This method is defined as
- * <a href="https://epsg.org/coord-operation-method_1026/Mercator-Spherical.html">
+ * <a
+ * href="https://epsg.org/coord-operation-method_1026/Mercator-Spherical.html">
  * EPSG:1026</a>.
  *
  * @param properties See \ref general_properties of the conversion. If the name

diff --git a/test/unit/test_operation.cpp b/test/unit/test_operation.cpp
@@ -1241,7 +1241,7 @@ TEST(operation, utm_export) {
               "    PARAMETER[\"False northing\",10000000,\n"
               "        LENGTHUNIT[\"metre\",1],\n"
               "        ID[\"EPSG\",8807]],\n"
-              "    ID[\"EPSG\",17001]]");
+              "    ID[\"EPSG\",16101]]");
 
     EXPECT_EQ(
         conv->exportToWKT(