遥感图像的辐射校正、辐射归一化

%辐射归一化 %输入参考图像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')