import arcpy, os arcpy.CheckOutExtension("Spatial") arcpy.env.workspace = r"Z:\BP\TST" Image16 = r"Z:\BP\TST\16.tif" Image17 = r"Z:\BP\TST\17.tif" Image16_red = arcpy.sa.Raster(Image16+"\\Band_1") Image16_green = arcpy.sa.Raster(Image16+"\\Band_2") Image16_blue = arcpy.sa.Raster(Image16+"\\Band_3") Image17_red = arcpy.sa.Raster(Image17+"\\Band_1") Image17_green = arcpy.sa.Raster(Image17+"\\Band_2") Image17_blue = arcpy.sa.Raster(Image17+"\\Band_3") if Image17_red == "205": if Image17_green == "255": if Image17_blue == "194": output_obj = (Image17_red = Image16_red) & (Image17_green = Image16_green) & (Image17_blue = Image16_blue) output_obj.save(r"Z:\BP\TST\output.tif")