|
读取NetCDF文件（***.nc） 回顾：【气候软件】Python读取气象数据 NetCDF文件（***.nc）
本例以一套CMIP6的BCC-CSM2-MR模式输出的全球tas（表面气温）变量1850-2014年的历史格点资料为例：（数据下载地址：www.52ocean.cn） 文件信息查看如下： . [1 m% F) z: v
可以看到，文件的坐标有时间time， 经度lon，纬度lat，高度height，变量有表面气温tas。 %读取nc文件用到xarray库包中.open_dataset（）函数； %提取某时间范围内、某经纬度范围内的数据用到.loc[]函数； 6 u' }9 p% p$ X
- 0 M# t+ q& M5 h7 E7 x( H
+ d8 H) }1 J2 O. O T4 [& S
1 q5 @$ C6 Z: |3 K0 t
( {9 w( i5 \- D' B) K8 K$ M9 h6 t/ |- 0 n9 s; N+ _ ~, t. x4 o
- 4 o% o& h8 O: n- u5 b% c7 R% D8 M
- & g" G( `9 H/ D- r& B1 O! b) f
5 F3 B6 |/ T; f
1 |* ^1 F1 n& G) L/ R; k% A% t9 W' L) y9 j
#1引库import xarray as xr#2数据（数据免费下载后台输入“BCC-CSM2-MR数据”回复获取）data = xr.open_dataset('tas_Amon_BCC-CSM2-MR_historical_r1i1p1f1_gn_185001-201412.nc', decode_times=False)print(data) #输出显示数据信息，便于了解tas = data['tas'] #提取变量tas数据 a = tas.loc[15.5:74.5] # 提取time在15.5至74.5区间内的数据b = tas.loc[15.5:74.5, 40:55, 115:135] #选取了lat在40°N-55°N，lon在115°E-135°E范围的数据% T' `, }0 T3 A; F& S2 [: h
7 H; B! X6 [7 L% g% D读取（time时间在15.5至74.5区间内，lat纬度在40°N-55°N，lon精度在115°E-135°E）范围结果显示：
7 Z) n, V" q/ G) r8 C1 F. ^2 t- R* h- t& {3 u
读取TXT文件 $ z2 \% w( l0 x
回顾：【气候软件】Python2：读取TXT文本格式的数据 站点数据是气象数据很重要的一部分，尤其是在观测数据方面，站点数据一般存放于.txt(.csv)文件中，首先就需要读取数据。 以一套中国气象数据网下载的某省的1979-2019年逐月气温数据为例，数据文件格式为TXT文本格式：（数据下载地址：www.52ocean.cn） TXT文本保存的数据如下图所示：
%读取txt文件用到pandas库包中.read_csv（）函数；
5 Y' B) G( ^" Q, F0 N% j- ! b0 B/ q. a* `9 B1 o
% z" `3 o+ L. C* t) \/ X
6 f- Y1 e- y; l. W- R0 D- 4 F6 A) A @3 s8 G" o0 c9 R/ F
4 C( x5 i3 t. j- g* I/ v- ( o3 [3 e7 S* A+ Q
- . A8 z5 {4 J, y& U& v& B( M7 _
7 w8 C3 h y7 q& B+ x' [) f- $ p5 u- s1 e# o6 Y! k
# U* f- |0 G' g; p6 G: @
#1引库import pandas as pd#2数据##skiprows跳过了前1行的文件头##sep=','为文件中数据以逗号间隔；若以空格为间隔修改为sep='\s+'#设置导入 数据 的列名称data = pd.read_csv("tem.txt", skiprows=1, sep=',', header=None) print(data). s \7 b! ~( j4 G: W: d
9 A, E' I" h' t- b8 }
读取结果显示：
可以看到，文件的数据有22列，492行，第0列是年份，第1列是月份，剩余列数是变量数据。 skiprows=1表示跳过了前1行的文件头（因为第一行的数据中含有非数字的字符，所以选择跳过不读取）。 如果其他文件中数据是以空格为间隔，则sep='\s+'定义了数据间隔为空格。 header：设置导入 数据 的列名称，默认为 "infer"，注意它与下面介绍的 names 参数的微妙关系。 names：
% `3 }( \9 L- `! H9 I2 U3 x9 @' }& o8 W- m
| 
