если сделать через минимум среднеквадратичного значения посчитанного сверху/снизу:
Код
horStDevData = Table[size[[2]] - Ordering[Table[StandardDeviation[Flatten[data[[1;; y + 1, x, col]]]]^2 + StandardDeviation[Flatten[data[[y;; size[[2]], x, col]]]]^2, {y, 1, size[[2]] - 1}], 1][[1]], {x, size[[1]]}];
то начинает цепляться за более резкие переходы (светлозеленые точки вдоль дороги, зелёный - фит):

можно наверное совместить оба метода, и считать среднеквадратичное с поправкой на цвет неба, тогда наверное совсем хорошо получится.
upd: тут косяк, только по одному цветовому каналу считается, исправлять уже лень.