fastICA_imag.zip_fastica图像_图像分离_混合图像分离

%%%%快速ICA算法在混合图像分离中的应用%%%%% %_________% clc clear all I1 = imread('rice.png'); I3 = imread('E:\李福\陈庆强\资料\lena256.png'); I2 = imread('E:\李福\陈庆强\资料\cameraman.png'); subplot(331),imshow(I1),subplot(332),imshow(I2),subplot(333),imshow(I3); %_________% s1 = reshape(I1,[1,256*256]); s2 = reshape(I2,[1,256*256]); s3 = reshape(I3,[1,256*256]); s = [s1;s2;s3]; sig = double(s); %Aorig = rand(3); Aorig = [ 0.7446 0.9334 0.8392 0.2679 0.6833 0.6288 0.4399 0.2126 0.1338 ]; mixedsig = Aorig*sig; omixedsig = Aorig*sig; %____显示混合图像_____% ms1 = reshape(mixedsig(1,:),[256,256]); ms2 = reshape(mixedsig(2,:),[256,256]); ms3 = reshape(mixedsig(3,:),[256,256]); MI1 = uint8(round(ms1));MI2 = uint8(round(ms2));MI3 = uint8(round(ms3)); subplot(334),imshow(MI1),subplot(335),imshow(MI2),subplot(336),imshow(MI3), %_____高通滤波_____% %[b,a]=butter(10,0.5,'high'); %mixedmighaha=colfilt(ms1,[256 256],'distinct','filter(P1,P2,x)',b,a); %Y=im2col(mixedmighaha,[256 256],'distinct'); %mixedsig(1,:)=Y'; %mixedmighaha=colfilt(ms2,[256 256],'distinct','filter(P1,P2,x)',b,a); %Y=im2col(mixedmighaha,[256 256],'distinct'); %mixedsig(2,:)=Y'; %mixedmighaha=colfilt(ms3,[256 256],'distinct','filter(P1,P2,x)',b,a); %Y=im2col(mixedmighaha,[256 256],'distinct'); %mixedsig(3,:)=Y'; %_____中心化_____% mixedsig = zeros(size(omixedsig)); meanValue = mean(omixedsig')'; mixedmean = meanValue; mixedsig = omixedsig - meanValue*ones(1,size(omixedsig,2)); %_____白化_____% covarianceMatrix = cov(mixedsig',1); [E,D] = eig(covarianceMatrix); eigenvalues = flipud(sort(diag(D))); whiteningMatrix = inv(sqrt(D))*E'; dewhiteningMatrix = E*sqrt(D); whitesig = whiteningMatrix*mixedsig; X = whitesig; [vectorSize,numSamples] = size(X); B = zeros(vectorSize); numOfIC = vectorSize; %_____迭代过程_____% for r = 1:numOfIC i = 1; maxNumIterations = 20; w = ones(vectorSize,1) ; w = w/norm(w); while i<=maxNumIterations +1 w = w - B*B'*w; w = w/norm(w); w = (X*((X'*w).^3))/numSamples - 3*w; w = w/norm(w); i = i + 1; end A(:,r) = dewhiteningMatrix*w; W(r,:) = w'*whiteningMatrix; B(:,r) = w; end icasig = W*mixedsig + (W*mixedmean)*ones(1,numSamples); icasig = abs(50*icasig); Is1 = reshape(icasig(1,:),[256,256]); Is2 = reshape(icasig(2,:),[256,256]); Is3 = reshape(icasig(3,:),[256,256]); II1 = uint8(round(Is1)); II2 = uint8(round(Is2)); II3 = uint8(round(Is3)); %______显示结果______% subplot(337),imshow(II1),subplot(338),imshow(II2),subplot(339),imshow(II3)