ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

日记

2022-08-31 15:32:44  阅读:186  来源: 互联网

标签:tmax 37 filename 62 import copy 日记


import xarray as xr
import numpy as np
import os
import matplotlib.pyplot as plt
from pandas.core.frame import DataFrame 
from scipy.stats.mstats import ttest_ind
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import cartopy.mpl.ticker as cticker
import matplotlib.patches as patches
filePath = '文件夹路径'
#获取文件夹下所有文件名
filename = os.listdir(filePath)
#文件名过滤
filename = list(filter(lambda x:x.startswith('文件名头'), filename))

tmax_all = np.zeros((9,37*62))
for i in range(9):
    f = xr.open_dataset('/home/mw/input/mysqzj8460/{}'.format(filename[i]))
    tmax = f.TMAX.loc[f.time.dt.month.isin([7,8])].loc['1979-07-01':'2015-08-31']
    tmax_all[i] = tmax.interpolate_na(dim="time", method="linear")
tmax_mean = np.nanmean(tmax_all,axis=0)
tmax_2d = tmax_mean.reshape((37,62))

tmax_copy = tmax_mean.copy()
for i in range(1,37*62-1,1):
    if ((tmax_copy[i]>=33)&(tmax_copy[i-1]<33)&(tmax_copy[i+1]<33)):
        tmax_copy[i]=32
tmax_copy_2d = tmax_copy.reshape((37,62))

 

标签:tmax,37,filename,62,import,copy,日记
来源: https://www.cnblogs.com/eroeg/p/16643241.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[[email protected]]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有