保姆级教程：用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 basemap

1.2 获取FY-2G卫星数据

中国风云二号(FY-2G)静止气象卫星的全圆盘数据可以从国家卫星气象中心获取。对于"利奇马"台风，我们需要以下两种文件：

1. HDF数据文件：包含红外通道的原始观测数据
2. 经纬度对照表：用于将卫星数据映射到地理坐标

典型的数据文件命名格式为：FY2G_FDI_ALL_NOM_YYYYMMDD_HHMM.hdf，其中包含了特定时间的观测数据。

2. 理解卫星数据与亮温(TBB)

气象卫星的红外通道测量的是地气系统的辐射亮度，通过定标可以转换为亮温(Brightness Temperature)。亮温是分析台风强度的重要参数。

2.1 HDF文件结构解析

FY-2G的HDF文件包含多个数据集，我们需要重点关注：

2.2 亮温计算原理

亮温的计算过程可以概括为：

1. 读取原始观测值(DN值)
2. 使用定标表将DN值转换为辐射亮度
3. 通过普朗克公式反算得到亮温

以下是关键的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, dtype=np.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 tbb

3. 地理坐标转换与地图投影

卫星数据的坐标系统与常规地理坐标不同，需要进行转换才能正确显示在地图上。

3.1 读取经纬度对照表

FY-2G卫星提供了专门的经纬度对照文件(.dat)，用于将图像像素映射到地理坐标：

def read_geolocation(dat_file): with open(dat_file, 'rb') as f: data = np.fromfile(f, dtype='float32') # 风云二号卫星数据特有的存储格式 data = data.reshape([2288, 2288, 2], order='F') longitude = data[:,:,0] + 104.5 # 卫星星下点经度修正 latitude = data[:,:,1] return longitude, latitude

3.2 地图投影设置

对于台风分析，我们通常使用圆柱投影(Plate Carrée)来保持区域形状：

m = Basemap(projection='cyl', llcrnrlat=15, urcrnrlat=45, # 纬度范围 llcrnrlon=105, urcrnrlon=145, # 经度范围 resolution='i')

4. 完整云图绘制流程

现在我们将所有步骤整合起来，绘制台风"利奇马"的卫星云图。

4.1 数据处理流程

1. 读取HDF文件获取亮温数据
2. 读取.dat文件获取地理坐标
3. 设置地图投影参数
4. 绘制海岸线、国界等地理要素
5. 将亮温数据投影到地图上
6. 添加色标、标题等辅助元素

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(projection='cyl', llcrnrlat=15, urcrnrlat=45, llcrnrlon=105, urcrnrlon=145, resolution='i') # 绘制地理要素 m.drawcoastlines(linewidth=0.5) m.drawcountries(linewidth=0.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, levels=np.linspace(180, 300, 121), cmap='jet') # 添加色标 cbar = m.colorbar(cs, location='right', pad="5%") cbar.set_label('Brightness Temperature (K)') # 添加标题 plt.title('Typhoon Lekima - FY2G IR1 Brightness Temperature\n2019-08-10 12:00 UTC', fontsize=14, pad=20) # 保存图像 plt.savefig(output_png, dpi=300, bbox_inches='tight') plt.close()

5. 台风分析进阶技巧

生成卫星云图只是第一步，更重要的是从中提取有用的气象信息。

5.1 识别台风眼

在亮温图像中，台风眼通常表现为：

• 小而圆的温暖区域(亮温较高)
• 周围被极低温的云墙包围
• 直径通常在30-60公里之间

我们可以通过寻找局部温度最大值来定位台风眼：

def find_eye_center(tbb, threshold=280): # 寻找高于阈值的区域 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_y

5.2 估算台风强度

虽然精确的强度估算需要更复杂的方法，但亮温可以提供初步判断：

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, levels=np.linspace(180, 300, 121), cmap='jet') m.drawcoastlines() anim = FuncAnimation(fig, update, frames=len(hdf_files), interval=200) anim.save('typhoon_animation.mp4', writer='ffmpeg')

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 地图投影问题

问题：数据在地图上的位置不正确
检查步骤：

1. 确认卫星星下点经度修正值是否正确
2. 检查地图范围是否包含数据区域
3. 验证经纬度数据的符号（东经/西经，北纬/南纬）

6.3 图像显示问题

问题：云图颜色分布不合理
调整方法：

• 修改contourf的levels参数
• 尝试不同的色彩映射(colormap)
• 检查数据范围是否异常

# 查看数据统计信息 print(f"亮温最小值：{np.nanmin(tbb)}") print(f"亮温最大值：{np.nanmax(tbb)}") print(f"亮温平均值：{np.nanmean(tbb)}")

7. 实际应用案例

以2019年台风"利奇马"为例，我们可以通过时间序列分析其发展过程：

1. 8月7日：台风在菲律宾以东形成，亮温图像显示结构松散
2. 8月9日：台风眼开始形成，中心亮温升高
3. 8月10日：达到巅峰强度，眼区温度达290K以上
4. 8月11日：登陆后结构逐渐破坏

通过对比不同时间的云图，可以清晰地看到台风结构的演变过程。在实际操作中，我发现在处理连续时次数据时，最好先统一所有图像的颜色标尺，这样更有利于比较。