tests/testthat/test-round-trip.R

In AustralianAntarcticDivision/rbgm: Tools for 'Box Geometry Model' (BGM) Files and Topology for the Atlantis Ecosystem Model

library(testthat)
context("round-trip of BGM links")

#devtools::install_github("mdsumner/rbgm")
library(rbgm)
library(dplyr)
bfile <- grep("antarctica_99", bgmfiles::bgmfiles(), value = TRUE)
#for (bfile in bgmfiles::bgmfiles()) {
  context(sprintf("reading file %s", basename(bfile)))
  bgm <- bgmfile(bfile)
  polybgm <- rbgm::boxSpatial(bgm)
  linebgm <- rbgm::faceSpatial(bgm)
  
#}
  
## run on some systems
  # library(bgmfiles)
  # library(rbgm)
  # bgm <- bgmfile(grep("NGOM", bgmfiles::bgmfiles(), value = TRUE))
  # 
  # boxes <- boxSpatial(bgm)
  # 
  # boxLL <- spTransform(boxes, "+proj=longlat +ellps=WGS84")
  # 
  # writeOGR(boxLL, ".", "NGOM_LL", "ESRI Shapefile", overwrite_layer = TRUE)
  # 
  # ## re-create our original bgm? 
  # ##system('java -jar bgmeriser-stripped.jar  -from "+proj=longlat +ellps=WGS84" -to  "+proj=utm +a=6378137.0 +es=0.006694380022900787 +lon_0=-81d00 +lat_0=0d00 +x_0=500000.0 +y_0=0.0 +k=0.9996 +zone=17" NGOM_LL.shp')
  #