简介：RClimDex（1.0）官方用户手册的中文资源，由加拿大环境部气候研究分支的Xuebin Zhang和Feng Yang开发维护，面向气候科学家、研究人员及气候变化监测工作者。手册系统说明如何基于R环境安装、加载RClimDex，并完成数据质量控制与极端气候指数计算；涵盖高温日数、低温日数、极端降雨量等27个核心指数，以及自定义阈值与基准期的方法。资源包仅包含1个doc文档，大小约600KB，内容结构完整。手册不仅列出已知bug与报告流程，附录还包含气候指数清单、输入数据格式、指数定义、阈值计算详解及Windows环境常见问题解答，便于用户逐项对照操作。目前已有3836人学习，适合气候领域学生、科研人员及业务工作者阅读；也可作为相关课程的辅助参考资料，帮助快速上手极端气候指数分析。

1. 装 RClimDex 之前：先搞懂它为什么必须依赖 R 和 TclTk

接手一套逐日气象数据，要算 ETCCDMI 专家组推荐的 27 个核心极端气候指数，最省事的方案就是 RClimDex。这个包 2004 年就发布了 1.0 版，开发背景很明确:早期有个 Excel 版 ClimDex，但 2003 年发现基于百分位的温度指数在 Excel 里算会产生序列不匀性，修复需要做 bootstrap 重采样，Excel 环境几乎实现不了，开发团队才把整套逻辑移植到 R 上。R 免费、跨平台、统计和绘图能力都够硬，于是就有了你手上这份 RClimDex 1.0 手册对应的工具。装这套东西的坑不在软件本身，而在前置环境:R 的版本选不对、TclTk 缺失、数据格式写错，都会让它在 QC 阶段直接罢工。这篇笔记按我实际装过的顺序来写，适合要跑 27 个核心指数、又被老 R 脚本兼容性问题卡住的人。

2. 装 R 环境：版本选择与 TclTk 依赖，Windows 和 Linux 分开走

RClimDex 1.0 手册里写得很直白：它需要 R 的 base 包和图形界面依赖 TclTk。先说结论，再讲版本，最后落到安装命令上。

2.1 版本怎么选：1.9.0 是 TclTk 的分水岭

手册里写的是“developed under R 1.84, should run with R 1.84 or a later version”，同时强调“TclTk is included in the default installation of R 1.9.0 or later versions”。也就是说，R 1.9.0 之前 TclTk 要手动单独装，1.9.0 之后默认就带上。这个细节放在 2004 年是常识，放到今天很多人会忽略，因为现在装 R 4.x 时大家默认勾选全套组件，TclTk 早就不是什么稀罕物。

但真正的坑在这里:RClimDex 1.0 是 2004 年的老脚本，我在 R 4.x 上新装过一次，source 进去之后 TclTk 相关的变量初始化偶尔会报错，表现为控制台刷一段 tcl 错误后主菜单出不来。这不是 R 的 bug，是老脚本对太新版本的接口变化没跟上。我一般的处理是:先拿手头最新的 R 试跑，不行就退到 R 3.x 系列，专给 RClimDex 建一个独立环境。手册要求的 1.84 实在太老，不推荐刻意装那么旧的版本，那会儿的 R 在 Windows 上的安装程序和现在的 UCRT 运行时完全不是一回事。

2.2 Windows 下安装 R：注意安装器里的 Tcl/Tk 组件

Windows 上装 R 没什么技术含量，去 CRAN 任一镜像下载 Windows 安装程序，双击运行。唯一要注意的是安装界面里有一个组件勾选列表，Tcl/Tk 这一项必须确保被选中。R 1.9.0 之后的默认安装确实包含 TclTk，但某些精简安装流程或自定义安装路径的人会顺手把组件剪掉，等 RClimDex 加载时才想起来，又得回头补。

装完之后桌面上会有 R 的快捷方式，双击进入 RGui。这里有个老生常谈的问题:手册附录 E 提到，部分 Windows 机器需要手动设置 HOME 环境变量，否则 R 可能起不来或行为异常。现代 R 一般会自动识别用户目录，我在 Windows 10/11 上没再手动设过 HOME，但如果你的机器是精简版系统、用户目录异常，就在系统环境变量里补一个 HOME，指向你放数据的目录既简单又稳。

2.3 Linux 下安装 R：二进制优先，源码编译留作备选

Linux 上装 R，我从不建议源码编译，除非发行版仓库里没有预编译包。Debian/Ubuntu 系直接一条命令：

