s050.m

来自「简单的数字图像处理程序」· M 代码 · 共 39 行

% 真彩图像分段线性变换  s050
clear;
I50=imread('1.jpg');
I500=I50(:,:,1);
I501=I50(:,:,2);
I502=I50(:,:,3);
I503=im2double(I500);
I504=im2double(I501);
I505=im2double(I502);
[m,n]=size(I500);
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 I503(i,j)<KL
     I506(i,j)=k1*I503(i,j);
     I507(i,j)=k1*I504(i,j);
     I508(i,j)=k1*I505(i,j);
     elseif I503(i,j)>KH
     I506(i,j)=(1-k3*(1-KH))+k3*(I503(i,j)-KH);
     I507(i,j)=(1-k3*(1-KH))+k3*(I504(i,j)-KH);
     I508(i,j)=(1-k3*(1-KH))+k3*(I505(i,j)-KH);
     else
     I506(i,j)=k1*KL+k2*(I503(i,j)-KL);
     I507(i,j)=k1*KL+k2*(I504(i,j)-KL);
     I508(i,j)=k1*KL+k2*(I505(i,j)-KL);
     end
  end
end
I509(:,:,1)=im2uint8(I506);
I509(:,:,2)=im2uint8(I507);
I509(:,:,3)=im2uint8(I508);
imshow(I509)