cover = imread('s3898.pgm'); %原始图像路径
stego = imread('ad3898.pgm');%修改图像像素后的图像路径
figure;
subplot(1, 2, 1); imshow(cover); title('cover');
subplot(1, 2, 2); imshow((double(stego) - double(cover) + 1)/2); title('embedding changes: +1 = white, -1 = black');
%fprintf('\n\nImage embedded , change rate: %.4f, distortion per pixel: %.6f\n',sum(cover(:)~=stego(:))/numel(cover));
fprintf('change bits %d \n', sum(cover(:)~=stego(:)))
输出对比
change bits 1991