Определитель матрицы hessian полутонового изображения слишком мал в matlab

Question

Я пытаюсь найти определитель матрицы hessian изображения в оттенках серого 50x50. Определитель матрицы, которую я получаю, является очень маленьким значением i.e 4.7612e-134. Я думаю, что чего-то не хватает. Мой код ниже. Благодаря

% computing second derivatives in each direction first 

[gx, gy] = gradient(double(sliceOfImageK2)); 

[gxx, gxy] = gradient(gx); 

[gyx, gyy] = gradient(gy); 

hessianMatrix = [gxx gxy; gxy gyy]; 

determinantHessianMatrix = det(hessianMatrix) 

Answer 1

Я не думаю, что вы должны собрать 100x100 матрицу, если вы хотите назвать это Hessian. Соберите вместо этого матрицу 2x2 на каждый из 50x50 (2500) пикселей, где вы отбираете свои производные.

Это 2500 гессианы, выраженное в 2500x4 матрицы:

H = [gxx(:) gxy(:) gyx(:) gyy(:)] 

Здесь выражается в виде 2500 2x2 матриц:

H_ = reshape(H', 2, 2, length(H)) 

и эти детерминанты каждой 2x2 матрицы:

D = H(:,1).*H(:,4) - H(:,2).*H(:,3) 

Здесь как 50x50 матрица с определителем Гесса в каждом пикселе, если это то, что вы после:

reshape(D, 50, 50)