Not Logged In

stormtracks 0.1.2

Warning, this project is in Alpha, it will be hard to get working and you can’t trust the documentation!

The main aim of this project is to develop an algorithm to detect and track tropical cyclones in the C20 Reanalysis Project. The algorithm will then be run against these data looking for trends in the time range of the data, from 1871 to 2013. Currently the algorithm uses vorticity maxima to locate potential candiates for cyclones, and tracks these maxima using a neareat neighbour approach from time frame to time frame. These tracks are then matched to best tracks from the IBTrACS catalogue, which it is hoped will allow for automated categorisation of these tracks (and the corresponding cyclones) through a machine learning technique: Support Vector Machines (SVMs).

Installing and Running

To run the analysis, you will first need to download some C20 Reanalysis data and the ibtracs data (download_2005.py):

importstormtracks.downloadasdl# Data will be saved to ~/stormtracks/data/dl.download_ibtracs()# N.B. one year is 4.2 GB of data! This will take a while.# It will download 3 files, two with the wind velocities at ~sea level (u9950/v9950)# and Pressure at Sea Level (prmsl).dl.download_full_c20(2005)

importdatetimeasdtimportpylabaspltfromstormtracks.c20dataimportC20Datafromstormtracks.ibtracsdataimportIbtracsDatafromstormtracks.detectimportGlobalCyclonesfromstormtracks.matchimportmatch2importstormtracks.plottingaspl# Create a wrapper for the C20 Reanalysis data.c20data=C20Data(2005)c20data.first_date()# Plot PSL for 1st of Jan 2005.pl.plot_on_earth(c20data.lon,c20data.lat,c20data.psl)plt.show()# Plot vorticity for 1st of Jan 2005.pl.plot_on_earth(c20data.lon,c20data.lat,c20data.vort)plt.show()# Load IBTrACS data for 2005.ibtracs=IbtracsData()best_tracks=ibtracs.load_ibtracks_year(2005)# Run the analysis for 2005.gdata=GlobalCyclones(c20data,ensemble_member=0)gdata.track_vort_maxima(dt.datetime(2005,6,1),dt.datetime(2005,7,1))# Match the generated tracks agains the best tracks.matches=match2(gdata.vort_tracks_by_date,best_tracks)good_matches=[maformainmatches.values()ifma.av_dist()<5andma.overlap>6]forgmingood_matches:pl.plot_match(gm,None)plt.show()