-
下载安装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)
'''
网友评论