SF.net SVN: matplotlib:[8808] branches/v1_0_maint

[<wea...@us...>]

Revision: 8808
          http://matplotlib.svn.sourceforge.net/matplotlib/?rev=8808&view=rev
Author:   weathergod
Date:     2010-11-22 16:05:17 +0000 (Mon, 22 Nov 2010)

Log Message:
-----------
Fixed a mistake with the Hammer projection in the geo module.
Thanks to Matthias Plum for reporting and Tobias Winchen for confirming.

Modified Paths:
--------------
    branches/v1_0_maint/CHANGELOG
    branches/v1_0_maint/lib/matplotlib/projections/geo.py

Modified: branches/v1_0_maint/CHANGELOG
===================================================================
--- branches/v1_0_maint/CHANGELOG	2010-11-22 14:01:21 UTC (rev 8807)
+++ branches/v1_0_maint/CHANGELOG	2010-11-22 16:05:17 UTC (rev 8808)
@@ -1,3 +1,5 @@
+2010-11-22 Fixed error with Hammer projection. - BVR
+
 2010-11-12 Fixed the placement and angle of axis labels in 3D plots. - BVR
 
 2010-11-07 New rc parameters examples.download and examples.directory

Modified: branches/v1_0_maint/lib/matplotlib/projections/geo.py
===================================================================
--- branches/v1_0_maint/lib/matplotlib/projections/geo.py	2010-11-22 14:01:21 UTC (rev 8807)
+++ branches/v1_0_maint/lib/matplotlib/projections/geo.py	2010-11-22 16:05:17 UTC (rev 8808)
@@ -346,7 +346,7 @@
             cos_latitude = np.cos(latitude)
             sqrt2 = np.sqrt(2.0)
 
-            alpha = 1.0 + cos_latitude * np.cos(half_long)
+            alpha = np.sqrt(1.0 + cos_latitude * np.cos(half_long))
             x = (2.0 * sqrt2) * (cos_latitude * np.sin(half_long)) / alpha
             y = (sqrt2 * np.sin(latitude)) / alpha
             return np.concatenate((x, y), 1)


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