5－1－1－2 分段线性变换的源程序：

%%%%   5－1－1－2 分段线性变换的源程序：
clear all;
X1=imread('E:\role0
 
 
%%%%   5－1－1－2 分段线性变换的源程序：
clear all;
X1=imread('E:\role0\001i.bmp');
figure,
imshow(X1)
f0=0;g0=0;
f1=20;g1=10;
f2=180;g2=230;
f3=255;g3=255;
figure,plot([f0,f1,f2,f3],[g0,g1,g2,g3])
axis tight,xlabel('f'),ylabel('g')
title('intensitytransformation')%绘制变换曲线
r1=(g1-g0)/(f1-f0);
b1=g0-r1*f0;
r2=(g2-g1)/(f2-f1);
b2=g1-r2*f1;
r3=(g3-g2)/(f3-f2);
b3=g2-r3*f2;
[m,n]=size(X1);
X2=double(X1);
for i=1:m
    for j=1:n
        f=X2(i,j);
        g(i,j)=0;
        if(f>=f1)&(f<=f2)
            g(i,j)=r1*f+b2;
        elseif(f>=f2)&(f<=f3)
            g(i,j)=r3*f+b3;
        end
    end
end
figure,imshow(mat2gray(g))
 
 
 
  
  

1i.bmp');
figure,
imshow(X1)
f0=0;g0=0;
f1=20;g1=10;
f2=180;g2=230;
f3=255;g3=255;
figure,plot([f0,f1,f2,f3],[g0,g1,g2,g3])
axis tight,xlabel('f'),ylabel('g')
title('intensitytransformation')%绘制变换曲线
r1=(g1-g0)/(f1-f0);
b1=g0-r1*f0;
r2=(g2-g1)/(f2-f1);
b2=g1-r2*f1;
r3=(g3-g2)/(f3-f2);
b3=g2-r3*f2;
[m,n]=size(X1);
X2=double(X1);
for i=1:m
    for j=1:n
        f=X2(i,j);
        g(i,j)=0;
        if(f>=f1)&(f<=f2)
            g(i,j)=r1*f+b2;
        elseif(f>=f2)&(f<=f3)
            g(i,j)=r3*f+b3;
        end
    end
end
figure,imshow(mat2gray(g))%%%%   5－1－1－2 分段线性变换的源程序：
clear all;
X1=imread