Rosenbrock 함수의 최소점을 찾는 예제를 돌려 봄.
% rosenbrock.m
function f = rosenbrock(x)
f = 100*(x(2) - x(1)^2)^2 + (1 - x(1))^2;
% unitdisk.m
function [c, ceq] = unitdisk(x)
c = x(1)^2 + x(2)^2 - 1;
ceq = [ ];
Below is the ICP registration pseudo code.
read model image
fill distance volume D (--> Jacobian 계산 가능할듯)
read input image --> make 3D vertices
Objective function
function e = objective(h)
for all image (3D) points,
(x y w)T = H(x y 1)T
e += D(x,y,z) (interpolated)