python - DataArray.mean 不保留坐标
问题描述
DataArray.mean
不保留取决于应用平均值的维度的坐标。
注意:XLAT 和 XLONG 不依赖于时间;但是,一些 netcdf 文件在这两个上具有时间坐标。
我有这个 netcdf 文件 wrfout_d03.nc,我在其中打开文件:
ds = xr.open_dataset('/Users/jacob/Desktop/wrfpy/wrfout_d03_may.nc')
然后给出一个 DataSet 对象:
<xarray.Dataset>
Dimensions: (Time: 193, bio_emissions_dimension_stag: 41, bottom_top: 50, bottom_top_stag: 51, klevs_for_dvel: 1, seed_dim_stag: 12, soil_layers_stag: 4, south_north: 115, south_north_stag: 116, west_east: 115, west_east_stag: 116)
Coordinates:
XLAT (Time, south_north, west_east) float32 ...
XLONG (Time, south_north, west_east) float32 ...
XTIME (Time) datetime64[ns] ...
XLAT_U (Time, south_north, west_east_stag) float32 ...
XLONG_U (Time, south_north, west_east_stag) float32 ...
XLAT_V (Time, south_north_stag, west_east) float32 ...
XLONG_V (Time, south_north_stag, west_east) float32 ...
Dimensions without coordinates: Time, bio_emissions_dimension_stag, bottom_top, bottom_top_stag, klevs_for_dvel, seed_dim_stag, soil_layers_stag, south_north, south_north_stag, west_east, west_east_stag
Data variables:
datavars...
然后我使用以下代码访问PM2_5_DRY变量:
pm25 = ds.PM2_5_DRY
生成的对象pm25具有以下维度和坐标:
<xarray.DataArray 'PM2_5_DRY' (Time: 193, bottom_top: 50, south_north: 115, west_east: 115)>
[127621250 values with dtype=float32]
Coordinates:
XLAT (Time, south_north, west_east) float32 ...
XLONG (Time, south_north, west_east) float32 ...
XTIME (Time) datetime64[ns] ...
Dimensions without coordinates: Time, bottom_top, south_north, west_east
Attributes:
FieldType: 104
MemoryOrder: XYZ
description: pm2.5 aerosol dry mass
units: ug m^-3
stagger:
然后我操纵pm25对象并通过以下方式获得时间维度的平均值:
pm25_mean = pm25.mean(dim='Time', keep_attrs = True)
结果对象是一个 DataArray,但没有坐标 XLAT 或 XLON。
<xarray.DataArray 'PM2_5_DRY' (bottom_top: 50, south_north: 115, west_east: 115)>
array([[[14.73083 , 14.756626 , 14.796355 , ..., 20.325712 ,
20.855696 , 21.381271 ],
[14.651459 , 14.34477 , 14.371858 , ..., 18.00389 ,
18.4109 , 21.337002 ],
[14.59026 , 14.257076 , 14.293012 , ..., 17.391146 ,
18.217058 , 20.882664 ],
...,
[27.356459 , 27.21468 , 27.757051 , ..., 8.084272 ,
8.010168 , 7.989942 ],
[27.185486 , 27.02623 , 27.776043 , ..., 7.944748 ,
7.8795266 , 7.8552976 ],
[26.926008 , 27.724253 , 28.427626 , ..., 7.8269224 ,
7.773637 , 7.741844 ]],
Dimensions without coordinates: bottom_top, south_north, west_east
Attributes:
FieldType: 104
MemoryOrder: XYZ
description: pm2.5 aerosol dry mass
units: ug m^-3
stagger:
以下要检查的代码给出:
pm25_mean.coords
Coordinates:
*empty*
我尝试查看 xarray 中 mean 函数的文档;但是,我找不到将坐标从前一个对象复制到新对象的任何选项。
关于如何解决这个问题的任何提示?我想我需要从文件中访问这些坐标,然后再次组合它们。但我不确定如何完成这个过程。
还有,跟这个有关系吗?
XLAT (Time, south_north, west_east) float32
XLAT 是一个多维坐标,它也取决于时间。因为我得到了时间维度的平均值,所以维度 pm25 的数量减少到 3 而不是 4。这是否也会以某种方式影响坐标?
我需要最终对象具有坐标 XLAT 和 XLONG,因为我将对此进行可视化。
谢谢你的帮助!
解决方案
您只需要将坐标变量从 3 维转换为 2 维。
d = xr.open_dataset('.../pm25_sample.nc')
d['XLAT'] = d.XLAT.mean(dim = 'Time')
d['XLONG'] = d.XLONG.mean(dim = 'Time')
d['PM2_5_DRY'].mean(dim = 'Time')
...,
[ 0.03839084, 0.03837739, 0.03835952, ..., 0.03414929,
0.03412561, 0.03410038],
[ 0.03837854, 0.03836632, 0.03834687, ..., 0.03414606,
0.0341224 , 0.03409675],
[ 0.03836945, 0.03835024, 0.03833132, ..., 0.03414177,
0.03411727, 0.03409337]]], dtype=float32)
Coordinates:
XLAT (south_north, west_east) float32 14.086891 14.086907 ... 15.111996
XLONG (south_north, west_east) float32 120.49799 120.50717 ... 121.55791
Dimensions without coordinates: bottom_top, south_north, west_east
推荐阅读
- java - SwingFXUtils 替代图像序列化(Javafx、Swing、Raspberry Pi)
- php - 如何为所有用户刷新页面?
- elasticsearch - 如何在弹性搜索聚合中搜索?
- angular - angular 6 ionic 4 默认组件未注册
- laravel-5.4 - 方法类别不存在。拉拉维尔
- javascript - ExtJs:如何在 html 中使用变量
- css - CSS3 改变动画中的内容
- apache-spark - Apache Spark 中的 While 循环
- node.js - node.js 服务器中的 Firebase 身份验证持久性
- google-apps-script - 如何为同一个电子表格上的不同工作表应用相同的脚本?