title: Python 读取和处理nc文件
date: 2021-07-22 23:45:52
tags:
import netCDF4 as nc
import pandas as pd
import numpy as np
file = 'HEMCO_diagnostics.201301010000.nc'
dataset = nc.Dataset(file)
print(dataset.variables.keys())
EmisLEVO_Anthro = dataset.variables['EmisLEVO_Anthro'][:][0,0,:,:]
AREA = dataset.variables['AREA'][:]
EmisLEVO_Anthro_load = EmisLEVO_Anthro*AREA
lon = dataset.variables['lon'][:]
lat = dataset.variables['lat'][:]
long,lati = np.meshgrid(lon,lat)
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')