Несколько графиков отличных кругов с базовой картой Python

python matplotlib-basemap great-circle

276 просмотра

1 ответ

183 Репутация автора

Трудно найти решение этой проблемы:

Построение нескольких путей большого круга из массива

сообщение об ошибке:

Traceback (последний вызов был последним): файл "example.py", строка 41, в файле eq_map.drawgreatcircle (y, x, y2, x2, linewidth = 6, color = 'b') File "/ Library / Frameworks / Python. framework / Версии / 3.5 / lib / python3.5 / site-packages / mpl_toolkits / basemap / init .py ", строка 2893, в drawgreatcircle npoints = int ((dist + 0.5 * 1000. * del_s) / (1000. * del_s )) TypeError: можно только объединить список (не "float") в список


Код Python:

import csv

# Open the data file.
filename = 'sample.csv'

# Create empty lists for the latitudes and longitudes.
lats, lons = [], []
lats2, lons2 = [], []

# Read through the entire file, skip the first line,
#  and pull out just the lats and lons.
with open(filename) as f:
    # Create a csv reader object.
    reader = csv.reader(f)

    # Ignore the header row.
    next(reader)

    # Store the latitudes and longitudes in the appropriate lists.
    for row in reader:
        lats.append(float(row[1]))
        lons.append(float(row[2]))
        lats2.append(float(row[4]))
        lons2.append(float(row[5]))        

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

eq_map = Basemap(projection='merc',llcrnrlat=-80,urcrnrlat=80,llcrnrlon=-180,urcrnrlon=180,lat_ts=20,resolution='c')

eq_map.drawcoastlines()
eq_map.drawcountries()
eq_map.fillcontinents(color = 'gray')
eq_map.drawmapboundary()
eq_map.drawmeridians(np.arange(0, 360, 30))
eq_map.drawparallels(np.arange(-90, 90, 30))

x,y = eq_map(lons, lats)
x2,y2 = eq_map(lons2,lats2)
eq_map.drawgreatcircle(y,x,y2,x2,linewidth=6,color='b')

plt.show()

Sample.csv

Origin,origLat,origLon,Dest,destLat,destLon
"jfk",40.641311,-73.778139,"lax",33.941589,-118.40853
"teb",40.849023,-74.062953,"mia",34.730283,136.508588
Автор: sn4ke Источник Размещён: 18.07.2016 05:13

Ответы (1)


0 плюса

1 Репутация автора

Ваши друзья и латы в данный момент находятся в списке. Вам нужно перебрать, чтобы получить координаты.

print(type(lons))

Вот что я придумал:

for x, y, x2, y2 in zip(lats[0::1], lons[0::1], lats2[0::1], lons2[0::1]):
    eq_map.drawgreatcircle(x, y, x2, y2, linewidth=2, color='b')

Если вы не хотите создавать сюжеты, то вам не нужно

x,y = eq_map(lons, lats)
x2,y2 = eq_map(lons2,lats2)

С участками

for lats, lons, lats2, lons2 in zip(lats[0::1], lons[0::1], lats2[0::1], lons2[0::1]):
    x,y = eq_map(lons, lats)
    x2,y2 = eq_map(lons2,lats2)
    eq_map.plot(x, y, marker='.', color='r', markersize=10)
    eq_map.plot(x2, y2, marker='.', color='r', markersize=10)
    eq_map.drawgreatcircle(lons, lats, lons2, lats2, linewidth=2, color='b')
Автор: Blacktwin Размещён: 12.12.2016 07:12
Вопросы из категории :
32x32