spgrass6: Interface between GRASS geographical information system and R

Description Details Author(s) Examples

Description

Interpreted interface between GRASS geographical information system, versions 6 and 7, and R, based on starting R from within the GRASS environment, or on running R stand-alone and creating a throw-away GRASS environment from within R. The interface uses classes defined in the sp package to hold spatial data.

Details

Index:

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
readRAST                read GRASS raster files
writeRAST               write GRASS raster files
readVECT                read GRASS vector object files
writeVECT               write GRASS vector object files
gmeta                   read GRASS metadata from the current LOCATION
getLocationProj         return a PROJ.4 string of projection information
gmeta2grd               create a GridTopology object from the GRASS region
vInfo                   return vector geometry information
vColumns                return vector database columns information
vDataCount              return count of vector database rows
vect2neigh              return area neighbours with shared boundary length

Further information may be found in the document doc/spgrass_0.3.pdf in the directory returned by system.file("", package="spgrass6").

Author(s)

Roger Bivand

Maintainer: Roger Bivand <Roger.Bivand@nhh.no>

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
if (nchar(Sys.getenv("GISRC")) > 0) {
  require(rgdal)
  soilsph <- readRAST("soils.ph", ignore.stderr=TRUE, plugin=FALSE)
  summary(soilsph)
  grd <- gmeta2grd(ignore.stderr=TRUE)
  grd
  set.seed(1)
  pts <- spsample(soilsph, 200, "random")
  smple <- SpatialPointsDataFrame(pts, data=over(pts, soilsph))
  summary(smple)
  writeVECT(smple, "sp_dem", v.in.ogr_flags="overwrite", ignore.stderr=TRUE)
  bugsDF <- readVECT("bugsites", ignore.stderr=TRUE, mapset="PERMANENT")
  summary(bugsDF)
  vInfo("streams", ignore.stderr=TRUE)
  vColumns("streams", ignore.stderr=TRUE)
  vDataCount("streams", ignore.stderr=TRUE)
  streams <- readVECT("streams", type="line,boundary",
    remove.duplicates=FALSE, ignore.stderr=TRUE, plugin=FALSE)
  summary(streams)
}

Example output

Loading required package: sp
Loading required package: XML
GRASS GIS interface loaded with GRASS version: (GRASS not running)

spgrass6 documentation built on May 2, 2019, 9:15 a.m.

Related to spgrass6 in spgrass6...