warpimage.py and Mollweide projection

Hi,

Using the basemap toolkit, I try to warp an input cylindrical projection image to a Mollweide projection, following the basemap example script warpimage.py.

Unfortunately, it seems there's a bug in the masking of points outside the elliptical limits of the Mollweide projection. This is illustrated on the attached plot, generated by warpimage.py hacked to use a 'moll' projection (instead of 'ortho' which works fine).

For your information, I use numpy 1.1.0, matplotlib 0.91.2 and basemap 0.9.9.

Cheers.

warp2moll.png

···

--
    .~. Yannick COPIN (o:>* Doctus cum libro
    /V\ Institut de physique nucleaire de Lyon (IN2P3 - France)
   // \\ Tel: (33/0) 472 431 968 AIM: YnCopin ICQ: 236931013
  /( )\ http://snovae.in2p3.fr/ycopin/
   ^`~'^