This code that I've written runs fine up until Extract by Attributes, then hangs (continues running, with no output). What did I do wrong?
import arcpy
arcpy.CheckOutExtension ("Spatial")
from arcpy import env
from arcpy.sa import *
#Allows overwriting
arcpy.env.overwriteOutput = True
##Set workspace; necessary for arcpy scripts
##note that following commands only work for the defined workspace: i.e., the ExtractbyAttributes command automatically convert all rasters in the geodatabase specified
arcpy.env.workspace = r"S:\FIRE\year1992to2000"
####iteratively defines each subfolder as a workspace; for each raster in each subfolder, converts to specified geodatabase
workspaces = arcpy.ListWorkspaces()
for workspace in workspaces:
arcpy.env.workspace = workspace
for raster in arcpy.ListRasters():
arcpy.RasterToGeodatabase_conversion (raster, r"S:\FIRE\t92to00.gdb")
#resets workspace to geodatabase mentioned above
#the "r" before the pathname tells python that this is an absolute pathname. Hard to explain, but some combinations of slashes and characters are read as commands (for example, /g means the end of a line) by arcgis and will cause it to parse the pathname incorrectly. The "r" prevents this.
arcpy.env.workspace = r"S:\FIRE\t92to00.gdb"
#adds table for final append operation
#table created manually--you could create one within the program if you wanted, using the CreateTable and AddField commands
outtable = r"S:\FIRE\73to80.gdb\y92to00"
#for each raster in geodatabase, extract by attributes to include only pixels with the correct values, in this case, 6 and 2
#rename resulting rasters; "af" refers to conversions from forest and ag only
for raster in arcpy.ListRasters():
outpoint = raster + "af"
extract = arcpy.sa.ExtractByAttributes (raster, " VALUE = 6 OR VALUE = 8")
extract.save (outpoint)
#"for" loop calls only rasters with "af" in the name
#RastertoPoint_conversion converts rasters to point files to enable use of Mean Center
#rename resulting shapefiles; keep original file information
for raster in arcpy.ListRasters ("af", "*"):
outpoint= raster[:8] + "p"
arcpy.RastertoPoint_Conversion (raster, outpoint)
#"for" loop calls only point feature classes in the workspace
#MeanCenter_stats creates a new point shapefile with a single point, the mean center of the points extracted in the previous step
#rename the resulting shapefiles
for fc in arcpy.ListFeatureClasses ("p", "Point"):
outmean= fc[:8] + "m"
arcpy.MeanCenter_stats (fc, outmean)
uid = fc[1:8]
arcpy.AddField_management(outmean, "UID", "TEXT")
arcpy.CalculateField_management(outmean, "UID", '"'+uid+'"')
#"for" loop calls only point feature classes with the phrase "mean" in their filename
#Near_analysis adds new columns to the mean center point file
#"infeat" is a single shapefile containing the borders of public lands in California
for fc in arcpy.ListFeatureClasses ("m", "Point"):
infeat= "S:\FIRE\ben_data\NPSFSBLM"
arcpy.Near_analysis (fc, infeat)
arcpy.MakeTableView (fc, "tv")
arcpy.AppendManagement ("tv", outtable, "NO_TEST", "", "")
