算法程序电子版.docx

一、一点校正程序程序如下：

%一点校正

%----------------------Part1：数据的读取和显示 %

fidlow=fopen(E:\di.dat,r); %打开文件fidhigh=fopen(E:\gao.dat,r);fidhand27=fopen(E:\yuan.dat,r);x1=fread(fidlow);

x2=fread(fidhigh);

x3=fread(fidhand27); %读取文件y1=reshape(x1,[200,200]);

y2=reshape(x2,[200,200]);

y3=reshape(x3,[200,200]); %矩阵转换200*200c1=uint8(y1);

c2=uint8(y2);

c3=uint8(y3); %转换为uint8型

fclose(all);

%----------------part2：数据的处理部分，一点校正，高温定标 %

C=0;

fori=1:1:200

forj=1:1:200

C=C+y2(i,j);

end

end

C1=C/40000; %对焦平面图像单元求平均

%对高温图像进行校正high_high=C1*E;

D=zeros(200); %生成200*200零矩阵，用以储存差值

fori=1:1:200

forj=1:1:200

D(i,j)=y2(i,j)-C1;

end

end

%对低温图像进行校正

fori=1:1:200

forj=1:1:200

Y10(i,j)=y1(i,j)-D(i,j);

end

end

%对于目标图像进行校正fori=1:1:200

forj=1:1:200

Y30(i,j)=y3(i,j)-D(i,j);

end

end

%图像生成

figure(1);subplot(2,3,1);

uint_high_high=uint8(high_high);imshow(uint_high_high);

title(一点校正高温图像（高温定标）);

subplot(2,3,2);a6=uint8(Y10);imshow(a6);

title(一点校正低温图像（高温定标）);

subplot(2,3,3);a7=uint8(Y30);imshow(a7);

title(一点校正后目标图像（高温定标）);

%三维响应图subplot(2,3,4);colormap(gray);

surf(high_high);

axis([0200,0200,0255]);

title(高温图像三维响应图（高温定标）);

subplot(2,3,5);colormap(gray);surf(Y10);

axis([0200,0200,0255]);

%低温

%类型转换

%类型转换

%类型转换

%高温

title(低温图像三维响应图（高温定标）);

subplot(2,3,6);colormap(gray);

surf(Y30); %目标三维响应图

axis([0200,0200,0255]);

title(目标图像三维响应图（高温定标）);

%----------------Part3：数据的处理部分，一点校正，低温定标

------------------%

A=0;

fori=1:1:200

forj=1:1:200

A=A+y1(i,j);

end

end

A1=A/40000; %对焦平面图像单元求平均

E=zeros(200)+1;

low_low=A1*E; %对低温图像进行校正

B=zeros(200); %生成200*200零矩阵，用以储存差值

fori=1:1:200

forj=1:1:200

B(i,j)=y1(i,j)-A1;

end

end %求差值

fori=1:1:200

forj=1:1:200

Y2(i,j)=y2(i,j)-B(i,j);%对高温图像进行校正

end

end

fori=1:1:200

forj=1:1:200

Y3(i,j)=y3(i,j)-B(i,j);

end

end

%对于目标图像进行校正

a4=uint8(Y2);

figure(2);

%类型转换

subplot(2,3,1);

imshow(a4);

title(一点校正后高温图像(低温定标）);

uint_low_low=uint8(low_low); %类型转换subplot(2,3,2);

imshow(uint_low_low);

title(一点校正后低温图像（低温定标）);

a5=uint8(Y3);subplot(2,3,3);imshow(a5);

title(一点校正后目标