留言板
- sunjiajia
- 指的是spi值，谢谢您的回复，让我明白很多，还要虚心学习
- lq_nuist
- #井号表示解释说明，在R中是不运行的，这句话是给使用者看的，我的站点纬度是44.25，你要根据你研究的站点的纬度位置修改file$PET <- thornthwaite(file$TMED, 纬度数值)最后面的数值。
sunjiajia: 谢谢，有几个问题，能帮忙解答一下吗？
（1）请问该方法是得到一个站点的spi数值，那么要计算一个地区的多个站点，是计算平均值还是总和？
（2）“44.25为纬度位 ...
3个月尺度的spi就是用前3个月的数据计算得到的，从第三个数据开始才有3个数值用来计算，前两个自然就没有数据，没有数据就是缺测，缺测就是NA。
有批处理的方法，具体要看你的原始数据和要算的时间尺度。基本思想是：程序不变，改变你的输入数据。比如你的原始数据是多年多个站点（小于等于12）的年总降水和年平均温度，可以把月份看做站点，时间尺度为1，得到的就是多年多个站点的年尺度spi值。
计算平均值还是总和你是指计算spi吗？我不太明白
- lq_nuist
- 可以，不过这也是我们上课实习的内容，我研究的也不深
sunjiajia: hello,看到您的帖子，"R中计算SPEI和SPI程序"http://bbs.06climate.com/forum.php?mod=viewthread&tid=58249&extra=&page=1想请教您一些问 ...