sudo apt update sudo apt install r-base r-base-dev tcl8.6 tcl8.6-dev tk8.6 tk8.6-dev

这条命令把 R 本体和 TclTk 的系统库一起装掉。r-base-dev 提供编译扩展组件时要用的头文件，虽然跑 RClimDex 用不上编译，但后面你自己写包、装其他依赖时一定会用到。tcl8.6 和 tk8.6 是 R 的 tcltk 包在 Linux 下的后端，缺了它们，R 里 library(tcltk) 会直接报找不到共享库。

CentOS/RHEL 系用 dnf install R R-devel tcl tcl-devel tk tk-devel，包名大同小异。装完后在终端敲 R 进入控制台，看到版本信息就说明 R 本体没问题。

2.4 装完先验证 TclTk：这一步不做，后面加载必翻车

不管 Windows 还是 Linux，装完 R 的第一件事不是急着下 RClimDex，而是验证 tcltk 是否真的可用。在 R 控制台里执行：

library(tcltk) capabilities("tcltk")

第一条把 tcltk 包载入，第二条返回逻辑值 TRUE 才说明 TclTk 环境完整。如果第二条返回 FALSE，先别怀疑 RClimDex，回头检查 Tcl/Tk 组件和系统库装没装全。这一步是我每次搭老 R 环境必做的体检项，能省掉后面一大半玄学问题。

3. 加载 RClimDex：source 命令的路径坑和控制台报错怎么读

R 环境就绪之后，目标就是让 RClimDex 的主菜单弹出来。RClimDex 不是一个传统意义上的 R 包，没有 install.packages 流程，它的全部代码在一个叫 rclimdex.r 的 R 脚本里，用 source 命令加载。

3.1 先把 rclimdex.r 放到工作目录：路径里别带中文和空格

从手册附带材料或你下载的资源里找到 rclimdex.r，我习惯把它单独放到一个工作目录，比如 D:\climdata 或 ~/climdata，和后面的站点数据分开。加载前先设置工作目录，Windows 下路径分隔符最容易写错：

setwd("D:/climdata") source("rclimdex.r")

注意 setwd 里的斜杠。R 在 Windows 下也接受正斜杠，反斜杠是转义符，写成 "D:\climdata" 会被 R 当成 \c 转义处理，直接报错。这是我最常见到的初级翻车点。另外路径里不要带中文和空格，老脚本对编码的处理很脆弱，路径一复杂，source 阶段就黑匣子给你看。

3.2 source 加载与 RGui 菜单加载：两种方式的结果一样

手册给了两种加载方式。第一种就是在 R 控制台直接敲 source("rclimdex.r")，脚本会在当前工作目录下找文件。第二种是针对 Windows 的 RGui 用户:菜单栏点 File，然后选 Source R code，在弹出的窗口里找到 rclimdex.r 文件双击。这个菜单操作本质上也是在执行 source，只是不用手动敲路径。

加载成功的标志是 R 控制台出现 RClimDex 主菜单，菜单里只有三个选项:Load Data and Run QC、Indices Calculation、Exit。看到这个菜单，说明脚本解析完成，TclTk 图形界面也初始化成功了。主菜单没出现之前，后面所有操作都无从谈起。

3.3 加载失败的三种常见表现：别急着改数据，先看环境

我遇到过的加载失败，九成是下面三种情况。第一种是 R 报 “cannot open file 'rclimdex.r': No such file or directory”，这是工作目录不对，文件不在这。第二种是报 tcl 相关的错误，比如 “TclError: can't find package”，这是 TclTk 环境没装齐，回到 2.4 的验证步骤查。第三种是脚本加载过程报一串类似 “object not found” 的错误，这通常发生在 R 4.x 上跑老脚本，某些过时的接口调用失败。第三种最麻烦，我的处理是换一个 R 3.x 版本单独建环境，不跟主 R 版本抢位置。

手册里还提过一种在线加载方式:直接 source 一个网址，从官方站点拉最新版 rclimdex.r。这个功能听着方便，实操我不推荐，一是那个 2004 年的地址早就不是稳定资源，二是网络源一旦挂了整个环境就瘫了。把脚本下载到本地再加载，才是能复现的部署方式。

4. 准备输入数据：stationname.txt 六列格式与 -99.9 缺失值约定

