-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathnc文件读取.py
68 lines (59 loc) · 1.9 KB
/
nc文件读取.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
# -*- coding: utf-8 -*-
from netCDF4 import Dataset
nc_obj = Dataset(r'C:\Users\DELL\Downloads\CMS_Flux_Fire_2010_v2.nc')
# 查看nc文件有些啥东东
print(nc_obj)
print('---------------------------------------')
# 查看nc文件中的变量,结果是:['lon', 'lat', 'lev', 'time', 'tem']
print(nc_obj.variables.keys())
for i in nc_obj.variables.keys():
print(i)
print('---------------------------------------')
# 查看每个变量的信息
print(nc_obj.variables['lat'])
print(nc_obj.variables['lon'])
print(nc_obj.variables['lev'])
print(nc_obj.variables['time'])
print(nc_obj.variables['tem'])
print('---------------------------------------')
# 查看每个变量的属性
# print(nc_obj.variables['lat'].ncattrs())
# print(nc_obj.variables['lon'].ncattrs())
# print(nc_obj.variables['lev'].ncattrs())
# print(nc_obj.variables['time'].ncattrs())
# print(nc_obj.variables['tem'].ncattrs())
#
# print(nc_obj.variables['lat'].units)
# print(nc_obj.variables['lon'].units)
# print('---------------------------------------')
# 读取数据值
lat = (nc_obj.variables['lat'][:])
lon = (nc_obj.variables['lon'][:])
lev = (nc_obj.variables['lev'][:])
# print(lev)
time = (nc_obj.variables['time'][:])
tem = (nc_obj.variables['tem'][:])
# lat = (nc_obj.variables['lat'][1:2])
# lon = (nc_obj.variables['lon'][1:2])
# lev = (nc_obj.variables['lev'][1:2])
# time1 = (nc_obj.variables['time'][1:2])
# tem = (nc_obj.variables['tem'][1:2])
# print(lat)
# print(len(lat))
# print(lon)
# print(lev)
# print(time1)
# print(len(time))
# print(tem)
# print('---------------******-------------------')
file = open('ppm_lat.txt', 'a')
file.write('lat,lon,time,tem' + '\n')
file.write('lat,lon' + '\n')
file.write('lat' + '\n')
for i in range(len(lat)):
file.write(str(lat[i]) + ',\n')
file.write(str(lon[i]) + ',')
file.write(str(lev[i]) + ',')
file.write(str(time[i]) + ',')
file.write(str(tem[i]) + ',\n')
file.close()