R/sourceModel.R

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

##' Application of Source model for dust deposited rate estimation to all receptor and source locations
##' 
##'
##' \code{sourceModel} Estimation of Dust Deposition Model (EDDM) established to estimate dust deposited rate at receptor
##' 
##' \code{sourceModel} require the source and receptor location in Easting and Northing format.
##' The \code{sourceActivity} for this model is the type of material processing activity in the quarry or
##' if the \code{sourceInput} is the quarries location, the \code{sourceActivity} is the name of the quarries.
##' 
##' @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 windInput A data frame containing \code{ws} for windspeed (m/s) and \code{wd} for wind direction
##' 
##' @param sourceActivity The pointsource name such as \sQuote{primaryCrusher} or \sQuote{point A} which is
##' according to the \code{souceInput} \code{sourceActivity}.
##' 
##' 
##' @export
##' 
##' @return The results will be the estimation of dust deposited rate at all the receptors from all point source.
##' The unit measurement for dust deposited rate is ug/m2/month.
##' 
##' @author Zul Fadhli & Dr. Izhar Abadi
##' 
##' @examples 
##' 
##' #demo
##' sourceModel((quarryInput = quarryInput, sourceInput = sourceInput, receptorInput = receptorInput, windInput = windInput, particleSize = "tsp")

sourceModel <- function(quarryInput = quarryInput,
                        sourceInput = sourceInput,
                        receptorInput = receptorInput,
                        windInput = windInput,
                        particleSize = "tsp"
                        ){

  size <- particleSize

quarryInput <- quarryInput
sourceInput <- sourceInput
receptorInput <- receptorInput
windInput <- windInput



a <- primaryCrusher(quarryInput, particleSize = size)

b <- secondaryCrusher(quarryInput, particleSize = size)

c <- tertiaryCrusher(quarryInput, particleSize = size)

d <- stockpile(quarryInput, particleSize = size)

e <- screening(quarryInput, particleSize = size, QuarryProduction = quarryInput$quarryProduction,
               totalScreen = quarryInput$nScreen)

f <- unpavedRoad(quarryInput, particleSize = size)

g <- pavedRoad(quarryInput, particleSize = size)

h <- blasting(quarryInput, particleSize = size, areaBlasted = quarryInput$areaBlasted,
              blastFreq = quarryInput$blastingFreq, blastHole.depth = quarryInput$blasthole, moistureContent = quarryInput$moisture)


ii <- drilling(quarryInput, particleSize = size,
               drillHole = quarryInput$drillhole)


results <- rep(NA,nrow(receptorInput))

for(i in 1:nrow(sourceInput)){

  results[i] <- list(


    source.model(sourceInput, receptorInput, sourceActivity =

                  sourceInput[sourceInput$sourceActivity == sourceInput[i, "sourceActivity"],"sourceActivity"]


                , windInput, vg = 5, AnemometerHeight = 16,
                dustGenerated = if(sourceInput[i, "type"]== "primary crusher"){
                  a
                }else if(sourceInput[i, "type"]== "secondary crusher"){
                  b
                }else if(sourceInput[i, "type"]== "tertiary crusher"){
                  c
                }else if(sourceInput[i, "type"]== "stockpile"){
                  d
                }else if(sourceInput[i, "type"]== "screening"){
                  e
                }else if(sourceInput[i, "type"]== "unpaved road"){
                  f
                }else if(sourceInput[i, "type"]== "paved road"){
                  g
                }else if(sourceInput[i, "type"]== "blasting"){
                  h
                }else if(sourceInput[i, "type"]== "drilling"){
                  ii
                }


    )


  )
}


myresults <- data.frame(matrix(unlist(results), nrow=NROW(results), ncol = nrow(receptorInput), byrow=T),stringsAsFactors=F)

receptorName <- receptorInput$receptor

myresults <- rbind(receptorName, myresults)

sourceSite <- c("sites" ,sourceInput$sourceActivity)


myresults <- cbind(sourceSite, myresults)

names(myresults) <- as.matrix(myresults[1, ])
myresults <- myresults[-1, ]
myresults[] <- lapply(myresults, function(x) type.convert(as.character(x)))

myresults



}