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

Improve default nearest resampling radius of influence #527

Open
djhoese opened this issue Sep 16, 2022 · 0 comments
Open

Improve default nearest resampling radius of influence #527

djhoese opened this issue Sep 16, 2022 · 0 comments
Assignees
Labels
component:remap Remapping or resampling enhancement

Comments

@djhoese
Copy link
Member

djhoese commented Sep 16, 2022

See #521. The KDTree-based resampler in Satpy is powered by pykdtree. If the user doesn't force the value for radius_of_influence (the geocentric distance to look for valid pixels until giving up and putting in a fill value), then pyresample checks the "geocentric resolution" estimates of the input area and output area and takes the maximum between the two. In the AMSR2 high resolution band case that means the 5000m hard-coded nadir resolution is multiplied by 2 to get an estimate for the resolution at the limb of the swath. For most high resolution output grids this value of 10000m will always be larger and therefore be the place where work should be done.

In pyresample/geometry.py there is a SwathDefinition.geocentric_resolution method. This method should be updated to check for a limb_resolution attribute and use that to make a better estimate of the limb resolution. Either the maximum of the nadir and limb or an average (could cause holes). The amsr2_l1b reader in Satpy should then also be updated to include a limb_resolution attribute for the updated pyresample code to use.

@djhoese djhoese added enhancement component:remap Remapping or resampling labels Sep 16, 2022
@djhoese djhoese added this to the Polar2Grid 3.1 milestone Sep 16, 2022
@djhoese djhoese self-assigned this Sep 16, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
component:remap Remapping or resampling enhancement
Projects
None yet
Development

No branches or pull requests

1 participant