Plot eine Linie entlang einer konstanten Breite in der Basiskarte, Python, wenn eine konische Projektion

Question

würde Ich mag eine Linie, das zwei Punkte der gleicher Breite ziehen auf ein basemap Beispiel eine konische Kartenprojektion (das heißt, wo Breite sind nicht gerade Linien).

Unabhängig davon, ob ich m.drawgreatcircle oder m.plot, ist die resultierende Linie ein gerade (-I denken ...?) Linie zwischen den beiden Punkten, wie auf einer Linie entgegengesetzt, die entlang einer konstanten Breite geht. Weiß jemand, wie man dieses Problem löst? Einige Beispielcode und das resultierende Bild ist unten. Ich würde diese gelbe gestrichelte Linie lieben, entlang der 55N Linie zu laufen.

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

#set up the map 
m = Basemap(resolution='l',area_thresh=1000.,projection='lcc',\ 
     lat_1=50.,lat_2=60,lat_0=57.5,lon_0=-92.5,\ 
     width=6000000,height=4500000) 


#add some basic map features 
m.drawmeridians(np.arange(-155,-5,10),\ 
     labels=[0,0,0,1],fontsize=8,linewidth=0.5) 
m.drawparallels(np.arange(30,85,5),\ 
     labels=[1,0,0,0],fontsize=8,linewidth=0.5) 
m.drawcoastlines(linewidth=0.5) 
m.drawcountries(linewidth=1) 
m.drawstates(linewidth=0.3) 

#plot some topography data 
m.etopo() 

#draw a line between two points of the same latitude 
m.drawgreatcircle(-120,55,-65,55,linewidth=1.5,\ 
     color='yellow',linestyle='--') 

Entschuldigt, wenn ich etwas sehr einfach fehle ...!

Answer 1

drawgreatcicle funktioniert offenbar nicht korrekt auf einer lcc Projektionskarte.

Sie können immer nur selbst eine Zeile erstellen, anstatt sich auf diese Hilfsfunktion zu verlassen. Zu diesem Zweck erstellen Sie Koordinaten entlang der Linie, projizieren sie und rufen Sie plot.

lon = np.linspace(-120,-65) 
lat = np.linspace(55,55) 
x,y = m(lon,lat) 
m.plot(x,y, linewidth=1.5, color='yellow',linestyle='--')