- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
K_L变换
K-L变换·定义
以矢量信号X的协方差矩阵Ф的归一化正交特征矢量q所构成的正交矩阵Q,来对该矢量信号X做正交变换Y=QX,则称此变换为K-L变换(K-LT或KLT),K-LT是Karhuner-Loeve变换的简称,有的文献资料也写作KLT。可见,要实现KLT,首先要从信号求出其协方差矩阵Ф,再由Ф求出正交矩阵Q。Ф的求法与自相关矩阵求法类似。
K-L变换的特性
(1)去相关特性。
K-L变换是变换后的矢量信号Y的分量互不相关。
能量集中性。
所谓能量集中性,是指对N维矢量信号进行K-L变换后,最大的方差见集中在前M个低次分量之中(MN)。
(3)最佳特性。
K-L变换是在均方误差测度下,失真最小的一种变换,其失真为被略去的各分量之和。由于这一特性,K-L变换被称为最佳变换。许多其他变换都将K-L变换作为性能上比较的参考标准。
(4)无快速 \o 算法与数据结构知识库 \t /aalbertini/article/details/_blank 算法,且变换矩阵随不同的信号样值集合而不同。
这是K-L变换的一个缺点,是K-L变换实际应用中的一个很大障碍。
K-L变换的原理
代码
% K_L 变换
X=imread(E:\beijingcut.tif); %读图像 X是一个波段图像图像的宽和高610*740
XX=reshape(X,[610*740,6]); %将图像转换化为一个波段的图像
pic1=XX(:,1);
pic1=reshape(pic1,[610,740]);
figure(1);
imshow(pic1,[]);
title(Original_Image);
XX=im2double(XX); %图像转化为双精度
covx = cov(XX); % 求图像的协方差
[pc,variances,explained] = pcacov(covx); % 求解K_L变换矩阵,即PCA
pc
myKL=XX*pc; %进行K_L变换
pic2=myKL(:,1);
pic2=reshape(pic2,[256,256]);
figure(2);
imshow(pic2,[]);
title(K_L_a);
figure(3);
pic3=myKL(:,2);
pic3=reshape(pic3,[256 256]);
imshow(pic3,[]);
title(K_L_b);
my=myKL*pc;
my1=my(:,1);
my1=reshape(my1,[256,256]);
figure(4);
imshow(my1,[]);
展示
文档评论(0)