import arcpy arcpy.CheckOutExtension("Spatial") input_raster = "D:\\temp\\test.img" band1_obj = arcpy.sa.Raster(input_raster+"\\Layer_1") band2_obj = arcpy.sa.Raster(input_raster+"\\Layer_2") band3_obj = arcpy.sa.Raster(input_raster+"\\Layer_3") output_obj = (band1_obj == 112) & (band2_obj == 115) & (band3_obj == 5) output_obj.save(r"D:\\temp\\output.img")