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

Я столкнулся с проблемой при попытке нанести прямоугольник изучаемой области на графику базовой карты. Я пробовал два разных варианта: 1. Используя Polygon. 2. Использование drawgreatcircle.

Это мой код с Polygon:

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

m = Basemap(width=6000000,height=5000000,rsphere=(6378137.00,6356752.3142),\
                resolution='l',projection='lcc',\
                lat_1=0.,lat_2=30,lat_0=15,lon_0=-75.)

lw = 1
dashes = [5,7]
graticules_color = 'black'

fig1 = plt.figure(figsize=(16,20))
fig1.patch.set_facecolor('#e6e8ec')
ax = fig1.add_axes([0.1,0.1,0.8,0.8])

m.etopo()

m.drawparallels(np.arange(-90.,90.,10.),labels=[1,0,0,0],fontsize=20,linewidth=lw, dashes=dashes, color=graticules_color)
m.drawmeridians(np.arange(-180.,180.,10.),labels=[0,0,0,1],fontsize=20,linewidth=lw, dashes=dashes, color=graticules_color)
m.drawmapboundary(fill_color='white')
m.drawcountries()

x1,y1 = m(-100,0)
x2,y2 = m(-100,30)
x3,y3 = m(-50,30)
x4,y4 = m(-50,0)

poly = Polygon([(x1,y1),(x2,y2),(x3,y3),(x4,y4)],facecolor='none',edgecolor='k',linewidth=1.5)
plt.gca().add_patch(poly)

Вот результат: Box с использованием набора инструментов Polygon введите здесь описание изображения< /a> По сути, прямоугольник не находится в координатах проекции, которые я использую на графике (или, по крайней мере, линии не следуют параллелям и меридианам, как это делает сетка).

Глядя на второй вариант (с функцией drawgreatcircle), он работает так, как мне нужно для правой, левой и нижней линий, но не для верхней. Дополнительный код:

lo1 = -100
lo2 = -50
la1 = 0
la2 = 30
#l1
m.drawgreatcircle(lo1, la1, lo1, la2,linewidth=2,color='r')
#l2
m.drawgreatcircle(lo1, la2, lo2, la2,linewidth=2,color='r')
#l3
m.drawgreatcircle(lo2, la2, lo2, la1,linewidth=2,color='r')
#l4
m.drawgreatcircle(lo2, la1, lo1, la1,linewidth=2,color='r') 

И рисунок: Коробка с нарисованным кружком введите здесь описание изображения Может ли кто-нибудь помочь мне решить этот вопрос?

Большое спасибо!

Анджела.


person Angela Maria Gomez Garcia    schedule 03.05.2018    source источник


Ответы (1)


Добавьте этот фрагмент кода к себе. Он построит линию, которую вы хотите.

lons = np.linspace(-100, -50, 15)
lats = np.ones(len(lons))*30
m.plot(lons, lats, color="r", linewidth=1.5, latlon=True)
person swatchai    schedule 03.05.2018
comment
Работает отлично!, я использовал самый сложный способ ^^ - person Angela Maria Gomez Garcia; 10.05.2018