Google Earth Engine（GEE）——缨帽变换

最新推荐文章于 2024-12-11 23:43:51 发布

我也想吃a

最新推荐文章于 2024-12-11 23:43:51 发布

阅读量6k

点赞数

分类专栏： GEE教程文章标签：开发语言 javascript

本文链接：https://blog.csdn.net/weixin_48048649/article/details/128994700

版权

GEE教程专栏收录该内容

19 篇文章

订阅专栏

缨帽（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");

展示