%辐射归一化
%输入参考图像A和待配准图像B
A=imread('G:\ENVI遥感图像处理方法数据盘\sub3处理\源图像\05-1.tif');
B=imread('G:\ENVI遥感图像处理方法数据盘\sub3处理\源图像\02-1.tif');
%分别求出A,B的均值和标准差
a=double(A);
A1=mean2(a)
A2=std2(a)
b=double(B);
B1=mean2(b)
B2=std2(b)
%辐射归一化
C=(b-B1)*A2/B2+A1;
%剔除不符合的点
s=size(C);
for i=1:s(1)
for j=1:s(2)
if C(i,j)<0
C(i,j)=0;
end
if C(i,j)>255
C(i,j)=255;
end
end
end
%像素四舍五入取整
D=round(C);
%求配准后的均值和方差
C1=mean2(D)
C2=std2(D)
%转化为单精度
E=uint8(D);
%显示输入图像的第一波段
figure;
imshow(A(:,:,1,1))
figure;
imshow(B(:,:,1,1))
figure;
imshow(E(:,:,1,1))
%保存图像
imwrite(E,'02-RadiometricCorrection.tif')
