| ArcGIS Geoprocessing Tools {RPyGeo} | R Documentation |
Wrappers for a small selection of ArcGIS geoprocessing
functions based on the rpygeo.geoprocessor.
rpygeo.EucDistance.sa(in.data, out.raster, maxdist = NULL,
cellsize = NULL, out.direction.raster = NULL,
env = rpygeo.env, ...)
rpygeo.Hillshade.sa(in.raster, out.raster, azimuth = 315,
altitude = 45, model.shadows = c("NO_SHADOWS", "SHADOWS"),
z.factor = 1, ...)
rpygeo.Slope.sa(in.raster, out.raster,
unit = c("DEGREE", "PERCENT_RISE"), z.factor = 1, ...)
rpygeo.Aspect.sa(in.raster, out.raster, ...)
rpygeo.Delete.management(in.data, data.type = NULL, ...)
in.raster, in.data, out.raster |
Names of ArcGIS raster or vector
datasets or feature classes in a geodatabase (relative to the
current workspace defined in a rpygeo.env environment).
Shapefiles must include the extension “.shp”. |
env |
A list defining an RPyGeo working environment as built
by rpygeo.build.env. |
maxdist, cellsize, out.direction.raster, |
|
azimuth, altitude, model.shadows, z.factor, |
|
unit, data.type |
Arguments to be passed to the Python geoprocessing function. See ArcGIS help files for information on the usage of scripting commands and their arguments. |
... |
Additional arguments to be passed to rpygeo.geoprocessor. |
These functions simply try to replicate the behaviour of the
ArcGIS/Python geoprocessing functions of the same name.
See rpygeo.geoprocessor for details on what happens behind
the scenes.
The function return NULL if no error occurred, otherwise
a character vector containing the error message.
Alexander Brenning
rpygeo.geoprocessor, rpygeo.build.env
# Allow ArcGIS to overwrite existing datasets:
## Not run: rpygeo.env$overwriteoutput = 1
# Calculate the slope of a DEM raster dataset
# in the current ArcGIS workspace:
## Not run: rpygeo.geoprocessor("Slope_sa",c("dem","slope"))
# Same:
## Not run: rpygeo.geoprocessor("Slope_sa('dem','slope')")
# Same, using the more convenient wrapper:
## Not run: rpygeo.Slope.sa("dem","slope")
# Three at a time or separately:
## Not run: date()
## Not run:
rpygeo.geoprocessor("Slope_sa('dem','slope')",
"Aspect_sa('dem','aspect')", "Hillshade_sa('dem','hshd')")
## End(Not run)
## Not run: date() # ~20 sec on my computer
## Not run: rpygeo.Slope.sa("dem","slope")
## Not run: rpygeo.Aspect.sa("dem","aspect")
## Not run: rpygeo.Hillshade.sa("dem","hshd")
## Not run: date() # ~50 sec
## Not run: rpygeo.Delete.management("slope")
## Not run: rpygeo.Delete.management("aspect")
## Not run: rpygeo.Delete.management("hshd")
# Calculate the Euclidian distance from railway lines
# up to a max. distance of 1000 map units:
## Not run:
rpygeo.geoprocessor("EucDistance_sa",
args=list("rail.shp","raildist",1000))
## End(Not run)
# Same:
## Not run: rpygeo.EucDistance.sa("rail.shp","raildist",maxdist=1000)
# Use MapAlgebra to calculate a distance-decay function:
## Not run:
rpygeo.geoprocessor("SingleOutputMapAlgebra_sa",
args=c("exp( raildist / -100 )","distdecay"))
## End(Not run)
# Or why not in just one step if you like MapAlgebra:
## Not run:
rpygeo.geoprocessor( "SingleOutputMapAlgebra_sa",
args=c("exp( EucDistance( rail.shp, \#, \#, 1000 ) / -100 )","distdecay") )
## End(Not run)