icacalcpi.m

基于ica的图像处理工具包 里面有大量的经典算法和用于处理的图像 有图形界面直接运行icalab.p文件

function [G, pi] = icacalcpi( W, Q, A );

% ICACALCPI.M
% Copyright (c) by Krzysztof Siwek

G = abs( W * Q * A ) ;
[Nrows, Ncols] = size( G );
N = Nrows;
[I, J] = max( G' );        
G = pinv(diag(I))*G;
G = G';
if Nrows ~= Ncols,
   pi = -1;
else
   pi = (sum(sum(G)-max(G))/(Nrows*Ncols-length(max(G))));
end
G = G';