>

[leekk8]

Це MATLAB файл обчислення орієнтації області пальців.Однак, я, можливо, буде потрібно математичні рівняння або теорії коду.Будь ласка, допоможіть, якщо хто-небудь зрозуміти.TQ

[Orientim функції, надійність] = ...
ridgeorient (IM, gradientsigma, blocksigma, orientsmoothsigma)

[Рядки, колонки] = розмір (мкм);

% Розрахувати зображень градієнтів.
Цзе виправити = (6 * gradientsigma), якщо ~ мода (Цзе, 2); Зі Зі = 1; кінець
F = fspecial ('гаусівської, Цзе, gradientsigma);% Створити фільтр Гауса.
[/ Х, ф] = градієнт (F);% градієнт Gausian.

Gx = фільтр2 (FX, мкм);% Градієнт зображення в х
Гр = фільтр2 (ФГ, КМ);% ...і у

Розрахунковий% орієнтацією на місцевий хребта в кожній точці знаходження
% Головною осі зміни в образі градієнтів.

Gxx = Gx. ^ 2;% коваріантність даних для зображень градієнтів
Gxy = Gx .* Гр;
Gyy = Гр. ^ 2;

% Тепер гладкою коваріаційний даних для здійснення зваженого підсумовування
% Даних.
Цзе виправити = (6 * blocksigma), якщо ~ мода (Цзе, 2); Зі Зі = 1; кінець
F = fspecial ('гаусівської, Цзе, blocksigma);
Gxx фільтр2 = (F, Gxx);
Gxy = 2 * фільтр2 (F, Gxy);
Gyy фільтр2 = (F, Gyy);

% Аналітичне рішення основних напрямки
знаменника по кожному показнику = SQRT (Gxy. ^ 2 (Gxx - Gyy). ^ 2) EPS;
sin2theta = Gxy. / знаменника по кожному показнику;% синус і косинус кута подвоїв
cos2theta = (Gxx-Gyy). / знаменника по кожному показнику;

Цзе виправити = (6 * orientsmoothsigma), якщо ~ мода (Цзе, 2); Зі Зі = 1; кінець
F = fspecial ('гаусівської, Цзе, orientsmoothsigma);
cos2theta фільтр2 = (F, cos2theta);% згладжених синус і косинус
sin2theta фільтр2 = (F, sin2theta);% подвоїлася кутів

orientim = р / 2 atan2 (sin2theta, cos2theta) / 2;