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

[Question] Resample of mesoscale data gives blank data #2749

Closed
agp-earth opened this issue Feb 22, 2024 · 1 comment
Closed

[Question] Resample of mesoscale data gives blank data #2749

agp-earth opened this issue Feb 22, 2024 · 1 comment

Comments

@agp-earth
Copy link

agp-earth commented Feb 22, 2024

Hi! I have a quick question! I am resampling a GOES 17 mesoscale L1B radiance image with an area I defined. I don't get any errors (so it does not have to do with coordinate overlap issues). However, when I check the data of the resampled image I have NaN values and plotting it gives me a blank image. I know it is reading the original file fine as that gives a nice image when plotting (plt.imshow(scn['C07']).

scn = Scene(filenames=glob.glob(myfile), reader='abi_l1b')
scn.load(['C07'],calibration='radiance')
projection = '+proj=utm +zone=3 +ellps=WGS84 +datum=WGS84 +units=m +no_defs'
proj_dict = utils.proj4_str_to_dict(projection)
area_def = geometry.AreaDefinition(area_id, description, proj_id, proj_dict, x_size, y_size, area_extent)
cropscn = scn.resample(my_area)
cropdata = cropscn['C07'].data
nan_check = np.isnan(np.array(cropdata.max()))
@djhoese
Copy link
Member

djhoese commented Feb 27, 2024

While I agree this question should be in the satpy repository, I didn't see this until I had already started responding in pytroll/pyresample#587. Closing as a duplicate of that one.

@djhoese djhoese closed this as completed Feb 27, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants