Python读取nc文件数据

---

title: Python 读取和处理nc文件
date: 2021-07-22 23:45:52
tags:

# 读取nc文件并对其中的数据进行处理
import netCDF4 as nc
import pandas as pd
import numpy as np
# 2013年1月1日
file = 'HEMCO_diagnostics.201301010000.nc'
dataset = nc.Dataset(file)

# 查看一些文件中都有哪一些变量
print(dataset.variables.keys())

# 读取人为排放的左旋葡萄糖的数据,单位为kg/m2/s，空间分辨率为5*5
# 数据是4维的，分别是时间，高度，纬度，经度
EmisLEVO_Anthro = dataset.variables['EmisLEVO_Anthro'][:][0,0,:,:]
# 读取网格面积数据units = "m2"
AREA = dataset.variables['AREA'][:]

# 网格面积*排放强度，单位为kg/s
EmisLEVO_Anthro_load = EmisLEVO_Anthro*AREA

# 读取经纬度数据lat=46, lon=72
lon = dataset.variables['lon'][:]
lat = dataset.variables['lat'][:]
long,lati = np.meshgrid(lon,lat)

# 下面将数据输出为csv数据
long = long.reshape(46*72)
lati = lati.reshape(46*72)
EmisLEVO_Anthro_load = EmisLEVO_Anthro_load.reshape(46*72)

total = {'lon':long,'lat':lati,'EmisLEVO_Anthro_load':EmisLEVO_Anthro_load}
total = pd.DataFrame(total)
total.to_csv('EmisLEVO_Anthro_load.csv')