经纬度和墨卡托坐标互相转换

import math
#经纬度转墨卡托坐标
def ls(lat,lng):
    # list=ll_wl.split(',')
    lat=float(lat)
    lng=float(lng)
    x =  lng* 20037508.34 / 180
    y = math.log(math.tan((90 + lat) * math.pi / 360)) / (math.pi / 180)
    y = y * 20037508.34 / 180
    print str(x),str(y)
#墨卡托坐标转经纬度坐标
def ll(x,y):
    x=float(x)
    y=float(y)
    x=x/20037508.34*180
    y=y/20037508.34*180
    y=180/math.pi*(2*math.atan(math.exp(y*math.pi/180))-math.pi/2)
    print x,y
if __name__ == '__main__':
    ls('30.697218','104.073694')
    ll('12985924.36','4789962.05')

python版本为2.7  如有不恰当  请指出  谢谢