Spatial Modeler Tutorials

Tasseled Cap Transformation for Landsat 8

by on - edited on by (7,357 Views)

Download model

Description:

The Tasseled Cap Transformation is a technique commonly used in land cover mapping or other classification projects. It takes the linear combination of satellite imagery bands and a specialized coefficient matrix to create a n-band image with the first 3 bands containing the majority of the useful information, similar to Principal Component Analysis. The first three bands created are generally held to represent Brightness, Greenness and Wetness. The coefficient matrix, which is unique to each imaging sensor, is based on image statistics and empirical observations. The matrix^[1] used in this model for the Landsat 8 OLI sensor is as follows (note that the first band of Landsat 8 OLI imagery, coastal aerosol, is effectively dropped):

	LS 8 Band 1	LS 8 Band 2	LS 8 Band 3	LS 8 Band 4	LS 8 Band 5	LS 8 Band 6	LS 8 Band 7
Brightness	0.0000	0.3029	0.2786	0.4733	0.5599	0.5080	0.1872
Greenness	0.0000	-0.2941	-0.2430	-0.5424	0.7276	0.0713	-0.1608
Wetness	0.0000	0.1511	0.1973	0.3283	0.3407	-0.7117	-0.4559
TCT4	0.0000	-0.8239	0.0849	0.4396	-0.0580	0.2013	-0.2773
TCT5	0.0000	-0.3294	0.0557	0.1056	0.1855	-0.4349	0.8085
TCT6	0.0000	0.1079	-0.9023	0.4119	0.0575	-0.0259	0.0252

Other Landsat sensors can be transformed to Tasseled Cap from the Raster tab / Resolution group / Spectral pulldown.

Assumptions: The Landsat 8 OLI imagery should have been previously corrected to at-satellite reflectance values. This greatly increases the accuracy of the tasseled-cap transform. Information on this procedure can be found here: http://landsat.usgs.gov/Landsat8_Using_Product.php

Below you can see how the resulting first 3 bands compare individually to the original image, which is displayed RGB to 4/3/2.

First Band of TCT, Brightness, compared to the original image on right:

Second Band of TCT, Greenness, compared to the original image on right::

Third Band of TCT, Wetness, compared to the original image on right::

TasseledCapLandsat8-v15-1-1.gmdx

Reference:

Coefficient matrix derived from:

Muhammad Hasan Ali Baig, Lifu Zhang, Tong Shuai & Qingxi Tong (2014)

Derivation of a tasselled cap transformation based on Landsat 8 at-satellite reflectance,
Remote Sensing Letters, 5:5, 423-431, DOI: 10.1080/2150704X.2014.915434

Rate this article:

Yes No

Comments

Here's another relevant spatial recipe for Landsat 8 imagery from our partner Sterling Geo: "Calculate the Top of the Atmosphere Reflectance on Landsat 8 Images" found at http://www.sterlinggeo.com/spatial-modeler-library-index/ Thanks @JohnnieJ for this tip!

Permalink

