绘制 WRF（Weather Research and Forecasting）模型的模拟区域或域（domain）通常涉及以下几个关键步骤： 1. **读取 WRF 输出文件**：WRF 模型输出通常为 NetCDF 格式文件，可以使用 `netCDF4` 或 `xarray` 读取。 2. **提取地理信息**：从 WRF 输出中获取经纬度信息，通常可以通过 `wrf-python` 提供的 `wrf.latlon_coords()` 函数提取。 3. **设置地图投影**：使用 `Cartopy` 或 `Basemap` 设置合适的地图投影，如 `PlateCarree`、`LambertConformal` 等。 4. **绘制地理范围**：使用 `matplotlib` 和地图库绘制 WRF 模拟域的边界、地形、海岸线等。 ### 使用 Cartopy 绘制 WRF 模拟域 Cartopy 是一个现代的地图投影和地理可视化库，推荐用于 WRF 数据的可视化。 ```python import matplotlib.pyplot as plt import cartopy.crs as ccrs import cartopy.feature as cfeature import xarray as xr from wrf import getvar, latlon_coords # 读取 WRF 输出文件 wrf_file = xr.open_dataset("wrfout_d01_2023-01-01_00:00:00") # 提取变量（例如温度）用于绘图 temp = getvar(wrf_file, "temp", timeidx=0) # 获取经纬度坐标 lats, lons = latlon_coords(temp) # 创建地图投影 projection = ccrs.LambertConformal(central_longitude=-95, central_latitude=35) # 创建绘图 fig = plt.figure(figsize=(10, 8)) ax = plt.axes(projection=projection) # 设置地图范围 ax.set_extent([-120, -70, 25, 50], crs=ccrs.PlateCarree()) # 添加地理要素 ax.add_feature(cfeature.COASTLINE) ax.add_feature(cfeature.BORDERS, linestyle=':') ax.add_feature(cfeature.STATES, linestyle='--') # 绘制数据 contour = ax.contourf(lons, lats, temp, transform=ccrs.PlateCarree(), cmap='viridis') plt.colorbar(contour, ax=ax, orientation='horizontal', pad=0.04) # 添加标题 ax.set_title('WRF Domain Simulation Area (Cartopy)') plt.show() ``` ### 使用 Basemap 绘制 WRF 模拟域 虽然 `Basemap` 已不再积极维护，但仍然广泛使用于 WRF 可视化。 ```python import matplotlib.pyplot as plt from mpl_toolkits.basemap import Basemap import xarray as xr from wrf import getvar, latlon_coords # 读取 WRF 输出文件 wrf_file = xr.open_dataset("wrfout_d01_2023-01-01_00:00:00") # 提取变量和经纬度 temp = getvar(wrf_file, "temp", timeidx=0) lats, lons = latlon_coords(temp) # 创建地图 m = Basemap(projection='lcc', resolution='l', llcrnrlat=lats.min(), urcrnrlat=lats.max(), llcrnrlon=lons.min(), urcrnrlon=lons.max(), lat_0=35, lon_0=-95) # 创建绘图 fig, ax = plt.subplots(figsize=(10, 8)) m.ax = ax # 绘制地理要素 m.drawcoastlines(linewidth=0.5) m.drawcountries(linewidth=0.5) m.drawstates(linewidth=0.3) # 绘制数据 x, y = m(lons, lats) contour = m.contourf(x, y, temp, cmap='viridis') plt.colorbar(contour, orientation='horizontal', pad=0.04) # 添加标题 ax.set_title('WRF Domain Simulation Area (Basemap)') plt.show() ``` ### 注意事项 - **WRF 域的层级**：在多域模拟中，可以通过读取 `wrfout_d01`, `wrfout_d02` 等文件来绘制不同层级的域。 - **投影设置**：确保地图投影与 WRF 模型使用的投影一致，例如 Lambert Conformal 投影通常用于 WRF 模拟。 - **性能优化**：处理大规模数据时，建议使用 `dask` 或 `xarray` 的并行功能以提高性能[^2]。