NetCDF数据的标准化处理（一）

大家好，欢迎来到IT知识分享网。

NetCDF数据的标准化处理（一）

  我们在业务应用中经常会碰到NetCDF(network Common Data Form)网络通用数据格式，相信大家对这种数据格式也比较熟悉了，它是一种自定义的多维科学数据集，在海洋、气象、环保、地球物理等领域应用非常广泛。

一、NetCDF数据结构

二、ArcGIS对NetCDF的数据管理

三、简单的NetCDF标准化处理

• 第一步：读取原始NetCDF文件的维度信息

oldNC=Dataset(oldncfilepath,'r',format='NETCDF4') len_lon=len(oldNC.dimensions["west_east"]) len_lat=len(oldNC.dimensions["south_north"]) len_z=len(oldNC.dimensions["bottom_top"]) len_time=len(oldNC.dimensions["Time"]) 

• 第二步：创建新NetCDF文件的维度

newNC=Dataset(outnewncfilepath,'w',format='NETCDF4') newNC.createDimension("lon",len_lon) newNC.createDimension("lat",len_lat) newNC.createDimension("z",len_z) newNC.createDimension("time",len_time) 

• 第三步：创建新NetCDF文件的各个变量

var_lon=newNC.createVariable("lon","f8",("lon",)) var_lat=newNC.createVariable("lat","f8",("lat",)) var_z=newNC.createVariable("z","f8",("z",)) var_time=newNC.createVariable("time","f8",("time",)) var_pm25=newNC.createVariable("pm25","f8",("time","z","lat","lon",),fill_value=9.96921E36) var_pm10=newNC.createVariable("pm10","f8",("time","z","lat","lon",),fill_value=9.96921E36) var_u10=newNC.createVariable("U10","f8",("time","z","lat","lon",),fill_value=9.96921E36) var_v10=newNC.createVariable("V10","f8",("time","z","lat","lon",),fill_value=9.96921E36) 

• 第四步：定义新NetCDF文件变量的属性信息

var_lon.units="degree_east" var_lon.standard_name = "longitude"; var_lat.units="degree_north" var_lat.standard_name = "latitude"; #var_time.units = 'hours since 2015-01-01 00:00:00' #var_time.calendar = "gregorian" 

• 第五步：给新NetCDF文件的变量赋值

var_lon[:]=oldNC.variables["west_east"][:] var_lat[:]=oldNC.variables["south_north"][:] var_z[:]=numpy.arange(0,len_z,1) var_time[:]=numpy.arange(0,len_time,1) var_pm25[:]=oldNC.variables["PM2_5_DRY"][:] var_pm10[:]=oldNC.variables["PM10"][:] var_u10[:]=oldNC.variables["U"][:][:,:,0:len_lat,0:len_lon] var_v10[:]=oldNC.variables["V"][:][:,:,0:len_lat,0:len_lon] 

• 第六步：关闭NetCDF文件

oldNC.close() newNC.close() 

  上述示例的完整代码已经上传至网盘，链接：https://pan.baidu.com/s/1zSMsgBVOT7RXo4TJrK3n2w
提取码：y4h6，欢迎大家留言，互相交流学习。

  想了解ArcGIS最新的技术动态和环保最新的应用，请关注微信公众号“环保GIS技术与应用”