netCDF文件后缀名是.nc,是一种二进制的网格文件,无法直接打开,下面给出基于python的netCDF文件生成的脚本。
输入数据说明:
数据为31°N-40°N,130°E-142°E区域表层的地震波速,网格分布大小不均匀。数据存储在txt文件中,分为三列,无表头。
部分数据记录见下图,第一列是经度,第二列是维度,第三列是波速扰动值%
过程与代码
为了将其生成均匀网格的二维nc文件,我们需要1.文件读取,2.数据插值,3.写入nc文件。
1 '''
2 program: ncdata.py3 goal: interpolate txt data and generate grd file4 input: vp1_02.txt5 output: vp1_03.nc6 author: Liang Xuran7 e-mail: xuranliang@hotmail.com8 time: 2020.08.239 '''
10 importnetCDF4 as nc11 importpandas as pd12 importnumpy as np13 importcopy14
15 file_txt = "vp1_02.txt"#txt文件位置
16 file_nc = "vp1_03.nc" #nc文件输出位置
17 npa = 37 #纬向格点数(uninpolate)
18 nra = 45 #经向格点数(uninpolate)
19 npx = 180 #