缨帽(Kauth-Thomas)变换,又称K-T变换,是一种基于图像物理特征的固定转换,变换后指向与地物类型和变化有密切关系的方向,特别是与植物生长过程和土壤有关,其转换系数对同一传感器是固定的。缨帽变换可应用于图像压缩、图像去噪、图像增强、判读解译等方面。
缨帽变换是一种特殊的主成分分析,与主成分分析不同的是其转换系数是固定的,因此它独立于单个图像,不同图像产生的土壤亮度和绿度可以互相比较。
代码
// 定义变换系数
var coefficients = ee.Array([
[0.3037, 0.2793, 0.4743, 0.5585, 0.5082, 0.1863],
[-0.2848, -0.2435, -0.5436, 0.7243, 0.0840, -0.1800],
[0.1509, 0.1973, 0.3279, 0.3406, -0.7112, -0.4572],
[-0.8242, 0.0849, 0.4392, -0.0580, 0.2012, -0.2768],
[-0.3280, 0.0549, 0.1075, 0.1855, -0.4357, 0.8085],
[0.1084, -0.9022, 0.4120, 0.0573, -0.0251, 0.0238]
]);
//影像转换为array(筛选影像对应波段,将波段数据转换成二维数组矩阵)
var image = ee.Image("LANDSAT/LT05/C01/T1_TOA/LT05_044034_20081011")
.select(['B1', 'B2', 'B3', 'B4', 'B5', 'B7']);
var arrayImage1D = image.toArray();
var arrayImage2D = arrayImage1D.toArray(1);
//相乘变换
//根据缨帽变换规则将系数矩阵和影像矩阵做矩阵运算,然后根据行做平展变换
var componentsImage = ee.Image(coefficients)
.matrixMultiply(arrayImage2D)
.arrayProject([0])
.arrayFlatten([[
'brightness', 'greenness', 'wetness',
'fourth', 'fifth', 'sixth'
]]);
var vizParams = {
"bands": ['brightness', 'greenness', 'wetness'],
"min": -0.1,
"max": [0.5, 0.1, 0.1]
};
Map.centerObject(image, 7);
Map.addLayer(componentsImage, vizParams, "componentsImage");
展示