两点校正高温或低温图片--红外图像校正

Matlab程序实现，图像处理，红外图像处理，两点校正高温或低温图片--红外图像校正%两点校正程序文件clear;clc;a1=fopen('handdat_34.dat','rb','b');%打开二进制数据信息pic1=fread(a1,[200200]);%将二进制数据信息读出到200*200的多维数组里面a2=fopen('hidat_62.dat','rb','b');pic2=fread(a2,[200200]);a3=fopen('low31dat_62.dat','rb','b');pic3=fread(a3,[200200]);%求解阵列焦平面高温定标点相应均值A1A1=mean(pic2(:));%求pic2的平均值将结果返回给b%A1=b1*ones(200);%K1=A1./pic2;%高温定点校正因子200*200矩阵数组K%求解阵列焦平面低温定标点相应均值A2A2=mean(pic3(:));%求pic2的平均值将结果返回给b%A2=b2*ones(200);%K2=A2./pic3;%计算校正因子200*200矩阵数组KGij=(A1-A2)./(pic2-pic3);%求增益校正系数GijOij=(A1.*pic3-A2.*pic2)./(pic2-pic3);%求偏置校正系数Oij%两点校正高温图片Yij_hi=Gij.*pic2-Oij;%两点校正低温图片Yij_low=Gij.*pic3-Oij;%两点校正目标成像图片Yij_handdat=Gij.*pic1-Oij;figure(1);subplot(131),imshow(Yij_hi,[]),title('两点校正高温图片');subplot(132),imshow(Yij_low,[]),title('两点校正低温图片');subplot(133),imshow(uint8(Yij_handdat),[]),title('两点校正目标成像图片');figure(2);subplot(131),mesh(1:200,1:200,roundn(Yij_hi,-3)),title('两点校正高温图片三维图');subplot(132),mesh(1:200,1:200,roundn(Yij_low,-3)),title('两点校正低温图片三维图');subplot(133),mesh(1:200,1:200,Yij_handdat),title('两点校正目标成像图片三维图');colormapgray;%选为灰度%colormapgray;%选为灰度Yij_handdat=uint8(Yij_handdat);save('Yij_handdat','Yij_handdat')