Check lat lon coordinates to see if each is approx. in general area of US Island Areas Guam, USVI, Amer Samoa or N Marianas

See islandareas

Usage

latlon_is.islandareas(lat, lon, exact_but_slow_islandareas = FALSE)

Arguments

lat

vector of latitudes

lon

vector of longitudes

exact_but_slow_islandareas

optional logical, set it to TRUE to check each point vs boundaries in states_shapefile to identify which ones are in Island Areas according to that shapefile. The default method here is much faster, but just checks if a point is within a bounding box that should approximate each of the Island Areas, found in the object islandareas.

Value

vector of TRUE / FALSE values indicating a given lat lon pair is approximately in one of the rough bounding boxes that includes the 4 Island Areas.

See also

is.island() latlon_is.usa() latlon_is.islandareas() latlon_is.available() latlon_is.possible()

Examples

  isles <- stateinfo2[ EJAM:::latlon_is.islandareas(lat = stateinfo2$lat, lon = stateinfo2$lon) & !is.na(stateinfo2$lat), ]
  mapfast(isles)