- 积分
- 2997
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-5-9
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 claire11 于 2022-1-19 10:00 编辑
这个东西改了我一周了。。。。救救孩子,再改不好要被自己气死了!
问题描述:
想修改经纬度的数字大小,但是各种尝试均无果。
图:
代码:
# -*- coding: utf-8 -*-
"""
Created on Tue Dec 21 15:28:28 2021
@author: 18112
"""
import xarray as xr
import numpy as np
import proplot as plot
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
file700 = "g:/data/hubei/700.nc"
ds1 = xr.open_dataset(file700)
file850 = "g:/data/hubei/850.nc"
ds2 = xr.open_dataset(file850)
file950 = "g:/data/hubei/950.nc"
ds3 = xr.open_dataset(file950)
fileland = "g:/data/hubei/land1.nc"
ds4 = xr.open_dataset(fileland)
file1000 = "g:/data/hubei/1000.nc"
ds5 = xr.open_dataset(file1000)
#省略我自己的变量计算部分
shr6=np.sqrt(delta1)
shr3=np.sqrt(delta2)
shr1=np.sqrt(delta3)
lon=ds1.longitude
lat=ds1.latitude
i=7
time=ds1.time
#创建画布
fig,axs=plot.subplots(proj='cyl',ncols=1,dpi=200)
ax=axs[0]
fig.format(lonlim=(108,117),latlim=(29,37),labels=True,grid= False,lonlocator=2,latlines=1)
plt.title('6km vertical wind shear and temperature \n2021-08-12T05:00',fontdict={'weight':'normal','size': 8}) #改变图标题字体
map = Basemap(llcrnrlon=108, llcrnrlat=29, urcrnrlon=117, urcrnrlat=34)
map.readshapefile(r'g:/China2020/bou2_4l','states',drawbounds = True) #这个是省 bount-line是县界
#设置底图以及变量绘图
m0 = ax.contourf(tem,globe=True,cmap='ice',levels= np.linspace(22,30,9))
q0 = ax.contour(shr6,linewidth=0.3,color='r',levels = np.arange(0.0, 10.0 , 1))
plt.clabel(q0,fontsize=3,fmt='%.1f',colors='r')
ax.colorbar(m0,fraction=0.05,pad=0.05)
#ax.set_xticklabels(lon, labelsize=2)
#ax.yticklabels(fontsize='8')
plt.ylabel('temperture (℃)', fontdict={'weight': 'normal', 'size': 8})#改变坐标轴标题字体
plt.tick_params(labelsize=8)
plt.savefig('g:/工作文件2021部分/812柳林/图/08120500shr6 and tem.png')
——已解决,原因在于basemap不能直接用我列出来的格式去修改刻度标记字体等,参照了csdn大佬的帖子已经成功修改:https://blog.csdn.net/weixin_44092702/article/details/99690821
不要靠近basemap,会变得不幸.我还是学学cartopy吧。。。
|
|