1.卫星图像绘制

import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
# 图像绘制
plt.figure(dpi=300)
m = Basemap(projection='cyl', llcrnrlat=11, llcrnrlon=105, urcrnrlat=35, urcrnrlon=135)raw_lat = raw_lat[490:1080, 655:1470]
raw_lon = raw_lon[490:1080, 655:1470]
data = CH11[490:1080, 655:1470] # 在保证图片完整性的前提下缩小画图范围，减少画图时间x, y = m(raw_lon, raw_lat)
# 绘制轮廓图
cf = m.pcolormesh(x, y, data, cmap='gray')cbar = m.colorbar(cf, location='right', size='5%', pad='2%')
font = {'family': 'serif','color': 'darkred','weight': 'normal','size': 8,}
m.drawmeridians(np.arange(105, 135, 3), linewidth=0.2, labels=[0, 0, 0, 1])
m.drawparallels(np.arange(11, 35, 3), linewidth=0.2, labels=[1, 0, 0, 0])

2.画出点坐标和坐标的图例

xlonT_1, xlatT_1 = m(centerLonT_1, centerLatT_1)
plt.plot(xlonT_1, xlatT_1, 'co', color='orange',markersize=5,label='t-1')  # 添加label参数即可加上图例
# plt.text(xlonT_1, xlatT_1, ' t-1', fontsize=12) # 显示文字标注xlon1, xlat1 = m(centerLonT1, centerLatT1)
plt.plot(xlon1, xlat1, 'co', color='blue',markersize=5,label='t+1')  # 添加label参数即可加上图例
# plt.text(xlon1, xlat1, ' t+1', fontsize=12) # 显示文字标注plt.legend() # 显示图例，这句必须要加不然不显示
plt.xticks(fontsize=8)
plt.yticks(fontsize=8)plt.show()

效果：

另外，参数‘co’可以指定点坐标的形状

co ------- 实心圆
c^ ------- 三角形
c* ------- 五角星