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

if (FALSE) { # \dontrun{
# this would require the testpoints_1000 data from the EJAM package:
  isles <- which(latlon_is.islandareas(lat = testpoints_1000$lat, lon = testpoints_1000$lon))
  mapfast(testpoints_1000[isles, ]) # c(213,785)
  which(!(latlon_is.usa(lat = testpoints_1000$lat, lon = testpoints_1000$lon)))
} # }