window 11 python3.7 读写 ERA5Land
2022-12-30 本文已影响0人
zyhthinking
-
下载安装wrf_python
https://www.lfd.uci.edu/~gohlke/pythonlibs/#wrf_python -
pip 安装
pip install ecmwflibs
pip install eccodes
pip install cfgrib
pip install xarray -
代码编写:
'''
nc_obj = xr.load_dataset("LAI_19981019.grib", engine="cfgrib")
lvdata=nc_obj["lai_lv"].data
lvdata1=lvdata[:,0:1800]
lvdata2=lvdata[:,1800:3600]
xy = np.hstack([lvdata2, lvdata1])
run=GRID()
'''
0:图像左上角的X坐标;
1:图像东西方向分辨率;
2:旋转角度,如果图像北方朝上,该值为0;
3:图像左上角的Y坐标;
4:旋转角度,如果图像北方朝上,该值为0;
5:图像南北方向分辨率;
'''
outadfGeoTransform=[-180,0.1,0,90,0,-0.1]
def write_img(self, filename, im_proj, im_geotrans, im_data)
run.write_img("lv3.tif", "", outadfGeoTransform,xy)
print(nc_obj)
'''