s049.m

简单的数字图像处理程序

% 灰度图像分段线性变换  s049
clear;
I49=imread('1.jpg');
I490=rgb2gray(I49);
I491=im2double(I490);
[m,n]=size(I490);
L=1;
KL=0.2*L;
KH=0.8*L; 
k1=0.5;
k3=0.5;
k2=(1-k3*(1-KH)-k1*KL)/(KH-KL);
for i=1:m
  for j=1:n
     if I491(i,j)<KL
     I492(i,j)=k1*I491(i,j);
     elseif I491(i,j)>KH
     I492(i,j)=(1-k3*(1-KH))+k3*(I491(i,j)-KH);
     else
     I492(i,j)=k1*KL+k2*(I491(i,j)-KL);
     end
  end
end
I493=im2uint8(I492);
imshow(I493)