pytroll / satpy

Python package for earth-observing satellite data processing
http://satpy.readthedocs.org/en/latest/
GNU General Public License v3.0
1.07k stars 295 forks source link

create rgb air mass by meteosat #1825

Open gisgeo3 opened 3 years ago

gisgeo3 commented 3 years ago

Describe the bug Different colors appear to me

A clear and concise description of what the bug is. Dear, I created an rgb air mass using arcgis and the colors appeared to me very different, I am a PhD student who is not specialized in programming and the python language and I do not know how to create an rgb air mass inside a satpy python environment. We seek knowledge and knowledge from each other. So please help me solve this problem. If I don't know how, my thesis topic will be cancelled.

To Reproduce

# Your code here
I haven't code.

Expected behavior A clear and concise description of what you expected to happen.

Actual results Text output of actual results or error messages including full tracebacks if applicable.

Screenshots If applicable, add screenshots to help explain your problem.

Environment Info:

Additional context Add any other context about the problem here.

djhoese commented 3 years ago

Have you read the Satpy documentation: https://satpy.readthedocs.io/en/stable/

There are a couple examples for loading different products but switching to airmass should be a small change: https://satpy.readthedocs.io/en/stable/examples/index.html

There's also my tutorial which should help introduce you to Satpy: https://github.com/pytroll/tutorial-satpy-half-day

Your basic code should look like:

from satpy import Scene
from glob import glob

scn = Scene(reader='seviri_l1b_hrit', filenames=glob('/path/to/your/files/*'))
scn.load(['airmass'])
scn.save_datasets(base_dir='/path/to/output/dir')

It is possible you will need resampling, but I don't work with SEVIRI data very often so I've forgotten. We can fix that if you hit any issues with this. If your SEVIRI data is not in the HRIT format than you will need to use a different reader (see the table of readers at the bottom of the Satpy docs home page).

gisgeo3 commented 3 years ago

Can I contact you personally.