correctCoordGPS: Correct the GPS coordinates

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

Description

This function build the most probable GPS coordinates of the plot corners from multiple GPS measurements.

Usage

correctCoordGPS(longlat = NULL, projCoord = NULL, coordRel, rangeX,
  rangeY, maxDist = 10, drawPlot = FALSE, rmOutliers = FALSE)

Arguments

longlat	(optionnal) data frame with the coordinate in longitude latitude (eg. cbind(longitude, latitude)).
projCoord	(optionnal) data frame with the projected coordinate in X Y
coordRel	data frame with the relative coordinate in the same order than the longlat or projCoord
rangeX	a vector of length 2 giving the range for plot relative X coordinates
rangeY	a vector of length 2 giving the range for plot relative Y coordinates
maxDist	a numeric giving the maximum distance above which GPS measurements should be considered as outliers (by default 10 m)
drawPlot	a logical if you want to display a graphical representation
rmOutliers	a logical if you want to remove the outliers from coordinates calculation

Details

You must give one parameter between longlat and projCoord

Value

If you there are no outliers or rmOutliers = TRUE, a list with:

• corner: a matrix with the coordinates of the corners

• polygon: a spatial polygon

• outliers: Coordinates lines considered as outliers, if any

• codeUTM: the code UTM for the coordinate if the parameter longlat is set

Author(s)

Arthur PERE, Maxime REJOU-MECHAIN

Arthur PERE, Maxime REJOU-MECHAIN

Examples

projCoord <- data.frame(
  X = c(
    runif(5, min = 9, max = 11), runif(5, min = 8, max = 12),
    runif(5, min = 80, max = 120), runif(5, min = 90, max = 110)
  ),
  Y = c(
    runif(5, min = 9, max = 11), runif(5, min = 80, max = 120),
    runif(5, min = 8, max = 12), runif(5, min = 90, max = 110)
  )
)
projCoord <- projCoord + 1000
coordRel <- data.frame(
  X = c(rep(0, 10), rep(100, 10)),
  Y = c(rep(c(rep(0, 5), rep(100, 5)), 2))
)

aa <- correctCoordGPS(
  projCoord = projCoord, coordRel = coordRel,
  rangeX = c(0, 100), rangeY = c(0, 100)
)
bb <- correctCoordGPS(
  projCoord = projCoord, coordRel = coordRel,
  rangeX = c(0, 100), rangeY = c(0, 100), rmOutliers = TRUE
)
## Not run: 
correctCoordGPS(
  projCoord = projCoord, coordRel = coordRel,
  rangeX = c(0, 100), rangeY = c(0, 100), drawPlot = TRUE
)

## End(Not run)

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