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

Cannot read AVHRR in AAPP format #179

Closed
howff opened this issue Feb 1, 2018 · 8 comments
Closed

Cannot read AVHRR in AAPP format #179

howff opened this issue Feb 1, 2018 · 8 comments
Labels

Comments

@howff
Copy link

howff commented Feb 1, 2018

satpy 0.8.1 pyresample 1.7.1

Trying to load and resample AVHRR in AAPP format I get an image like the attached - the projection is wrong and there are vertical stripes (probably both the same cause).

2018_1_20_1524_unknown_4_s1_grid

@djhoese
Copy link
Member

djhoese commented Feb 1, 2018

Could you please provide the commands you are using to run this? Which satellite is this data from? @adybbroe Aren't there some AVHRRs that have degraded data quality?

Edit: Next time, please follow the github issue template as much as possible and when it makes sense. It exists specifically so I don't have to ask the above questions.

@adybbroe
Copy link
Contributor

adybbroe commented Feb 1, 2018

@davidh-ssec Yes, NOAA-15 in particular, but this manifests through drop out of scan lines. This I believe is something else

@howff
Copy link
Author

howff commented Feb 1, 2018

global_scene = Scene(reader='aapp_l1b', filenames=['hrpt_noaa19_20180120_1524_46109.l1b'])
global_scene.load(['4'])
local_scene = global_scene.resample('euro4')
local_scene.save_dataset('4', filename='ch4.png')

@adybbroe
Copy link
Contributor

adybbroe commented Feb 1, 2018

@howff I can't reproduce your issue.
I did this:

import os
from satpy import Scene
from satpy.utils import debug_on
debug_on()

BASEDIR = "/home/a000680/data/pps/import/PPS_data/source/metop01_20170710_0830_24957"
FILES = os.path.join(BASEDIR, "hrpt_metop01_20170710_0830_24957.l1b")

if __name__ == '__main__':

    global_scene = Scene(reader='aapp_l1b', filenames=[FILES])
    global_scene.load(['4'])
    local_scene = global_scene.resample('euro4')
    local_scene.save_dataset('4', filename='ch4.png',
                             overlay={'coast_dir': '/home/a000680/data/shapes/'})

Using pyresample master branch (version 1.7.1)
Pyproj 1.9.5.1
Numpy 1.13.1
Satpy master 0.8.1
pykdtree 1.2.1

Could you check your versions of the above?

@howff
Copy link
Author

howff commented Feb 2, 2018

I've got
pyproj - same
numpy - newer 1.14.0
satpy - same
pykdtree newer 1.2.2

Data is https://wetransfer.com/downloads/c7bf44ab69ef40f8283ca9ac96d5f14a20180202090725/bede6aa5aa5eb1e2637f2c73620d277520180202090725/f315fb

@adybbroe
Copy link
Contributor

adybbroe commented Feb 2, 2018

@howff Ok. Where do you get the pykdtree version from. The latest version is 1.2.1 I thought!?
Numpy 1.14.0 is buggy. There is a really annoying bug concerning the writing of npz files, which pyspectral use for the LUTs.

@howff
Copy link
Author

howff commented Feb 2, 2018

Argh, I don't believe it, the image viewer has a bug displaying certain PNGs!!
Sorry for wasting time.
BTW, pykdtree 1.2.2 is from pypi.
Just before closing, is the bug in numpy 1.14.0 a deal-breaker?

@mraspaud
Copy link
Member

mraspaud commented Feb 2, 2018

@howff if you use the savez function, then yes it is, otherwise it should be fine I guess

@howff howff closed this as completed Feb 2, 2018
@djhoese djhoese added the bug label Feb 19, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

4 participants