Базовая карта (Matlibplot) не отображает точные координаты GPS WGS84

У меня есть список координат широты и долготы с фотографий Flickr, сделанных в Сиднее (см. Здесь)

Когда я нанесу эти точки на карты Google, их местоположение будет правильным (см. Здесь)

Однако, когда я использую свою Базовую карту, результаты не верны. Это особенно заметно на береговых линиях: ,

Я безуспешно пытался изменить Базовую карту на множество кодов EPSG, безуспешно пытался преобразовать координаты в разные коды EPSG.

Я не совсем уверен, куда идти отсюда. Сейчас я просто пробую что-то, но, похоже, ничего не работает, поэтому любая помощь будет принята с благодарностью! ура

Основной код:

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np
import math
import pyproj


location_data = get_location_data("sydney")[:500]
bmap = create_map()

for lat, lon in location_data:
    plot_on_map(bmap, float(lat), float(lon))

show_map()

методы:

#
#   Create a Basemap and return it
#
def create_map():
    lower_left_lat = -34.126127
    lower_left_lon = 150.951113
    upper_right_lat = -33.751265
    upper_right_lon = 151.425672


    bmap = Basemap(projection='merc', # doesn't need to be mercator
                  llcrnrlat=lower_left_lat,
                  urcrnrlat=upper_right_lat,
                  llcrnrlon=lower_left_lon,
                  urcrnrlon=upper_right_lon,
                  resolution='f',
                  epsg='4326')

    bmap.drawcoastlines()
    bmap.drawcountries()
    bmap.fillcontinents(lake_color='aqua')
    bmap.drawmapboundary()

    return bmap

#
#   Plot coordinates lat, lon on the Basemap map
#
def plot_on_map(bmap, lat, lon):
    #gda94 = pyproj.Proj(init='epsg:4283')
    #mgaz56 = pyproj.Proj(init='epsg:4326')
    #x, y = pyproj.transform(gda94, mgaz56, lon, lat)

    x,y = bmap(lon, lat)
    bmap.plot(x, y, 'bo', markersize=3)

def show_map():
    plt.show()

0 ответов

Другие вопросы по тегам