Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import arcpy
- import os
- import glob
- from arcpy import env
- from arcpy.sa import *
- arcpy.CheckOutExtension("Spatial")
- env.workspace = r"File_pathway"
- arcpy.env.extent = r"File_pathway2"
- indef_file = "File_pathway3"
- inRasters = glob.glob(r"File_pathway4\\*.tif")
- inMaskData = "Quads"
- extraction_area = "OUTSIDE"
- analysis_extent = "Quads"
- for inRaster in inRasters:
- outExtractByMask = ExtractByMask(inRaster, inMaskData)
- classifiedraster = ClassifyRaster(outExtractByMask, indef_file)
- classifiedraster.save = (r"File_pathway")
- newname = inRaster + "CF"
- arcpy.management.Rename(classifiedraster,newname)
- outCellstats = CellStatistics([classifiedraster], "SUM", "DATA")
- outCellstats.save=(r"file_pathway")
- newname2 = inRaster + "CS"a
- arcpy.management.Rename(outCellstats, newname2)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement