如何对多边形内的所有网格单元求和以创建每日总计Netcdf

问题描述

我对python和netcdf文件非常陌生,无法将它们包裹住。我已将存储在netcdf中的每日平均降水量数据网格化。我有一个多边形,我曾使用它来裁剪netcdf文件以供使用:

AWRA_Daily = xarray.open_dataarray('qtot.nc')
AWRA_Daily.rio.set_spatial_dims(x_dim="longitude",y_dim="latitude",inplace=True)
AWRA_Daily.rio.write_crs("epsg:4283",inplace=True)
Narran_Shape = geopandas.read_file('shapefile.shp',crs="epsg:4283")                                                        
clipped = AWRA_Daily.rio.clip(Narran_Shape.geometry.apply(mapping),Narran_Shape.crs,drop=False)

我遇到的问题是,我需要每天汇总所有网格单元,以针对该总降水量区域创建时间序列。因为在每个单元格中,我需要每天绘制的区域都需要为总“降雨”做出贡献。任何帮助将不胜感激。

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)