保姆级教程：手把手教你准备DSSAT、WOFOST、APSIM三大作物模型的气象数据文件（附Python脚本）

三大作物模型气象数据实战：从原始数据到DSSAT/WOFOST/APSIM标准格式

第一次接触作物模型的研究者，往往会在气象数据准备阶段卡壳——手头的CSV数据怎么变成模型要求的特殊格式？单位不统一怎么办？关键参数如何计算？本文将用真实气象站数据演示完整处理流程，提供可直接复用的Python代码模板。

1. 气象数据预处理基础

无论使用哪种作物模型，原始气象数据都需要经过清洗和标准化处理。假设我们从气象站获取了包含以下字段的CSV文件：

• 日期（YYYY-MM-DD）
• 日最高温（℃）
• 日最低温（℃）
• 太阳辐射（MJ/m²/d）
• 降水量（mm）
• 风速（m/s）
• 相对湿度（%）

常见数据问题处理：

import pandas as pd # 读取原始数据示例 raw_data = pd.read_csv('weather_station.csv', parse_dates=['日期']) raw_data = raw_data.dropna() # 删除缺失值 raw_data = raw_data[raw_data['太阳辐射'] >= 0] # 剔除负值

单位换算对照表：

特别注意：WOFOST的辐射量需要乘以1000转换单位，这是新手最容易出错的地方

2. DSSAT气象文件(.WTH)生成全流程

2.1 文件命名规范

DSSAT要求严格的8字符文件名格式：

• 前4位：气象站代码（大写字母）
• 5-6位：年份后两位
• 7-8位：站点编号（通常为01）

示例代码生成文件名：

station_code = "UFGA" year_suffix = "22" site_number = "01" filename = f"{station_code}{year_suffix}{site_number}.WTH"

2.2 关键参数计算

AMP（年温度振幅）计算：

# 计算月平均温度 monthly_avg = raw_data.resample('M', on='日期').mean() # 计算各月(Tmax+Tmin)/2 monthly_mean = (monthly_avg['日最高温'] + monthly_avg['日最低温']) / 2 # AMP = 最热月平均 - 最冷月平均 AMP = monthly_mean.max() - monthly_mean.min()

2.3 完整文件生成

使用pandas生成标准格式：

header = f""" *WEATHER DATA : {station_code} @ INSI LAT LONG ELEV TAV AMP REFHT WNDHT {station_code} 34.58 -103.20 1348 {TAV:.1f} {AMP:.1f} -99.0 -99.0 @DATE SRAD TMAX TMIN RAIN """ with open(filename, 'w') as f: f.write(header) for _, row in raw_data.iterrows(): day_num = row['日期'].dayofyear f.write(f"{row['日期'].year%1000}{day_num:03d} {row['太阳辐射']:.1f} {row['日最高温']:.1f} {row['日最低温']:.1f} {row['降水量']:.1f}\n")

3. WOFOST气象文件(.XXX)制作要点

3.1 文件结构特点

• 文件扩展名为3位年份数字（如2022年→022）
• 每文件只包含单年数据
• 辐射量单位必须为KJ/m²/d

单位转换示例：

# MJ/m²/d → KJ/m²/d raw_data['太阳辐射'] = raw_data['太阳辐射'] * 1000

3.2 文件内容生成

典型WOFOST气象文件结构：

站点名称 纬度 经度 海拔 年 日序 辐射 最高温 最低温 水汽压 风速 降水量 2022 1 12500 12.5 3.2 1.2 3.0 0.0 2022 2 9800 10.8 2.5 1.1 2.8 2.4

Python实现代码：

wofrost_header = f"{station_name} {latitude} {longitude} {elevation}" with open(f"{station_code}022.XXX", 'w') as f: f.write(wofrost_header + "\n") for _, row in raw_data.iterrows(): day_num = row['日期'].dayofyear line = f"{row['日期'].year} {day_num} {int(row['太阳辐射'])} {row['日最高温']:.1f} {row['日最低温']:.1f} {row['水汽压']:.1f} {row['风速']:.1f} {row['降水量']:.1f}" f.write(line + "\n")

4. APSIM气象文件(.met)高效生成方案

4.1 Excel模板法

1. 下载官方模板（通常位于APSIM安装目录/Examples/WeatherFiles/MetData.xlsx）
2. 按模板格式填写数据
3. 另存为"文本（制表符分隔）(*.prn)"
4. 手动修改文件扩展名为.met

4.2 Python直接生成

apsim_header = f""" [weather.met.weather] stationname = {station_name} latitude = {latitude} longitude = {longitude} tav = {TAV} amp = {AMP} year day radn maxt mint rain """ with open(f"{station_name}.met", 'w') as f: f.write(apsim_header) for _, row in raw_data.iterrows(): day_num = row['日期'].dayofyear f.write(f"{row['日期'].year} {day_num} {row['太阳辐射']:.1f} {row['日最高温']:.1f} {row['日最低温']:.1f} {row['降水量']:.1f}\n")

5. 常见问题解决方案

问题1：DSSAT运行提示找不到气象文件

• 确保文件放在Weather子目录
• 检查文件名是否为8字符（不含扩展名）
• 确认文件扩展名大写.WTH

问题2：WOFOST模拟产量异常低

• 首先检查辐射单位是否为KJ/m²/d
• 验证温度数据范围是否合理
• 检查降水数据是否有负值

问题3：APSIM无法读取.met文件

• 确认文件是纯文本格式
• 检查表头是否符合标准格式
• 尝试用文本编辑器查看特殊字符

三个模型气象参数对比表：

实际项目中，我通常会先准备DSSAT格式文件，然后通过单位转换生成WOFOST版本，最后用Python脚本转为APSIM格式。这种方法可以确保三个模型使用完全一致的基础数据，避免因数据预处理差异导致的模拟结果偏差。