R/inside.reg.bc1.R

In geo: Draw and Annotate Maps, Especially Charts of the North Atlantic

#' Determine which bormicon (or gadget) region data belong to.
#' 
#' Determine which bormicon (or approximately gadget) region data belong to.
#' 
#' 
#' @param data Data set with coordinates in components \code{lat, lon}.
#' @return Returns original data with area code as an added component
#' \code{area}.
#' @note Needs further elaboration ?
#' @seealso Calls \code{\link{geoinside}}, called by
#' \code{\link{inside.reg.bc}} and \code{\link{Reitur2Svaedi1to10}}. Data set
#' \code{\link{reg.bc}} with bormicon-area outlines is used.
#' @keywords manip
#' @export inside.reg.bc1
inside.reg.bc1 <-
function(data)
{
        if(nrow(data) > 1)
                tmpdata <- data[, c("lat", "lon")]
        else tmpdata <- as.data.frame(data[, c("lat", "lon")])
        tmpdata$area <- rep(0, nrow(tmpdata))
        i <- 1
        ind <- geoinside(tmpdata, reg = geo::reg.bc[[i]], option = 0, robust = F)
        if(length(ind) > 0)
                tmpdata[ind, "area"] <- i
        i <- 2
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 3
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 4
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 5
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 6
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 7
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 8
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 9
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 10
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 11
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 12
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 13
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 14
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 15
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        i <- 16
        j <- tmpdata$area == 0
        j1 <- c(1:length(j))
        j1 <- j1[j == T]
        if(length(j1) > 0) {
                ind <- geoinside(tmpdata[j1,  ], reg = geo::reg.bc[[i]], option = 0,
                        robust = F)
                if(length(ind) > 0)
                        tmpdata[j1[ind], "area"] <- i
        }
        data$area <- tmpdata$area
        return(data)
}