centroid function

Centroid of spherical polygons

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).

centroid(x, ...)

Arguments

• x: SpatialPolygons* object, or a 2-column matrix or data.frame reprenting a single polgyon (longitude/latitude)
• ...: Additional arguments. None implemented

Note

For multi-part polygons, the centroid of the largest part is returned.

Returns

A matrix (longitude/latitude)

See Also

area, perimeter

Author(s)

Robert J. Hijmans

Examples

pol <- rbind(c(-180,-20), c(-160,5), c(-60, 0), c(-160,-60), c(-180,-20))
centroid(pol)