import math


def latlng2tile_TD(lat_deg, lng_deg, zoom):
    """
    convert latitude, longitude and zoom into tile in x and y axis
    referencing http://www.cnblogs.com/Tangf/archive/2012/04/07/2435545.html
    Keyword arguments:
    lat_deg -- latitude in degree
    lng_deg -- longitude in degree
    zoom    -- map scale (0-18)
    Return two parameters as tile numbers in x axis and y axis
    """
    if lat_deg >= 85.05112877980659 or lat_deg <= -85.05112877980659:
        raise Exception('wmts latitude error lat')
    n = math.pow(2, int(zoom + 1))
    reg = 360.0 / n
    x = (lng_deg + 180.0) // reg
    y = (90.0-lat_deg) // reg
    return x, y