kt4

From DavinciWiki
Revision as of 12:29, 2 May 2013 by Cedwards (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Description

Generates the "Tasseled Cap" KT4 derived imagery from a given set of LandSat Thematic Mapper imagery.

Arguments and Return Values

Arguments: 6 Thematic Mapper bands (1-5 and 7), and optionally a mode parameter

Return Value: A float greyscale image

Usage

Syntax: kt4(tm1, tm2, tm3, tm4, tm5, tm7 [, mode = INT])

This function (and the related functions kt1, kt2, kt3, kt5, and kt6) generate imagery derived from LandSat Thematic Mapper data using the "Tasseled Cap" method (similar to pca, but with a fixed matrix of coefficients).

The optional mode parameter specifies the format of the input data: 0 = LandSat 4 reflectance data, 2 = LandSat 7 reflectance data, 3 = LandSat 4 DN ("digital number") data, 4 = LandSat 5 DN ("digital number") data. If this parameter is not specified, then it defaults to 0.

The coefficients for this function are based on the tables in the Arizona Remote Sensing Center, Office of Arid Lands Studies' documentation.


Examples

dv> tm1
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> tm2
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> tm3
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> tm4
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> tm5
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> tm7
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> result = kt4(tm1, tm2, tm3, tm4, tm5, tm7)
Generating a KT4 (Tasseled Cap - Fourth Component) image.
2000x1000x1 array of float, bsq format [8,000,000 bytes]

dv> max(result)
0.1139493585

dv> min(result)
-0.2797231674

Images

Input bands:

NYC-Band Comparison.PNG

Result of KT4 processing:

NYC-KT4.PNG

False colour image where red = KT4, green = KT5, and blue = KT6:

NYC-KT4-KT5-KT6.PNG

False colour image where red = KT4, green = KT2, and blue = KT6:

NYC-KT4-KT2-KT6.PNG

DavinciWiki Mini-Nav Bar

Contents


Contact Developers

  • davinci-dev [AT] mars.asu.edu

All other topics

  • See navigation on the left

Recent Library Changes

Created On: 07-16-2012
Modified On: 05-02-2013

Personal tools