code: def _deg2num(self, lat_deg, lon_deg, zoom, leave_float=False):
          
          lat_rad = mod_math.radians(lat_deg)
          n = 2.0 ** zoom
          xtile = (lon_deg + 180.0) / 360.0 * n
          ytile = (1.0 - mod_math.log(mod_math.tan(lat_rad) + (1 / mod_math.cos(lat_rad))) / mod_math.pi) / 2.0 * n
          if not leave_float:
              xtile = int(xtile)
              ytile = int(ytile)
          return TileInfo(xtile, ytile, zoom)
docs: Taken from http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Python