保姆级教程：用Python和Basemap绘制台风‘利奇马’期间的卫星云图（附完整代码）

用Python和Basemap绘制台风卫星云图实战指南去年夏天当超强台风利奇马逼近华东沿海时气象部门发布了大量卫星云图数据。这些看似简单的彩色图片背后其实隐藏着台风的强度、移动路径等关键信息。作为气象爱好者或相关领域的研究者你是否想过自己动手处理这些数据生成专业的卫星云图本文将带你从零开始使用Python和Basemap库处理FY-2G气象卫星的HDF数据最终生成可用于台风分析的亮温(TBB)云图。1. 环境准备与数据获取在开始之前我们需要准备好Python环境和必要的数据。不同于普通的图像处理气象卫星数据有其特殊的格式和处理方式。1.1 安装必要的Python库处理气象卫星数据需要几个关键库pip install numpy h5py matplotlib basemap注意Basemap库在某些系统上可能需要通过conda安装conda install -c anaconda basemap1.2 获取FY-2G卫星数据中国风云二号(FY-2G)静止气象卫星的全圆盘数据可以从国家卫星气象中心获取。对于利奇马台风我们需要以下两种文件HDF数据文件包含红外通道的原始观测数据经纬度对照表用于将卫星数据映射到地理坐标典型的数据文件命名格式为FY2G_FDI_ALL_NOM_YYYYMMDD_HHMM.hdf其中包含了特定时间的观测数据。2. 理解卫星数据与亮温(TBB)气象卫星的红外通道测量的是地气系统的辐射亮度通过定标可以转换为亮温(Brightness Temperature)。亮温是分析台风强度的重要参数。2.1 HDF文件结构解析FY-2G的HDF文件包含多个数据集我们需要重点关注数据集路径描述/CALChannelIR1红外通道1的定标表/NOMChannelIR1红外通道1的原始观测值/NomFileInfo文件元信息包括卫星位置2.2 亮温计算原理亮温的计算过程可以概括为读取原始观测值(DN值)使用定标表将DN值转换为辐射亮度通过普朗克公式反算得到亮温以下是关键的Python实现代码def calculate_brightness_temperature(hdf_file): with h5py.File(hdf_file, r) as hdf: calibration hdf[/CALChannelIR1][:] observation hdf[/NOMChannelIR1][:] # 初始化亮温数组 tbb np.zeros_like(observation, dtypenp.float32) # 逐像素计算亮温 for i in range(observation.shape[0]): for j in range(observation.shape[1]): dn observation[i,j] if dn not in [65534, 65535]: # 排除无效值 tbb[i,j] calibration[0, dn] return tbb3. 地理坐标转换与地图投影卫星数据的坐标系统与常规地理坐标不同需要进行转换才能正确显示在地图上。3.1 读取经纬度对照表FY-2G卫星提供了专门的经纬度对照文件(.dat)用于将图像像素映射到地理坐标def read_geolocation(dat_file): with open(dat_file, rb) as f: data np.fromfile(f, dtypefloat32) # 风云二号卫星数据特有的存储格式 data data.reshape([2288, 2288, 2], orderF) longitude data[:,:,0] 104.5 # 卫星星下点经度修正 latitude data[:,:,1] return longitude, latitude3.2 地图投影设置对于台风分析我们通常使用圆柱投影(Plate Carrée)来保持区域形状m Basemap(projectioncyl, llcrnrlat15, urcrnrlat45, # 纬度范围 llcrnrlon105, urcrnrlon145, # 经度范围 resolutioni)4. 完整云图绘制流程现在我们将所有步骤整合起来绘制台风利奇马的卫星云图。4.1 数据处理流程读取HDF文件获取亮温数据读取.dat文件获取地理坐标设置地图投影参数绘制海岸线、国界等地理要素将亮温数据投影到地图上添加色标、标题等辅助元素4.2 完整实现代码import numpy as np import h5py from mpl_toolkits.basemap import Basemap import matplotlib.pyplot as plt def plot_typhoon_tbb(hdf_path, dat_path, output_png): # 读取亮温数据 tbb calculate_brightness_temperature(hdf_path) # 读取地理坐标 lon, lat read_geolocation(dat_path) # 创建地图 plt.figure(figsize(12, 10)) m Basemap(projectioncyl, llcrnrlat15, urcrnrlat45, llcrnrlon105, urcrnrlon145, resolutioni) # 绘制地理要素 m.drawcoastlines(linewidth0.5) m.drawcountries(linewidth0.5) m.drawparallels(np.arange(15, 46, 5), labels[1,0,0,0]) m.drawmeridians(np.arange(105, 146, 5), labels[0,0,0,1]) # 投影亮温数据 x, y m(lon, lat) cs m.contourf(x, y, tbb, levelsnp.linspace(180, 300, 121), cmapjet) # 添加色标 cbar m.colorbar(cs, locationright, pad5%) cbar.set_label(Brightness Temperature (K)) # 添加标题 plt.title(Typhoon Lekima - FY2G IR1 Brightness Temperature\n2019-08-10 12:00 UTC, fontsize14, pad20) # 保存图像 plt.savefig(output_png, dpi300, bbox_inchestight) plt.close()5. 台风分析进阶技巧生成卫星云图只是第一步更重要的是从中提取有用的气象信息。5.1 识别台风眼在亮温图像中台风眼通常表现为小而圆的温暖区域(亮温较高)周围被极低温的云墙包围直径通常在30-60公里之间我们可以通过寻找局部温度最大值来定位台风眼def find_eye_center(tbb, threshold280): # 寻找高于阈值的区域 warm_area (tbb threshold) # 计算质心 y, x np.indices(tbb.shape) center_x np.mean(x[warm_area]) center_y np.mean(y[warm_area]) return center_x, center_y5.2 估算台风强度虽然精确的强度估算需要更复杂的方法但亮温可以提供初步判断亮温范围(K)对应的云顶高度可能的对流强度 200极高极强200-220高强220-240中高中等 240低弱5.3 多时次动画制作分析台风移动和强度变化时动画比单张图片更有说服力。我们可以使用Matplotlib的FuncAnimation创建动画from matplotlib.animation import FuncAnimation def create_animation(hdf_files, dat_file): fig plt.figure(figsize(10, 8)) m setup_basemap() def update(frame): tbb calculate_brightness_temperature(hdf_files[frame]) lon, lat read_geolocation(dat_file) x, y m(lon, lat) plt.clf() m.contourf(x, y, tbb, levelsnp.linspace(180, 300, 121), cmapjet) m.drawcoastlines() anim FuncAnimation(fig, update, frameslen(hdf_files), interval200) anim.save(typhoon_animation.mp4, writerffmpeg)6. 常见问题与解决方案在实际操作中可能会遇到各种问题。以下是几个常见问题及其解决方法6.1 数据读取问题问题HDF文件无法打开或数据读取错误可能原因文件路径错误文件损坏库版本不兼容解决方案try: with h5py.File(hdf_path, r) as hdf: print(文件打开成功包含的数据集, list(hdf.keys())) except Exception as e: print(f文件读取失败{str(e)})6.2 地图投影问题问题数据在地图上的位置不正确检查步骤确认卫星星下点经度修正值是否正确检查地图范围是否包含数据区域验证经纬度数据的符号东经/西经北纬/南纬6.3 图像显示问题问题云图颜色分布不合理调整方法修改contourf的levels参数尝试不同的色彩映射(colormap)检查数据范围是否异常# 查看数据统计信息 print(f亮温最小值{np.nanmin(tbb)}) print(f亮温最大值{np.nanmax(tbb)}) print(f亮温平均值{np.nanmean(tbb)})7. 实际应用案例以2019年台风利奇马为例我们可以通过时间序列分析其发展过程8月7日台风在菲律宾以东形成亮温图像显示结构松散8月9日台风眼开始形成中心亮温升高8月10日达到巅峰强度眼区温度达290K以上8月11日登陆后结构逐渐破坏通过对比不同时间的云图可以清晰地看到台风结构的演变过程。在实际操作中我发现在处理连续时次数据时最好先统一所有图像的颜色标尺这样更有利于比较。