IA = rgb2gray (im2double (A));
c_diag = corrcoef(ia(1:1:1,1:end-1),ia(2:结束,2:结束));
c_vert = corrcoef(IA(1:end-1,:), IA(2:end,:));
c_horz = corrcoef(ia(:,1:结束-1,:),ia(:,2:结束));
数字;
子图(3,3,1)
plot(IA(1:end-1, 1:end-1), IA(2:end, 2:end),'b'),title(sprintf('%s Horizontal original image',b));
子图(3,3,2)
图(1
子图(3,3,3)
绘图(IA(:,1:结束-1,:),ia(:,2:结束),'b'),标题(Sprintf('%s对角线映像',b));
IAA = rgb2gray (im2double (im12));
C_Diage = CorrCoef(IAA(1:1:1,1:结束-1),IAA(2:结束,2:结束)))
c_verte = corrcoef(iaa(1:end-1,:),iaa(2:结束,:))
c_horze = corrcoef (IAA (: 1: end-1,:), IAA(:, 2:结束)
子图(3,3,4)
图(IAA(1:结束-1,1:结束-1),IAA(2:结束,2:结束),'r'),标题(Sprintf('%s水平挤压图像',b));
子图(3,3,5)
绘图(IAA(1:1 :),IA(2:结束,:),'R'),标题(Sprintf('%s垂直呼气图像',b));
子图(3,3,6)
绘图(IAA(:,1:结束-1,:),ia(:,2:结束),'r'),标题(Sprintf('%s对角暗图像',b));