bi 0.1324, example: outCon = Con(inRas1 < 45,1, Con((inRas1 >= 45) & (inRas1 < 47),2, Con((inRas1 >= 47) & (inRas1 < 49),3, Con(inRas1 >= 49,4)))) -------------- ndvi_flaash2 = (Float("B5_NIR_reflectance_FLAASH.tif")-Float("B4_red_reflectance_FLAASH.tif"))/(Float("B5_NIR_reflectance_FLAASH.tif")+Float("B4_red_reflectance_FLAASH.tif")) Pv_flaash2 = ("ndvi_flaash2"-0.2)/(0.5-0.2) emissivity_sobrino_sjtsk30m = Con("ndvi_flaash2"<0.2,0.973,Con(("ndvi_flaash2">=0.2) & ("ndvi_flaash2"<0.5),"Pv_flaash2"*0.004 + 0.986, Con("ndvi_flaash2">=0.5,0.99))) BTs_sobrino=("L8_B10_toa_radiance_sjtsk30m"- 1.82 - 0.78*(1-"emissivity_sobrino_sjtsk30m")*3)/(0.78*"emissivity_sobrino_sjtsk30m") L8_20160630_LST_sobrino = 1321.08 /Ln(774.89/"BTs_sobrino"+1)