Questions tagged [satpy]

Satpy is a GPL-licensed Python library for processing Earth observation satellite data (see also ), with a major focus on data from operational meteorological satellite instruments such as Meteosat SEVIRI, GOES ABI, HIMAWARI AHI, AVHRR, VIIRS, or many others. It is part of the wider pytroll framework that also includes , trollimage, pyorbital, and many other packages.

11 questions
2
votes
1 answer

SatPy load all channels at once

I have a scene object, I would like to load all channels into a numpy array of shape (24,24,3). Where 3 is the number of channels. scene_xybox = scn.crop(xy_bbox=box) I have to select each channel: channel= scene_xybox['VIS006'].values repeat, and…
1
vote
1 answer

Using satpy, Blending multiple satellite image

I want blend multiple satellite images. but an error occured. I followed the example on the satpy document. here is code and netcdf file is here : https://drive.google.com/drive/folders/1zp6EBVfjuh41LDRRZo4PJoeGGGn13AKy?usp=sharing from glob import…
singsung
  • 27
  • 3
1
vote
1 answer

Why do I using satpy on Himawari-8 standard data failed?

When I use satpy to read Himawari-8 standard datas, the terminal always says " 'filenames' was provided but is empty. " The Himawari-8 has 16 observation bands, and it take one full disk picture every 10 minutes so for one folder of observation…
1
vote
0 answers

Plotting Sentinel-3 Scene with cartopy

I'm following along the Pyresample example for plotting a Sentinel-3 Scene along with coastline. However, it seems as if pyresample.kd_tree.resample_nearest and pyresample.save_quicklook expect a 2D array as data, rather than the 3D (3-band)…
sluque
  • 495
  • 4
  • 12
1
vote
1 answer

Spatial merge/combine of multiple netcdf with xarray or satpy

I have two spatial dataset in netcdf format. They have same time, dimensions, coordinates, and data variable. But they are for different spatial coordinates. In below I tried to show my two dataset by a polygon: import glob import xarray as xr…
HMadadi
  • 391
  • 5
  • 22
0
votes
1 answer

Resample and plot satpy scene data from geos to mercator using pyresample

I have a scene object from eumetsat satellite. My area of interest: (xmin_corner, ymin_corner, xmax_corner, ymax_corner) [40,-10,105,40] degrees looks like this area definition looks like this: Area ID: msg_seviri_unknown_3km Description: MSG…
Rex
  • 117
  • 8
0
votes
1 answer

How to change the shape of a custom color bar using matplotlib?

My data has integer values mapped to categorical variables and I created a custom color bar to reflect these variables. When I've done this previously, I needed to add an extra color in my list and add the 0 as an integer to balance out the small…
0
votes
0 answers

Remove black values in a modis merge image

the following image shows my actuall result of a merged modis image. I used pyhdf.SD to open and manipulate and plot the data in a loop. It started on the right upper corner. The plot is ok... put on some areas overlaped. The easiest way is to…
S.Kociok
  • 168
  • 1
  • 14
0
votes
1 answer

Could not load HDF Modis data with satpy

one maybe silly question. I want to load MODIS hdf data with python. The datasets are available at Earthdata. First of all a generated a list of data, which match the…
S.Kociok
  • 168
  • 1
  • 14
0
votes
1 answer

Merging GOES17, EUMETSAT, GK-2A Meteorological satellite image

I want to generate global weather satellite image using GOES17, EUMETSAT, and GK-2A. I want make it Plate carree coordinate. (GOES 17 netcdf file convert to Plate Carree) First, using Satpy, I made plate carree image. from satpy import Scene from…
singsung
  • 27
  • 3
0
votes
1 answer

GOES 17 netcdf file convert to Plate Carree

Sample file can download from https://noaa-goes17.s3.amazonaws.com/ABI-L1b-RadF/2021/213/00/OR_ABI-L1b-RadF-M6C13_G17_s20212130000319_e20212130009396_c20212130009445.nc I am trying to convert a projection to a plate carree. I loaded netcdf using…
singsung
  • 27
  • 3