geosphere (version 1.5-10)

centroid: Centroid of spherical polygons

Description

Compute the centroid of longitude/latitude polygons. Unlike other functions in this package, there is no spherical trigonometry involved in the implementation of this function. Instead, the function projects the polygon to the (conformal) Mercator coordinate reference system, computes the centroid, and then inversely projects it to longitude and latitude. This approach fails for polygons that include one of the poles (and is rather biased for anything close to the poles). The function should work for polygons that cross the -180/180 meridian (date line).

Usage

centroid(x, ...)

Arguments

x

SpatialPolygons* object, or a 2-column matrix or data.frame reprenting a single polgyon (longitude/latitude)

...

Additional arguments. None implemented

Value

A matrix (longitude/latitude)

See Also

area, perimeter

Examples

# NOT RUN {
pol <- rbind(c(-180,-20), c(-160,5), c(-60, 0), c(-160,-60), c(-180,-20))
centroid(pol) 
# }