Determinant of Jacobian of a displacement field

最新推荐文章于 2021-11-29 17:13:46 发布

juliosun

最新推荐文章于 2021-11-29 17:13:46 发布

阅读量556

点赞数

分类专栏：学习笔记

学习笔记专栏收录该内容

56 篇文章

订阅专栏

本文介绍了一个计算位移场雅可比行列式的MATLAB函数。该函数通过计算位移场的梯度并加上单位矩阵来构建雅可比矩阵，最后计算其行列式。此方法适用于图像配准等领域。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

% J = jacobian(S)
%   Determinant of Jacobian of a displacement field
%
% Herve Lombaert, Jan. 8th, 2013
%
function det_J = jacobian(sx,sy,sz)

    % Gradients
    [gx_y,gx_x,gx_z] = gradient(sx);
    [gy_y,gy_x,gy_z] = gradient(sy);
    [gz_y,gz_x,gz_z] = gradient(sz);
    
    % Add identity
    gx_x = gx_x + 1;
    gy_y = gy_y + 1;
    gz_z = gz_z + 1;
    
    % Determinant
    det_J = gx_x.*gy_y.*gz_z + ...
            gy_x.*gz_y.*gx_z + ...
            gz_x.*gx_y.*gy_z - ...
            gz_x.*gy_y.*gx_z - ...
            gy_x.*gx_y.*gz_z - ...
            gx_x.*gz_y.*gy_z;
end