points.on.land: Return coordinates over land Returns the set of input...

In JeffWeinell/misc.wrappers: Miscellaneous Wrapper Functions for Spatial Population Genomics

Return coordinates over land Returns the set of input coordinates that intersect with land. Land is considered to incude all areas covered by a feature in the rnaturalearth::ne_countries 10 meter resolution dataset. Input coordinates are are assumed to be WGS84 (EPSG:4326)

Description

Return coordinates over land

Returns the set of input coordinates that intersect with land. Land is considered to incude all areas covered by a feature in the rnaturalearth::ne_countries 10 meter resolution dataset. Input coordinates are are assumed to be WGS84 (EPSG:4326)

Usage

## S3 method for class 'on.land'
points(x, return.as = "data.frame")

Arguments

x	Matrix or data frame with longitude and latitude decimal degree in the first two columns, respectively.
return.as	Character string with class to use for object returned. Default "data.frame". Can also be "matrix" or "SP" (SpatialPoints).

Value

An object with class equal to the value of 'return.as', which contains the subset of input points that occur on land, or more specifically, not in an ocean.

Examples

library(misc.wrappers)
# Sample 50 points from 10-degree radius area with center located on land somewhere between -50 and 50 degrees latitude.
coords50 <- rcoords(r=10,size=50,limits=c(-180,180,-50,50))
# Use method 1 to create the outer set of points and generate a figure with points on map. Method 1 may not work well for points on continents vs. island archipelagos.
coords50_outer1 <- create.outer(coords=coords50, method=1, output.path="coords50_outer1.txt", plot.output.path="coords50_outer1.pdf")
# Use method 2
coords50_outer2 <- create.outer(coords=coords50, method=2, output.path="coords50_outer2.txt", plot.output.path="coords50_outer2.pdf")
# Use method 3
coords50_outer3 <- create.outer(coords=coords50, method=3, output.path="coords50_outer3.txt", plot.output.path="coords50_outer3.pdf")

JeffWeinell/misc.wrappers documentation built on Sept. 20, 2023, 12:42 p.m.