import xarray as xr import matplotlib.pyplot as plt if __name__=='__main__': ds=xr.open_dataset(r'D:\pythonsx\zl\sst.mnmean.nc') sst=ds['sst'] sstsel=sst.loc[:,5:-5,190:240] #5°S-5°N,170°W-120°W enso=sstsel.mean(dim=['lat', 'lon'])#区域平均修改这段代码,计算enso的距平
时间: 2023-06-13 07:09:18 浏览: 52
import xarray as xr
import matplotlib.pyplot as plt
if __name__=='__main__':
ds=xr.open_dataset(r'D:\pythonsx\zl\sst.mnmean.nc')
sst=ds['sst']
sstsel=sst.loc[:,5:-5,190:240] #5°S-5°N,170°W-120°W
enso_mean=sstsel.mean(dim=['lat', 'lon']) #区域平均
enso_anom= enso_mean - enso_mean.mean() #计算enso的距平
相关问题
import xarray as xr import matplotlib.pyplot as plt if name=='main': ds=xr.open_dataset(r'D:\pythonsx\zl\sst.mnmean.nc') sst=ds['sst'] sstsel=sst.loc[:,5:-5,190:240] #5°S-5°N,170°W-120°W enso=sstsel.mean(dim=['lat', 'lon'])#区域平均 ensom=enso.mean('time') #print(ensom) ensoa=enso-ensom #print(ensoa)修改这段代码,筛选出大于0.5的值
import xarray as xr
import matplotlib.pyplot as plt
if __name__ == '__main__':
ds = xr.open_dataset(r'D:\pythonsx\zl\sst.mnmean.nc')
sst = ds['sst']
sstsel = sst.loc[:, 5:-5, 190:240] # 5°S-5°N,170°W-120°W
enso = sstsel.mean(dim=['lat', 'lon']) # 区域平均
ensom = enso.mean('time')
ensoa = enso - ensom
# 选取大于0.5的值
ensoa_selected = ensoa.where(ensoa > 0.5, drop=True)
print(ensoa_selected)
import xarray as xr import matplotlib.pyplot as plt from netCDF4 import Dataset import cartopy.crs as ccrs#投影方式 import cartopy.feature as cfeature import cartopy.io.shapereader as shpreader import numpy as np import pandas as pd from matplotlib.font_manager import fontManager da=pd.read_csv(r'E:\python11\STATION_58237.txt',sep='\s+') print(da) tem=da['TEM'] #print(tem.shape) rhu=da['RHU'] tem_ave=[] for i in range(0,24): tem_data=tem[i:147:24] print(tem_data) tem_ave[i]=np.mean(tem_data) print(tem_ave) plt.plot(tem_ave)帮我改一下
import xarray as xr
import matplotlib.pyplot as plt
from netCDF4 import Dataset
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import cartopy.io.shapereader as shpreader
import numpy as np
import pandas as pd
from matplotlib.font_manager import fontManager
da = pd.read_csv(r'E:\python11\STATION_58237.txt', sep='\s+')
print(da)
tem = da['TEM']
rhu = da['RHU']
tem_ave = []
for i in range(0, 24):
tem_data = tem[i:147:24]
print(tem_data)
tem_ave.append(np.mean(tem_data)) # 将结果添加到列表中
print(tem_ave)
plt.plot(tem_ave)
plt.show() # 显示图形
相关推荐
![gz](https://img-home.csdnimg.cn/images/20210720083447.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)