mat7labs的个人空间 http://blog.eetop.cn/1631436 [收藏] [复制] [分享] [RSS]

日志

LVQ编码

已有 597 次阅读2019-10-1 20:57

function y=lvqcode(x,v)
v=double(v);
xorig=x;
nx=size(x,1);
ny=size(x,2);
x=double(x);
s=sqrt(size(v,2));
%take blockwise
%--------------------------------
x1=im2col(x,[s s],'distinct');
x=x1';
%---------------------------------

c=size(v,1);
n=size(x,1);

%-------------find nearest code--------------------
for i=1:n
mini=(((repmat(x(i,:),[c 1])-v(:,:)).^2))';
if(not (size(x,2)==1))
    mini=(sum(mini))';
end
mini=sqrt(mini);

m=min(mini);
temp=find(mini==m);
if(length(temp)==1)
    index(i)=temp;
else
index(i)=temp(1);
end
clear temp;
end
%------------------------------------------------------
%y=col2im(index,[1 1],[nx/s ny/s],'distinct');
y=reshape(index,nx/s,ny/s);

评论 (0 个评论)

facelist

您需要登录后才可以评论 登录 | 注册

关闭

站长推荐上一条 /3 下一条

小黑屋| 关于我们| 联系我们| 在线咨询| ET创芯网 ( 京ICP备:10050787号 京公网安备:11010502037710 )

GMT+8, 2019-12-6 22:46 , Processed in 0.026327 second(s), 5 queries , Gzip On, Redis On.

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

返回顶部