Если внимательно посмотреть на последние заметки нашего блога, то можно подумать, что он превратился в чисто фотографический. Однако недавно мы хорошенько обдумали перспективы развитя наших сайтов, не побоялись и поставили всё с ног на голову. Поскольку фотографическая тематика постов отныне будет жить только на , этот сайт вновь станет выглядеть так, как он выглядел в самом начале. Да-да, начиналось всё не с фотографических обзоров, а с заметок, типа:
Вязание крючком игрушки Бендера (Футурама)
Взрослая загадка на смекалку — расшифровать послание инопланетянам
Работа с GMT
Изготовление трафарета
и т.д.
Сегодня мы вернёмся к теме картирования. Задача стоит следующая.
В отдельных точках акватории проводились измерения океанологических хараткеристик. Для каждой точки нам известны географические координаты (широта и долгота). Необходимо нанести в виде кружков положение точек на карте. Делать всё мы будем в Python (я обычно при этом использую ipython notebook, так удобнее).
Для начала импортируем всё, что нам понадобится:
%pylab inline import matplotlib.pylab as plt from mpl_toolkits.basemap import Basemap
затем: создадим две переменные с широтами и долготами. В данном случае проще и быстрее было создать их вручную:
lat = [66.297,66.299,66.298,66.295,66.301,66.304,66.288,66.289,66.286,66.289] lon = [33.640,33.660,33.690,33.747,33.829,33.908,33.891,33.839,33.781,33.740]
Ну и теперь само построение карты:
m = Basemap(projection='merc',llcrnrlat=66.22,urcrnrlat=66.37,\
llcrnrlon=33.60,urcrnrlon=34,resolution='f')
figsize(10,15)
x, y = m(lon,lat)
m.drawcoastlines()
m.fillcontinents(color='gray',lake_color='white')
m.drawmapboundary(fill_color='white')
m.drawparallels(np.arange(66.22,66.37,.04),labels=[1,0,0,0], fontsize =14)
m.drawmeridians(np.arange(33.60,34.,.1),labels=[0,0,0,1], fontsize =14)
m.scatter(x,y,20,marker='o',color='k')
plt.title("Location of the measurement points", fontsize =14)
plt.show()
Некоторые пояснения к тому, что мы сделали (я не стал это делать в коде, чтобы не перегружать его): Читать далее «Построение точек на карте по заданным координатам при помощи Basemap в Python.»