Learn R Programming

rnoaa (version 0.6.0)

storm_shp: Get NOAA wind storm tabular data, metadata, or shp files from IBTrACS

Description

Get NOAA wind storm tabular data, metadata, or shp files from IBTrACS

Usage

storm_shp(basin = NULL, storm = NULL, year = NULL, type = "points", path = "~/.rnoaa/storms", overwrite = TRUE)
storm_shp_read(x)
storm_data(basin = NULL, storm = NULL, year = NULL, path = "~/.rnoaa/storms", overwrite = TRUE, ...)
storm_meta(what = "storm_columns")

Arguments

basin
(character) A basin name, one of EP, NA, NI, SA, SI, SP, or WP.
storm
(character) A storm serial number of the form YYYYJJJHTTNNN. See Details.
year
(numeric) One of the years from 1842 to 2014
type
(character) One of points or lines. This gives shp files with points, or with lines.
path
(character) A path to store the files, Default: ~/.rnoaa/storms
overwrite
(logical) To overwrite the path to store files in or not, Default: TRUE.
x
Output from storm_shp, a path to shp file to read in.
...
Curl options passed on to GET (optional)
what
(character) One of storm_columns or storm_names.

Details

International Best Track Archive for Climate Stewardship (IBTrACS)

Details for storm serial numbers:

  • YYYY is the corresponding year of the first recorded observation of the storm
  • JJJ is the day of year of the first recorded observation of the storm
  • H is the hemisphere of the storm: N=Northern, S=Southern
  • TT is the absolute value of the rounded latitude of the first recorded observation of the storm (range 0-90, if basin=SA or SH, then TT in reality is negative)
  • NNN is the rounded longitude of the first recorded observation of the storm (range 0-359)

For example: 1970143N19091 is a storm in the North Atlantic which started on May 23, 1970 near 19N 91E

See http://www.ncdc.noaa.gov/ibtracs/index.php?name=numbering for more.

The datasets included in the package storm_names, and storm_columns may help in using these storm functions.

References

http://www.ncdc.noaa.gov/ibtracs/index.php?name=wmo-data

Examples

Run this code
## Not run: 
# # Metadata
# head( storm_meta() )
# head( storm_meta("storm_columns") )
# head( storm_meta("storm_names") )
# 
# # Tabular data
# ## Get tabular data for basins, storms, or years
# storm_data(basin='WP')
# storm_data(storm='1970143N19091')
# storm_data(year=1940)
# storm_data(year=1941)
# storm_data(year=2010)
# 
# # shp files
# ## storm_shp downloads data and gives a path back
# ## to read in, use storm_shp_read
# res <- storm_shp(basin='EP')
# storm_shp_read(res)
# 
# ## Get shp file for a storm
# (res2 <- storm_shp(storm='1970143N19091'))
# 
# ## Plot shp file data, we'll need sp library
# library('sp')
# 
# ### for year 1940, points
# (res3 <- storm_shp(year=1940))
# res3shp <- storm_shp_read(res3)
# plot(res3shp)
# 
# ### for year 1940, lines
# (res3_lines <- storm_shp(year=1940, type="lines"))
# res3_linesshp <- storm_shp_read(x=res3_lines)
# plot(res3_linesshp)
# 
# ### for year 2010, points
# (res4 <- storm_shp(year=2010))
# res4shp <- storm_shp_read(res4)
# plot(res4shp)
# ## End(Not run)

Run the code above in your browser using DataLab