spAndBndConversion: Convert sp and bnd format into each other
In datacamp/R2BayesX: Estimate Structured Additive Regression Models with 'BayesX'

Description Usage Arguments Value Author(s) See Also Examples

Convert geographical information objects of class "SpatialPolygons" (or specializations) from R-package sp to objects of class "bnd" from R-package R2BayesX and vice versa.

1 2	sp2bnd(spObject, regionNames, height2width, epsilon) bnd2sp(bndObject)

`spObject`	object of class `"SpatialPolygons"` (or specializations).
`regionNames`	character vector of region names (parallel to the Polygons list in spObject), defaults to the IDs.
`height2width`	ratio of total height to width, defaults to the bounding box values.
`epsilon`	how much can two polygons differ (in maximum squared Euclidean distance) and still match each other?, defaults to machine precision.
`bndObject`	object of class `"bnd"`.

Equivalent object in the other format.

Daniel Sabanes Bove.

nb2gra, gra2nb for conversion between the neighborhood structure formats and read.bnd, write.bnd for the interface to the R2BayesX files.

## bnd to sp:
file <- file.path(find.package("R2BayesX"), "examples", "Germany.bnd")
germany <- read.bnd(file)
spGermany <- bnd2sp(germany)

## plot the result together with the neighborhood graph
if(require("sp") & require("spdep")) {
  plot(spGermany)
  nbGermany <- poly2nb(spGermany)
  plot(nbGermany, coords = coordinates(spGermany), add = TRUE)

  ## example with one region inside another
  spExample <- spGermany[c("7231", "7235"), ]
  plot(spExample)
  plot(poly2nb(spExample), coords = coordinates(spExample), add = TRUE)

  ## now back from sp to bnd:
  bndGermany <- sp2bnd(spGermany)
  plotmap(bndGermany)

  ## compare names and number of polygons
  stopifnot(
    identical(names(bndGermany), names(germany)),
    identical(length(bndGermany), length(germany))
  )
}