Id : 表示 Rn 的恒等映射

Question

我想知道Identity mapping Id的像素值；它通常用于图像配准。说h是变形；然后，---对于 Id+h--- 他们称之为扭曲的图像。 但是，我不知道如何将 Id 嵌入到 matlab；它的价值是什么？它是单位矩阵吗？ 最好的

Answer 1

a=0.1;
stopcr=eps
I1=dicomread('1.dcm');
I2=dicomread('2.dcm');
maxiter=1000;
normdiv=inf;
while( (abs(b)>= stopcr) &(iter <=maxiter))

div1= (I1-I2).* uint16(imgradient(I2));

lap= 0.6 * (imgradient(imgradient(I2)));
lap=uint16(lap);
b=(div1+lap);
div= a* b;
div=uint16(div);
normdiv=abs(div);
solnew= w +div;

w=solnew;
end
w;
b
imcontrast(imshow(w))

这就是我所做的；我通过为其分配固定图像（I1）来初始化 w。此代码旨在将 I2 注册到 I1 中。 我不知道评估我的代码。另一方面，Id 表示身份映射。有人能分享一下他/她的想法吗，