Revision: 6722
          http://matplotlib.svn.sourceforge.net/matplotlib/?rev=6722&view=rev
Author:   jswhit
Date:     2008-12-31 17:00:34 +0000 (Wed, 31 Dec 2008)

Log Message:
-----------
add some comments

Modified Paths:
--------------
    trunk/toolkits/basemap/lib/mpl_toolkits/basemap/proj.py

Modified: trunk/toolkits/basemap/lib/mpl_toolkits/basemap/proj.py
===================================================================
--- trunk/toolkits/basemap/lib/mpl_toolkits/basemap/proj.py     2008-12-31 
16:29:37 UTC (rev 6721)
+++ trunk/toolkits/basemap/lib/mpl_toolkits/basemap/proj.py     2008-12-31 
17:00:34 UTC (rev 6722)
@@ -96,13 +96,17 @@
              urcrnrlat == 90):
             self._fulldisk = True
             self._proj4 = pyproj.Proj(projparams)
+            # raise an exception for ellipsoids - there appears to be a bug
+            # in proj4 that causes the inverse transform to fail for points
+            # more than 90 degrees of arc away from center point for ellipsoids
+            # (works fine for spheres)
             if self.ellipsoid:
                 msg = dedent("""
                 full disk (whole world) Azimuthal Equidistant projection can
                 only be drawn for a perfect sphere""")
                 raise ValueError(msg)
-            llcrnrx = -np.pi*self.rminor
-            llcrnry = -np.pi*self.rminor
+            llcrnrx = -np.pi*self.rmajor
+            llcrnry = -np.pi*self.rmajor
             self._width = -llcrnrx
             self._height = -llcrnry
             urcrnrx = -llcrnrx


This was sent by the SourceForge.net collaborative development platform, the 
world's largest Open Source development site.

------------------------------------------------------------------------------
_______________________________________________
Matplotlib-checkins mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/matplotlib-checkins

Reply via email to