Python用Basemap画位势高度距平图和5880线

最新推荐文章于 2025-05-04 14:00:00 发布

账号经常没登陆，私信评论忘回见谅

最新推荐文章于 2025-05-04 14:00:00 发布

阅读量3.4k

点赞数 4

CC 4.0 BY-SA版权

分类专栏：气象绘图文章标签： python

本文链接：https://blog.youkuaiyun.com/mohenxiang/article/details/125210902

import xarray as xr
import matplotlib.pyplot as plt
import numpy as np
from mpl_toolkits.basemap import Basemap
plt.rcParams['font.sans-serif'] = ['KaiTi','SimHei']
plt.rcParams['font.size']= 12
plt.rcParams['axes.unicode_minus']=False
hgt=r"D:\论文\论文数据\nc文件1998hgt.mon.mean.nc"  #读取1998年目标文件
def DrawSeaLPJuly():
    Hgt = xr.open_dataset(hgt)["hgt"]
    time_label = (Hgt["time"].dt.month == [6,7,8])    #选定6，7，8三个月
    Hgt_main=Hgt[time_label].loc[dict(level=500)].mean("time").loc[60:0,70:180]  # 选定气压为500hpa，选定经纬度为0~60N，70~180E
    lat = np.array(Hgt_main['lat'])  # 设置好经纬度
    lon=np.array(Hgt_main['lon'])  # 设置好经纬度
   #  算气候状态（1981~2010）30年：
    i=1981
    hgt1=r"D:\论文\论文数据\nc文件1998hgt.mon.mean.nc"
    Hgt1= xr.open_dataset(hgt1)["hgt"]
    time_label = (Hgt1["time"].dt.month == [6, 7, 8])
    Hgt_total = Hgt1[time_label].loc[dict(level=500)].mean("time").loc[60:0,70:180]
    i=i+1
    while i<=2010:
        hgt1=r"D:\论文\论文数据\nc文件"+str(i)+"hgt.mon.m