I've recently implemented a simple equal area (Lambert) projection in
mpl, using two polar Axes (to draw each hemisphere), and this
currently works fine to draw lines and points, since I'm transforming
the points in my own code, from 3d-polar to 2d-projected-polar
co-ordinates, and then plotting those.
Now I have a need to also be able to draw coloured intensity maps on
the Axes using the same transform, and I'm wondering what the best way
to go about this would be. When I started the first implementation, I
looked at basemap and its transforms, but they didn't seem to do what
I wanted any more easily than just implementing it myself. I'm also
vaguely aware of some transformation capabilities (that are relatively
new?) in mpl itself, but I'm not sure whether they're useful to me for
polar co-ordinate systems and plot contents, rather than plot
decorations. Ideally, the transform would apply to the interpolation
too, but I think it would be sufficient to linearly interpolate
between the points: I can simply provide them at high enough density
to minimize the discrepancy.
Any suggestions on the best way to do this would be gratefully received.
Post-doctoral research fellow
Neurobiology, University of Pittsburgh