Codecov Report

✅ All modified and coverable lines are covered by tests.
✅ Project coverage is 72.2%. Comparing base (556c088) to head (6a91e9a).
⚠️ Report is 25 commits behind head on develop.

@@             Coverage Diff             @@
##           develop    #2565      +/-   ##
===========================================
+ Coverage     55.5%    72.2%   +16.7%     
===========================================
  Files          644      642       -2     
  Lines       124135   122025    -2110     
===========================================
+ Hits         68947    88179   +19232     
+ Misses       55188    33846   -21342     

I see different results with method="cubic", since this method handles the internal masked option differently. Perhaps only apply the scale option (only on L533) when method in ("linear", "nearest")?

@mwtoews did you mean, keep the second rescale on 583 unchanged?

@mwtoews did you mean, keep the second rescale on 583 unchanged?

Probably fine with rescale=True, but I'm not familiar with the extrapolate_edges option. And I should point out that I didn't see specific test coverage for it.