对气象文件的数据提取与处理，包含.nc文件与.csv文件

 

import netCDF4 as nc
import sys
import glob
import csv

#读取t文件夹下的nc文件的北京温度值  t 气温

def t_getcsv(path):
    sites1_t_data = glob.glob(path + 'data/t2015/*.nc')
    sites1_t_data.sort(key=lambda x: float(x[-14:-3]))
    print('t',sites1_t_data, len(sites1_t_data))
    beij1_t = []
    for i in sites1_t_data:
        nc_obj = nc.Dataset(i)
        lev = list(nc_obj.variables['level'][:])
        print(lev)
        UU = nc_obj.variables['lon'][151:160]
        NN = nc_obj.variables['lat'][63:72]
        #print(UU,'\n',NN)

        lev700 = lev.index(700)
        lev850 = lev.index(850)
        lev925 = lev.index(925)

        for c in [lev700, lev850, lev925]:
            cc = []
            hour = int(i[-5:-3])+7
            if hour > 24:
                q = sites1_t_data.index(i)
                q1 = i
                if q < len(sites1_t_data)-1:
                    q1 = sites1_t_data[q + 1]

                cc.append(q1[-14:-5]+'01')
            elif hour ==7:
                cc.append(i[-14:-5] +'07')
            else:
                cc.append(i[-14:-5] + str(hour))

            c = int(c)
            cc.append(lev[c])
            for m in range(151, 161, 1):
                #mm = nc_obj.variables['lon'][m]
                m = int(m)
                for n in range(63, 73, 1):
                    # nn = nc_obj.variables['lat'][n]
                    n = int(n)
                    # cc.extend([m, n])

                    # print(m,n,c)
                    hh = nc_obj.variables['t'][c][n][m]
                    cc.append(hh)
            print('t','\n', cc,)
            beij1_t.append(cc)

            '''if c == lev700:
                data700.extend(cc[2:])
            elif c == lev850:
                data850.extend(cc[2:])
            else: