Create a boundary polygon by tracking around coordinates stored in a RasterStack.
romsboundary(cds) databoundary(x, mask = NULL, ...)
cds | two-layer Raster |
---|---|
x | a thing |
mask | logical |
... | dots |
The first layer in the stack is treated as the X coordinate, second as Y.
## ok in this one #library(rgdal) # proj4string(ice_grid_boundary) <- CRS("+init=epsg:4326") # pweird <- "+proj=laea +lon_0=147 +lat_0=-42 +ellps=WGS84" # laea_world <- spTransform(antarctica, pweird) # plot(extent(laea_world) + 8e6, type = "n", asp = 1) # plot(laea_world, add = TRUE) # plot(spTransform(ice_grid_boundary, pweird), add = TRUE, border = "darkgrey")