Using Python in GIS: Part 7 'The Finale!' Land Surface Temperature

Показать описание

This final video in our series Using Python in GIS looks at performing a Land Surface Temperature analysis.
This episode uses intensive Raster Calculations and even gives a tip on extracting information from a previously made NDVI layer.
Please be sure to Subscribe, share and hit that notification Bell!
Thank you for sticking with me during this long series! You all are great!

Twitter: @my_physical
Instagram: @MyPhysicalWorld

Code used in this video:

# ======================================================== Land Surface Temperature ============================================================
print("Creating the LST image.....")
# ========================== At Satellite Temp for each band = K2/Ln(K1/RAD - 1) - 273.15 =======================================================
print("Performing Cell Statistics for averaging of Band 10 and Band 11 At Satellite Temperature......")
# =========================== Emissivity = PV = SQ(NDVI - NDVImin / NDVImax - NDVImin) = 0.004(PV) + 0.986 ======================================
print("Calculating Proportion of Vegetation.......")
inRastN = "NDVI.img"
NDVImin = arcpy.GetRasterProperties_management(inRastN, "MINIMUM")
print(NDVImin)
NDVImax = arcpy.GetRasterProperties_management(inRastN, "MAXIMUM")
print(NDVImax)
PvNum = Raster(inRastN) - NDVImin
PvDenom = NDVImax - NDVImin
PvSq = Square(PvDiv)
# Pv = Square((inRastN - NDVImin) / (NDVImax - NDVImin))
print("Calculating Emissivity......")
# ======================================= Land Surface Temperature Final Calculation ==============================================================
# ----- Formula is: AtSat / 1 + (wavelength of emitted rad (11.5)) * (AtSat / 14380) * ln(LSE)
print("Calculating Final Land Surface Temperature......")
LSE = Raster("LSE.img")
LST = ((AtSat / 1) + (11.5 * (AtSat / 14380)) * Ln(LSE))
LST.save("LST.img")