gdalwarp -crop_to_cutline: relax epsilon to avoid extremely long loop (…

…fixes #4826)
OSGeo · Nov 14, 2021 · 784a426 · 784a426
1 parent 3c11b98
commit 784a426
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/apps/gdalwarp_lib.cpp b/apps/gdalwarp_lib.cpp
@@ -451,7 +451,9 @@ static CPLErr CropToCutline( OGRGeometryH hCutline, char** papszTO,
         if( hCTCutlineToSrc != nullptr )
             OGR_G_Transform( hGeomInSrcSRS, hCTCutlineToSrc );
 
-        const double epsilon = std::numeric_limits<double>::epsilon();
+        // Do not use a smaller epsilon, otherwise it could cause useless
+        // segmentization (https://github.com/OSGeo/gdal/issues/4826)
+        constexpr double epsilon = 1e-10;
         for(int nIter=0;nIter<10;nIter++)
         {
             OGR_G_DestroyGeometry(hTransformedGeom);