我想知道Identity映射Id的像素值；它通常用于图像配准。假设h是变形；那么，对于Id+h，他们称之为扭曲图像。但是，我不知道如何将Id嵌入到matlab中；它的价值是什么？它是什么恒等矩阵吗？最佳
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))
这是我所做的；我通过给w分配固定图像（I1）来初始化它。此代码旨在将I2注册到I1中。我不知道如何评估我的代码。另一方面，Id表示身份映射。你有人能分享他/她的想法吗，