Heres what you do
Compute the gradient at each point in
the image
Create the H matrix from the entries in the gradient
Compute the eigenvalues.
Find points with large response (l- > threshold)
Choose those points where l- is a local maximum as features