|
读取NetCDF文件（***.nc） 回顾：【气候软件】Python读取气象数据 NetCDF文件（***.nc）
本例以一套CMIP6的BCC-CSM2-MR模式输出的全球tas（表面气温）变量1850-2014年的历史格点资料为例：（数据下载地址：www.52ocean.cn） 文件信息查看如下： 5 h8 d* B* V! ?- i
可以看到，文件的坐标有时间time， 经度lon，纬度lat，高度height，变量有表面气温tas。 %读取nc文件用到xarray库包中.open_dataset（）函数； %提取某时间范围内、某经纬度范围内的数据用到.loc[]函数； 7 @1 p! [& A$ i' r7 `( j% N
0 k" [, I, H6 y: Z# l( J `4 q- ' u. g) b4 E" ^3 b4 p+ _
- : J& S" c- v3 E# d
- . W8 V! ~6 h3 Z3 l
& G' m. W2 ? C- & ~/ k6 r. K" ^8 S$ i
- * Z5 v6 ]! x+ G0 S/ L
- ! s- F; ~: B/ n9 s$ T& L
- 6 U9 \' u# }3 e2 O/ ~+ E
7 T$ ^5 Y6 @8 W" c a c V
#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范围的数据% ^/ ?, b% s" k9 [" T- |! g3 z
4 ] s. p; O) a2 e1 c
读取（time时间在15.5至74.5区间内，lat纬度在40°N-55°N，lon精度在115°E-135°E）范围结果显示：
- G- W! W' z2 e: S
% G6 L0 ~- C2 f" c9 L! v读取TXT文件
7 |( S% K, N: H7 t, P回顾：【气候软件】Python2：读取TXT文本格式的数据 站点数据是气象数据很重要的一部分，尤其是在观测数据方面，站点数据一般存放于.txt(.csv)文件中，首先就需要读取数据。 以一套中国气象数据网下载的某省的1979-2019年逐月气温数据为例，数据文件格式为TXT文本格式：（数据下载地址：www.52ocean.cn） TXT文本保存的数据如下图所示：
%读取txt文件用到pandas库包中.read_csv（）函数； 3 C* }" P8 c( E2 n8 W4 A/ W4 Q. \
- - d6 S. i+ C$ h7 S4 {
- 1 j, H/ H& }9 q2 ^' V
( v+ g: e9 t3 [3 D& A5 {& E% ?
) N* C5 s/ E6 q" @6 H8 h- # @2 R" H1 n) n
- 6 w' b, F" g6 w! f7 v! I6 T W
& e9 _/ Y. B- y* A5 w% k) U- 4 t3 p, o3 {3 Q
8 P: n r3 n/ H4 j
' o; U: s+ M, @) N- v( T8 r
#1引库import pandas as pd#2数据##skiprows跳过了前1行的文件头##sep=','为文件中数据以逗号间隔；若以空格为间隔修改为sep='\s+'#设置导入 数据 的列名称data = pd.read_csv("tem.txt", skiprows=1, sep=',', header=None) print(data)1 W* c7 [) y8 L
1 h' \6 E) e7 d# L5 V( |7 H& ]* D读取结果显示：
可以看到，文件的数据有22列，492行，第0列是年份，第1列是月份，剩余列数是变量数据。 skiprows=1表示跳过了前1行的文件头（因为第一行的数据中含有非数字的字符，所以选择跳过不读取）。 如果其他文件中数据是以空格为间隔，则sep='\s+'定义了数据间隔为空格。 header：设置导入 数据 的列名称，默认为 "infer"，注意它与下面介绍的 names 参数的微妙关系。 names： 5 f9 G" f( y4 a# a! V0 }
3 E4 O: W( a: L' d! {7 i* `/ |$ p | 
