独立成分分析作为图像

时间:2013-11-10 17:47:08

标签: image matlab plot

我使用FastICA对面部图像集应用了独立分量分析。我已成功检索到独立组件和混合矩阵。独立组件的值是双倍的,我想将这些组件显示为Web上可用的图像,例如http://research.ics.aalto.fi/ica/imageica/。我不知道如何表明这一点。

2 个答案:

答案 0 :(得分:1)

使用montage功能可以更自动地执行此操作的一种方法。或者你可以查看subimage。另请参阅this StackOverflow question and answer

答案 1 :(得分:0)

尝试以下功能

function [XX,fh]=dispImgs(X,cols,gap,ihw,fh)
% Courtesy A. Leonardis, D. Skocaj
% see http://vicos.fri.uni-lj.si/danijels/downloads

[M,N]=size(X);
if nargin<2 cols=floor(sqrt(N)); end;
if nargin<3 gap=0; end;
if nargin<4 ihw=[sqrt(M),sqrt(M)]; end;
if nargin<5 fh = figure; end; % new figure

ih=ihw(1);iw=ihw(2);
maxv=max(X(:));
rows=floor(N/cols);
XX=zeros((rows*ih)+(rows-1)*gap,(cols*iw)+(cols-1)*gap)+maxv;

for i=1:N
   a=(iw+gap)*mod(i-1,cols)+1;
   b=(iw+gap)*mod(i-1,cols)+iw;
   c=(ih+gap)*(floor((i-1)/cols))+1;
   d=(ih+gap)*(floor((i-1)/cols))+ih;
   XX(c:d,a:b)=reshape(X(:,i)',ih,iw);
end;

xxmax=max(XX(:));
xxmin=min(XX(:));

fh = figure(fh);
imshow((XX-xxmin)/(xxmax-xxmin));
axis off;

示例:

X: imsize by N
dispImgs( X, 8, 4, imsize );%show all N images in 8 columns with Gap=4