import geomorph_routines_module
class LicenseError(Exception):
pass
try:
#Check for spatial analyst license
if arcpy.CheckExtension("Spatial") == "Available":
arcpy.CheckOutExtension("Spatial")
else:
raise LicenseError
#
inR = arcpy.GetParameterAsText(0)
r = geomorph_routines_module.checkExt(inR)
rOne = arcpy.GetParameterAsText(0)
rasterOne = geomorph_routines_module.checkExt(rOne)
rasterOneDesc = arcpy.Describe(rasterOne)
rTwo = arcpy.GetParameterAsText(1)
rasterTwo = geomorph_routines_module.checkExt(rTwo)
rasterTwoDesc = arcpy.Describe(rasterTwo)
# Set overwrite option
env.overwriteOutput = True
analysisWindow = arcpy.GetParameterAsText(3)
#Set message about running
arcpy.AddMessage("Calculating Statistics......")
import geomorph_routines_module
class LicenseError(Exception):
pass
try:
#Check for spatial analyst license
if arcpy.CheckExtension("Spatial") == "Available":
arcpy.CheckOutExtension("Spatial")
else:
raise LicenseError
inR = arcpy.GetParameterAsText(0)
r = geomorph_routines_module.checkExt(inR)
inZUnits = arcpy.GetParameterAsText(1)
zFactor = geomorph_routines_module.getZFactor(r, inZUnits)
# Set overwrite option
env.overwriteOutput = True
# Getting info on raster
dscRaster = arcpy.Describe(r)
inSpRef = dscRaster.spatialReference
spRefType = inSpRef.type
cellSize = dscRaster.meanCellHeight
if spRefType == "Geographic":
cellSize = cellSize / zFactor
slope = Slope(r, "DEGREE", zFactor)
aspect = Aspect(r)
import geomorph_routines_module
class LicenseError(Exception):
pass
try:
#Check for spatial analyst license
if arcpy.CheckExtension("Spatial") == "Available":
arcpy.CheckOutExtension("Spatial")
else:
raise LicenseError
r = arcpy.GetParameterAsText(0)
demRaster = geomorph_routines_module.checkExt(r)
dscRaster = arcpy.Describe(demRaster)
zUnits = arcpy.GetParameterAsText(1)
zFactor = geomorph_routines_module.getZFactor(demRaster, zUnits)
t = arcpy.GetParameterAsText(2)
s = arcpy.GetParameterAsText(3)
# Set overwrite option
env.overwriteOutput = True
#Need to have appropriate extent, will not work in ArcMap if spRefs differ
#Set message about running
arcpy.AddMessage("Running Slope Impedance......")
tmp1 = Slope(demRaster, "PERCENT_RISE", zFactor)
