python 绘制nc格式源清单

最新推荐文章于 2024-06-24 17:59:43 发布

跺碎那片云

最新推荐文章于 2024-06-24 17:59:43 发布

阅读量2.1k

点赞数

分类专栏： Python

本文链接：https://blog.csdn.net/shaw_1989/article/details/78682965

版权

该博客通过Python使用matplotlib和netCDF4库展示了如何读取nc格式文件，并绘制NO2排放率的地图。代码中首先定义了地图范围，然后从文件中提取NO2数据，使用Basemap创建地图并进行颜色填充，最后添加了颜色条、海岸线、州和国家边界等元素。

摘要由CSDN通过智能技术生成

#!/usr/bin/python
#coding=utf-8

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap, cm
from netCDF4 import Dataset
import sys
import cmaps

lon_0 =  116.75
lat_0 =  37.25

filename = './emiss_27.ncf'
print(filename)
readfile = Dataset(filename, mode='r', open=True)
no2    = np.squeeze(readfile.variables['NO2'][:][:][:])
no2_units = readfile.variables['NO2'].units
#no2 = readfile.variables['NO2'][:]

m = Basemap(lat_0=lat_0, lon_0=lon_0, llcrnrlat=24.75 ,urcrnrlat=49.75,\
            llcrnrlon=104.25,urcrnrlon=129.25,\
            rsphere=6371200.,resolution='l')

lons = np.arange(104.25,129.25,0.25)
lats = np.arange(24.75,49.75,0.25)

lon, lat = np.meshgrid(lons, lats)
xi, yi = m(lon, lat)

######  提取NO 网格信息COL,ROW
#for i in [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12,