finding critical points with Hessian Matrix



I'm trying to find morse critical points in an image by using the
eigenvalues of hessian matrix. For some reason i get erroneous
results, think it must be my coding.

here is an example:

img=peaks(250);

[imgX,imgY] = gradient(img);
[imgxx,imgyx]=gradient(imgX);
[imgxy,imgyy]=gradient(imgY);

Hess=cell(length(img(:,1)));
for row=1:length(img(:,1));
for col=1:length(img(1,:));
Hess(row,col)={eig([imgxx(row,col) imgxy(row,col); imgyx(row,col)
imgyy(row,col)])};
end
end

i get nothing like a should get, most of the results are non-zero and
i cannot then pin point the max, min or saddle points this way. Am i
missing something?
.



Relevant Pages

  • Re: want to save memory in sparse matrix
    ... is (PARPACK and ARPACK) answer for me. ... basically we construct a hessian matrix. ... then we diagonalise it and get ... the eigenvalues. ...
    (comp.lang.fortran)
  • Re: ARPACK to get all eigenvalues
    ... i have a code with a very large hessian matrix. ... i have to diagonalise it to find the eigenvalues for the calculation of normal mode analysis. ... after searching on net i came across ARPACK and the parallel version PARPACK. ... RWTH - Aachen University ...
    (sci.math.num-analysis)
  • Re: random matrix eigenvalues
    ... outside the unit circle, creating curious exceptions. ... title(['scaled eigenvalues of random matrix of rank ' ... Gershgorin disks give you the reason why it ... Note that the same thing holds for UNIFORM ...
    (comp.soft-sys.matlab)
  • Re: why does professor david c ullrich have to put people down to feel good about himself?
    ... I already know about eigenvalues. ... reason") and v an eigenvector of the matrix. ... invertible matrices A and B, ... Curiouser and curiouser. ...
    (sci.math)
  • ARPACK to get all eigenvalues
    ... hessian matrix. ... i have to diagonalise it to find the ... eigenvalues for the calculation of normal mode analysis. ... eigenvalues with ARPACK? ...
    (sci.math.num-analysis)

Loading