R/mapasp.R

Defines functions mapasp

Documented in mapasp

mapasp <- function(data, xlim = bbox(data)[1,], ylim = bbox(data)[2,]) {
	# calculates aspect ratio for levelplot of geographic data,
	# using proportial units (compare eqscplot)
	if (!is(data, "Spatial"))
		stop("cannot extract coordinates bounding box from data")
	if (!(is.na(slot(slot(data, "proj4string"), "projargs")) || is.projected(data)))
		return( (diff(ylim)/diff(xlim)) / cos((mean(ylim) * pi)/180))
	if (is.R() && version$major >= 2)
		return("iso")
	else
		return(diff(ylim)/diff(xlim))
}

Try the sp package in your browser

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

sp documentation built on Nov. 27, 2023, 1:08 a.m.