|
读取NetCDF文件(***.nc) 回顾:【气候软件】Python读取气象数据 NetCDF文件(***.nc)
本例以一套CMIP6的BCC-CSM2-MR模式输出的全球tas(表面气温)变量1850-2014年的历史格点资料为例:(数据下载地址:www.52ocean.cn) 文件信息查看如下: ; p6 {8 P* Q' h. I" N/ c( m
可以看到,文件的坐标有时间time, 经度lon,纬度lat,高度height,变量有表面气温tas。 %读取nc文件用到xarray库包中.open_dataset()函数; %提取某时间范围内、某经纬度范围内的数据用到.loc[]函数; * p8 d- G1 q/ f2 U$ `# n6 T1 ]# L5 j
) _7 v; T+ r% {9 j
# M+ S+ `( h3 x9 V- A/ g' C' L+ t
- - ]: ^, R) G3 a6 G' L. l
/ `9 F3 a3 M% Q" n' q- k$ U- / ^$ r* r# J0 W4 X
- ( E' D- C" J4 ]2 @ d8 P
- S/ g! O% n7 R, m1 q
8 @8 a5 P/ Y ]3 v$ h) X
/ g, M) H' F: U! b B! {
#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范围的数据
: H4 g# P9 x* L; W; a4 C $ t) N" J/ _; c- ]7 f& [" B0 \" y
读取(time时间在15.5至74.5区间内,lat纬度在40°N-55°N,lon精度在115°E-135°E)范围结果显示:
) D* H. I5 o4 p; B
9 d$ y/ F( |; g: l5 w
读取TXT文件 ' d. G# b8 c7 i
回顾:【气候软件】Python2:读取TXT文本格式的数据 站点数据是气象数据很重要的一部分,尤其是在观测数据方面,站点数据一般存放于.txt(.csv)文件中,首先就需要读取数据。 以一套中国气象数据网下载的某省的1979-2019年逐月气温数据为例,数据文件格式为TXT文本格式:(数据下载地址:www.52ocean.cn) TXT文本保存的数据如下图所示:
%读取txt文件用到pandas库包中.read_csv()函数; + e" ]5 e2 k9 ?" Z3 z& P" g, z* n
- : p$ K V( y2 E, L- S
- 9 `) U( w f" B5 v: l& |1 ~
- # c) G7 a; e9 u: v
: j0 x! _' S- K/ e' x7 h- ; v2 c( }: t* x: |5 _( U! \
- * C# v, o5 W1 O8 v! N; _! j
3 V' Y' Q2 p# B. i5 f
" x( Q( R, u7 _; k6 o$ u) E# R
- I: A7 D/ ?5 M; }& x
, N) X; Y& q; \) R. P
#1引库import pandas as pd#2数据##skiprows跳过了前1行的文件头##sep=','为文件中数据以逗号间隔;若以空格为间隔修改为sep='\s+'#设置导入 数据 的列名称data = pd.read_csv("tem.txt", skiprows=1, sep=',', header=None) print(data)
' q0 l. K1 a2 G8 M ' i; O+ u9 P* y7 P; G# G; [
读取结果显示:
可以看到,文件的数据有22列,492行,第0列是年份,第1列是月份,剩余列数是变量数据。 skiprows=1表示跳过了前1行的文件头(因为第一行的数据中含有非数字的字符,所以选择跳过不读取)。 如果其他文件中数据是以空格为间隔,则sep='\s+'定义了数据间隔为空格。 header:设置导入 数据 的列名称,默认为 "infer",注意它与下面介绍的 names 参数的微妙关系。 names: ( |1 S* m) x# l- l
8 O9 S$ E+ n* R& @, p2 M! ~
|