topic ExtractByMask doesn&acute;t work (AG 10 SP5) in Python Questions https://community.esri.com/t5/python-questions/extractbymask-doesn-acute-t-work-ag-10-sp5/m-p/239398#M18653 <HTML><HEAD></HEAD><BODY><SPAN>Hi i have one problem which drive me crazy. i used ExtractByMask which worked fine (SP4) but now i dont know why it doesn´t. any idea? thanks a lot&nbsp; (see code please&nbsp; line: pThiessen = ExtractByMask(pThiessenTemp, pMask)&nbsp;&nbsp; err: 000582 Error occured during execution)</SPAN><BR /><BR /><PRE class="plain" name="code"> import arcpy from arcpy.sa import * def main(): &nbsp;&nbsp;&nbsp; try: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.CheckOutExtension("Spatial") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.overwriteOutput = True &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inPoints = arcpy.GetParameterAsText(0) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inField = arcpy.GetParameterAsText(1) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inType = arcpy.GetParameterAsText(2) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inMask = arcpy.GetParameterAsText(3) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inWorkSpace = arcpy.GetParameterAsText(4) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pMask = Raster(inMask) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.snapRaster = inMask &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; #Nastavenie union extentu body a povodie &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; desc = arcpy.Describe(inPoints) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmin = min(desc.extent.XMin,pMask.extent.XMin) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ymin = min(desc.extent.YMin,pMask.extent.YMin) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmax = max(desc.extent.XMax,pMask.extent.XMax) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ymax = max(desc.extent.YMax,pMask.extent.YMax) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.extent = arcpy.Extent(xmin, ymin, xmax, ymax) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; #Koniec nastavenia extentu &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.mask = inMask &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.cellSize = inMask &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pCellsize = (pMask.meanCellHeight + pMask.meanCellWidth) / 2 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.workspace = inWorkSpace &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (inWorkSpace.find(".gdb") == -1) | (inWorkSpace.find(".gdb") == -1): &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; db = ".shp" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; else: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; db = "" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if inField == "FID": &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inField = "OBJECTID" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.CreateThiessenPolygons_analysis(inPoints, "ThTemp1", "ALL") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.FeatureToRaster_conversion("ThTemp1" + db, inField, "thtemp2", pCellsize) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.Delete_management("ThTemp1"+ db) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessenTemp = Raster(inWorkSpace + "\\thtemp2") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.ClearEnvironment("extent") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.extent = pMask.extent &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessen = ExtractByMask(pThiessenTemp, pMask) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.Delete_management(inWorkSpace +"\\thtemp2") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if inType == "precipitation": &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessen.save(inWorkSpace + "\\thiessen_p") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Name = "thiessen_p" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; elif inType == "evapotranspiration": &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessen.save(inWorkSpace + "\\thiessen_e") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Name = "thiessen_e" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; else: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessen.save(inWorkSpace + "\\thiessen_t") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Name = "thiessen_t" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; del pMask &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; del pThiessenTemp &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; del pThiessen &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.snapRaster = None &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.extent = None &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.mask = None &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.cellSize = None &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; AddRasterLayerToMap(inWorkSpace, Name) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.AddMessage("Script Run Succesfully") &nbsp;&nbsp;&nbsp; except Exception as e: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.AddError(e.message) </PRE><DIV style="display:none;"> </DIV></BODY></HTML> Wed, 24 Oct 2012 07:13:40 GMT martinmaretta 2012-10-24T07:13:40Z ExtractByMask doesn´t work (AG 10 SP5) https://community.esri.com/t5/python-questions/extractbymask-doesn-acute-t-work-ag-10-sp5/m-p/239398#M18653 <HTML><HEAD></HEAD><BODY><SPAN>Hi i have one problem which drive me crazy. i used ExtractByMask which worked fine (SP4) but now i dont know why it doesn´t. any idea? thanks a lot&nbsp; (see code please&nbsp; line: pThiessen = ExtractByMask(pThiessenTemp, pMask)&nbsp;&nbsp; err: 000582 Error occured during execution)</SPAN><BR /><BR /><PRE class="plain" name="code"> import arcpy from arcpy.sa import * def main(): &nbsp;&nbsp;&nbsp; try: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.CheckOutExtension("Spatial") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.overwriteOutput = True &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inPoints = arcpy.GetParameterAsText(0) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inField = arcpy.GetParameterAsText(1) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inType = arcpy.GetParameterAsText(2) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inMask = arcpy.GetParameterAsText(3) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inWorkSpace = arcpy.GetParameterAsText(4) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pMask = Raster(inMask) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.snapRaster = inMask &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; #Nastavenie union extentu body a povodie &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; desc = arcpy.Describe(inPoints) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmin = min(desc.extent.XMin,pMask.extent.XMin) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ymin = min(desc.extent.YMin,pMask.extent.YMin) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; xmax = max(desc.extent.XMax,pMask.extent.XMax) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ymax = max(desc.extent.YMax,pMask.extent.YMax) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.extent = arcpy.Extent(xmin, ymin, xmax, ymax) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; #Koniec nastavenia extentu &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.mask = inMask &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.cellSize = inMask &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pCellsize = (pMask.meanCellHeight + pMask.meanCellWidth) / 2 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.workspace = inWorkSpace &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (inWorkSpace.find(".gdb") == -1) | (inWorkSpace.find(".gdb") == -1): &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; db = ".shp" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; else: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; db = "" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if inField == "FID": &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; inField = "OBJECTID" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.CreateThiessenPolygons_analysis(inPoints, "ThTemp1", "ALL") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.FeatureToRaster_conversion("ThTemp1" + db, inField, "thtemp2", pCellsize) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.Delete_management("ThTemp1"+ db) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessenTemp = Raster(inWorkSpace + "\\thtemp2") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.ClearEnvironment("extent") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.extent = pMask.extent &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessen = ExtractByMask(pThiessenTemp, pMask) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.Delete_management(inWorkSpace +"\\thtemp2") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if inType == "precipitation": &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessen.save(inWorkSpace + "\\thiessen_p") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Name = "thiessen_p" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; elif inType == "evapotranspiration": &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessen.save(inWorkSpace + "\\thiessen_e") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Name = "thiessen_e" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; else: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pThiessen.save(inWorkSpace + "\\thiessen_t") &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Name = "thiessen_t" &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; del pMask &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; del pThiessenTemp &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; del pThiessen &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.snapRaster = None &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.extent = None &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.mask = None &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.env.cellSize = None &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; AddRasterLayerToMap(inWorkSpace, Name) &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.AddMessage("Script Run Succesfully") &nbsp;&nbsp;&nbsp; except Exception as e: &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy.AddError(e.message) </PRE><DIV style="display:none;"> </DIV></BODY></HTML> Wed, 24 Oct 2012 07:13:40 GMT https://community.esri.com/t5/python-questions/extractbymask-doesn-acute-t-work-ag-10-sp5/m-p/239398#M18653 martinmaretta 2012-10-24T07:13:40Z