如何将从netCDF4获得的plt图保存为地理参考geoTiff？

问题描述

以下代码来自http://www.acgeospatial.co.uk/sentinel-5p-and-python/ 这是一个关于如何处理 Sentinel 5p netCDF 数据的非常好的教程。 现在我尝试将获得的 plt 图像保存为地理参考 Tiff (geoTiff)。如何做到这一点？

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from netCDF4 import Dataset
import numpy as np  

file = 'S5P_OFFL_L2__O3_____20201025T093552_20201025T111721_15720_01_020103_20201027T030032.nc'
fh = Dataset(file,mode='r')

  
lons = fh.groups['PRODUCT'].variables['longitude'][:][0,:,:]
lats = fh.groups['PRODUCT'].variables['latitude'][:][0,:]
o3 = fh.groups['PRODUCT'].variables['ozone_total_vertical_column'][0,:]
o3_units = fh.groups['PRODUCT'].variables['ozone_total_vertical_column_precision'].units
#o3_units = o3_units

m = Basemap(width=5000000,height=3500000,\
            llcrnrlon=-15.,llcrnrlat=30.,urcrnrlon=80.,urcrnrlat=80.,\
            resolution='c',projection='merc',\
            #lat_ts=40,lat_0=lat_0,lon_0=lon_0)
            lat_ts=40,lat_0=50,lon_0=20)

xi,yi = m(lons,lats)

# Plot Data
cs = m.pcolor(xi,yi,np.squeeze(o3*2241.15),cmap='jet')

#plt.axis('off')
plt.show()

如您所见，结果以 plt.show 方法呈现。但我想将它保存到 geoTiff。谢谢！

解决方法

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

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

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