how to gamma transform a matrix cancel inf

I would like to know if its possible to use the transformation gamma on a matrix of images directly. Problem is with the eigen function. I have to change the format to a double i did something like this

dbx=[];    
ddbx=[];
for i=1:M
temp=double(S(:,i));
dbx=[dbx temp];
gammaimage= gamma(dbx(:,i));
ddbx=[ddbx gammaimage];
end

i am using the eigen function eig() and get

Error using eig
Input to EIG must not contain NaN or Inf.

not sure how to go about this. Will i have to abandon the idea of using gamma transformation?

Answer 1

[Walter Roberson]

gamma() is the Gamma function, the extension of the factorial function. gamma(171) = factorial(170) is your limit before you hit infinity.

I would suggest that you are looking for Gamma Correction rather than gamma()