Filtering don't preserve value in the edge

Q&A
9

This could be done with two grdclip operations but we would need a another grdmath operator. One that would do A if A != NaN else B.
Trivial to implement in Julia/Matlab.