diff --git a/rtengine/histmatching.cc b/rtengine/histmatching.cc
index 0faa8c357..5db504816 100644
--- a/rtengine/histmatching.cc
+++ b/rtengine/histmatching.cc
@@ -100,7 +100,6 @@ void mappingToCurve(const std::vector<int> &mapping, std::vector<double> &curve)
 {
     curve.clear();
     
-    const int npoints = 8;
     int idx = 15;
     for (; idx < int(mapping.size()); ++idx) {
         if (mapping[idx] >= idx) {
@@ -114,7 +113,6 @@ void mappingToCurve(const std::vector<int> &mapping, std::vector<double> &curve)
             }
         }
     }
-    int step = std::max(int(mapping.size())/npoints, 1);
 
     auto coord = [](int v) -> double { return double(v)/255.0; };
     auto doit =
@@ -147,9 +145,17 @@ void mappingToCurve(const std::vector<int> &mapping, std::vector<double> &curve)
     while (start < idx && (mapping[start] < 0 || start < idx / 2)) {
         ++start;
     }
-    
-    doit(start, idx, idx > step ? step : idx / 2, true);
-    doit(idx, int(mapping.size()), step, idx - step > step / 2 && std::abs(curve[curve.size()-2] - coord(idx)) > 0.01);
+
+    const int npoints = 8;
+    int step = std::max(int(mapping.size())/npoints, 1);
+    int end = mapping.size();
+    if (idx <= end / 3) {
+        doit(start, idx, idx / 2, true);
+        doit(idx, end, (end - idx) / 3, false);
+    } else {
+        doit(start, idx, idx > step ? step : idx / 2, true);
+        doit(idx, int(mapping.size()), step, idx - step > step / 2 && std::abs(curve[curve.size()-2] - coord(idx)) > 0.01);
+    }
     
     if (curve.size() > 2 && (1 - curve[curve.size()-2] <= step / (256.0 * 3))) {
         curve.pop_back();