R/helper.R

Defines functions .isPolygon .normalizeLonRad .normalizeLonDeg

# Author: Robert J. Hijmans
# April 2010
# version 1
# license GPL3

.normalizeLonDeg <- function(x) {
	(x + 180) %% 360 - 180
}

.normalizeLonRad <- function(x) {
	(x + pi) %% (2*pi) - pi 
}


.isPolygon <- function(x, fix=FALSE) {
	x <- stats::na.omit(x)
	if (nrow(x) < 4) {
		stop('this is not a polygon (insufficent number of vertices)')
	}
	if (length(unique(x[,1]))==1) {
		stop('All longitudes are the same (not a polygon)')
	}
	if (length(unique(x[,2]))==1) {
		stop('All latitudes are the same (not a polygon)')
	}
	if (! all(!(is.na(x))) ) {
		stop('polygon has NA values)')
	}
	if (! isTRUE(all.equal(x[1,], x[nrow(x),]))) {
		stop('this is not a valid (closed) polygon. The first vertex is not equal to the last vertex')	
	}
	return(x)
}

Try the geosphere package in your browser

Any scripts or data that you put into this service are public.

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