加载成功只是开始，真正的门槛在数据格式。RClimDex 对输入文件的格式有严格约定，手册附录 B 专门讲这个，格式错了 QC 阶段会弹一堆错误。这份数据格式是所有环节里最不能自由发挥的地方。

4.1 六列 ASCII 结构：年、月、日、降水、Tmax、Tmin

输入文件是纯文本 ASCII，文件名按站点名命名，手册里叫 stationname.txt，示例用的是一个叫 21946.txt 的站点文件。每一行是一条逐日记录，六列，按顺序是年份、月份、日、日降水量、日最高温、日最低温。手册在讲某个 bug 时给出了一个标准行格式示例:

1960 12 31 -99.9 -99.9 -99.9

这一行就是基准期 bug 的解决方案:在数据最前面补一条基准期开始前一天的缺失记录。它也顺带把六列结构交代清楚了——第四列降水，第五列最高温，第六列最低温。实际数据里各列用空格或制表符分隔，R 读取时默认按空白符切分，关键是列顺序不能错。

我做数据预处理时的一般习惯是:年份、月份、日期用整数，数值列保留一位小数，缺失统一填 -99.9。单位方面，手册没有强制规定，但气候学通用做法是降水用毫米、温度用摄氏度，你自己心里有数就行，指数计算不关心单位，只关心数值。

4.2 缺失值 -99.9 与不合理值的分水岭：别混淆两者

-99.9 是 RClimDex 约定的缺失值编码，QC 阶段会被统一替换成 R 内部的 NA。除此之外，还有两类不合理值会被 QC 自动替换成 NA:日降水量小于 0 的负值，以及日最高温小于日最低温的记录。这两类在真实数据里并不罕见，尤其是降水负值，很多老观测站的整理数据里会把微量降水记成负值，格式上看着像有效数，逻辑上却站不住脚。

这里有个很容易被忽略的点:被标记的异常温度和负降水会进 QC 报告，但 QC 输出的 indcal.csv 里只把缺失值和不合理值替换成 NA，异常值只是被标记、不会被自动修改。人工复核完异常值后如果决定修正原始数据，必须重新跑一遍 QC，指数才会基于修正后的数据重新计算。

4.3 用脚本把日常表格转成 RClimDex 格式：多留一个心眼

很多人的原始数据是 Excel 表格，直接另存为 txt 往往带着制表符和多余表头，RClimDex 读不进去。我一般用一段简短脚本转格式，先读原始 CSV，再按六列写出。下面这段不要照抄，按你的实际列名改：

import csv input_file = "raw_station.csv" output_file = "station.txt" with open(input_file, newline="", encoding="utf-8") as fin, \ open(output_file, "w", newline="", encoding="utf-8") as fout: reader = csv.reader(fin) next(reader) # 跳过表头 rows = [] for row in reader: y, m, d = row[0], int(row[1]), int(row[2]) prcp, tmax, tmin = row[3], row[4], row[5] # 缺失值统一改成 -99.9 prcp = prcp if prcp not in ("", "NA", "NaN") else "-99.9" tmax = tmax if tmax not in ("", "NA", "NaN") else "-99.9" tmin = tmin if tmin not in ("", "NA", "NaN") else "-99.9" rows.append(f"{y} {m:02d} {d:02d} {prcp} {tmax} {tmin}") with open(output_file, "w", encoding="utf-8") as fout: fout.write("\n".join(rows) + "\n")

这段代码做三件事:字段映射、缺失值统一替换、日期补零对齐。月份和日期补成两位数是为了让输出文件规格统一。写成空格分隔、不带任何表头，输出文件就能直接进 RClimDex。转换完我永远会人工抽查几行，看看年份是否连续、日期有没有跳变，这一步能拦下大部分后续 QC 报错。

5. 运行数据质量控制与常见问题排查：从 -99.9 到 NA 再到 log 目录

数据文件就绪后，在主菜单选 Load Data and Run QC，弹窗里选择站点 txt 文件。这一步完成的是三件事:把 -99.9 全替换成 NA、把不合理值替换成 NA、标记可能的温度异常值。QC 是计算指数前的强制前置步骤，原始数据文件不会被改动。

5.1 QC 内部到底做了什么：缺失值替换、不合理值清理、异常值标记

