s048.m

简单的数字图像处理程序

% 真彩图像对数变换  s048
I48=imread('1.jpg');
I480=im2double(I48);
a0=I480(:,:,1);
a1=I480(:,:,2);
a2=I480(:,:,3);
[m,n]=size(I480(:,:,1));
q=1/log10(50);  
for i=1:m
    for j=1:n
       a1(i,j)=q*log10(1+49*I480(:,:,1));
       a2(i,j)=q*log10(1+49*I480(:,:,2));
       a3(i,j)=q*log10(1+49*I480(:,:,3));
    end
end
b1=im2uint8(a1);
b2=im2uint8(a2);
b3=im2uint8(a3);
imshow(cat(3,b1,b2,b3))