TP-probleme-inverse-3D/TP6/delta_z.m
2023-06-25 16:38:01 +02:00

11 lines
244 B
Matlab

function [d_z] = delta_z(I, voisins, distances)
I_max = max(I, [], "all");
abs_grad_z = sqrt( (I_max ./ I).^2 - 1 );
p1 = voisins(:, 1);
p2 = voisins(:, 2);
d_z = distances .* ( abs_grad_z(p1) + abs_grad_z(p2) ) / 2;
end