Learn R Programming

rnoaa (version 0.8.4)

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",
  overwrite = TRUE)

storm_shp_read(x)

storm_data(basin = NULL, storm = NULL, year = NULL, 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.

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 HttpClient (optional)

what

(character) One of storm_columns or storm_names.

Value

A tibble (data.frame)

File storage

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.

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)
# }

Run the code above in your browser using DataLab