centroid: Centroid of spherical polygons

In geosphere: Spherical Trigonometry

Description

Compute the centroid of longitude/latitude polygons. Unlike other functions in this package, there is no spherical trigonomery 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	a 2-column matrix (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) 

Example output

           lon      lat
[1,] -133.3333 -23.8934

geosphere documentation built on May 2, 2019, 5:16 p.m.