Python 使用netCDF4读写nc文件以及截取指定经纬度范围内的数据生成新的nc文件

使用Python中的netCDF4库可以很方便地读写nc文件，并且可以根据需要截取指定经纬度范围内的数据生成新的nc文件。以下是一个示例代码：

from netCDF4 import Dataset
import numpy as np

# 读取nc文件
nc_file = Dataset("file_path.nc", mode="r")

# 获取经纬度和变量数据
lons = nc_file.variables["lon"][:]
lats = nc_file.variables["lat"][:]
data = nc_file.variables["variable_name"][:]

# 获取经纬度范围
lon_range = [lon_min, lon_max]  # 经度范围
lat_range = [lat_min, lat_max]  # 纬度范围

# 截取指定经纬度范围内的数据
lon_index = np.where((lons >= lon_range[0]) & (lons <= lon_range[1]))[0]
lat_index = np.where((lats >= lat_range[0]) & (lats <= lat_range[1]))[0]
data_subset = data[lat_index[0]:lat_index[-1]+1, lon_index[0]:lon_index[-1]+1]

# 新建nc文件并写入数据
new_nc_file = Dataset("new_file_path.nc", mode="w", format="NETCDF4")
new_nc_file.createDimension("lon", len(lon_index))
new_nc_file.createDimension("lat", len(lat_index))
new_lon = new_nc_file.createVariable("lon", np.float32, ("lon",))
new_lon[:] = lons[lon_index]
new_lat = new_nc_file.createVariable("lat", np.float32, ("lat",))
new_lat[:] = lats[lat_index]
new_data = new_nc_file.createVariable("variable_name", np.float32, ("lat", "lon"))
new_data[:] = data_subset
new_nc_file.close()

# 关闭nc文件
nc_file.close()

在这里，`nc_file`是原始nc文件。然后，代码获取了nc文件中的经纬度和变量数据。接着，通过指定的经纬度范围截取需要的数据，并新建一个nc文件并写入数据。新建nc文件的过程中，需要先创建新的经纬度维度，并将截取到的经纬度数据写入到新的nc文件中。最后，将截取到的变量数据写入到新的nc文件。注意，在使用完nc文件后，需要将其关闭。