R/distance.R

In zf-ibrahim/myqdmi: Malaysia Quarry Dust Monitoring Artificial Intelligent

#Distance



#can use this package
#library(geosphere)
#distm (c(lon1, lat1), c(lon2, lat2), fun = distHaversine)

##' Distance
##' 
##' A function to calculate distance between receptor point to source location.
##'
##' @param sourceInput A data frame containing \code{sourceActivity} is a name for the pointsource, \code{type} (the type of material processing),
##' \code{x} for Easting, \code{y} for Northing and \code{z} for elevation height of the location.
##' 
##' @param receptorInput A data frame containing \code{receptor} the name/label for receptor point,
##' \code{x} for Easting, \code{y} for Northing and \code{z} for elevation height of the location.
##' 
##' @param sourceActivity 
##'
##' @return distance between receptor location to source location
##' @export
##' 
##' @author Zul Fadhli & Dr Izhar Abadi
##'
##' @examples
##' #demo
##' distance(sourceInput, receptorInput, sourceActivity)
distance <- function(sourceInput, receptorInput, sourceActivity ){

  sourceInput <- sourceInput
  receptorInput <- receptorInput
  sourceActivity <- sourceActivity


  distance.storage <- rep(NA,nrow(receptorInput))

  for (i in 1:nrow(receptorInput)){
 distance.storage[i] <- sqrt(
      ((sourceInput[sourceInput$sourceActivity == sourceActivity,"x"] - receptorInput[i,"x"])^2) +
      ((sourceInput[sourceInput$sourceActivity == sourceActivity,"y"] - receptorInput[i,"y"])^2) +
      ((sourceInput[sourceInput$sourceActivity == sourceActivity,"z"] - receptorInput[i,"z"])^2)

     )

  }
  
  #x = SQRT((PcoorX-RcoorX)^2+(PcoorY-RcoorY)^2+(PcoorZ-RcoorZ)^2)  
  distance.storage
}