python读取nc文件并简单绘制co2浓度等值线图

zhuyidong

1. import matplotlib
   
2. from netCDF4 import Dataset,num2date
   
3. import datetime
   
4. from mpl_toolkits.basemap import Basemap
   
5. import matplotlib.pyplot as plt
   
6. from matplotlib import colors
   
7. import numpy as np
   
8. #让图中中文字体显示方式
   
9. matplotlib.rcParams['font.sans-serif']=['SimHei']   # 用黑体显示中文
   
10. matplotlib.rcParams['axes.unicode_minus']=False     # 正常显示负号
    
11. 
    
12. # 自定义时间段
    
13. month_B=7;day_B=1;month_E=8;day_E=31  #这里选取盛夏
    
14. # 自定义纬度范围
    
15. lat_B=10;lat_E=90;lon_B=30;lon_E=190  #这里选取欧亚
    
16. # 自定义等压面层次
    
17. level=500
    
18. 
    
19. # 打开文件
    
20. filename='E:\PycharmProjects\yiyue\GOSATTFTS2017010120170131_4BP01MV0206.nc'
    
21. fh=Dataset(filename,mode='r')
    
22. # 对绘图变量进行赋值
    
23. lons = fh.variables['lon'][:]
    
24. lats = fh.variables['lat'][:]
    
25. conc = (fh.variables['conc'][1,1,:,:])*100000000
    
26. print(conc)
    
27. # 新建地图并设定经纬度范围
    
28. m = Basemap(llcrnrlat = lat_B, urcrnrlat = lat_E, llcrnrlon = lon_B, urcrnrlon = lon_E)
    
29. # 网格化经纬度并形成坐标矩阵
    
30. lon, lat = np.meshgrid(lons, lats)
    
31. xi, yi = m(lon, lat)
    
32. print(len(xi))
    
33. # 设置等值线范围及间隔
    
34. levels=range(40000,43000,400)
    
35. # 创建填色图层，alpha为透明度0-1，Vmax表示大于该值使用同样的填色以突出副热带高压主体
    
36. cs = m.contourf(xi, yi,conc,levels=levels,vmin=40000,vmax=43000,alpha=1,cmap="RdBu_r")
    
37. # 创建等值线图层
    
38. cs2 = m.contour(xi, yi, conc,levels=levels,colors='black',linewidths=0.5)
    
39. # 绘制经纬线背景网格 label为经纬度标签位置[左,右,上,下]
    
40. m.drawparallels(np.arange(-90., 91., 20.), labels=[1,0,0,0], fontsize=10)
    
41. m.drawmeridians(np.arange(-180., 181., 40.), labels=[0,0,0,1], fontsize=10)
    
42. # 添加海岸线国界等
    
43. m.drawcoastlines()
    
44. m.drawstates()
    
45. m.drawcountries()
    
46. cbar = m.colorbar(cs, location='bottom', pad="10%")
    
47. plt.title('co2浓度',size=20)
    
48. # 执行可视化(弹出窗口)
    
49. plt.savefig('u-shaded.png')
    
50. plt.show()
    
51. 
    

复制代码

qinhantang10

楼主，这数据在哪下载？

zhuyidong

qinhantang10 发表于 2021-4-21 22:17
楼主，这数据在哪下载？

https://ladsweb.modaps.eosdis.nasa.gov/

qinhantang10

能告诉是在哪一栏里下载吗？

熏风习习

你好数据名称叫 ：GOSAT？

zhuyidong

熏风习习 发表于 2021-7-6 17:23
你好数据名称叫 ：GOSAT？

是的、、、、、、、、、、、、、

qinhantang10

熏风习习 发表于 2021-7-6 17:23
你好数据名称叫 ：GOSAT？

能下载  co2吗？

helei0229

zhuyidong 发表于 2021-4-23 08:43
https://ladsweb.modaps.eosdis.nasa.gov/

你好，gosat数据不是h5文件嘛，这个怎么转成nc

zhuyidong

helei0229 发表于 2022-2-3 23:29
你好，gosat数据不是h5文件嘛，这个怎么转成nc

我用python直接读