RClimDex (1.0)极端气候指数计算软件用户手册张学斌 Feng Yang加拿大环境部气候研究中心2004年9月10日南京信息工程大学遥感学院陈昌春译注2013.8作者致谢RClimDex 由Xuebin Zhang(张学斌)and Feng Yang(加拿大气象局气候研究部)开发与维护,最初的开发由加拿大国际发展办事处通过《加中气候协作项目,C5》资助。
Lisa Alexander, Francis Zwiers, Byron Gleason, David Stephenson, Albert Klan Tank, Mark New, Lucie Vincent与Tom Peterson对R包的开发与测试作出了重要贡献。
CCl/CLIVAR ETCCDMI的有关研讨会也对RClimDex的改进提供了宝贵的意见。
.译者的话原英文说明中所介绍的下载网址链接已无效,新网址包括http://www.pcic.uvic.ca/tools-and-data/climdex/software.shtmlRClimdex可计算极端气候指数27项,以下摘录来自一硕士论文《内蒙古地区极端气候事件时空变化及其与NDVI的相关性》(使用RClimdex软件)的15项指数名称翻译及一段简要说明。
1.指数名称、解释、单位FD0 霜日一年中日最低温<0℃的日数天SU25 夏日日数日最高气温>25℃的日数天GSL 作物生长期连续6 日>5℃或<5℃的时间跨度天TN10p 冷夜日数日最低气温<10%分位值的日数天TN90p 暖夜日数日最低气温>90%分位值的日数天TX10p 冷昼日数日最高温<10%分位值的日数天TX90p 暖昼日数日最高温>90%分位值的日数天WSDI 热持续指数连续6 日最高温在90%分位值日数天CSDI 冷持续指数连续6 日最低温在10%分位值日数天RX5day 5 日最大降水量每月内连续五日的最大降水量 mmCDD 持续干燥指数日降水量<1mm 的最长连续日数天CWD 持续湿润指数日降水量≥mm 的最大持续日数天SDII 普通日降水强度降水量≥1mm 的总量与日数之比 mmR10 强降水日数每年日降水量>=10mm 的总日数天R95pTOT 强降水量 95%分位值强降水之和 mm2.简要说明在应用RClimDex 处理数据之前,必须确保每个站点的数据以文本格式储存,并且储存的气象数据必须按照年、月、日、24 小时日降水量、日最高气温、日最低气温等顺序排列,各记录项之间通过空格将其隔开。
由于研究的气象记录年限跨度很大,难免有些记录值缺测,因此程序要求所有的缺测值都统一设为-99.9。
启动R 编辑器,读入RClimDex 程序,在桌面窗口的引导下读入某站点气象数据的文本文件,程序将对文件数据进行检验处理,用于提高所用数据的质量。
检验处理主要涉及几个方面的错误记录检查:(1)日最低温>日最高温(2)日降水量<0mm(3)记录值严重偏离本地区气象实际情况,即超出3 倍标准差的值定义为出界值,通过人工检查同相邻站点的记录进行比对,合理的保留,不合理的按缺测进行处理。
进行完质量控制处理后就可以对数据进行计算处理,得出一系列的excel 格式的文件结果。
目录1.简介2.R语言安装2.1How to install R2.2How to run R3.如何使用RClimDex?3.1RClimDex的加载与使用3.2数据质量控制3.3指数计算4.已知的bugs5.Bug报告Appendix A: 气候指数列表Appendix B: 输入数据格式Appendix C: 指数定义Appendix D: Threshold and in-base period temperature indices calculation Appendix E: R for Windows FAQ1. 介绍加拿大气象研究中心的Xuebin Zhang(张学斌,华人)与Feng Yang 开发的RClimDex(1.0)是对ClimDex软件(由加拿大国家气候数据中心的Byron Gleason 开发)的R语言改写与改进。
RClimDex(1.0)是基于R 编辑器开发的用于计算多种极端气候指数。
ClimDex是以Excel为基础的程序,意在提供容易使用的软件包用于计算监测气候变化的气候极值。
它由加拿大国家气候数据中心的Byron Gleason开发,已经用于CCl/CLIVAR workshops on climate indices from 2001.使用者可以在窗口界面的提示和引导下将研究区域的气象记录值文件进行处理和计算,结果将以excel 文件的形式输出。
ClimDex项目软件开发之初就强调运行环境不依赖于特定的操作系统。
因此使用R语言就就非常自然,因为R语言是健壮的、强大的统计分析与绘图软件。
它既可运行于Windows 系统,也可运行于Unix系统。
2003年发现,在ClimDex及其他程序中,计算百分位基础的温度指数在指数序列组中出现了不一致性,修正工作需要的引导(又称“自展”)工作在EXCEL中难以进行,这使开发ClimDex的R语言包显得更加迫切。
RClimDex (1.0)有着比较友好的操作界面,计算由CCl/CLIVAR推荐及可由用户定义门槛值的重要极端气候指数27 种,包括16 个极端温度指数和11 个极端降水指数,这27 个极端气候指数都是世界气象组织气候委员会推荐使用的核心极端气候指数,已包括了ClimDex (Version 1.3)中的几乎全部指数。
RClimDex (1.0)在R 1.84下进行开发,它应该可以运行于R 1.84及其后的版本。
开发极端气候指数软件的目的之一是用于气候监测研究,这要求各项指数必须均一化(一致化)。
数据均一化在本版中已经设想并未实施。
当前的RClimDex版只包括一个简单的数据质量控制程序,由ClimDex计划(译者注:ClimDex是一个专门的计划,有专门的网站。
当前分别有EXCEL版本、R语言版本、FORTRAN 语言版本的ClimDex软件)提供。
对于ClimDex的用户,我们希望,在计算极端指数前已进行数据质量控制。
本用户手册提供逐步操作的说明。
2. R语言的安装与运行(译者注:很方便,几乎不需要设置,一路点击即可)R is a language and environment for statistical computing and graphics. It is a implementation of the S language developed by John Chambers and colleagues Laboratories (formerly AT&T, now Lucent Technologies). S-plus provides a commercial implementation of the S language.2.1 How to install RRClimDex需要R包与图形用户界面TclTk。
R的安装过程非常简单:1) 登录网站, 2) Follow the links to download the most recent version of R for your computer operating system from any mirror site of CRAN.对于Microsoft Windows (95, 98, 2000, and XP)系统,下载Windows版安装程序。
进行安装后,R将出现在你的机器中,并且在你的桌面增加快捷图标。
TclTk 程序已包含在R 1.9.0及以后版本中。
对于Linux系统, download proper precompiled binaries and follow the instruction to install R. For other unix systems, you many need to download source code and compile it yourself.2.2 How to run R在 Windows系统环境下,双击桌面的R快捷图标,或者在“开始”菜单中点击, R用户界面将会出现。
第一次使用时你可能需要配置称做“HOME”的环境变量。
详情可参见附录E中用户常问的问题(FAQ)。
在unix 环境下,直接在 R程序界面运行。
欲在R界面退出,输入q()即可。
在Windows系统下,还可以选择“File”菜单然后点击“Exit”选项。
3. 如何使用RClimDex3.1 RClimDex的加载1. 直接在R窗口中操作,RClimDex可以通过下拉菜单加载R包。
在RGui菜单下选“File” ,再选“Source R code”。
这将弹出一个窗口,你应指向你保存RClimDex的位置。
2.在R语言命令提示符“>”下,输入源文件(R包)名称:source(“rclimdex.r”)。
如文件不在缺省目录(什么目录?)中,在rclimdex.r名称之前,你可能需要用全路径名,此时在提示符下输入:source("f:/rclimdex/rclimdex.r"(注意:全路径名之间的间隔符为unix下采用的“/”而不是windows操作系统命令行(即以前的DOS操作系统)下的“\”)。
这将加载RClimDex。
一旦源程序加载成功,RClimDex主菜单就会出现:译者注:英文原文中,3.1中1、2方法出现次序相反。
为便利用户使用作了调换。
3.2. 加载数据,运行质量控制模块(QC)数据质量控制模块是计算极端指数的预备性处理工作。
RClimDex质量控制模块执行如下功能: 1) 用R语言认可的数据格式代替缺失值(假设原数据以-99.9表示),如用“NA”表示“不可得”;2) 用“NA”代替所有的不合理值。
这些值包括 a) 小于0的日降水量;b) 小于日最小气温的日最大气温。
此外,质量控制模块还能识别日最大气温与日最大气温的离群值(不正常值)。
这些离群值是由用户定义的区域之外的日特征值。
当前,这个区域被定义为均值±日特征值的n倍标准差(翻译待审核!),即[mean –n*std, mean+n*std]。
std表示日特征值的标准差,n 表示一个用户的输入,mean(均值)表示气候特征的日值。
在RClimDex主菜单中选择“Load Data and Run QC” ,打开的窗口形式如下。
用户能够选取用于计算极端指数的气候数据文件。