# PCA OF AN IMAGE...

21 visualizzazioni (ultimi 30 giorni)
divya il 23 Mag 2013
Risposto: Shaveta Arora il 1 Feb 2016
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 = 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;
##### 0 CommentiMostra -2 commenti meno recentiNascondi -2 commenti meno recenti

Accedi per commentare.

### Risposta accettata

Image Analyst il 24 Mag 2013
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 CommentiMostra 2 commenti meno recentiNascondi 2 commenti meno recenti
divya il 24 Mag 2013
Undefined function 'eig' for input arguments of type 'uint8'.
Error in pca (line 46) [PC, V] = eig(covariance);
still getting this error.....:(
Image Analyst il 24 Mag 2013
So cast it to double, or do
which eig

Accedi per commentare.

### Più risposte (3)

yagnesh il 24 Mag 2013
try data = double(data)- double(repmat(mn,1,N));
##### 0 CommentiMostra -2 commenti meno recentiNascondi -2 commenti meno recenti

Accedi per commentare.

farahnaz il 22 Mar 2014
[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 CommentiMostra 2 commenti meno recentiNascondi 2 commenti meno recenti
Shaveta Arora il 31 Gen 2016
pls elaborate about the pcimage variable. how to get PC image as 2D
Image Analyst il 31 Gen 2016
I don't see how he's getting PC images. There aren't any. You can see my attached demo if you want.

Accedi per commentare.

Shaveta Arora il 1 Feb 2016
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.
##### 0 CommentiMostra -2 commenti meno recentiNascondi -2 commenti meno recenti

Accedi per commentare.

### Categorie

Scopri di più su Dimensionality Reduction and Feature Extraction in Help Center e File Exchange

### Tag

Non è stata ancora inserito alcun tag.

### Community Treasure Hunt

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

Start Hunting!

Translated by