Шейп-файл, используемый для вычисления часового пояса, не поддерживается больше.
Я только что столкнулся с той же проблемой сегодня, и я не уверен, насколько уместен мой ответ по прошествии всего этого времени, но я просто написал функцию Python, которая делает то, что вы хотите. Вы можете найти это здесь.
https://github.com/cstich/gpstotz
Редактировать:
Как уже упоминалось в комментариях, я также должен опубликовать код. Код основан на шейп-файле часовых поясов Эрика Мюллера, который вы можете получить здесь - http://efele.net/maps/tz/world/.
Редактировать 2:
Как оказалось, шейп-файлы имеют несколько архаичное определение внешних и внутренних колец (в основном, внешние кольца используют правило правой руки, в то время как внутренние кольца используют правило левой руки). В любом случае Фиона, кажется, позаботилась об этом, и я соответствующим образом обновил код.
from rtree import index # requires libspatialindex-c3.deb
from shapely.geometry import Polygon
from shapely.geometry import Point
import os
import fiona
''' Read the world timezone shapefile '''
tzshpFN = os.path.join(os.path.dirname(__file__),
'resources/world/tz_world.shp')
''' Build the geo-index '''
idx = index.Index()
with fiona.open(tzshpFN) as shapes:
for i, shape in enumerate(shapes):
assert shape['geometry']['type'] == 'Polygon'
exterior = shape['geometry']['coordinates'][0]
interior = shape['geometry']['coordinates'][1:]
record = shape['properties']['TZID']
poly = Polygon(exterior, interior)
idx.insert(i, poly.bounds, obj=(i, record, poly))
def gpsToTimezone(lat, lon):
'''
For a pair of lat, lon coordiantes returns the appropriate timezone info.
If a point is on a timezone boundary, then this point is not within the
timezone as it is on the boundary. Does not deal with maritime points.
For a discussion of those see here:
http://efele.net/maps/tz/world/
@lat: latitude
@lon: longitude
@return: Timezone info string
'''
query = [n.object for n in idx.intersection((lon, lat, lon, lat),
objects=True)]
queryPoint = Point(lon, lat)
result = [q[1] for q in query
if q[2].contains(queryPoint)]
if len(result) > 0:
return result[0]
else:
return None
if __name__ == "__main__":
''' Tests '''
assert gpsToTimezone(0, 0) is None # In the ocean somewhere
assert gpsToTimezone(51.50, 0.12) == 'Europe/London'