%讀入兩幅圖像
imgA=imread（‘A.jpg‘）;
imgB=imread（‘B.jpg‘）;
x=rgb2gray（imgA）;
y=rgb2gray（imgB）;

%精度轉換
imgA=double（x）;
imgB=double（y）;

%健壯性判斷
%%[imgArowimgAcol]=size（imgA）;
%%[imgBrowimgBcol]=size（imgB）;
%%[imgArowimgAcol]
%%[imgBrowimgBcol]
%%if?imgArow<1||imgAcol<1||imgBrow