Я пытаюсь построить контурную диаграмму из совокупной суммы ударов молнии по региону Пуэрто-Рико, однако данные очерчиваются в правом нижнем углу графика.Все остальное работает правильно, за исключением этого расположения графика, которое должно появиться на вершине Пуэрто-Рико около верхнего левого угла.Я сбит с толку, почему это происходит из-за того, что на графиках рассеяния я сделал те же данные по широте и долготе, которые идеально позиционируют его по региону, только контур ниже доставляет мне проблемы.
Из-за исследования подобных проблемПохоже, что решение заключается в том, чтобы как-то настроить сетку, но когда я изменяю ее, график становится больше или меньше, но не перемещает позиции.Возможно, что-то не так с соотношением размеров графика по сравнению с базовой картой?
Вот что у меня есть:
import numpy as n
import pandas as pd
import matplotlib.pyplot as p
import datetime as datetime
from scipy import stats
import matplotlib.pyplot as p
from mpl_toolkits.basemap import Basemap
from matplotlib.colors import Normalize
data = n.genfromtxt("/home/008918949/eni_tl.txt", usecols=
(0,1,2,3,4,5,6,7,8,9,10,11))
month= data[:,0]
day= data[:,1]
year = data[:,2]
hour= data[:,3]
minute = data[:,4]
second = data[:,5]
latitude = data[:,6]
longitude = data[:,7]
amplitude = data[:,8]
l_type = data[:,10]
n_strikes = data[:,11]
n_strikes=n.cumsum(n_strike) #cumulative sum of the strikes
#Plotting Basemap:
norm = Normalize()
p.figure(figsize=(17,7))
lllat=17.4
ulat=18.6
lllon=-67.6
ulon=-64.4
my_map = Basemap(projection='mill', lat_0=18, lon_0=-66,
resolution = 'h', area_thresh = 50,
llcrnrlon=lllon, llcrnrlat=lllat,
urcrnrlon=ulon, urcrnrlat=ulat)
my_map.drawcoastlines()
my_map.drawcountries()
my_map.fillcontinents(color='lightgreen')
my_map.drawmapboundary(fill_color = 'white')
my_map.drawmeridians(n.arange(lllon, ulon, 1.),
color = '0.25', linewidth = 0.5,
labels=[False, False, False, True])
my_map.drawparallels(n.arange(lllat, ulat, 1.),
color = 'black', linewidth = 0.1,
labels=[True, False, False, False])
my_map.drawmapscale(18 + 0.08, -66 + 0.015,
18, -66,
10.,
barstyle='fancy', labelstyle='simple',
fillcolor1='w', fillcolor2='#555555',
fontcolor='#555555',
zorder=5)
#Transform lon / lat coordinates from data to map projection
lon2012=map(int, longitude) #making sure all data is in integers
lat2012=map(int, latitude)
lon2012, lat2012 = my_map(*(lon2012, lat2012))
#Grid data:(These are the values for the minumum and maximum latitude and longitude of data in map coordinates):
lon2012min= 66716.924570808071 #meters
lon2012max= 400301.54742485163
lat2012min= -45792.225173141342
lat2012max= 68807.774434396997
xi = n.linspace(lon2012min, lon2012max)
yi = n.linspace(lat2012min, lat2012max)
xi, yi = n.meshgrid(xi, yi)
#Interpolate
x, y, z = lon2012, lat2012, n_strikes
from scipy.interpolate import griddata as gd
zi = gd((lon2012, lat2012),n_strikes,(xi, yi),
method='linear')
zimin=n.min(zi)
zimax=n.max(zi)
#Contour plot alpha 0.6
con = my_map.contourf(xi, yi, zi, zorder=4, alpha=0.6, cmap='YlOrRd')
#Add color bar and title
#Add color bar, title, and scale
cbar = p.colorbar(con, orientation='horizontal', fraction=.057, pad=0.05)
cbar.set_label("Cumulative Sum of Lightning Strikes (Counts)")
my_map.drawmapscale(
-67.4, 19.5, 28., -13,
100,
units='km', fontsize=10,
yoffset=None,
barstyle='fancy', labelstyle='simple',
fillcolor1='w', fillcolor2='#000000',
fontcolor='#000000',
zorder=5)
p.title("Cumulative Lightning Strikes Over Puerto Rico and Virgin Islands
Contour Plot (2012-2017)",
fontsize=18)
p.savefig("contourtest.png", dpi=100)
p.show()
Вот каков результат:
Также, если есть какие-то вопросы, минимальные и максимальные значения для lat2012 и lon2012 в сетке верны, потому что я выполнил min () и max ()функционировать на них.Я просто выписал действительные числовые значения в качестве проверки, но в любом случае дает тот же результат.
Мои извинения, если это легко исправить, но я потратил много времени на эту проблему, поэтомулюбая помощь приветствуется.