Угловые координаты плитки статической карты Google

12

Как рассчитать угловые координаты (ограничивающий прямоугольник) статического изображения карты Google, полученного, например, с помощью http://maps.googleapis.com/maps/api/staticmap?center=40.714728,-73.998672&zoom=12&size=400x400&sensor=false ?

Boocko
источник
Можете ли вы объяснить, зачем вам нужны угловые координаты каждого изображения?
Mapperz
Я не ОП, но есть аналогичная необходимость. У меня есть GPS-координаты, которые я хочу построить, и мне нужно фоновое изображение для отображения данных. Чтобы графики работали, мне нужно убедиться, что у меня есть точные координаты ограничительной рамки из изображения карты Google.
Джун

Ответы:

9

Я думаю, что это не так сложно решить проблему, но у меня есть сомнения по поводу того, является ли точность абсолютно правильной. Прежде всего, вам нужно конвертировать ваши центральные широты / долготы в пиксели с кодами gdal2tiles. Если я найду время и при желании смогу преобразовать его в стабильный код для нахождения угловых координат.

Это код Python:

tileSize = 256
initialResolution = 2 * math.pi * 6378137 / tileSize
# 156543.03392804062 for tileSize 256 pixels
originShift = 2 * math.pi * 6378137 / 2.0
# 20037508.342789244

def LatLonToMeters( lat, lon ):
        "Converts given lat/lon in WGS84 Datum to XY in Spherical Mercator EPSG:900913"

        mx = lon * originShift / 180.0
        my = math.log( math.tan((90 + lat) * math.pi / 360.0 )) / (math.pi / 180.0)

        my = my * originShift / 180.0
        return mx, my


def MetersToPixels( mx, my, zoom):
        "Converts EPSG:900913 to pyramid pixel coordinates in given zoom level"

        res = Resolution( zoom )
        px = (mx + originShift) / res
        py = (my + originShift) / res
        return px, py


# Dont forget you have to convert your projection to EPSG:900913
mx = -8237494.4864285 #-73.998672
my = 4970354.7325767 # 40.714728
zoom = 12

pixel_x, pixel_y = LatLonToMeters(MetersToPixels( mx, my, zoom))

Тогда вы можете использовать сложение или вычитание, посмотрев следующее изображение:

МАТЕМАТИКА

Если вы хотите найти точку А:

x = pixel_x - 200
y = pixel_y + 200

или вы хотите найти точку B:

x = pixel_x + 200
y = pixel_y + 200

и последнее, что вам нужно сделать, это конвертировать ваши пиксели в широты и долготы.

def PixelsToMeters( px, py, zoom):
    "Converts pixel coordinates in given zoom level of pyramid to EPSG:900913"

    res = Resolution(zoom)
    mx = px * res - originShift
    my = py * res - originShift
    return mx, my

def MetersToLatLon( mx, my ):
    "Converts XY point from Spherical Mercator EPSG:900913 to lat/lon in WGS84 Datum"

    lon = (mx / originShift) * 180.0
    lat = (my / originShift) * 180.0

    lat = 180 / math.pi * (2 * math.atan(math.exp(lat * math.pi / 180.0)) - math.pi / 2.0)
    return lat, lon




#Result

llx, lly = MetersToLatLon( PixelsToMeters( x, y, zoom) )

итак, результат у меня есть:

point A - UpperLeftLatLon = 40.7667530977 -74.0673365509
point B - UpperRightLatLon = 40.7667530977 -73.9300074493
point C - LowerRightLatLon = 40.6626622172 -73.9300074493
point D - LowerLeftLatLon = 40.6626622172 -74.0673365509

Я надеюсь, что это поможет вам ....

Арагон
источник
Спасибо. Точность в порядке. Любой простой способ рассчитать Resolution(zoom)? Это известно для OSM, но я не могу найти его для карт Google.
Boocko
попробуй использовать initialResolution / (2**zoom).
Арагон
У меня есть сомнения относительно последней части. Как преобразовать проекцию (-73.998672, 40.714728) в EPSG: 900913 (-8237494.4864285, 4970354.7325767)?
norman784
Здесь намекают, что разрешение (масштабирование) зависит от широты: gis.stackexchange.com/questions/108373/…
arivero