Fix special handling of Azimuth parameter of Krovak

2 files changed, 4 insertions, 2 deletions

diff --git a/src/coordinateoperation.cpp b/src/coordinateoperation.cpp
index d187f2cc..ced0ab9f 100644
--- a/src/coordinateoperation.cpp
+++ b/src/coordinateoperation.cpp
@@ -5177,7 +5177,9 @@ void Conversion::_exportToPROJString(
         double latitudePseudoStandardParallel = parameterValueNumeric(
             EPSG_CODE_PARAMETER_LATITUDE_PSEUDO_STANDARD_PARALLEL,
             common::UnitOfMeasure::DEGREE);
-        if (std::fabs(colatitude - 30.28813972222222) > 1e-8) {
+        // 30deg 17' 17.30311'' = 30.28813975277777776
+        // 30deg 17' 17.303''   = 30.288139722222223 as used in GDAL WKT1
+        if (std::fabs(colatitude - 30.2881397) > 1e-7) {
             throw io::FormattingException(
                 std::string("Unsupported value for ") +
                 EPSG_NAME_PARAMETER_COLATITUDE_CONE_AXIS);
diff --git a/src/io.cpp b/src/io.cpp
index f396f1df..ec35bdc5 100644
--- a/src/io.cpp
+++ b/src/io.cpp
@@ -6674,7 +6674,7 @@ CRSNNPtr PROJStringParser::Private::buildProjectedCRS(
                 } else if (step.name == "krovak") {
                     if (param->epsg_code ==
                         EPSG_CODE_PARAMETER_COLATITUDE_CONE_AXIS) {
-                        value = 30.2881397222222;
+                        value = 30.28813975277777776;
                     } else if (
                         param->epsg_code ==
                         EPSG_CODE_PARAMETER_LATITUDE_PSEUDO_STANDARD_PARALLEL) {