Merge pull request #281 from micahcochran/patch-1

Add diagnostics functions: proj4 ver. & Inv. Hammer

Author: WeatherGod

lib/mpl_toolkits/basemap/diagnostic.py

@@ -2,6 +2,23 @@
 These are diagnostic and debugging functions for basemap.
 """
 
+def proj4_version():
+    """
+    Gives the proj.4 library's version number. (requires pyproj to be installed)
+
+    returns string, so proj.4 version 4.9.3 will return "4.9.3"
+    """
+    import pyproj
+    
+    # Get PROJ4 version in a floating point number
+    proj4_ver_num = pyproj.Proj(proj='latlong').proj_version
+    
+    # reformats floating point number into string (4.90 becomes '4.9.0')
+    # Exploits single number version numbers for proj4,
+    # This will need likely need to be change at some point as proj.4 version 4.10.0???
+    return '.'.join( str(int(proj4_ver_num*100)) )
+    
+    
 def package_versions():
     """
     Gives version information for dependent packages.
@@ -13,7 +30,7 @@ def package_versions():
 
     from matplotlib import __version__ as matplotlib_version
     from numpy import __version__ as numpy_version
-    import pyproj
+    from pyproj import __version__ as pyproj_version
     from shapefile import __version__ as pyshp_version
 
     import _geoslib
@@ -35,10 +52,6 @@ def package_versions():
         pil_version = 'not installed'
         pillow_version = 'not installed'
 
-    # Get PROJ.4 version info in a floating point number
-    proj_ver_num = pyproj.Proj(init='epsg:4326').proj_version
-    # reformats floating point number into string (4.90 becomes '4.9.0')
-    proj4_version = '.'.join(list(str(int(proj_ver_num*100))))
 
     BasemapPackageVersions = namedtuple(
                                'BasemapPackageVersions',
@@ -51,11 +64,45 @@ def package_versions():
                    basemap = basemap_version,
                    matplotlib = matplotlib_version,
                    numpy = numpy_version,
-                   pyproj = pyproj.__version__,
+                   pyproj = pyproj_version,
                    pyshp = pyshp_version,
-                   PROJ4 = proj4_version,
+                   PROJ4 = proj4_version(),
                    GEOS = _geoslib.__geos_version__,
                    # optional dependencies below
                    OWSLib = OWSLib_version,
                    PIL = pil_version,
                    Pillow = pillow_version)
+
+def check_proj_inv_hammer(segfault_protection=True):
+    """
+    Check if the inverse of the hammer projection is supported by installed
+    version of PROJ4.
+    
+    segfault_protection   True (default) - test while protecting from segfault
+                          False -  testing that might cause Python to segfault.
+                                   BE CAREFUL setting this flag to False!
+                                   If it segfaults, this the inverse hammer is not supported.
+
+    returns True      - inverse hammer is supported
+            False     - inverse hammer is not supported
+            "Unknown" - support is Unknown
+    """
+    from distutils.version import LooseVersion
+    from pyproj import __version__ as pyproj_version
+    
+    if LooseVersion(proj4_version()) > LooseVersion('4.9.2'):
+        return True
+    
+    if LooseVersion(pyproj_version) > LooseVersion('1.9.5.1') \
+            or segfault_protection is False:
+        from pyproj import Proj
+        hammer = Proj(proj='hammer')
+        
+        x, y = hammer(-30.0, 40.0)
+        try:
+            lon, lat = hammer(x, y, inverse=True)
+            return True
+        except RuntimeError:            
+            return False
+    
+    return 'Unknown'