We haven’t always moved to DST on the second Sunday in March. Query = [n.object for n in idx.intersection((lon, lat, lon, lat),Īssert gpsToTimezone(0, 0) is None # In the ocean somewhereĪssert gpsToTimezone(51.50, 0.I’m drafting this article on the morning of March 10th, 2019, which is the morning when the United States collectively moves our clocks forward by one hour and makes the switch to daylight saving time (DST). Does not deal with maritime points.įor a discussion of those see Timezone info string If a point is on a timezone boundary, then this point is not within the Idx.insert(i, poly.bounds, obj=(i, record, poly))įor a pair of lat, lon coordiantes returns the appropriate timezone info. TzshpFN = os.path.join(os.path.dirname(_file_),Īssert shape = 'Polygon'Įxterior = shape ''' Read the world timezone shapefile ''' In any case fiona seems to take care of that and I updated the code accordingly. Īs it turns out shapefiles have a somewhat archaic definition of exterior and interior rings (basically exterior rings are using the right hand rule, while interior rings are using the left hand rule). The code is based on Eric Muller's shapefile of timezones, which you can get here. You can find it here.Īs mentioned in the comments I should also post code. I just faced the same issue today, and I am not sure how relevant my answer is after all this time, but I basically just wrote a Python function that does what you want. The shapefile used to compute the timezone is not maintained anymore.
0 kommentar(er)
