centroid: Centroid of spherical polygons

In geosphere: Spherical Trigonometry

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)

Note

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

Author(s)

Robert J. Hijmans

See Also

area, perimeter

Examples

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

geosphere documentation built on Nov. 16, 2022, 1:06 a.m.