import sys import pyproj extents = 122115.0, 485775.0 # lower left coordinate of science park map cellsize = 5, 5 # cell size of science park map src = pyproj.Proj(proj='latlong',datum='WGS84') # regular GPS coordinate system (used in google maps) dst = pyproj.Proj(init='EPSG:28992') # coordinate system of NL/science park map reverse = False if sys.argv[1] == '--reverse': reverse = True del sys.argv[1] if not reverse: lat = float(sys.argv[1]) lng = float(sys.argv[2]) mx, my = pyproj.transform(src, dst, lng, lat) x = int((mx - extents[0]) / cellsize[0]) y = int((extents[1] - my) / cellsize[1]) print 'Lat/Lng:', lat, lng print 'Map X, Y:', mx, my print 'X, Y:', x, y else: x = float(sys.argv[1]) y = float(sys.argv[2]) mx = (x * cellsize[0]) + extents[0] my = (y * cellsize[1]) + extents[1] lng, lat = pyproj.transform(dst, src, mx, my) print 'X, Y:', x, y print 'Map X, Y:', mx, y print 'Lat/Lng:', lat, lng