H60_howto

Project information

The objective is to create contour or scatter maps from the H60 products and to clip the original dataset from a geographical extent.

H60 products are delivered in the original geostationary coordinate frame and the corresponding longitude/latitude coordinates are delivered in a separated file.

Note that imshow is the most efficient way to display the map and it does not require the corresponding longitude/latitude grid, just the correct CRS.

But, the final objective is to collocate the precipitation products with third party missions. It is thus necessary to be able to compute the corresponding (x,y) and (longitude,latitude) grid.

There is an additional difficulty here since the AreaDefinition (in the sense of pyresample) has to be manually defined in order to correctly display the dataset using the global attributes defining the projection.

Python xarray satpy pyresemple cartopy