彩色平衡处理 matlab程序

% balance.m % 彩色平衡 clc; clear; close all; im=double(imread('h5.tif')); [m,n,p]=size(im); F1=im(1,1,:); F2=im(1,2,:); F1_(1,1,1)=F1(:,:,2); F1_(1,1,2)=F1(:,:,2); F1_(1,1,3)=F1(:,:,2); F2_(1,1,1)=F2(:,:,2); F2_(1,1,2)=F2(:,:,2); F2_(1,1,3)=F2(:,:,2); K1=(F1_(1,1,1)-F2_(1,1,1))/(F1(1,1,1)-F2(1,1,1)); K2=F1_(1,1,1)-K1*F1(1,1,1); L1=(F1_(1,1,3)-F2_(1,1,3))/(F1(1,1,3)-F2(1,1,3)); L2=F1_(1,1,3)-L1*F1(1,1,3); for i=1:m for j=1:n new(i,j,1)=K1*im(i,j,1)+K2; new(i,j,2)=im(i,j,2); new(i,j,3)=L1*im(i,j,3)+L2; end end im=uint8(im); new=uint8(new); figure; imshow(im); title('原始图'); figure; imshow(new); title('平衡后');