# -*- coding: utf-8 -*-import sys
import os
import netCDF4 as nc
import numpy as np
"""
NC数据文件读写操作
"""# 打开nc文件defNCopen(ncfile):# load in the netCDF4 file
ncdata = nc.Dataset(ncfile)return ncdata
# read the variables from the netCDF4 file and assign #获取数据数组1defReadNcDataset1(ncdata,datasetPath):if datasetPath in ncdata.variables.keys():
dataset = ncdata.variables[datasetPath]
dataset1 = np.array(dataset)return dataset1
#获取数据数组2defReadNcDataset2(ncdata,datasetPath):if datasetPath in ncdata.variables.keys():
dataset = ncdata.variables[datasetPath][:,:]
dataset1 = np.array(dataset)return dataset1
#获取数据数组3defReadNcDataset3(ncdata,datasetPath):if datasetPath in((ncdata.groups)["PRODUCT"]).variables.keys():
dataset =((ncdata.groups)["PRODUCT"]).variables[datasetPath]
dataset1 = np.array(dataset)return dataset1
# 关闭nc文件defNCclose(ncdata):
ncdata.close()
# -*- coding: utf-8 -*-import os
import time
import glob
import numpy as np
from osgeo import gdal, osr
from scipy import ndimage as nd
import NC_Util
def