arcpy.MakeFeatureLayer_management(infc_gdb, "fc_lyr") path, fc = os.path.split(infc_gdb) clippath, fc_clip = os.path.split(clipfc) print "Starting Clip for {0} with ".format(fc, fc_clip) outfc = out_working + os.sep + fc + "_LandClip" if arcpy.Exists(outfc): print "{0} Already clipped".format(outfc) else: arcpy.env.workspace = out_working arcpy.env.overwriteOutput = True arcpy.Clip_analysis("fc_lyr", clipfc, outfc) else: if MutliGDB: for v in infc_gdb: ingdb = v for fc in functions.fcs_in_workspace(ingdb): start_loop = datetime.datetime.now() path, gdb = os.path.split(ingdb) inloopfc = ingdb + os.sep + fc arcpy.Delete_management("fc_lyr") arcpy.MakeFeatureLayer_management(inloopfc, "fc_lyr") clippath, fc_clip = os.path.split(clipfc) outgdbpath = out_working + os.sep + gdb CreateGDB(out_working, gdb, outgdbpath) print "Starting Clip for {0} with ".format(fc, fc_clip) outfc = outgdbpath + os.sep + fc + "_LandClip" if arcpy.Exists(outfc): print "{0} Already clipped".format(outfc) else: arcpy.env.workspace = out_working arcpy.env.overwriteOutput = True
for line in inputFile: line = line.split(',') finalGDB = str(line[0]) finalGDB = finalGDB.strip('\n') use = str(line[1]) use = use.strip('\n') outtable = finalGDB + os.sep + "Topo_results_" + use if arcpy.Exists(outtable): intable = outtable else: intable = topolist print finalGDB for fc in functions.fcs_in_workspace(finalGDB): #print fc current_update = {} CHchecklist = fc.split("_") CHcheck = CHchecklist[0] if CHcheck != 'Topo': continue elif str(CHchecklist[1]) == 'T': continue elif str(CHchecklist[1]) == 'R': entlist = fc.split("_") ent = str(entlist[2]) filename = str(entlist[1]) + "_" + str(entlist[2]) #print filename else:
for k, v in globals().items(): if k.startswith('huc'): num = k.replace("huc", "") HUCdict[num] = v start = time.time() start_times(start) # Create Workspaces create_directory(File_dir) create_directory(LYR_dir) create_directory(DBF_dir) # Execute loop that will check each input comp fc to see if a species in the species list is present, and if so it will # extract all streams with in the range and save to a dbf for fc in functions.fcs_in_workspace(MasterSpeFC): check_prj = fc.split('_') try: if 'AlbersUSA' == check_prj[4]: continue except: pass print fc input_fc = MasterSpeFC + os.sep + fc stream_crosswalk(input_fc, HUC2_lwr48, LYR_dir, DBF_dir, HUCdict, noNHD, csvpath) print noNHD done = time.time() end_times(done, start)
createdirectory(DBF_dir) createdirectory(CSV_dir) counter = 0 if MultiFolder: FolderList = os.listdir(inFolder) print FolderList for value in FolderList: FinalGDB = inFolder + os.sep + value + os.sep + 'Final.gdb' # print FinalGDB outfolder = DBF_dir + os.sep + value createdirectory(outfolder) if fc: for fc in functions.fcs_in_workspace(FinalGDB): print fc counter += 1 table = "{0}.dbf".format(fc) outtable = outfolder + os.sep + table tableview = "tbl_view_" + str(counter) arcpy.Delete_management("fc_lyr") if arcpy.Exists(outtable): continue else: arcpy.MakeFeatureLayer_management(fc, "fc_lyr") arcpy.MakeTableView_management("fc_lyr", tableview) arcpy.TableToTable_conversion(tableview, outfolder, table) else: for raster in functions.rasters_in_workspace(FinalGDB):
def CreateDirectory(DBF_dir): if not os.path.exists(DBF_dir): os.mkdir(DBF_dir) def CreateGDB(OutFolder, OutName, outpath): if not arcpy.Exists(outpath): arcpy.CreateFileGDB_management(OutFolder, OutName, "CURRENT") start_script = datetime.datetime.now() print "Started: {0}".format(start_script) for fc in functions.fcs_in_workspace(inGDB): compFC = inGDB + os.sep + str(fc) group = fc.split('_') group = str(group[2]) fish = fishnet outpath_working = outpath + os.sep + "Working.gdb" outpath_final = outpath + "GAP_" + group + '.gdb' resultfolder, workinggdb = os.path.split(outpath_working) finalgdb = os.path.basename(outpath_final) lyr = resultfolder + os.sep + 'lyr' CreateDirectory(resultfolder) CreateDirectory(lyr) CreateGDB(resultfolder, workinggdb, (resultfolder + os.sep + workinggdb)) CreateGDB(resultfolder, finalgdb, (resultfolder + os.sep + finalgdb))
# import system modules import os import csv import time import datetime import functions import arcpy # Title : Checks projection for all fc infolder to make sure they are in the STD geo projection before moving forward Infolder = 'J:\Workspace\ESA_Species\CriticalHabitat\NAD_Final' # ##################################################################################################################### ############################################################################################################## start_script = datetime.datetime.now() print "Script started at {0}".format(start_script) # Prints Elapse clock prj_list = [] for fc in functions.fcs_in_workspace(Infolder): # set local variables ORGdsc = arcpy.Describe(fc) ORGsr = ORGdsc.spatialReference ORGprj = ORGsr.name.lower() if ORGprj not in prj_list: prj_list.append(ORGprj) print "{0} {1}".format(fc, ORGsr.name) else: pass
symbologyLayer = "J:\Workspace\ESA_Species\Step3\ZonalHis_GAP\GAP.lyr" arcpy.ApplySymbologyFromLayer_management( "rdlayer", symbologyLayer) functions.ZonalHist(lyr, zoneField, "rdlayer", dbf, outpath_final, scratchpath) del lyr print "Loop completed in: {0}".format(datetime.datetime.now() - start_loop) start_script = datetime.datetime.now() print "Script started at: {0}".format(start_script) if inlocation[-3:] == 'gdb': for fc in functions.fcs_in_workspace(inlocation): infc = inlocation + os.sep + fc outpath_final = results for v in pathlist: rasterLocation = v for raster in rasters_in_workspace(rasterLocation): ZonalHist(rasterLocation, raster, outpath_final, infc, FullFile) else: infc = inlocation outpath_final = results for v in pathlist: rasterLocation = v for raster in rasters_in_workspace(rasterLocation): ZonalHist(rasterLocation, raster, outpath_final, infc, FullFile)