Just a heads up, I was having trouble with this model, I then worked out that I needed to change the input raster type to a float rather than the default integer. My input was a float as this was the output of from the TOA conversion model (https://community.hexagongeospatial.com/t5/Spatial-Recipes/Landsat-8-Top-of-Atmosphere-Reflectance-C...)

tct input.PNG

Permalink

Thanks Lisa,

I'm not sure why the original author hardwired the Raster Input operator to set the raster steam to Integer. If your input data happens to be Reflectance values that are still in the 0 - 1 floating point range, they'll get truncated to 0s or 1s. They may have set it this way because a common practice for Reflectance data is to multiply up by 10,000 or 20,000 and truncate to 16-bit integer to save a bit of disk space. If the input had been left its default "Native" setting it should work with both Float and Integer style Reflectance inputs.

I'll fix the model tomorrow and re-post a new version to the original article.

Thanks for pointing this out!

Cheers

Permalink

Hi, sorry I can´t open the link: http://www.sterlinggeo.com/spatial-modeler-library-index/

is there another way?

TIA

Permalink

Sure - just use the Download Model link at the top of the article (you will need to be registered and logged in to do so).

Cheers

Permalink

[{"id":"c_2","title":"Hexagon Geospatial Recipes","items":[{"id":"26453","title":"Report Factor Loadings (Correlation) between Principal Components and Input Bands","view_href":"/t5/Spatial-Recipes/Report-Factor-Loadings-Correlation-between-Principal-Components/ta-p/26453","isArticle":true,"isUnordered":false,"items":[]},{"id":"22616","title":"Machine Learning Example: Classifying Features","view_href":"/t5/Spatial-Recipes/Machine-Learning-Example-Classifying-Features/ta-p/22616","isArticle":true,"isUnordered":false,"items":[]},{"id":"1609","title":"Tasseled Cap Transformation for Landsat 8","view_href":"/t5/Spatial-Recipes/Tasseled-Cap-Transformation-for-Landsat-8/ta-p/1609","isArticle":true,"isUnordered":false,"items":[]},{"id":"1766","title":"Simple Cloud Mask for Landsat Imagery","view_href":"/t5/Spatial-Recipes/Simple-Cloud-Mask-for-Landsat-Imagery/ta-p/1766","isArticle":true,"isUnordered":false,"items":[]},{"id":"1774","title":"Stretch Image through Saved Lookup Table (LUT)","view_href":"/t5/Spatial-Recipes/Stretch-Image-through-Saved-Lookup-Table-LUT/ta-p/1774","isArticle":true,"isUnordered":false,"items":[]},{"id":"739","title":"Urban Change Detection","view_href":"/t5/Spatial-Recipes/Urban-Change-Detection/ta-p/739","isArticle":true,"isUnordered":false,"items":[]},{"id":"1868","title":"Saving a No Stretch LUT into the header of an existing image","view_href":"/t5/Spatial-Recipes/Saving-a-No-Stretch-LUT-into-the-header-of-an-existing-image/ta-p/1868","isArticle":true,"isUnordered":false,"items":[]},{"id":"774","title":"ASTER Band Ratios","view_href":"/t5/Spatial-Recipes/ASTER-Band-Ratios/ta-p/774","isArticle":true,"isUnordered":false,"items":[]},{"id":"772","title":"Attach Attributes to Raster","view_href":"/t5/Spatial-Recipes/Attach-Attributes-to-Raster/ta-p/772","isArticle":true,"isUnordered":false,"items":[]},{"id":"770","title":"Blended Filters","view_href":"/t5/Spatial-Recipes/Blended-Filters/ta-p/770","isArticle":true,"isUnordered":false,"items":[]},{"id":"768","title":"Clump with Colors","view_href":"/t5/Spatial-Recipes/Clump-with-Colors/ta-p/768","isArticle":true,"isUnordered":false,"items":[]},{"id":"776","title":"PC with Percentage Contribution","view_href":"/t5/Spatial-Recipes/PC-with-Percentage-Contribution/ta-p/776","isArticle":true,"isUnordered":false,"items":[]},{"id":"758","title":"Fill Zeros","view_href":"/t5/Spatial-Recipes/Fill-Zeros/ta-p/758","isArticle":true,"isUnordered":false,"items":[]},{"id":"756","title":"Helicopter Landing Zones","view_href":"/t5/Spatial-Recipes/Helicopter-Landing-Zones/ta-p/756","isArticle":true,"isUnordered":false,"items":[]},{"id":"778","title":"Integer to String","view_href":"/t5/Spatial-Recipes/Integer-to-String/ta-p/778","isArticle":true,"isUnordered":false,"items":[]},{"id":"748","title":"Pan Sharpening: Several Methods","view_href":"/t5/Spatial-Recipes/Pan-Sharpening-Several-Methods/ta-p/748","isArticle":true,"isUnordered":false,"items":[]},{"id":"744","title":"Color Shaded Relief","view_href":"/t5/Spatial-Recipes/Color-Shaded-Relief/ta-p/744","isArticle":true,"isUnordered":false,"items":[]}]},{"id":"c_3","title":"Partner Spatial Recipes","items":[{"id":"28115","title":"Orthorectify using Ground Control Points (GCPs)","view_href":"/t5/Spatial-Modeler-Tutorials/Orthorectify-using-Ground-Control-Points-GCPs/ta-p/28115","isArticle":true,"isUnordered":false,"items":[]},{"id":"27448","title":"Extract Coastline from WorldView-2 & -3 satellite imagery","view_href":"/t5/Spatial-Modeler-Tutorials/Extract-Coastline-from-WorldView-2-amp-3-satellite-imagery/ta-p/27448","isArticle":true,"isUnordered":false,"items":[]},{"id":"10322","title":"Sentinel-2 Stack / Import","view_href":"/t5/Spatial-Modeler-Tutorials/Sentinel-2-Stack-Import/ta-p/10322","isArticle":true,"isUnordered":false,"items":[]},{"id":"1336","title":"IMAGINE UAV","view_href":"/t5/Spatial-Recipes/IMAGINE-UAV/ta-p/1336","isArticle":true,"isUnordered":false,"items":[]},{"id":"1338","title":"Fuze Go™ MS Sharp","view_href":"/t5/Spatial-Recipes/Fuze-Go-MS-Sharp/ta-p/1338","isArticle":true,"isUnordered":false,"items":[]}]}]

Courses

Spatial Modeler Tutorials

Tasseled Cap Transformation for Landsat 8

Custom Matrix Input

Landsat 8

Linear Combination

Spatial Modeler

Tasseled Cap

Download model

Description:

Reference: