|
读取NetCDF文件（***.nc） 回顾：【气候软件】Python读取气象数据 NetCDF文件（***.nc）
本例以一套CMIP6的BCC-CSM2-MR模式输出的全球tas（表面气温）变量1850-2014年的历史格点资料为例：（数据下载地址：www.52ocean.cn） 文件信息查看如下： 4 ^4 g. n/ H$ ?9 q/ o# l: _
可以看到，文件的坐标有时间time， 经度lon，纬度lat，高度height，变量有表面气温tas。 %读取nc文件用到xarray库包中.open_dataset（）函数； %提取某时间范围内、某经纬度范围内的数据用到.loc[]函数；
" ?4 q/ H b: F. A
' K/ {% |- T) ^% `
( v6 A. g# M' n1 T: c0 A
' I6 j! g% _2 m. |6 q; k# G9 b% H: s# o. C
5 j) R; G' G/ |2 O& w
5 C: U4 I4 y, |1 K) s) L- + S2 l) C' F( R5 F- G! o6 N
- # S# S4 J2 J4 [; r1 _0 I1 Z
2 R0 }0 Q& V( B f
( h: g9 X7 t) C4 @/ j1 i
' V/ m8 K) d- k3 ?8 K* g
#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范围的数据
" w+ W. H9 Q2 Z0 t/ c- |. ~5 s + E5 h; ^. o& [; C: D+ \2 l' t8 [8 m
读取（time时间在15.5至74.5区间内，lat纬度在40°N-55°N，lon精度在115°E-135°E）范围结果显示：
' g+ Z1 {: T5 d# c
8 m( e( F2 ]5 ]$ F' \: I6 r( W( H读取TXT文件
6 h3 i; R8 S: f5 G& D2 W回顾：【气候软件】Python2：读取TXT文本格式的数据 站点数据是气象数据很重要的一部分，尤其是在观测数据方面，站点数据一般存放于.txt(.csv)文件中，首先就需要读取数据。 以一套中国气象数据网下载的某省的1979-2019年逐月气温数据为例，数据文件格式为TXT文本格式：（数据下载地址：www.52ocean.cn） TXT文本保存的数据如下图所示：
%读取txt文件用到pandas库包中.read_csv（）函数；
, U% x9 Z. Y! u9 B z+ H7 h
* z8 \. O5 I1 T% d- P) {* Z( S% p/ K; X1 o
+ R8 F4 G- x- r; S9 q3 j- ; y2 M1 X. G8 g; m* t q6 W
7 g- H+ L) x% |6 n n/ @
! \' C% e5 P' ]
$ O" o3 b9 Y. A% {. v- : G4 L' n m5 N+ l
- 6 n- j" k v: j. i7 B
" g v# C3 n9 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)
; l* m, z1 H1 p, }/ R y
, G8 e8 p A$ B9 a读取结果显示：
可以看到，文件的数据有22列，492行，第0列是年份，第1列是月份，剩余列数是变量数据。 skiprows=1表示跳过了前1行的文件头（因为第一行的数据中含有非数字的字符，所以选择跳过不读取）。 如果其他文件中数据是以空格为间隔，则sep='\s+'定义了数据间隔为空格。 header：设置导入 数据 的列名称，默认为 "infer"，注意它与下面介绍的 names 参数的微妙关系。 names：
. ]. c, P1 F# X, b) N% e8 V$ D( n
! L1 E7 y" b8 M6 Q7 M | 
