Hi, the pyproj package seems to cause a problem in the polarmap example of the basemap toolkit.
Thanks,
david
huardda@…559…:~/svnrepos/toolkits/basemap/examples$ python polarmaps.py
min/max etopo20 data:
-9026.625 6228.8125
plotting North Polar Lambert Azimuthal Equal Area basemap …
plotting North Polar Stereographic basemap …
plotting North Polar Azimuthal Equidistant basemap …
Traceback (most recent call last):
File “polarmaps.py”, line 51, in ?
resolution=‘c’,area_thresh=10000.,lat_0=lat_0,lon_0=lon_0_ortho)
File “/usr/local/lib/python2.4/site-packages/matplotlib/toolkits/basemap/basemap.py”, line 1111, in init
az1,alpha21,dist=gc.inv(lon_0,lat_0,math.radians(lonsnew[0]),math.radians(latsnew[0]),radians=True)
File “/usr/local/lib/python2.4/site-packages/matplotlib/toolkits/basemap/pyproj.py”, line 478, in inv
_Geod._inv(self, inx, iny, inz, ind, radians=radians)
File “_geod.pyx”, line 123, in _geod.Geod._inv
ValueError: undefined inverse geodesic (may be an antipodal point)
This is from a fresh SVN version of both matplotlib and basemap. Linux, Ubuntu edgy, Xeon-64.