blockgroups.py 文件源码

def LatLonToMeters(lat, lon ):
    """"Converts given lat/lon in WGS84 Datum to XY in Spherical Mercator EPSG:900913
    We Use Lambert Azimuthal Equal Area projection.
    """

    proj = pyproj.Proj(proj='laea')
    mx,my = proj(lon,lat)

    return mx, my