|
读取NetCDF文件（***.nc） 回顾：【气候软件】Python读取气象数据 NetCDF文件（***.nc）
本例以一套CMIP6的BCC-CSM2-MR模式输出的全球tas（表面气温）变量1850-2014年的历史格点资料为例：（数据下载地址：www.52ocean.cn） 文件信息查看如下： 0 K% P! L8 c! W
可以看到，文件的坐标有时间time， 经度lon，纬度lat，高度height，变量有表面气温tas。 %读取nc文件用到xarray库包中.open_dataset（）函数； %提取某时间范围内、某经纬度范围内的数据用到.loc[]函数；
( }8 O; ?, ~9 v
, W! i( i/ b. [9 J1 Z. P& P8 E/ V- 2 k1 K- b% m; ~
; w$ `. V9 h" v$ Z9 w9 B4 L- ' C' p$ n& x) y$ `& o( i. W+ ?3 E
- % R; U2 Q8 m! a8 i9 W `
/ c0 z4 \; k5 G# ~/ C7 l
, w1 {, P" Q1 C4 J! _( G, y- + ?; \- m$ Q7 C, a' u; l3 S% h: P+ d
- - f6 n/ r) q( t: q6 y/ Q
8 ~) t! k4 x6 |0 q3 {
#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范围的数据- x3 }1 h: x) C/ S
6 O+ Q; c0 N0 ]# ]! z读取（time时间在15.5至74.5区间内，lat纬度在40°N-55°N，lon精度在115°E-135°E）范围结果显示：
% M B( |- F" w1 @* @; B' I2 E, {( H8 O. m9 {5 J" Z
读取TXT文件
: b4 w1 x7 ^# i& P1 \回顾：【气候软件】Python2：读取TXT文本格式的数据 站点数据是气象数据很重要的一部分，尤其是在观测数据方面，站点数据一般存放于.txt(.csv)文件中，首先就需要读取数据。 以一套中国气象数据网下载的某省的1979-2019年逐月气温数据为例，数据文件格式为TXT文本格式：（数据下载地址：www.52ocean.cn） TXT文本保存的数据如下图所示：
%读取txt文件用到pandas库包中.read_csv（）函数；
9 H" e; ~8 ?* M5 b6 s- y$ {, T- 6 L) |2 o* n5 \* Z" [! M
- : U1 ?! {4 T$ h9 i/ D1 n3 S
- 1 Q7 N8 K/ m( J9 a
+ I4 T, k& F6 {5 }* k- * K- B" i; J6 l. J* h
s5 p3 n: o# E) H5 Z" A5 J: s$ ^% J
3 k3 y/ A! V. I& g+ Y5 w; n Q- R
( M; ^- H/ y1 z k, I5 v- 7 t1 U9 y, \+ \, e5 z0 o& w
: o0 g: I8 k& L7 K3 q
#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. B9 G }: n( Q v
% F- }# F- q% t6 m, s( k
读取结果显示：
可以看到，文件的数据有22列，492行，第0列是年份，第1列是月份，剩余列数是变量数据。 skiprows=1表示跳过了前1行的文件头（因为第一行的数据中含有非数字的字符，所以选择跳过不读取）。 如果其他文件中数据是以空格为间隔，则sep='\s+'定义了数据间隔为空格。 header：设置导入 数据 的列名称，默认为 "infer"，注意它与下面介绍的 names 参数的微妙关系。 names： 3 \% k7 K9 K+ A& p* A
$ b8 Y3 v. I2 r) N# d | 
