Get NOAA wind storm tabular data, metadata, or shp files from IBTrACS
storm_shp(basin = NULL, storm = NULL, year = NULL, type = "points",
overwrite = TRUE)storm_shp_read(x)
storm_data(basin = NULL, storm = NULL, year = NULL,
overwrite = TRUE, ...)
storm_meta(what = "storm_columns")
(character) A basin name, one of EP, NA, NI, SA, SI, SP, or WP.
(character) A storm serial number of the form YYYYJJJHTTNNN. See Details.
(numeric) One of the years from 1842 to 2014
(character) One of points or lines. This gives shp files with points, or with lines.
(logical) To overwrite the path to store files in or not,
Default: TRUE
Output from storm_shp
, a path to shp file to read in.
Curl options passed on to HttpClient
(optional)
(character) One of storm_columns or storm_names.
A tibble (data.frame)
We use rappdirs to store files, see
user_cache_dir
for how
we determine the directory on your machine to save files to, and run
rappdirs::user_cache_dir("rnoaa/storms")
to get that directory.
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.
# 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)
# }
Run the code above in your browser using DataLab