如何读取多个nc文件并提取出里面的经度和纬度绘制到地图上

在Python中，你可以使用NumPy库处理NetCDF (Network Common Data Form) 文件，并结合Matplotlib或Cartopy库来绘制地图。以下是一个基本步骤：

1. 安装必要的库：

   pip install numpy xarray matplotlib cartopy netcdf4

2. 导入所需模块：

import numpy as np
import xarray as xr
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER

3. 读取多个nc文件：

nc_files = ['file1.nc', 'file2.nc', 'file3.nc'] # 将你的文件名替换这里
datasets = [xr.open_dataarray(fn) for fn in nc_files]
all_coordinates = [(ds.longitude.values, ds.latitude.values) for ds in datasets]

4. 提取经纬度数据：

longitudes = np.concatenate([coord[0] for coord in all_coordinates])
latitudes = np.concatenate([coord[1] for coord in all_coordinates])

5. 绘制地图：

fig, ax = plt.subplots(subplot_kw={'projection': ccrs.PlateCarree()})
ax.set_global()
ax.coastlines(resolution='110m')
ax.stock_img()

# 使用经纬度数据创建点图
ax.scatter(longitudes, latitudes, transform=ccrs.PlateCarree(), s=5)

# 添加网格线标签
gl = ax.gridlines(draw_labels=True,
                  linewidth=0.5, color='gray', alpha=0.5, linestyle='--')
gl.xlabels_top = gl.ylabels_right = False
gl.xformatter = LONGITUDE_FORMATTER
gl.yformatter = LATITUDE_FORMATTER

plt.show()

在这个例子中，我们首先打开每个nc文件并将经纬度数组合并。然后，在地图上使用`scatter()`函数将坐标点绘制成散点图。