data-raw/rawdata.R

In dariomasante/bnspatial: Spatial Implementation of Bayesian Networks and Mapping

## Build the .RData from raw
## Make an empty blank R session (or better restart the R session in vanilla mode)
rm(list=ls()) 

## Land use change example
library(raster)
library(bnspatial)

ConwyData <- list()

rawPath <- system.file("extdata/ConwyStatus.tif", package = "bnspatial")
ConwyData$ConwyStatus <- readAll(raster(rawPath))

rawPath <- system.file("extdata/ConwySlope.tif", package = "bnspatial")
ConwyData$ConwySlope <- readAll(raster(rawPath))

rawPath <- system.file("extdata/ConwyLU.tif", package = "bnspatial")
ConwyData$ConwyLU <- readAll(raster(rawPath))

rawPath <- system.file("extdata/LUclasses.txt", package = "bnspatial")
ConwyData$LUclasses <- importClasses(rawPath)

rawPath <- system.file("extdata/LandUseChange.net", package = "bnspatial")
ConwyData$LandUseChange <- loadNetwork(rawPath)

spDataLst <- linkMultiple(c(ConwyData$ConwyLU, ConwyData$ConwySlope, ConwyData$ConwyStatus), 
                          ConwyData$LandUseChange, ConwyData$LUclasses, verbose = FALSE)
coord <- aoi(ConwyData$ConwyLU, xy=TRUE)
ConwyData$evidence <- bulkDiscretize(spDataLst, coord)

rm(rawPath)
rm(coord)
rm(spDataLst)
rm(.Random.seed)
#save('ConwyData', file=system.file("data/ConwyData.RData", package = "bnspatial"))