netCDF中nc
# 读取nc文件并对其中的数据进行处理
import netCDF4 as nc
import pandas as pd
import numpy as np
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
# 2013年1月1日
file = 'NIRv.GPP.198205.v1.nc'
dataset = nc.Dataset(file)
# 查看一些文件中都有哪一些变量
print(dataset.variables.keys())
# 读取经纬度数据lat=46, lon=72
lon = dataset.variables['longitude'][:].data
lat = dataset.variables['latitude'][:].data
GPP = dataset.variables['GPP'][:].data
其中,GPP数据读取时报错invalid literal for int() with base 10: '-9999.0'
GPPdataset的属性为
Out[20]:
<class 'netCDF4._netCDF4.Variable'>
int16 GPP(longitude, latitude)
units: gC m-2 m-1
long_name: Gross Primary Productivity
scaling_factor: 0.001
missing_value: -9999.0
unlimited dimensions:
current shape = (7200, 3600)
filling on, default _FillValue of -32767 used