edge cases in matlab handles nan, and in python i have to cast them to 0

[iancharest]

this creates edge differences

[kendrickkay]

To elaborate, the issue is that the Matlab implementation explicitly sets to "badval" any location that is outside of valid locations. The Python implementation may be more forgiving than this