Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Bilinear interpolation doesn't work with StackedAreaDefinitions #735

Closed
pnuu opened this issue Apr 26, 2019 · 0 comments · Fixed by pytroll/pyresample#208
Closed

Bilinear interpolation doesn't work with StackedAreaDefinitions #735

pnuu opened this issue Apr 26, 2019 · 0 comments · Fixed by pytroll/pyresample#208

Comments

@pnuu
Copy link
Member

pnuu commented Apr 26, 2019

Describe the bug
Bilinear interpolation doesn't work when the area of the input data is StackedAreaDefinition. An example of this case is MSG/SEVIRI with a missing segment from the middle of the disk, or the HRV channel when there are segments available for both of the imaging windows.

To Reproduce

import glob
from satpy import Scene

fnames = glob.glob("H*__")
glbl = Scene(reader='seviri_l1b_hrit', filenames=fnames)
glbl.load(['HRV'])
lcl = glbl.resample('euro4', resampler='bilinear')
lcl.save_datasets(base_dir='/tmp')

Expected behavior
Image of HRV channel resampled to euro4 area.

Actual results
Empty image of euro4 area.

Environment Info:

  • Any recent (xarray/dask based) Satpy
  • PyResample Version: 0.12.0 (also earlier versions seem to have the same issue)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Development

Successfully merging a pull request may close this issue.

1 participant