... to nodata, NaNs, "folded" neighbours (when left and right are a the same
position...)

To avoid clashes with otbIteratorHelpers from OTB

Departs from output parameters to use the more concise `y = f(x)` form.