{"id":1446,"date":"2026-10-25T09:18:25","date_gmt":"2026-10-25T09:18:25","guid":{"rendered":"https:\/\/helpwithmygisassignment.com\/blog\/?p=1446"},"modified":"2026-10-25T09:18:41","modified_gmt":"2026-10-25T09:18:41","slug":"urban-heat-island-mapping-landsat-lst-tutorial","status":"publish","type":"post","link":"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/","title":{"rendered":"Urban Heat Island Mapping: Calculating Land Surface Temperature (LST) from Landsat TIR Bands"},"content":{"rendered":"<div id=\"ez-toc-container\" class=\"ez-toc-v2_0_89 counter-hierarchy ez-toc-counter ez-toc-grey ez-toc-container-direction\">\n<div class=\"ez-toc-title-container\">\n<p class=\"ez-toc-title\" style=\"cursor:inherit\">Table of Contents<\/p>\n<span class=\"ez-toc-title-toggle\"><a href=\"#\" class=\"ez-toc-pull-right ez-toc-btn ez-toc-btn-xs ez-toc-btn-default ez-toc-toggle\" aria-label=\"Toggle Table of Content\"><span class=\"ez-toc-js-icon-con\"><span class=\"\"><span class=\"eztoc-hide\" style=\"display:none;\">Toggle<\/span><span class=\"ez-toc-icon-toggle-span\"><svg style=\"fill: #999;color:#999\" xmlns=\"http:\/\/www.w3.org\/2000\/svg\" class=\"list-377408\" width=\"20px\" height=\"20px\" viewBox=\"0 0 24 24\" fill=\"none\"><path d=\"M6 6H4v2h2V6zm14 0H8v2h12V6zM4 11h2v2H4v-2zm16 0H8v2h12v-2zM4 16h2v2H4v-2zm16 0H8v2h12v-2z\" fill=\"currentColor\"><\/path><\/svg><svg style=\"fill: #999;color:#999\" class=\"arrow-unsorted-368013\" xmlns=\"http:\/\/www.w3.org\/2000\/svg\" width=\"10px\" height=\"10px\" viewBox=\"0 0 24 24\" version=\"1.2\" baseProfile=\"tiny\"><path d=\"M18.2 9.3l-6.2-6.3-6.2 6.3c-.2.2-.3.4-.3.7s.1.5.3.7c.2.2.4.3.7.3h11c.3 0 .5-.1.7-.3.2-.2.3-.5.3-.7s-.1-.5-.3-.7zM5.8 14.7l6.2 6.3 6.2-6.3c.2-.2.3-.5.3-.7s-.1-.5-.3-.7c-.2-.2-.4-.3-.7-.3h-11c-.3 0-.5.1-.7.3-.2.2-.3.5-.3.7s.1.5.3.7z\"\/><\/svg><\/span><\/span><\/span><\/a><\/span><\/div>\n<nav><ul class='ez-toc-list ez-toc-list-level-1 ' ><li class='ez-toc-page-1 ez-toc-heading-level-1'><a class=\"ez-toc-link ez-toc-heading-1\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Urban_Heat_Island_Mapping_Calculating_Land_Surface_Temperature_LST_from_Landsat_TIR_Bands\" >Urban Heat Island Mapping: Calculating Land Surface Temperature (LST) from Landsat TIR Bands<\/a><ul class='ez-toc-list-level-2' ><li class='ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-2\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#1_Introduction_Mapping_Urban_Heat_Islands_with_Landsat\" >1. Introduction: Mapping Urban Heat Islands with Landsat<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-3\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#2_Data_Acquisition_and_Requirements\" >2. Data Acquisition and Requirements<\/a><ul class='ez-toc-list-level-3' ><li class='ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-4\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Landsat_89_bands_and_spatial_resolution\" >Landsat 8\/9 bands and spatial resolution<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-5\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Download_the_imagery\" >Download the imagery<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-6\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Software_requirements\" >Software requirements<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-7\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#3_Step-by-Step_Workflow_in_QGIS\" >3. Step-by-Step Workflow in QGIS<\/a><ul class='ez-toc-list-level-3' ><li class='ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-8\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Step_1_Convert_thermal_Digital_Numbers_to_TOA_spectral_radiance\" >Step 1: Convert thermal Digital Numbers to TOA spectral radiance<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-9\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Step_2_Convert_spectral_radiance_to_brightness_temperature\" >Step 2: Convert spectral radiance to brightness temperature<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-10\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Step_3_Calculate_NDVI_and_fractional_vegetation_cover\" >Step 3: Calculate NDVI and fractional vegetation cover<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-11\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Step_4_Estimate_land_surface_emissivity\" >Step 4: Estimate land surface emissivity<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-12\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Step_5_Calculate_land_surface_temperature\" >Step 5: Calculate land surface temperature<\/a><ul class='ez-toc-list-level-4' ><li class='ez-toc-heading-level-4'><a class=\"ez-toc-link ez-toc-heading-13\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Option_A_Simplified_emissivity-corrected_single-channel_estimate\" >Option A: Simplified emissivity-corrected single-channel estimate<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-4'><a class=\"ez-toc-link ez-toc-heading-14\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Option_B_Full_mono-window_algorithm\" >Option B: Full mono-window algorithm<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-4'><a class=\"ez-toc-link ez-toc-heading-15\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#What_about_the_split-window_algorithm\" >What about the split-window algorithm?<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-16\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Step_6_The_recommended_shortcut_%E2%80%94_use_USGS_Collection_2_Level-2_surface_temperature\" >Step 6: The recommended shortcut \u2014 use USGS Collection 2 Level-2 surface temperature<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-17\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#4_Applied_Analysis_Classifying_Heat_Zones\" >4. Applied Analysis: Classifying Heat Zones<\/a><ul class='ez-toc-list-level-3' ><li class='ez-toc-heading-level-3'><a class=\"ez-toc-link ez-toc-heading-18\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#Link_hotspots_to_vegetation_and_built_surfaces\" >Link hotspots to vegetation and built surfaces<\/a><\/li><\/ul><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-19\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#5_Quality_Control_and_Interpretation\" >5. Quality Control and Interpretation<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-20\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#6_SEO_Best_Practices_and_Related_Tutorials\" >6. SEO Best Practices and Related Tutorials<\/a><\/li><li class='ez-toc-page-1 ez-toc-heading-level-2'><a class=\"ez-toc-link ez-toc-heading-21\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/#7_Conclusion_and_Key_Takeaways\" >7. Conclusion and Key Takeaways<\/a><\/li><\/ul><\/li><\/ul><\/nav><\/div>\n<h1><span class=\"ez-toc-section\" id=\"Urban_Heat_Island_Mapping_Calculating_Land_Surface_Temperature_LST_from_Landsat_TIR_Bands\"><\/span>Urban Heat Island Mapping: Calculating Land Surface Temperature (LST) from Landsat TIR Bands<span class=\"ez-toc-section-end\"><\/span><\/h1>\n<h2><span class=\"ez-toc-section\" id=\"1_Introduction_Mapping_Urban_Heat_Islands_with_Landsat\"><\/span>1. Introduction: Mapping Urban Heat Islands with Landsat<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-1605\" src=\"https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_duhopnduhopnduho-1024x559.jpg\" alt=\"Urban Heat Island Mapping Landsat LST Tutorial\" width=\"1024\" height=\"559\" srcset=\"https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_duhopnduhopnduho-1024x559.jpg 1024w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_duhopnduhopnduho-300x164.jpg 300w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_duhopnduhopnduho-768x419.jpg 768w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_duhopnduhopnduho.jpg 1408w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/p>\n<p>This <b data-path-to-node=\"2,0\" data-index-in-node=\"6\">Urban Heat Island Mapping Landsat LST Tutorial<\/b> explores how built-up areas exhibit higher temperatures than their surrounding rural or less-developed environments. Urban Heat Islands (UHIs) occur when built-up areas exhibit higher temperatures than their surrounding rural or less-developed environments. Asphalt roads, concrete buildings, rooftops, and industrial surfaces absorb and store incoming solar energy, while reduced vegetation limits evaporative cooling. These conditions can create localized <a href=\"https:\/\/helpwithmygisassignment.com\/blog\/idw-vs-kriging-gis\/\">thermal hotspots<\/a> that influence energy demand, outdoor comfort, and urban environmental planning.<\/p>\n<p>Remote sensing provides an effective way to map these patterns across an entire city. Landsat 8 and Landsat 9 carry the Thermal Infrared Sensor (TIRS), which measures emitted thermal radiation in two bands. These observations can be processed to estimate <strong>Land Surface Temperature (LST)<\/strong>, revealing spatial variations between dense urban development, parks, water bodies, and vegetated outskirts.<\/p>\n<p>A critical distinction is that LST is not the same as air temperature. Air temperature describes the temperature of the surrounding atmosphere, commonly measured in a shaded, ventilated screen. LST represents the radiometric temperature of the land surface observed by the satellite. A hot roof may have a substantially higher LST than the air temperature measured nearby.<\/p>\n<p>Consequently, a Landsat-derived UHI map is primarily a map of the <em>surface urban heat island<\/em>, not a direct map of pedestrian-level air temperature or human heat stress.<\/p>\n<p>This tutorial demonstrates how to calculate LST in QGIS using Landsat thermal infrared bands, radiometric calibration, vegetation indices, and emissivity correction. It also explains when a manually implemented algorithm is appropriate and when using the USGS-preprocessed surface-temperature product is the more defensible choice.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"2_Data_Acquisition_and_Requirements\"><\/span>2. Data Acquisition and Requirements<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<h3><span class=\"ez-toc-section\" id=\"Landsat_89_bands_and_spatial_resolution\"><\/span>Landsat 8\/9 bands and spatial resolution<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Landsat 8 and Landsat 9 combine the Operational Land Imager (OLI), which records reflected solar radiation, with TIRS, which records thermal infrared radiation.<\/p>\n<table>\n<tbody>\n<tr>\n<th>Band<\/th>\n<th>Sensor<\/th>\n<th>Main purpose<\/th>\n<th>Native spatial resolution<\/th>\n<\/tr>\n<tr>\n<td>Band 4<\/td>\n<td>OLI<\/td>\n<td>Red reflectance for NDVI<\/td>\n<td>30 m<\/td>\n<\/tr>\n<tr>\n<td>Band 5<\/td>\n<td>OLI<\/td>\n<td>Near-infrared reflectance for NDVI<\/td>\n<td>30 m<\/td>\n<\/tr>\n<tr>\n<td>Band 10<\/td>\n<td>TIRS<\/td>\n<td>Primary thermal input for single-channel LST<\/td>\n<td>100 m<\/td>\n<\/tr>\n<tr>\n<td>Band 11<\/td>\n<td>TIRS<\/td>\n<td>Second thermal channel for specialized retrievals<\/td>\n<td>100 m<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Although Landsat thermal data are distributed on a 30 m grid in standard products, the native thermal spatial resolution is approximately 100 m. Resampling the raster to 30 m does not create new thermal detail. This distinction matters when interpreting small roofs, narrow roads, or individual urban parcels. (<a href=\"https:\/\/www.usgs.gov\/faqs\/what-are-band-designations-landsat-satellites\" target=\"_blank\" rel=\"noopener\">USGS: Landsat band designations<\/a>)<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Download_the_imagery\"><\/span>Download the imagery<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Use <a href=\"https:\/\/earthexplorer.usgs.gov\/\" target=\"_blank\" rel=\"noopener\">USGS EarthExplorer<\/a> or the <a href=\"https:\/\/www.usgs.gov\/landsat-missions\/landsat-collection-2\" target=\"_blank\" rel=\"noopener\">USGS Landsat data access portal<\/a> to obtain Landsat Collection 2 data.<\/p>\n<p>For a manual calculation, download a matching Level-1 scene containing Band 10, the optical bands, and the MTL metadata file. For NDVI, obtain compatible optical surface-reflectance data from the same acquisition.<\/p>\n<p>For most applied projects, downloading the corresponding <strong>Collection 2 Level-2 Surface Temperature product<\/strong> is preferable. It includes the <code>ST_B10<\/code> layer, which has already undergone the USGS surface-temperature retrieval workflow.<\/p>\n<p>The Copernicus Data Space Ecosystem is particularly useful when incorporating Sentinel imagery as a complementary data source, but USGS is the recommended primary archive for Landsat scenes.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Software_requirements\"><\/span>Software requirements<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Use a recent QGIS release with the built-in Raster Calculator and GDAL processing tools. Extract the downloaded Landsat files into a working directory, retaining the original band files and metadata.<\/p>\n<p>Before processing, verify that all rasters cover the same scene, use compatible grids, and have the expected coordinate reference system. Do not mix scenes from different dates simply because their filenames refer to the same city.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"3_Step-by-Step_Workflow_in_QGIS\"><\/span>3. Step-by-Step Workflow in QGIS<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p>There are two practical approaches.<\/p>\n<ul>\n<li><strong>Recommended operational method:<\/strong> Convert the existing Collection 2 Level-2 <code>ST_B10<\/code> values to degrees Celsius.<\/li>\n<li><strong>Manual retrieval method:<\/strong> Start with Level-1 thermal data, calculate radiance and brightness temperature, estimate emissivity, and apply a documented LST algorithm.<\/li>\n<\/ul>\n<p>The first method is easier to reproduce and generally preferable for operational heat mapping. The second is useful for understanding thermal remote sensing and evaluating alternative retrieval algorithms.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_1_Convert_thermal_Digital_Numbers_to_TOA_spectral_radiance\"><\/span>Step 1: Convert thermal Digital Numbers to TOA spectral radiance<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Landsat Level-1 thermal images contain calibrated digital numbers (DNs). These values must be converted into top-of-atmosphere spectral radiance before brightness temperature can be calculated.<\/p>\n<p>The relationship is:<\/p>\n<p>[<br \/>\nL_\\lambda=M_LQ_{\\mathrm{cal}}+A_L<br \/>\n]<\/p>\n<p>Where:<\/p>\n<ul>\n<li>(L_\\lambda) is spectral radiance in (\\mathrm{W,m^{-2},sr^{-1},\\mu m^{-1}}).<\/li>\n<li>(M_L) is the multiplicative radiance scaling factor.<\/li>\n<li>(Q_{\\mathrm{cal}}) is the observed DN.<\/li>\n<li>(A_L) is the additive radiance scaling factor.<\/li>\n<\/ul>\n<p>Both coefficients are supplied in the scene&#8217;s MTL metadata file. For Band 10, locate <code>RADIANCE_MULT_BAND_10<\/code> and <code>RADIANCE_ADD_BAND_10<\/code>.<\/p>\n<p>In QGIS:<\/p>\n<ol start=\"1\">\n<li>Load the Band 10 raster.<\/li>\n<li>Open <strong>Raster \u2192 Raster Calculator<\/strong>.<\/li>\n<li>Enter an expression using the coefficients from your own metadata file.<\/li>\n<\/ol>\n<pre><code>(\"B10@1\" * RADIANCE_MULT_BAND_10)\r\n+ RADIANCE_ADD_BAND_10<\/code><\/pre>\n<p>Replace the placeholders with numerical coefficient values. For example, if your metadata lists a multiplicative factor of <code>0.0003342<\/code> and an additive factor of <code>0.1<\/code>, the expression becomes:<\/p>\n<pre><code>(\"B10@1\" * 0.0003342) + 0.1<\/code><\/pre>\n<p>The example coefficients are illustrative; always use the values associated with your actual scene.<\/p>\n<p>Save the output as <code>TOA_Radiance_B10.tif<\/code>, keeping its extent, alignment, and CRS consistent with the source.<\/p>\n<p>The USGS explains this conversion and the relevant metadata coefficients in its <a href=\"https:\/\/www.usgs.gov\/landsat-missions\/using-usgs-landsat-level-1-data-product\" target=\"_blank\" rel=\"noopener\">Level-1 data processing documentation<\/a>.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_2_Convert_spectral_radiance_to_brightness_temperature\"><\/span>Step 2: Convert spectral radiance to brightness temperature<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-1606\" src=\"https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_cen03lcen03lcen0-1024x559.jpg\" alt=\"Processing thermal bands and raster calculations in QGIS for land surface temperature\" width=\"1024\" height=\"559\" srcset=\"https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_cen03lcen03lcen0-1024x559.jpg 1024w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_cen03lcen03lcen0-300x164.jpg 300w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_cen03lcen03lcen0-768x419.jpg 768w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_cen03lcen03lcen0.jpg 1408w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/p>\n<p>Brightness temperature represents the blackbody-equivalent temperature associated with the observed radiance. It is an intermediate quantity rather than the final land surface temperature.<\/p>\n<p>The conversion uses the inverse Planck relationship:<\/p>\n<p>[<br \/>\nBT=\\frac{K_2}<br \/>\n{\\ln\\left(\\frac{K_1}{L_\\lambda}+1\\right)}<br \/>\n]<\/p>\n<p>Where (K_1) and (K_2) are the band-specific thermal calibration constants from the MTL file.<\/p>\n<p>For Band 10, retrieve <code>K1_CONSTANT_BAND_10<\/code> and <code>K2_CONSTANT_BAND_10<\/code>.<\/p>\n<p>In the QGIS Raster Calculator:<\/p>\n<pre><code>K2 \/ ln((K1 \/ \"TOA_Radiance_B10@1\") + 1)<\/code><\/pre>\n<p>Substitute the actual values for <code>K1<\/code> and <code>K2<\/code>. Save the raster as <code>Brightness_Temperature_K.tif<\/code>.<\/p>\n<p>Convert the brightness temperature from Kelvin to Celsius with:<\/p>\n<pre><code>\"Brightness_Temperature_K@1\" - 273.15<\/code><\/pre>\n<p>This produces top-of-atmosphere brightness temperature in degrees Celsius. It is not yet a complete atmospheric- and emissivity-corrected LST product.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_3_Calculate_NDVI_and_fractional_vegetation_cover\"><\/span>Step 3: Calculate NDVI and fractional vegetation cover<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Vegetation affects surface heating through shading and evapotranspiration. NDVI provides a useful indicator of vegetation cover:<\/p>\n<p>[<br \/>\nNDVI=\\frac{NIR-Red}{NIR+Red}<br \/>\n]<\/p>\n<p>For Landsat 8\/9, use Band 5 for near-infrared and Band 4 for red reflectance.<\/p>\n<p>If using Collection 2 Level-2 surface reflectance, first apply the published scaling relationship:<\/p>\n<p>[<br \/>\n\\rho=DN(0.0000275)-0.2<br \/>\n]<\/p>\n<p>For example, calculate the scaled near-infrared raster in QGIS:<\/p>\n<pre><code>(\"SR_B5@1\" * 0.0000275) - 0.2<\/code><\/pre>\n<p>Then calculate scaled red reflectance:<\/p>\n<pre><code>(\"SR_B4@1\" * 0.0000275) - 0.2<\/code><\/pre>\n<p>Finally, calculate NDVI using the resulting reflectance rasters:<\/p>\n<pre><code>(\"NIR_Reflectance@1\" - \"Red_Reflectance@1\")\r\n\/\r\n(\"NIR_Reflectance@1\" + \"Red_Reflectance@1\")<\/code><\/pre>\n<p>Use the actual output layer names in your project. Mask fill values, clouds, cloud shadows, and other unsuitable pixels before calculating the index. Landsat&#8217;s <code>QA_PIXEL<\/code> band contains bit-packed quality information for this purpose. The USGS documents these flags in its <a href=\"https:\/\/www.usgs.gov\/landsat-missions\/landsat-collection-2-quality-assessment-bands\" target=\"_blank\" rel=\"noopener\">Collection 2 QA band reference<\/a>.<\/p>\n<p>Next, estimate fractional vegetation cover ((P_v)):<\/p>\n<p>[<br \/>\nP_v=<br \/>\n\\left(<br \/>\n\\frac{NDVI-NDVI_{\\min}}<br \/>\n{NDVI_{\\max}-NDVI_{\\min}}<br \/>\n\\right)^2<br \/>\n]<\/p>\n<p>A frequently used illustrative threshold pair is (NDVI_{\\min}=0.2) for sparse vegetation or bare surfaces and (NDVI_{\\max}=0.5) for dense vegetation.<\/p>\n<p>The thresholds should be adjusted to the scene&#8217;s vegetation distribution and land-cover characteristics. Avoid accepting fractions below zero or above one.<\/p>\n<p>In QGIS, use the Raster Calculator:<\/p>\n<pre><code>if(\r\n  \"NDVI@1\" &lt;= 0.2,\r\n  0,\r\n  if(\r\n    \"NDVI@1\" &gt;= 0.5,\r\n    1,\r\n    ((\"NDVI@1\" - 0.2) \/ 0.3)^2\r\n  )\r\n)<\/code><\/pre>\n<p>Save this as <code>Fractional_Vegetation.tif<\/code>. This is an empirical fractional-cover estimate rather than a direct measurement of the physical fraction of vegetation within every pixel.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_4_Estimate_land_surface_emissivity\"><\/span>Step 4: Estimate land surface emissivity<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Land surface emissivity ((\\epsilon)) describes how effectively a surface emits thermal radiation relative to a blackbody at the same temperature. Different urban materials and natural surfaces can have different emissivities, influencing the relationship between measured radiance and actual surface temperature.<\/p>\n<p>A commonly used NDVI-based approximation is:<\/p>\n<p>[<br \/>\n\\epsilon=0.004P_v+0.986<br \/>\n]<\/p>\n<p>This approximation assigns emissivity of approximately 0.986 to bare or sparsely vegetated pixels and increases it with fractional vegetation cover. It is a simplified empirical assumption, not a universal emissivity model for all roofs, asphalt mixtures, water, soil types, and building materials.<\/p>\n<p>In the QGIS Raster Calculator:<\/p>\n<pre><code>0.004 * \"Fractional_Vegetation@1\" + 0.986<\/code><\/pre>\n<p>Save the result as <code>Surface_Emissivity.tif<\/code>.<\/p>\n<p>For a research-grade thermal study, consider a more appropriate emissivity dataset or land-cover-specific emissivity assumptions. The USGS Collection 2 retrieval uses ASTER Global Emissivity Database information alongside NDVI and atmospheric inputs, illustrating why a single NDVI-based formula should be treated as an approximation.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_5_Calculate_land_surface_temperature\"><\/span>Step 5: Calculate land surface temperature<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>There is no single universal formula that can produce a fully corrected LST raster from brightness temperature and vegetation cover alone. The appropriate method depends on the atmospheric information available and the desired accuracy.<\/p>\n<h4><span class=\"ez-toc-section\" id=\"Option_A_Simplified_emissivity-corrected_single-channel_estimate\"><\/span>Option A: Simplified emissivity-corrected single-channel estimate<span class=\"ez-toc-section-end\"><\/span><\/h4>\n<p>For an exploratory workflow, a commonly used emissivity correction is:<\/p>\n<p>[<br \/>\nLST_K=<br \/>\n\\frac{BT_K}<br \/>\n{1+\\left(\\frac{\\lambda BT_K}{\\rho}\\right)<br \/>\n\\ln(\\epsilon)}<br \/>\n]<\/p>\n<p>Where:<\/p>\n<ul>\n<li>(BT_K) is brightness temperature in Kelvin.<\/li>\n<li>(\\lambda) is the effective wavelength of Band 10, approximately (10.895,\\mu m).<\/li>\n<li>(\\rho=hc\/k) is approximately (14380,\\mu m,K).<\/li>\n<li>(\\epsilon) is estimated surface emissivity.<\/li>\n<\/ul>\n<p>In the Raster Calculator:<\/p>\n<pre><code>\"Brightness_Temperature_K@1\"\r\n\/\r\n(\r\n  1\r\n  +\r\n  (\r\n    10.895 * \"Brightness_Temperature_K@1\" \/ 14380\r\n  )\r\n  * ln(\"Surface_Emissivity@1\")\r\n)<\/code><\/pre>\n<p>Convert the result to degrees Celsius:<\/p>\n<pre><code>\"LST_Emissivity_Corrected_K@1\" - 273.15<\/code><\/pre>\n<p>This approach corrects for emissivity under simplifying assumptions, but <strong>it does not explicitly correct for atmospheric absorption and emission<\/strong>. It should therefore be described as an approximate emissivity-corrected temperature estimate, not as a complete atmospheric-corrected surface-temperature retrieval.<\/p>\n<h4><span class=\"ez-toc-section\" id=\"Option_B_Full_mono-window_algorithm\"><\/span>Option B: Full mono-window algorithm<span class=\"ez-toc-section-end\"><\/span><\/h4>\n<p>A more rigorous single-channel approach explicitly incorporates atmospheric transmittance and an effective mean atmospheric temperature.<\/p>\n<p>A widely used form is:<\/p>\n<p>[<br \/>\nT_s=<br \/>\n\\frac{<br \/>\na(1-C-D)+<br \/>\n[b(1-C-D)+C+D]T_b-DT_a<br \/>\n}{C}<br \/>\n]<\/p>\n<p>Where:<\/p>\n<p>[<br \/>\nC=\\epsilon\\tau<br \/>\n]<\/p>\n<p>[<br \/>\nD=(1-\\tau)[1+(1-\\epsilon)\\tau]<br \/>\n]<\/p>\n<p>Here, (T_s) is the estimated surface temperature in Kelvin, (T_b) is brightness temperature, (\\epsilon) is emissivity, (\\tau) is atmospheric transmittance, and (T_a) is effective mean atmospheric temperature. The coefficients (a) and (b) depend on the selected Landsat algorithm and its parameterization.<\/p>\n<p>To implement this method in QGIS, prepare aligned rasters or constant rasters for the atmospheric parameters, calculate (C) and (D), and apply the equation in the Raster Calculator. Choose coefficients from a documented algorithm validated for the sensor and applicable temperature range; do not copy coefficients from a different sensor or study without checking their assumptions.<\/p>\n<p>The mono-window method is more demanding than the simplified emissivity correction because accurate atmospheric parameters must be obtained from suitable profiles or a defensible retrieval method.<\/p>\n<h4><span class=\"ez-toc-section\" id=\"What_about_the_split-window_algorithm\"><\/span>What about the split-window algorithm?<span class=\"ez-toc-section-end\"><\/span><\/h4>\n<p>A split-window approach uses both thermal channels, typically Bands 10 and 11 for Landsat 8\/9, to reduce atmospheric effects using the different absorption characteristics of the two bands. However, this requires an appropriate sensor-specific algorithm, calibration assumptions, and often atmospheric water-vapor information. Band 11 also has a documented calibration history that warrants care when selecting published algorithms.<\/p>\n<p>For general urban heat mapping, use a validated single-channel method or the supplied Level-2 surface-temperature product rather than assuming that a generic Band 10\/Band 11 expression will automatically improve accuracy.<\/p>\n<p><strong>[Image 1 placeholder: Landsat Band 10 raw scene vs resulting Land Surface Temperature map]<\/strong><\/p>\n<p><em>Alt text: &#8220;Choropleth comparison map of raw Landsat Band 10 thermal radiance versus converted Land Surface Temperature in degrees Celsius showing urban heat island hotspots.&#8221;<\/em><\/p>\n<p><strong>[Image 2 placeholder: Raster Calculator workflow flow chart for LST conversion]<\/strong><\/p>\n<p><em>Alt text: &#8220;Flowchart diagram illustrating mathematical steps to convert Landsat thermal digital numbers into brightness temperature, emissivity, and final land surface temperature.&#8221;<\/em><\/p>\n<h3><span class=\"ez-toc-section\" id=\"Step_6_The_recommended_shortcut_%E2%80%94_use_USGS_Collection_2_Level-2_surface_temperature\"><\/span>Step 6: The recommended shortcut \u2014 use USGS Collection 2 Level-2 surface temperature<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>For many applied GIS projects, the most defensible starting point is the official <code>ST_B10<\/code> raster. This product already incorporates a more comprehensive retrieval procedure using thermal observations, emissivity information, and atmospheric inputs. See the <a href=\"https:\/\/www.usgs.gov\/landsat-missions\/landsat-collection-2-level-2-science-products\" target=\"_blank\" rel=\"noopener\">USGS Collection 2 Level-2 Science Products documentation<\/a>.<\/p>\n<p>Load <code>ST_B10<\/code> into QGIS and first mask its fill pixels. Collection 2 Level-2 surface temperature has a fill value of zero. Convert valid pixels to Kelvin using:<\/p>\n<p>[<br \/>\nT_K=DN(0.00341802)+149<br \/>\n]<\/p>\n<p>Then convert Kelvin to Celsius:<\/p>\n<p>[<br \/>\nT_C=DN(0.00341802)+149-273.15<br \/>\n]<\/p>\n<p>In the Raster Calculator:<\/p>\n<pre><code>(\"ST_B10@1\" * 0.00341802) + 149 - 273.15<\/code><\/pre>\n<p>Ensure that zero-valued fill pixels are excluded rather than interpreted as physical temperatures. The <a href=\"https:\/\/www.usgs.gov\/faqs\/how-do-i-use-a-scale-factor-landsat-level-2-science-products\" target=\"_blank\" rel=\"noopener\">USGS scaling guide<\/a> documents these coefficients.<\/p>\n<p>The resulting raster is a more appropriate operational basis for UHI mapping than an emissivity-only correction, although it still requires examination of cloud contamination, quality flags, retrieval gaps, and known product limitations.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"4_Applied_Analysis_Classifying_Heat_Zones\"><\/span>4. Applied Analysis: Classifying Heat Zones<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-1607\" src=\"https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_nc4ooknc4ooknc4o-1024x559.jpg\" alt=\"Classifying urban thermal heat zones and hotspots from Landsat LST raster data\" width=\"1024\" height=\"559\" srcset=\"https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_nc4ooknc4ooknc4o-1024x559.jpg 1024w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_nc4ooknc4ooknc4o-300x164.jpg 300w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_nc4ooknc4ooknc4o-768x419.jpg 768w, https:\/\/helpwithmygisassignment.com\/blog\/wp-content\/uploads\/2026\/10\/Gemini_Generated_Image_nc4ooknc4ooknc4o.jpg 1408w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/p>\n<p>Once the LST raster is ready, style it in QGIS using <strong>Singleband pseudocolor<\/strong>. Choose a perceptually clear ramp that progresses from cool to warm colors and use consistent ranges when comparing different dates.<\/p>\n<p>For a single-scene analysis, classify temperature into relative thermal zones. One possible standard-deviation approach is:<\/p>\n<table>\n<tbody>\n<tr>\n<th>Relative class<\/th>\n<th>Example definition<\/th>\n<\/tr>\n<tr>\n<td>Cooler than typical<\/td>\n<td>LST below (\\mu-\\sigma)<\/td>\n<\/tr>\n<tr>\n<td>Moderate<\/td>\n<td>LST from (\\mu-\\sigma) to (\\mu)<\/td>\n<\/tr>\n<tr>\n<td>Warm<\/td>\n<td>LST from (\\mu) to (\\mu+\\sigma)<\/td>\n<\/tr>\n<tr>\n<td>Hotspot<\/td>\n<td>LST above (\\mu+\\sigma)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Here, (\\mu) is the mean LST for the valid study area and (\\sigma) is its standard deviation. These are relative statistical classes, not universal health-risk thresholds. They also depend on the area included in the statistics, so document whether the summary covers the urban core, the entire scene, or a defined study boundary.<\/p>\n<p>For a more direct surface urban heat-island indicator, compare average LST in the built-up area with an appropriate rural or vegetated reference zone:<\/p>\n<p>\\overline{T}_{\\mathrm{reference}}<br \/>\n]<\/p>\n<p>A positive value indicates that the urban surface is warmer than the chosen reference for that acquisition. It does not, by itself, establish a corresponding air-temperature difference.<\/p>\n<h3><span class=\"ez-toc-section\" id=\"Link_hotspots_to_vegetation_and_built_surfaces\"><\/span>Link hotspots to vegetation and built surfaces<span class=\"ez-toc-section-end\"><\/span><\/h3>\n<p>Overlay the LST map with NDVI, land-cover classifications, or canopy-cover data to investigate whether hotter pixels tend to have lower vegetation cover.<\/p>\n<p>You can also summarize temperature by administrative wards or land-use categories using zonal statistics. This allows an analyst to compare average LST across residential, industrial, commercial, and park areas.<\/p>\n<p>Use comparable acquisition dates, cloud masks, spatial resolutions, and processing methods when examining patterns over time. Differences in season, antecedent rainfall, surface moisture, and acquisition time can affect LST as much as changes in urban form.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"5_Quality_Control_and_Interpretation\"><\/span>5. Quality Control and Interpretation<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p>Before producing a final UHI map, check for the following:<\/p>\n<p><strong>Clouds and shadows:<\/strong> Mask them using quality flags. Cloud-contaminated thermal values should not be treated as valid surface temperatures.<\/p>\n<p><strong>Spatial resolution:<\/strong> Do not interpret resampled 30 m thermal pixels as independent 30 m thermal measurements. The native TIRS resolution is approximately 100 m.<\/p>\n<p><strong>Atmospheric correction:<\/strong> Distinguish brightness temperature, emissivity-adjusted estimates, and fully corrected surface-temperature products in your methodology.<\/p>\n<p><strong>Surface emissivity:<\/strong> NDVI-based emissivity is an approximation and can be unsuitable for certain urban materials. State the formula and thresholds used.<\/p>\n<p><strong>Comparison dates:<\/strong> Use consistent seasons and similar atmospheric conditions for change detection wherever possible.<\/p>\n<p><strong>Validation:<\/strong> Where available, compare the results with reliable ground measurements or independently processed temperature products, while recognizing that satellite overpass time and measurement footprints may differ from field observations.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"6_SEO_Best_Practices_and_Related_Tutorials\"><\/span>6. SEO Best Practices and Related Tutorials<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p>For a connected GIS learning resource, link this article to complementary technical guides:<\/p>\n<ul>\n<li><a href=\"https:\/\/chatgpt.com\/tutorials\/ndvi-calculation-landsat-qgis\/\" target=\"_blank\" rel=\"noopener\">Calculating NDVI from Landsat imagery<\/a><\/li>\n<li><a href=\"https:\/\/chatgpt.com\/tutorials\/remote-sensing-band-combinations\/\" target=\"_blank\" rel=\"noopener\">Remote sensing band combinations for land-cover interpretation<\/a><\/li>\n<li><a href=\"https:\/\/chatgpt.com\/tutorials\/spatial-statistics-hotspot-analysis\/\" target=\"_blank\" rel=\"noopener\">Spatial statistics and hotspot analysis in QGIS<\/a><\/li>\n<\/ul>\n<p>Replace these example paths with the canonical URLs on your website. Use descriptive anchor text and retain the main keyword in the title, introduction, and page metadata without repeating it unnaturally.<\/p>\n<p>For authoritative technical references, consult the <a href=\"https:\/\/www.usgs.gov\/landsat-missions\/landsat-8-oli-and-tirs-calibration-notices\" target=\"_blank\" rel=\"noopener\">USGS Landsat TIRS and calibration documentation<\/a>, the <a href=\"https:\/\/www.usgs.gov\/landsat-missions\/landsat-collection-2-surface-temperature\" target=\"_blank\" rel=\"noopener\">Collection 2 Surface Temperature documentation<\/a>, and the <a href=\"https:\/\/www.usgs.gov\/landsat-missions\/using-usgs-landsat-level-1-data-product\" target=\"_blank\" rel=\"noopener\">USGS Landsat Level-1 conversion guide<\/a>.<\/p>\n<h2><span class=\"ez-toc-section\" id=\"7_Conclusion_and_Key_Takeaways\"><\/span>7. Conclusion and Key Takeaways<span class=\"ez-toc-section-end\"><\/span><\/h2>\n<p>Landsat-based urban heat-island mapping combines thermal infrared observations with surface emissivity, vegetation information, and careful GIS processing. QGIS provides the raster tools needed to calculate radiance, brightness temperature, NDVI, fractional vegetation cover, and temperature maps, while Landsat Collection 2 Level-2 products provide a practical, preprocessed alternative.<\/p>\n<p>The most important distinction is between a simplified emissivity correction and a fully corrected LST retrieval. A scientifically defensible workflow should use a documented algorithm, appropriate atmospheric information, cloud and quality masking, and clear reporting of spatial-resolution limitations.<\/p>\n<p>For routine urban heat analysis, start with the USGS <code>ST_B10<\/code> surface-temperature product, convert it using the documented scale factor and offset, and classify thermal zones relative to a clearly defined study area. Use the manual Band 10 workflow when you need to understand, reproduce, or compare the retrieval process.<\/p>\n<p>A well-documented LST map can help identify potential surface heat hotspots, evaluate associations with vegetation and built-up cover, and support urban greening and climate-adaptation planning\u2014without confusing satellite-measured surface temperature with air temperature or direct human heat exposure.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Urban Heat Island Mapping: Calculating Land Surface Temperature (LST) from Landsat TIR Bands 1. Introduction: Mapping Urban Heat Islands with Landsat This Urban Heat Island Mapping Landsat LST Tutorial explores how built-up areas exhibit higher temperatures than their surrounding rural or less-developed environments. Urban Heat Islands (UHIs) occur when built-up areas exhibit higher temperatures than &#8230; <a title=\"Urban Heat Island Mapping: Calculating Land Surface Temperature (LST) from Landsat TIR Bands\" class=\"read-more\" href=\"https:\/\/helpwithmygisassignment.com\/blog\/urban-heat-island-mapping-landsat-lst-tutorial\/\" aria-label=\"Read more about Urban Heat Island Mapping: Calculating Land Surface Temperature (LST) from Landsat TIR Bands\">Read more<\/a><\/p>\n","protected":false},"author":1,"featured_media":1605,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1,4],"tags":[],"class_list":["post-1446","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized","category-gis-articles"],"_links":{"self":[{"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/posts\/1446","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/comments?post=1446"}],"version-history":[{"count":3,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/posts\/1446\/revisions"}],"predecessor-version":[{"id":1609,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/posts\/1446\/revisions\/1609"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/media\/1605"}],"wp:attachment":[{"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/media?parent=1446"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/categories?post=1446"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/helpwithmygisassignment.com\/blog\/wp-json\/wp\/v2\/tags?post=1446"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}