QC 的逻辑分三层。第一层，所有 -99.9 转成 R 能识别的 NA，这个没有商量余地。第二层，降水负值、Tmax 小于 Tmin 的记录全部替换成 NA，这是程序自动完成的。第三层，温度异常值识别，程序按某一天的多年平均值为中心，加上一个标准差区间，落在区间外部的点被标记为 potential outlier。这个区间的宽度由用户给一个参数 n 决定，默认 3，手册明确建议可以先设成 4，因为 n=3 时误标率很高，会把大量正常天气事件当成异常。

QC 结束后 log 子目录会生成四个 CSV 和四个 PDF。四个 CSV 分别是:温度不合理值记录、降水不合理值记录、温度潜在异常值完整列表、QC 后可用于指数计算的数据文件。四个 PDF 是逐日降水、逐日最高温、逐日最低温、逐日温度较差的时间序列图，缺失值在图上用红点标出。

5.2 异常值参数 n 设多少：默认 3 偏严，我一般从 4 开始

参数窗口里的 Criteria 默认值是 3，这个值代表 mean ± n*std 的区间倍数。n 越小、区间越窄、被标成异常的点越多。夏季 35℃ 的最高温在一个平均 28℃、标准差 3℃ 的站点上，n=3 时区间上限才 37℃，35℃ 虽然不会被标，但 38℃ 的极端热浪就会触发标记。这类真极端事件本来就是你做极端气候指数要研究的目标，所以标记不等于删数据，它只是提醒你人工复核这一刻是不是观测错误。

我一般先把 n 设为 4，跑完看 tepstdQC.csv 里到底标了多少点，如果单个站点每年标出几十上百个点，明显异常值占比过高，再逐步降到 3 或 3.5。关键是记住:这一步只标记，不改数，最终修正权在你手里。

5.3 log 目录里的四个文件和四个图怎么配合使用

log 目录是 QC 阶段唯一的产物输出地。四个文件的命名是站点名后接功能后缀，以 21946 为例，就是 21946tempQC.csv、21946prcpQC.csv、21946tepstdQC.csv、21946indcal.csv。前两个是水和温度的不合理值记录，用来追查原始数据问题;第三个是潜在异常值清单，按日期列出来供人工判断;第四个 indcal.csv 是 QC 后数据，指数计算的唯一数据源。

复核异常值时我一般只看 tepstdQC.csv，把每个标记日期对应回原始数据，看当天前后几天的数值是不是连续。如果前一天 30℃ 后一天 31℃，中间突然冒出一个 12℃，那基本可以断定是观测或录入错误，回原始文件修掉后重跑 QC。如果当天确实是一场极端天气过程，数值是合理记录，那就不动它，保留真实值参与指数计算。

5.4 五个常见问题排查：现象、原因、解决

第一个坑是 QC 慢得像死机。现象:点 OK 后 R 窗口长时间无响应，进度条都没有。原因:QC 本身涉及逐日气候态计算和异常值扫描，老脚本在慢 PC 上跑几百万次循环是正常的。解决:别急着强制关闭窗口，手册明确说这一步在慢机器上可能要几分钟。我会先等足 3 到 5 分钟，实在不动再考虑数据量是否过大。

第二个坑是加载数据时报格式错误。现象:R 控制台弹出一行或几行错误信息，QC 窗口没出来。原因:输入文件列顺序不对、有多余表头、分隔符不统一，字符编码异常也可能触发。解决:拿手册示例 21946.txt 的格式逐列对照，重点看第一行和第二行是否同构。我最快的一次排错，就只是发现文件里混了几个全角空格。

第三个坑是首年与基准期相同导致程序停跑。现象:指数计算阶段程序直接中断，没有任何明确报错。原因:这是手册里明确记载的已知 bug，某些计算需要基准期边界以外的数据。解决:在输入文件最前面补一条基准期开始前一天的缺失记录。基准期 1961-1990、数据也从 1961 年开始，就在第一行加 1960 12 31 -99.9 -99.9 -99.9，把时间轴往前垫一天，程序就有边界外的数据可用了。

第四个坑是改完原始数据后指数没变化。现象:修正了原文件的错误值，重跑指数，结果和之前完全一样。原因:指数计算用的不是原文件，而是 QC 阶段生成的 indcal.csv，原文件从来没被改过。解决:改完原数据必须重新执行 Load Data and Run QC，让 indcal.csv 重新生成，指数才会反映修正后的数据。

第五个坑是 Windows 下 R 起不来。现象:双击 R 图标没反应，或启动后立即退出。原因:老版本 Windows 下 HOME 环境变量未设置，R 找不到用户目录配置。解决:按手册附录 E 的做法，在系统环境变量里手动加 HOME。现代 R 一般自动处理，但这个老问题在精简系统和某些远程桌面环境里依然会出现。

