Goodman, Alexander (398J-Affiliate)
June 24, 2013 6:27 PM
for the quick reply. However I am a little lost with your suggested solution. The issue here is that if datain is a masked array, dataout will not be a masked array when order=3 is used unless masked=True is set, but this just results in none of the elements (besides those outside the domain of xin and yin) being masked. So are you saying to try setting the mask to what it would be after running interp() with order=0? I would appreciate a more concrete example.
After pondering this a bit more, I don’t think you can use order=3 with
masked arrays. Any reason why you can’t using a combo of bilinear and nearest neighbor as I described in my previous email?