From: Angus McMorland <amcmorl@gm...>  20090727 22:23:21

Hi all, 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 3dpolar to 2dprojectedpolar coordinates, 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 coordinate 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. Thanks, Angus.  AJC McMorland Postdoctoral research fellow Neurobiology, University of Pittsburgh 