I = imread('cameraman.tif');
J = uint8(filter2(fspecial('gaussian'), I));
K = imabsdiff(I,J);
%数据相减后取绝对值
%或者用指令K = imsubtract(I,J);数据相减后,小于0的用0处理。
%类似矩阵的减法运算,但特别注意类型的处理。
imshow(K,[])
I = imread('moon.tif');
I16 = uint16(I);
J = immultiply(I16,I16);
%类似矩阵点乘运算,但注意类型的处理
subplot(121),imshow(I)
subplot(122),imshow(J)
I = imread('cameraman.tif');
J = uint8(filter2(fspecial('gaussian'), I));
K = imlincomb(1,I,-1,J,128);
% 运算关系式K(r,c) = I(r,c) - J(r,c) + 128
imshow(K)
MATLAB图像处理—孙兆林编著—清华大学出版社