python 创建netcdf_python写入二维netCDF文件（txt转nc文件）

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 #