红外图像--分段线性变换

红外图像，分段线性变换，matlab实现

%分段线性变换

clc;clear;

% X1=imread('liang_dian_jiaozheng.png');

load Yij_handdat.mat;

X1 = Yij_handdat;

% X1= rgb2gray(X1) ;

f0=0;g0=0;%对图像灰度进行分段的点

f1=123;g1=70;

f2=240;g2=150;

f3=256;g3=220;

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;

subplot(131),imshow(X1),title('两点校正后的图像');

subplot(132),imshow(g, []),title('分段线性变换后的图像');

subplot(133),mesh(roundn(g,-3)),title('分段线性变换后图像三维图');