爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 20266|回复: 21

[源代码] 科罗拉多大学可在线运行的小波分析程序(fortran)

[复制链接]

新浪微博达人勋

发表于 2013-6-1 22:31:55 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 hillside 于 2013-6-1 22:51 编辑

      科罗拉多大学的Christopher Torrence and Gilbert Compo博士在网页提供了fortran、idl、matlab三个不同版本的小波分析程序。在线运行版本我没有测试。

wave_fortran.zip(见本帖附件)
Includes all of the following files...
Sample output
wavetest.f
wavelet.f
cfftpack.f
chisqr.f
sst_nino3.dat
   
   由http://www.tjporter.webs.com/linksdownloads.html网页Trevor J. Porter博士推荐, Trevor J. Porter介绍 :Drs. Christopher Torrence and Gilbert Compo provide a tool for wavelet analysis. Note that you can also perform your wavelet analysis online at Torrence and Compo's site.
    以下内容见于Trevor J. Porter博士网页:
Climate data, proxy data, analysis tools, etc.
International Tree-Ring Databank - A repository for >3000 tree-ring chronologies from around the world, all available for download free of charge. Don't forget to acknowledge the data contributors when using these data. [LINK]

Lamont-Doherty Earth Observatory Tree-Ring Lab - Aside from being one of the world's premier tree-ring labs, LDEO-TRL distributes many useful dendro programs written by Paul Krusic and Ed Cook including the widely used detrending program ARSTAN. [LINK]

The Science of Tree-Rings - Created and maintained by Henri Grissino-Mayer, this website is probably the most valuable dendro resource on the Web. There you can find info on pretty much all things related to dendrochronology (e.g., conference, software, job postings, etc.). [LINK]

NOAA Paleoclimatology Program - To help satisfy the objectives of the International Geosphere-Biosphere Programme, NOAA is host to numerous proxy-based paleoclimate and paleoenvironmental reconstructions. Visit their site to download these reconstructions free of charge. [LINK]

KMNI Climate Explorer - Developed by Dr. Geert Jan van Oldenborgh, this is by far my favourite online tool for climate analysis. Here you can explore spatial correlations between climate time-series (input your own datasets if you want) and a large selection of pre-loaded climate fields. [LINK]

NCEP/NCAR Reanalysis 1 - A gridded dataset developed from many global meteorological observations (land surface stations, ships, rawinsonde, aircraft, satellite, and other sources) which have been 're-analysed' by a Numerical Weather Prediction model in order to provide a spatially-continuous meteorological dataset covering 1948-present. [LINK]

HYSPLIT - NOAA's Hybrid Single Particle Lagrangian Integrated Trajectory Model allows you to calculate simple air parcel trajectories or more complex dispersion and particulate deposition. HYSPLIT works with the NCEP/NCAR Reanalysis 1 dataset and can be run on a PC, MAC, or online. [LINK]


Useful Matlab links
Matlab is a powerful, user-friendly programming language that can do amazing things for your research. If you haven't tried it, you ought to! Thanks to a number of Matlab-savvy scientists, ready-made functions and scripts for many research applications are available for download. The following are webpages of a few researchers whose Matlab code (available online) I recommend...

Dave Meko - Dr. Meko provides a comprehensive suite of Matlab functions for dendrochronological and general time-series analyses. I highly recommend his toolbox, lecture notes, and assignments! [LINK]

Torrence and Compo - Drs. Christopher Torrence and Gilbert Compo provide a tool for wavelet analysis. Note that you can also perform your wavelet analysis online at Torrence and Compo's site. [LINK]

Aslak Grinsted - Dr. Grinsted provides a useful tool for cross-wavelet and wavelet coherence analysis.[LINK]

Tom O'Haver - Dr. O'Haver's Interactive Fourier Filter is a great tool to help with low-, high-, band-pass-, and band-stop filtering, and best of all, you can view the effect of filtering parameters on your time-series as you change them! [LINK]


My Matlab functions
Here I've made a few simple functions available that I've developed over the years in hopes that others may find them useful too. Feel free to modify or use them as is, but please read the license before each download. If you found these functions helpful, let me know.

STATISTICAL FUNCTIONS

[nancorr] - Takes an input array with 2+ columns of time-series data and outputs: a correlation matrix; a matrix of corresponding significance (p) values; and a matrix indicating the number of data-points available for calculating each correlation coefficient. No special Matlab toolboxes are required.

[rbareps] - Takes an input array with 2+ columns of time-series data and outputs the running RBAR (mean series intercorrelation) and EPS (expressed population signal) statistics based on a window-length and lag specified by the user. No special Matlab toolboxes are required.

[zscores] - Takes an input array with 2+ columns of time-series data (the first column must be an equal-interval time-vector, e.g., years) and outputs an equivalent array converted to z-scores relative to a period of time specified by the user. No special Matlab toolboxes are required.

FILTERING FUNCTIONS

[splfilter] - Takes an input time-series (e.g., tree-ring chronology) and filters it with a n-year cubic smoothing spline with a 50% frequency cut-off; the wavelength of the filter is user-specified; the filter can be low- or high-pass filter as specified by the user; the input array must contain an equal-interval time-vector (e.g., years) as the first column; subsequent columns contain the time-series. Note, this function calls the 'csaps' function which is built into Matlab's Spline Toolbox. If your version of Matlab doesn't have the Spline Toolbox, you can purchase it from Mathworks online.

[rmfilter] - Takes an input time-series (e.g., tree-ring chronology) and filters it with a n-year running mean; the length of the filter is user-specified; the filter can be low-pass or high-pass as specified by the user; the input array must contain an equal-interval time-vector (e.g., years) as the first column; subsequent columns contain the time-series. No special Matlab toolboxes are required.

CLIMATE DATA FUNCTIONS

[hadcrut3extract] - Extracts gridded temperature data from the CRUTEM3, CRUTEM3v, HadCRUT3, and HadCRUT3v datasets (.dat) and outputs a monthly formatted time-series of all available data corresponding to the latitude/longitude coordinates specified by the user; the function will provide individual grid cell averages or regional averages. No special Matlab toolboxes are required.

PLOTTING FUNCTIONS

[ciplot] - Plots a single time-series with the area between its 95% confidence intervals shaded grey. No special Mathworks toolboxes are required to run this function.








wave_fortran.zip

33.79 KB, 下载次数: 45, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2013-6-2 09:13:02 | 显示全部楼层
这个貌似我介绍过 嘿嘿
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-6-2 21:01:57 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-6-11 19:36:06 | 显示全部楼层
顶一个哈哦
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-7-25 10:48:59 | 显示全部楼层
哎 英语太差了  找不到IDL版本的哎
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-7-26 11:10:26 | 显示全部楼层
牛牛牛牛牛
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 10:40:11 | 显示全部楼层
您好,请问fortran程序中wavetest在调用wavelet时是怎么调用的啊?为什么我每次运行都提示链接错误呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-7 10:45:49 | 显示全部楼层
好神奇啊!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-12-2 11:29:29 | 显示全部楼层
楼主 你测试过了吗 把图贴出来看看啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-12-3 21:17:41 | 显示全部楼层
麦田_smile 发表于 2013-12-2 11:29
楼主 你测试过了吗 把图贴出来看看啊

我在本帖第二行已指出,我没有测试过。有兴趣的网友可以试试。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表