numberCorner: Get the UTM coordinates with the corner of the plot

In ArthurPERE/biomass: Estimating Aboveground Biomass and Its Uncertainty in Tropical Forests

Description

Get the UTM coordinates from the latitude and longitude of the corners of a plot. The function also assign a number to the corners in a clockwise or counterclockwise way, with the number 1 for the XY origin. Corner numbering is done as followed:

• axis X: the corner 1 to the corner 2

• axis Y: the corner 1 to the corner 4

Usage

numberCorner(longlat = NULL, projCoord = NULL, plot, origin, clockWise)

Arguments

longlat	(optionnal) data frame with the coordinates in longitude latitude (eg. cbind(longitude, latitude)).
projCoord	(optionnal) data frame with the projected coordinates in X Y
plot	A vector of codes (names) of the plots
origin	A logical vector with TRUE corresponding of the origin of the axis of each plot.
clockWise	A logical, whether the numbering should be done in a clockwise (TRUE) or counterclockwise (FALSE) way.

Value

A data frame with:

• plot: The code of the plot

• X: The coordinates X in UTM

• Y: The coordinates Y in UTM

• corner: The corner numbers

Author(s)

Arthur PERE, Maxime REJOU-MECHAIN

Examples

coord <- data.frame(X = c(0, 200, 0, 200), Y = c(0, 0, 200, 200)) + 5000
plot <- rep("plot1", 4)
origin <- c(FALSE, FALSE, TRUE, FALSE)

# if you turn clock wise
corner <- numberCorner(projCoord = coord, plot = plot, origin = origin, clockWise = TRUE)

# Plot the plot
plot(coord, asp = 1)
text(coord, labels = corner$corner, pos = 1)


# Using a counterclockwise way
corner <- numberCorner(projCoord = coord, plot = plot, origin = origin, clockWise = FALSE)

# Plot the plot
plot(coord, asp = 1)
text(coord, labels = corner$corner, pos = 1)

ArthurPERE/biomass documentation built on May 18, 2019, 2:33 a.m.