cs2cs.py 文件源码

def geod2utm (zn, datum, lat, lon, elev) :
    '''   Convert geodetic coordinates to UTM   '''
    if zn == None :
        zn = lon2zone (lon)

    p = Proj (proj='utm', zone=zn, ellps=datum)

    X, Y = p (lon, lat)

    #   Return Y, X, Z
    return Y, X, elev