★ FEATURED ARTICLE 简介RClimDex极端气候指数计算软件用户手册面向气候科学研究者、气象数据分析人员及高校相关专业师生用于在 R 环境中完成极端气候指数的计算与变化检测。手册系统说明软件安装与运行、数据质量控制、指数计算流程及已知问题反馈并附有可用气候指数清单、输入数据格式、各指数定义与算法说明以及 Windows 下 R 常见问题解答便于读者按附录逐项对照上手。压缩包内共 1 个 doc 文档约 600KB即该手册的完整英文原版适合边操作边查阅。目前已有 3839 人学习下载是入门 RClimDex 与规范处理气象序列数据的实用参考。1. 极端气候指数计算软件不是点一下出图RClimDex 1.0 先解决数据入口问题极端气候指数计算软件 RClimDex 1.0 的入门门槛不在公式而在数据入口一份逐日气象文本如果列顺序、缺测编码或基准期没对齐后面的 27 个指数全都会歪。它适合手里已经有逐日降水、最高温、最低温序列需要按统一口径算 FD0、SU25、TX90p、R95p 这类指标的人不适合拿月值或年值来硬套。我见过最容易翻车的一步是把它当成独立安装的图形软件双击就想出结果。实际拿到的是一套 R 脚本加图形菜单先有 R 环境再 source 脚本再走 QC 和指数计算两条流程顺序错了后面全是无用功。这一章先把边界说清它做质量控制、核心指数计算和出图不做数据均一化也不替你判断异常值该不该改。适合做单站复现、课程实验和批量站点前的模板验证。2. 把 R 环境和 rclimdex.r 跑通从 source 命令到主菜单出现的四步RClimDex 1.0 的运行方式很朴素R 控制台加载rclimdex.r脚本把图形菜单挂出来。手册对应 R 1.84 及以后版本较新的 R 默认带 TclTk早期版本可能要单独装。很多人卡在“脚本加载了但菜单没出来”多半不是脚本坏而是工作目录、路径空格或图形组件没就位。下面按我平时复现的顺序走先保证 R 能启动再保证脚本能被找到最后才谈数据。2.1 R 版本与 TclTk 依赖为什么老手册总强调图形界面RClimDex 的交互入口是 TclTk 图形窗口不是纯命令行参数。R 环境装好之后先确认两件事R 能正常启动图形界面组件能弹窗。较新的 R 安装包通常默认包含 TclTk如果用的是很早的版本或者安装时裁掉了图形组件执行加载后可能只在控制台报错菜单不出现。我的做法是先开一个干净的 R 会话不加载任何额外包执行一条最简单的弹窗测试。能弹窗说明图形组件可用不能弹窗先修 R 安装不要急着改rclimdex.r。这一步看似笨但能省掉后面“到底是脚本问题还是环境问题”的拉扯。# 先确认当前 R 会话能正常启动 R.version.string # 测试 TclTk 图形组件是否能弹窗 # 能弹出窗口说明图形环境可用报错则先修 R 安装或图形组件 library(tcltk) tk_messageBox(type ok, message TclTk OK)这段代码只做环境体检。R.version.string用来看当前 R 版本后面排查“老脚本在新 R 上行为不一致”时有对照价值。library(tcltk)是加载图形组件tk_messageBox是弹窗测试如果这里就失败说明问题在 R 安装层不在 RClimDex 脚本层。参数上没有需要调的重点看报错信息里是找不到包还是图形显示不可用。2.2 本地 source 与菜单加载两条路径与工作目录设置加载脚本有两条常见路径控制台执行source(rclimdex.r)或者从图形菜单里选“Source R code”再挑文件。两条路径本质一样差别在路径解析。控制台路径相对当前工作目录菜单路径一般按你选中的文件走。我一般先把工作目录切到脚本所在目录再用相对路径加载避免脚本内部再去找别的文件时跑偏。工作目录里不要有中文、空格和特殊符号。不是脚本一定不支持而是老脚本在不同平台上的路径处理很容易出玄学问题。一个稳妥习惯是建一个纯英文目录把rclimdex.r和待处理数据分开放脚本目录放脚本数据目录放stationA.txt这类输入文件日志和结果让脚本自己建子目录。# 把工作目录切到脚本所在目录路径按自己的目录替换 setwd(D:/climate_work/rclimdex_demo) # 确认当前目录和文件是否存在 getwd() file.exists(rclimdex.r) # 加载脚本成功后应出现 RClimDex 主菜单 source(rclimdex.r)setwd()的参数是脚本所在目录建议用正斜杠或双反斜杠。getwd()用来确认没有切错盘符。file.exists(rclimdex.r)返回TRUE再执行source比直接 source 后看一屏报错更清楚。source()的参数是脚本文件名如果脚本不在当前目录就写完整路径但完整路径里尽量别带空格。加载成功后主菜单会出现没出现就回到控制台看第一条报错不要从最后一行倒着猜。2.3 加载失败的排查顺序路径、编码、包冲突加载失败我按固定顺序查先看工作目录再看文件名大小写再看脚本编码最后看是否有同名函数或旧会话残留。R 对文件名大小写敏感Rclimdex.r和rclimdex.r在两个平台上可能一个找一个不找。脚本编码如果是本地编码保存的在某些 R 会话里读进去会出乱码菜单文字可能变成问号甚至解析中断。旧会话残留也常见。上一次加载到一半报错变量还留在环境里第二次再 source 就可能出现函数覆盖或重复定义。我的习惯是每次重来先清空对象再重新 source。不要一边改脚本一边在同一个会话里反复 source出了问题很难判断是新代码还是旧对象。# 清空当前会话对象避免旧脚本残留干扰 rm(list ls()) # 重新确认目录和脚本 getwd() file.exists(rclimdex.r) # 重新加载 source(rclimdex.r, encoding UTF-8)rm(list ls())清空当前环境对象不影响磁盘文件。encoding UTF-8是给脚本读取指定编码如果脚本本身不是 UTF-8改成实际编码或者用编辑器另存为 UTF-8 再试。这个参数不是万能药路径错误时它救不了所以前面两步必须先过。加载成功后不要急着点数据先看主菜单是否完整菜单项缺失往往还是图形组件或脚本版本不匹配。注意RClimDex 1.0 是脚本加载式工具不是双击安装的独立程序先跑通 R 和 TclTk再谈数据格式。3. 数据质量控制-99.9、负降水与 TXTN 在计算前怎么被拦下RClimDex 的 QC 不是可选项手册明确要求先做质量控制再算指数。它处理三类问题缺测值-99.9转成 R 认识的NA不合理值转成NA包括负降水和日最高温小于日最低温温度异常值按“当日气候平均加减 n 倍标准差”标记。这里最容易误解的是“标记”和“修改”不是一回事异常值只被标出来不会被自动改掉QC 后的数据文件里仍然保留原值。下面把输入格式、n 值、日志文件和重跑规则拆开。3.1 输入格式与缺测编码stationname.txt 的列结构与日期边界输入文件名一般用stationname.txt内容按手册附录 B 的列顺序组织。常见逐日格式是年、月、日、降水、日最高温、日最低温文件里不放表头缺测统一写-99.9。不同版本的示例数据可能在列数上有差异但核心是让脚本按固定位置读列而不是按列名猜。你拿到的数据如果来自表格软件先另存为纯文本再检查有没有逗号、制表符混用。日期边界要特别注意。QC 和指数计算都依赖连续日序缺一天不会自动补月末日期写错会直接把后面列值错位。我的做法是先用 R 读一遍前几行和尾几行确认列数、缺测编码和日期范围再进图形菜单。不要小看这一步列错位导致的“异常值”弹窗十有八九不是气候异常而是数据入口歪了。# 读取输入文件前几行确认列数、缺测编码和日期范围 # 这里假设无表头列顺序为 YEAR MONTH DAY PRCP TMAX TMIN raw - read.table(stationA.txt, header FALSE, fill TRUE) head(raw, 5) tail(raw, 5) dim(raw)read.table()的header FALSE对应无表头输入fill TRUE只是防止列数不齐时直接中断方便先看问题。head()和tail()看首尾日期dim()看总行数和列数。参数上如果你的文件是逗号分隔要改用sep ,如果是制表符用sep \t。这一步不替代正式 QC只是提前发现列错位、空行和编码问题减少后面弹窗的误判。3.2 异常值窗口里的 n 值默认 3 为什么要先试 4温度异常值判定用“当日气候平均 ± n 倍标准差”n 默认是 3。手册里给了一个很实用的提醒3 可能标出大量值用户可以先用 4 起步。原因不复杂逐日温度序列在样本少、季节波动大或缺口多的时候3 倍标准差会显得很紧边缘年份和极端天气容易被大量标记。标出来的值不是让你全改而是让你抽查。我一般先填 4 跑一遍 QC看日志里异常值清单的长度和分布。如果清单集中在某些年份或某些月份先查数据源和录入问题如果清单很分散且数值本身合理再考虑是否收紧到 3。n 不是越小越严格越好它只是筛查尺度。真正要改数据还是回到原始文件改再重新跑 QC不要在 QC 输出文件上随手改完就算完。注意n 值影响的是“标记范围”不是自动清洗强度异常值是否修改由你根据原始记录判断。3.3 log 目录四类输出QC 后数据、异常清单与四张 PDFQC 跑完后log子目录里会留下几类文件温度不合理值清单、降水不合理值清单、温度异常值标记清单以及供指数计算用的 QC 后数据。手册里对应的是tempQC.csv、prcpQC.csv、tepstdQC.csv、indcal.csv这一类命名实际文件名会带站点前缀。前两个清单帮你定位负降水、TXTN 这类硬错误第三个清单列出被标记的异常值及日期最后一个才是后续指数计算的输入。同时还会出四张 PDF 时间序列图缺测值用红点标出分别看降水、最高温、最低温、日温差。PDF 图不是装饰它最适合看缺口和断裂某些年份整段空白或者某个季节突然密集红点往往比表格更直观。我的习惯是先翻四张图再打开异常值清单最后才决定要不要动原始数据。不要一上来就改indcal.csv因为原始输入文件没有被改动你改输出文件只会让下一次重跑对不上。输出类型主要看什么处理方式温度不合理值清单最高温小于最低温、缺测编码错位回原始文件改再重跑 QC降水不合理值清单负降水、非数值字符回原始文件改再重跑 QC温度异常值标记清单被 n 倍标准差标出的日期和值抽查原始记录决定是否修改QC 后数据文件缺失值和不合理值已转NA异常值仍保留作为指数计算输入不随手改四张 PDF 图缺测红点、序列断裂、季节缺口先看整体再看局部年份这张表里最容易被忽略的是最后一行“异常值仍保留”。很多人看到异常值清单里有记录就以为 QC 后数据已经被修正实际不是。QC 只把缺失值和不合理值转成NA被标记的异常值原样留在indcal.csv里。要不要改得回到原始文件改完再跑一次 QC新的结果才会进入指数计算。# 读取 QC 日志先看异常值标记清单的规模和日期分布 tepstd - read.csv(log/stationAtepstdQC.csv) head(tepstd) nrow(tepstd) # 按年份统计被标记的异常值数量快速看是否集中在某段 tepstd$year - substr(tepstd[[1]], 1, 4) table(tepstd$year)这段代码里read.csv()读的是log目录下的异常值清单文件名按你实际站点前缀替换。head()看列结构nrow()看标记总量最后用第一列截取年份做频数统计。参数上要注意第一列是不是日期列不同输出文件的列名可能不一样先head()看清再改。这个统计不替代人工判断只是让你快速知道异常值是全年分散还是集中在某几年后者更像数据源或录入问题。3.4 原始数据回改后的重跑规则QC 有一个很重要的边界指数计算用的是 QC 后数据原始输入文件不会被改动。你如果在原始文件里修正了某些值必须重新执行“Load Data and Run QC”让新的 QC 后数据覆盖旧的indcal.csv否则指数计算还在读旧结果。这个规则听起来简单实际很容易踩坑改完原始文件直接点“Indices Calculation”结果和没改一样。我的固定动作是改原始文件后先看文件修改时间再重跑 QC再看log里的清单有没有减少最后才进指数计算。批量站点时更要注意每个站点的log和indices输出都可能混在一起目录一定要按站点分开。不要在一个目录里反复覆盖不然出了问题连是哪一版数据都说不清。4. 27 个核心指数怎么算基准期、纬度符号与五个自定义阈值QC 过了之后才进入指数计算。RClimDex 1.0 能算手册附录 A 里的 27 个核心指数也可以只选需要的。计算前会弹参数窗口要求填基准期首末年份、站点纬度、日降水阈值 P以及四个自定义温度阈值。这里的参数不是随便填基准期决定百分位指数的 in-base 口径纬度符号决定南北半球的 GSL 计算P 和四个温度阈值直接决定 Rnn、SUmm、FDmm、TRmm、IDmm 的数值。下面按参数顺序拆。4.1 基准期与 in-base 百分位为什么不能拿整段序列直接算百分位类指数比如冷夜、暖夜、冷日、暖日依赖基准期内的日值分布。手册在附录 D 专门讲阈值和 in-base 期温度指数计算核心意思是百分位阈值要在基准期内算再应用到整段序列。你不能拿整段序列的分布去定阈值否则基准期外的极端年份会反过来影响阈值指数序列的可比性就没了。基准期首末年份要填完整不是只填一个起点。常见的坑是数据从基准期首年开始程序在计算百分位时需要基准期边界外的信息已知 bug 会让程序停下。手册给的绕行办法是在基准期开始前补一条缺测记录让数据首年早于基准期首年。这个补记录不是改真实数据只是给脚本一个边界缓冲后面 QC 会把它当缺测处理。注意基准期不是越长越好也不是随便截一段它要和你的研究口径一致填进去之后不要中途换。4.2 纬度符号与 GSL 的南北半球口径纬度参数用来判断站点在哪个半球南半球填负值。这个符号影响 GSL 生长季长度的计算窗口北半球按 1 月到 12 月找跨越 6 天的窗口南半球按 7 月到次年 6 月的口径走。符号填反了GSL 不会报错但结果会偏到另一个季节窗口这种错误最隐蔽。我的习惯是填完纬度先看一眼绝对值是否在合理范围再看正负号是否符合站点位置。有些用户会把纬度写成度分秒或者带方向字母脚本只认十进制度数值。南半球用负数不要写S。如果你的站点在赤道附近纬度接近 0符号影响很小但也不要省掉保持参数完整后面批量复现时不容易乱。4.3 P、SUmm、FDmm、TRmm、IDmm 五个自定义阈值怎么填参数窗口里有一个日降水阈值 P用来算日降水量超过 P 的天数也就是 Rnn 这一类指标。单位是毫米填你研究口径需要的阈值不要填成厘米。四个自定义温度阈值分别对应日最高温和日最低温的上限、下限手册把它们映射到 SUmm、FDmm、TRmm、IDmm其中mm是你填的阈值。填的时候不要只看标签要先确认它用的是 TX 还是 TN。参数位置作用对象常见映射填写要点日降水阈值 P日降水Rnn单位毫米按研究口径填日最高温上限TXSUmm日最高温超过该值的天数日最高温下限TXFDmm日最高温低于该值的天数日最低温上限TNTRmm日最低温超过该值的天数日最低温下限TNIDmm日最低温低于该值的天数这张表是填参数时最该放在手边的对照。实际窗口文字可能因版本略有差异但判断逻辑不变先看它问的是日最高温还是日最低温再看是上限还是下限。填错一个后面指数名对得上数值含义却变了。批量站点时我一般把五个阈值写在同一个配置表里逐站核对避免手填串行。4.4 指数选择与输出目录indices 与 plots 的文件命名参数填完窗口会列出可计算的指数默认全选。只算需要的就取消不需要的选项越多跑得越久。计算完成后指数序列放在indices子目录文件格式是逗号分隔的表格文件名一般是站点前缀加指数名图放在plots子目录JPEG 格式文件名规则相同。图里会画年序列、线性最小二乘趋势线和局部加权回归线并在图上给出线性趋势拟合统计量。这里要提醒一句图上的趋势线是可视化辅助不是自动结论。序列有缺口、异常值没处理、基准期口径变了趋势线都可能跟着变。我的做法是先看indices里的数值序列再对照plots里的图不要只看图就下判断。# 检查指数输出目录和文件是否生成 list.files(indices, pattern \\.csv$, full.names FALSE) list.files(plots, pattern \\.jpg$, full.names FALSE) # 读一个指数序列确认列名和年份范围 tx90p - read.csv(indices/stationA_TX90p.csv) head(tx90p) tail(tx90p)list.files()用正则筛选 CSV 和 JPGfull.names FALSE只看文件名方便核对命名规则。读指数序列时文件名按实际站点和指数名替换head()和tail()看年份范围是否覆盖你要求的起止年。如果目录为空先回控制台看有没有报错再确认指数选择窗口里是不是全取消了。参数上pattern只是过滤不会改变文件本身文件名对不上时先看站点前缀是否和输入文件一致。5. 常见问题排查基准期首年撞车、异常值误报和慢机等待RClimDex 的报错不算多但几个典型问题会反复出现。下面按“现象、原因、解决”记录我踩过的坑尽量把排查顺序写清楚。遇到问题时不要先改脚本先看控制台第一条报错和输出目录大多数问题都能从数据和路径上找到线索。5.1 基准期首年与数据首年相同已知 bug 的绕行办法现象参数填完一点计算程序停住或直接退出控制台没有特别清楚的业务报错。原因手册明确记录了一个已知 bug当可用数据的第一年与基准期第一年相同时百分位类指数计算需要基准期边界外的数据程序会停止。解决在输入文件最前面补一条基准期首年减一的缺测记录让数据首年早于基准期首年。补的记录用-99.9标记缺测QC 会把它转成NA不会进入真实统计。# 在输入文件首行补一条缺测记录年份取基准期首年减一 # 列顺序按 stationA.txtYEAR MONTH DAY PRCP TMAX TMIN # 下面只演示写入逻辑实际文件名和年份按你的数据替换 line - 1980 12 31 -99.9 -99.9 -99.9 writeLines(c(line, readLines(stationA.txt)), stationA_fixed.txt)这段代码把一条缺测记录拼到原文件前面生成新文件不动原始文件。line里的年份要按你的基准期首年减一填月日填基准期开始前最后一天后面三个-99.9对应降水、最高温、最低温。writeLines()会覆盖输出文件所以输出文件名不要和输入同名。补完之后用新文件重新走 QC不要直接拿旧 QC 结果算指数。5.2 弹出“不合理值”但数据看着正常先查编码与日序现象QC 弹窗提示存在不合理值打开原始文件看降水没有负数最高温也不小于最低温。原因常见是列错位、缺测编码不统一、日期格式错位或者小数点分隔符不对。比如把-99.9写成了-999脚本不会把它当缺测可能当成一个极端低值再比如某行少了一列后面的温度列整体前移TXTN 就出现了。解决先用 R 读前几行和尾几行检查列数和每列取值范围再查是否有空行、制表符和逗号混用。排查时不要只看弹窗汇总要打开log里的不合理值清单看它标的是哪一天、哪个值。如果错误集中在某一段连续日期优先怀疑录入或拼接如果分散在全序列优先怀疑缺测编码和列顺序。改完原始文件后重新 QC确认清单减少或消失再继续。5.3 计算跑很久或没有输出看控制台、看目录、看选择项现象点了指数计算界面像卡住几分钟没有结果或者等很久后indices目录还是空的。原因指数选得太多、机器较慢、输出目录没有写权限或者控制台其实已经报了错但被窗口挡住。手册也提醒过某些较慢的机器上 QC 可能要几分钟指数计算根据选择项也可能跑一阵。解决先把 R 控制台窗口调到前面看最后几行再确认indices和plots目录是否生成最后检查是不是全选了不需要的指数。我的习惯是第一次只选三到五个指数做冒烟测试确认输出正常再全选或按需批量。目录权限问题在受控电脑上很常见工作目录不要放在只读盘或同步盘里。计算过程中不要反复点按钮重复触发会让输出文件互相覆盖后面很难判断哪一版是有效结果。5.4 异常值被标了但没改QC 输出文件的边界要认清现象异常值清单里有记录打开 QC 后数据发现原值还在以为 QC 没生效。原因QC 只把缺测和不合理值转成NA被标记的异常值不会自动修改这是手册写明的行为。解决把异常值清单当抽查列表回原始文件判断是否修正修正后重新跑 QC让新的indcal.csv覆盖旧文件。不要直接在 QC 后数据上改完就算指数因为下一次重跑 QC 会把你的手改覆盖掉。这个边界也影响结果解释如果异常值没处理指数序列里会保留这些极端值趋势线可能被拉偏。我的做法是先在日志里标出异常值集中的年份再决定是回改原始数据还是在结果说明里注明未修正。无论哪种都要保持口径一致不能这个站改、那个站不改。6. 把 QC 日志当黑匣子读单站复现的五个验收点与结果抽样QC 日志不是跑完就丢的中间文件它更像黑匣子数据入口有没有歪、异常值集中在哪、缺测红点是不是成片都能从日志和四张 PDF 里倒推出来。单站复现时我一般不看“跑完了没有”而看五个验收点输入列数、缺测编码、异常值清单长度、QC 后数据年份范围、指数输出文件是否齐全。验收点看什么不通过时的动作输入列数年、月、日、降水、TX、TN 是否对齐回原始文件查空行和分隔符缺测编码-99.9是否统一统一编码后重跑 QC异常值清单标记是否集中在某几年抽查原始记录决定是否回改QC 后数据年份范围是否覆盖基准期和计算期补记录或修日期边界指数输出indices与plots是否成对生成看控制台报错和选择项抽样对照我一般做两步先挑一个指数比如 TX90p看首尾年份数值是否在合理范围再回原始日值里抽几天手工核对“超过阈值的天数”能不能对上。这个动作不复杂但能抓住大部分列错位和阈值填错。批量站点之前单站必须走完这一遍不然一个站点的参数错误会被复制到几十个站。从那以后我每次加载rclimdex.r都先跑环境体检再按“输入列数、缺测编码、QC 日志、基准期首年、指数抽样”的顺序强制走一遍不再凭感觉点窗口。需要原始手册附录表时按资源包内目录找附录 A 到 E先对列顺序和参数口径再谈结果。希望帮到你。本文还有配套的精品资源点击获取
阅读完成 · 觉得有帮助？
