I am working with (geostationary)
satellite data, and one of the things I want to do is plot a map (coastlines)
on top of the satellite image. As an IDL user I know how to do this in IDL (although
in IDL it is a bit of a hassle), but I don’t seem to be able to figure
out how this could be done in python.
The satellite data consists
of a rectangular field of N by N pixels. The data is already converted to
bitmap RGB values. I have figured out how this can be written to a bitmap image
like PNG or JPG.
I further figured out how to
do the satellite projection in Basemap, and how to plot Basemap coastlines. However,
the Basemap appears to have a non-transparent background, which overplots all
bitmap data if I plot the Bitmap data first. I had hoped that there would be
some transparency setting in Basemap, but alas.
(ps. Keep in mind that I do not want to use some
contour filling routine for plotting the satellite data. I want to retain the
original N x N pixels and image size)