我需要在netCDF CF文件中使用xarray的加权纬度计算全球平均值，然后将其转换为Pandas
回答此问题可获得 20 贡献值，回答如果被采纳可获得 50 分。
<p>我需要从netCDF cf data 3D（time，lat，lon）文件计算一个全局时间序列（time），然后将其转换为pandas/dataframe。我需要用cos（lat）来称量纬度。我一直在使用numpy进行平均，但是将pandas/dataframe转换为加权数组不起作用</p>
<pre><code>ds=xr.open_dataset('sample_data.nc')
data=ds.tas
start_time='1980-01-01'
end_time='2018-12-31'
time_slice = slice(start_time, end_time)
nrows=len(data.lat.values)
ncols=len(data.lon.values)
t=len(data.time.values)
weights=np.zeros([len(data.lat.values)])
latsr = np.deg2rad(data.lat.values).reshape((nrows,1))
weight_matrix=np.repeat(np.cos(latsr),ncols,axis=1)
wghtpr=np.zeros_like(data)
for i in range (0,t):
wghtpr[i,:,:]=data[i,:,:]*weight_matrix
new_data=wghtpr
wtdata=np.average(new_data,axis=1)
da=np.average(wtdata,axis=1)
</code></pre>
<p>这将导致一个没有“名称”的numpy数组</p>
<p>如果我做ds，我会得到：</p>
<pre><code><xarray.Dataset>
Dimensions: (bnds: 2, lat: 361, lon: 576, time: 477)
Coordinates:
* time (time) datetime64[ns] 1980-01-16T12:00:00 ... 2019-09-16
* lat (lat) float64 -90.0 -89.5 -89.0 -88.5 ... 88.5 89.0 89.5 90.0
* lon (lon) float64 0.0 0.625 1.25 1.875 ... 357.5 358.1 358.8 359.4
height float64 ...
Dimensions without coordinates: bnds
Data variables:
time_bnds (time, bnds) datetime64[ns] ...
lat_bnds (lat, bnds) float64 ...
lon_bnds (lon, bnds) float64 ...
tas (time, lat, lon) float32 244.15399 244.15399 ... 267.52875
Attributes:
institution: Global Modeling and Assimilation Office, NASA Goddard Sp...
institute_id: NASA-GMAO
experiment_id: MERRA-2
source: MERRA-2 Monthly tavgM_2d_slv_Nx
model_id: GEOS-5
references: http://gmao.gsfc.nasa.gov/research/merra/, http://gmao.g...
tracking_id: e77fd4de-19c2-45ad-afe2-ce3f6c1eb148
mip_specs: CMIP5
source_id: MERRA-2
product: reanalysis
frequency: mon
creation_date: 2015-10-11T23:12:34Z
history: 2015-10-11T23:12:34Z CMOR rewrote data to comply with CF...
Conventions: CF-1.4
project_id: CREATE-IP
table_id: Table Amon_ana (10 March 2011) c3ffdce87438d8df0839620ee...
title: Reanalysis output prepared for CREATE-IP.
modeling_realm: atmos
cmor_version: 2.9.1
doi: http://dx.doi.org/10.5067/AP1B0BA5PD2K
contact: MERRA-2, Steven Pawson (steven.pawson-1@nasa.gov)
#
</code></pre>
请先 登录 后评论
针对此问题: 更多的回答
- 89 关注
- 1 收藏，216 浏览
- 网友 提问于 2天前
相关Python问题
- 尽管Python中的所有内容都是引用，为什么Python导师在没有指针的列表中绘制字符串和整数？ 8 回答
- 尽管python中的表达式为false，但循环仍在运行 8 回答
- 尽管python代码正确，但从nifi ExecuteScript处理器获取语法错误 3 回答
- 尽管Python在Neovim中工作得很好，但插件不能识别Neovim中的Python主机 4 回答
- 尽管python字典包含了大量的条目，但它并没有增长 7 回答
- 尽管python说模块存在，为什么我会得到这个消息？ 2 回答
- 尽管setuptools和控制盘是最新的，但无法识别singleversionexternallymanaged 2 回答
- 尽管stdout和stderr重定向，但未捕获错误消息 10 回答
- 尽管Tensorboard的事件太大，但Tensorboard的步骤太少了 4 回答
- 尽管tkinter上的变量已更改，但显示未更改 10 回答
- 尽管try/except使用Python进行单元测试时出现断言错误 6 回答
- 尽管URL是sam，但仍会抛出“达到最大重定向” 1 回答
- 尽管url有效，Pandas仍读取url的\u csv错误 1 回答
- 尽管while中存在时间延迟，但LINUX线程的CPU利用率为100%（1） 6 回答
- 尽管x0在范围内，Scipy优化仍会引发ValueError 7 回答
- 尽管xpath正确，但使用selenium单击链接仍不起作用 2 回答
- 尽管下载了ffmpeg并设置了路径变量python，但没有后端错误 4 回答
- 尽管下载了i，但找不到型号“fr” 9 回答
- 尽管下载了plotnine包，但未获取名为“plotnine”的模块时出错 10 回答
- 尽管为所有行指定了权重，网格（0）仍不起作用 4 回答
