python 绘制nc格式源清单

最新推荐文章于 2025-05-26 15:24:48 发布

跺碎那片云

最新推荐文章于 2025-05-26 15:24:48 发布

阅读量2.2k

点赞数

CC 4.0 BY-SA版权

分类专栏： Python

本文链接：https://blog.youkuaiyun.com/shaw_1989/article/details/78682965

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

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#!/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, 13, 14, 15]:
for i in [0,]:
    print(u"第一个是", i)
    no2_0 = no2[i:i+1, 0:1, :, :]
    no2_0 = np.squeeze(no2_0)
    clevs = np.arange(0,1.3,