读取NetCDF文件(***.nc) 回顾:【气候软件】Python读取气象数据 NetCDF文件(***.nc)
本例以一套CMIP6的BCC-CSM2-MR模式输出的全球tas(表面气温)变量1850-2014年的历史格点资料为例:(数据下载地址:www.52ocean.cn) 文件信息查看如下:
$ H! D4 J- L7 N% ?( F Y
可以看到,文件的坐标有时间time, 经度lon,纬度lat,高度height,变量有表面气温tas。 %读取nc文件用到xarray库包中.open_dataset()函数; %提取某时间范围内、某经纬度范围内的数据用到.loc[]函数; 1 ?7 {# W9 X9 a; \: @. D
- 6 V! d. D" q; q- D
- : i7 w5 v5 n& q$ p$ r* Z' _# B9 v7 G
( D2 j$ u* k) I- k1 P
0 ^' a% x# }0 n9 C% Z! c" a- & ^& r% R* }$ H% U! j( U
- ( ~ b5 t1 D" v( C
- " Y/ ]7 a, Y6 ~9 N
; r6 s* H7 S, b* ]8 D- n- : j# y& F/ t3 ]) g# G5 c, B
& J; I% l3 M3 R4 ~6 z7 Y
#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范围的数据) ]2 d) A6 V) _1 e
, |2 w( k. F! i, ?& z) J* i9 k
读取(time时间在15.5至74.5区间内,lat纬度在40°N-55°N,lon精度在115°E-135°E)范围结果显示:
( h& P% y0 {% H, j
: y( Y3 V( x5 A1 l7 @; [# z# ^读取TXT文件
' [2 p" s0 f: ~2 ~1 M( T- k1 t! v) \8 j回顾:【气候软件】Python2:读取TXT文本格式的数据 站点数据是气象数据很重要的一部分,尤其是在观测数据方面,站点数据一般存放于.txt(.csv)文件中,首先就需要读取数据。 以一套中国气象数据网下载的某省的1979-2019年逐月气温数据为例,数据文件格式为TXT文本格式:(数据下载地址:www.52ocean.cn) TXT文本保存的数据如下图所示:
%读取txt文件用到pandas库包中.read_csv()函数;
8 B' M( V4 c! z- P2 `) R- 5 }+ _6 S( T% Z" I6 m; I
0 k7 f! Z0 ^% f9 [ g! e
& b2 m! C$ \3 h# D7 _ d% j- ' [8 @2 T3 t% @5 P: E. | F9 s
- + {2 p1 l- d5 P$ o
- 2 b0 b/ ?! e4 M5 R3 E t0 x3 ~
7 o4 [( @0 [. P+ x8 x5 V- 1 X9 i& L6 Z' H! I* ?) @
) a e! I4 G" O& Z
4 B$ N; r1 O& F$ b( L9 c+ I# Z& w
#1引库import pandas as pd#2数据##skiprows跳过了前1行的文件头##sep=','为文件中数据以逗号间隔;若以空格为间隔修改为sep='\s+'#设置导入 数据 的列名称data = pd.read_csv("tem.txt", skiprows=1, sep=',', header=None) print(data)
0 k+ B$ \/ K3 q5 _+ S/ s 1 x. C( }# q1 D* P) \( L
读取结果显示:
可以看到,文件的数据有22列,492行,第0列是年份,第1列是月份,剩余列数是变量数据。 skiprows=1表示跳过了前1行的文件头(因为第一行的数据中含有非数字的字符,所以选择跳过不读取)。 如果其他文件中数据是以空格为间隔,则sep='\s+'定义了数据间隔为空格。 header:设置导入 数据 的列名称,默认为 "infer",注意它与下面介绍的 names 参数的微妙关系。 names:
8 g% U* H3 V p) H8 I1 `) b# [
! T; a3 J+ Y2 N. M* i4 f |