6. 指数计算与验证跑通：从基准期参数到 plots 里的 JPEG 趋势图

QC 通过后，主菜单选 Indices Calculation，先弹参数设置窗口，再选指数，最后看输出。这一步是整个 RClimDex 的收口，我也是用它来验证环境是否真的装通。

6.1 参数窗口：基准期、纬度、降水阈值与四个温度阈值

参数窗口第一个要填的是基准期的起始年份和结束年份，这是百分位温度指数的计算基础。第二个是站点纬度，南半球填负数，程序用它判断生长季跨年计算的半球差异。第三个是自定义降水阈值 P，单位毫米，用来算超过 P 的降水日数(对应 Rnn 指标)。剩下四个是自定义温度阈值，分别对应:日最高温超过某值的天数、日最高温低于某值的天数、日最低温超过某值的天数、日最低温低于某值的天数。这四个阈值产生的指标在命名上会带上具体数值，比如阈值为 30℃ 时，日最高温超过 30℃ 的天数指标就是 SU30。

这里要注意，手动填写的只是一部分参数，真正输出的 27 个核心指数是程序预设的完整集合。手册附录 A 列了核心指标清单，以温度类为例，有霜冻日数 FD0(日最低温低于 0℃ 的天数)、夏日 SU25(日最高温超过 25℃ 的天数)、冰日 ID0、热带夜 TR20、生长季长度 GSL，以及月尺度极值 TXx、TNx、TXn、TNn 和百分位类 TN10p、TX10p、TN90p、TX90p 等。降水类则围绕逐日降水极值、连续干湿日和强降水贡献展开。计算全部 27 个指数是默认行为，你可以在接下来的窗口里按需反选。

6.2 指数选择与输出目录：indices 的 CSV 和 plots 的 JPEG

参数设置完成后进入指数选择窗口，默认全部勾选。如果你的研究只需要温度相关指标，就手动把降水类指数全部反选，能明显缩短计算时间。点 OK 后计算开始，同样可能出现几十秒到几分钟的等待。计算完成弹出提示窗口，输出按两个子目录组织:indices 目录存放每个指数一个 CSV 文件，文件名规则是站点名_指数名.csv，列之间用逗号分隔;plots 目录存放每个指数一张 JPEG 趋势图。

plots 里的图不是简单的时间序列，图上叠加了两条拟合线:实线是线性最小二乘回归趋势，虚线是局部加权回归趋势，线性趋势的拟合统计量也会直接显示在图上。这是查看某个极端指数是否在显著变化的最直观入口。

6.3 装完怎么验证跑通：用一份示例数据走完加载、QC、出图

装完这一整套环境，我会强制自己用示例数据走一遍完整流程，而不是直接上真实数据。首选就是手册自带的 21946.txt 这类站点文件，它格式标准、数据量适中，适合做验证样本。完整验证按四步走:

第一步加载。source("rclimdex.r") 后确认主菜单出现，三个选项齐全，这一步验证的是 R 与 TclTk 环境。第二步 QC。Load Data and Run QC，选择示例站点文件,把 n 设为 4，点 OK，等待 QC 完成弹窗。检查 log 目录是否生成了四个 CSV 和四个 PDF。第三步参数与计算。Indices Calculation 里填好基准期，注意避开数据首年与基准期首年相同的情况，纬度按站点实际填写，默认全选指数，点 OK 等待计算完成。第四步产物核对，用 shell 命令确认输出文件存在:

ls log/ # 应有 4 个 csv 和 4 个 pdf ls indices/ # 应有与所选指数数量一致的 csv ls plots/ # 应有与所选指数数量一致的 jpg head -5 indices/21946_TXx.csv # 查看指数序列开头几行

一个站点正常跑完，indices 和 plots 目录里的文件数应该等于你选择的指数个数。打开一张 JPEG 图，确认上面有时序曲线、红色缺失值标记点、实线和虚线两条拟合线，以及趋势统计量文本。看到这些，这套 RClimDex 环境就算真正跑通了。

从那以后我每次把 R 环境装到一台新机器上，都会先拿示例数据强制走一遍 QC 加指数计算，确认四类输出文件齐全才敢放真实数据进来。环境问题永远要比数据问题先解决，希望帮到你。

本文还有配套的精品资源，点击获取