PCA OF AN IMAGE...

6 次查看(过去 30 天)
divya
divya 2013-5-23
I want to find PCA of an image... but when i run the code i get following error..
Error using - Integers can only be combined with integers of the same class, or scalar doubles. Error in pca (line 26) data = data - repmat(mn,1,N);
here is the code....'
function [signals,PC,V] = pca()
I = imread('C:\Users\div\Pictures\Picture1.png');
I = rgb2gray( I);
[irow, icol] = size(I);
data = reshape(I',irow*icol,1);
% [sR ,sC ,eR ,eC] = deal(1, 3, 2, 4);
% Compute the sum over the region using the integral image.
% PCA1: Perform PCA using covariance.
% data - MxN matrix of input data
% (M dimensions, N trials)
% signals - MxN matrix of projected data
% PC - each column is a PC
% V - Mx1 matrix of variances
[M,N] = size(data);
% subtract off the mean for each dimension
mn = mean(data,2);
data = data - repmat(mn,1,N);
% calculate the covariance matrix
covariance = 1 / (N-1) * data * data;
% find the eigenvectors and eigenvalues
[PC, V] = eig(covariance);
% extract diagonal of matrix as vector
V = diag(V);
% sort the variances in decreasing order
[junk, rindices] = sort(-1*V);
V = V(rindices); PC = PC(:,rindices);
% project the original data set
signals = PC * data;

采纳的回答

Image Analyst
Image Analyst 2013-5-24
That is a dangerous way that you used size(). See Steve's blog: http://blogs.mathworks.com/steve/2011/03/29/even-more-information-about-the-size-function/
A common way to deal with your error is to just convert everything to double. But if you do then you need to use [] to see your images or else they will show up as all white because double images are expected to be in the range 0-1.
imshow(doubleImage, []);
The [] will adjust the display to handle whatever range your array may have, and not require it to be in the 0-1 range.
  4 个评论
divya
divya 2013-5-24
Undefined function 'eig' for input arguments of type 'uint8'.
Error in pca (line 46) [PC, V] = eig(covariance);
still getting this error.....:(
Image Analyst
Image Analyst 2013-5-24
So cast it to double, or do
which eig

请先登录,再进行评论。

更多回答(3 个)

yagnesh
yagnesh 2013-5-24
try data = double(data)- double(repmat(mn,1,N));

farahnaz
farahnaz 2014-3-22
please use this cod :
function [signals,PC,V] = pca() I=imread('G:\Users\farahnaz\Downloads\018.jpg'); b= im2double(I); imshow(b, []); [irow, icol] = size(b); data = reshape(b',irow*icol,1);
[M,N] = size(data);
mn = mean(data,2);
data = data - repmat(mn,1,N);
covariance = cov(data)';
[PC, V] = eig(covariance);
V = diag(V);
[junk, rindices] = sort(-1*V);
V = V(rindices); PC = PC(:,rindices);
  4 个评论
Shaveta Arora
Shaveta Arora 2016-1-31
pls elaborate about the pcimage variable. how to get PC image as 2D
Image Analyst
Image Analyst 2016-1-31
I don't see how he's getting PC images. There aren't any. You can see my attached demo if you want.

请先登录,再进行评论。


Shaveta Arora
Shaveta Arora 2016-2-1
actually I want to recover the image using principal components after the following instruction:
% project the original data set
signals = PC * data;
pls help me in recovering the image.

类别

Help CenterFile Exchange 中查找有关 Dimensionality Reduction and Feature Extraction 的更多信息

标签

尚未输入任何标签。

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by