MATLAB图像处理：18 SAR图像和光学图像的配准算法.zip

function Flag = RSRADFilter(P_Class,p); SIZE = size(P_Class); SIZEX = SIZE(1,1); SIZEY = SIZE(1,2); landa = 8.8;rou = 2; n = 12;q0 = 0.33333;%n = 20 P_Class0 = P_Class; %初始化数组 I_gradR = zeros(SIZEX,SIZEY,2);I_gradL = zeros(SIZEX,SIZEY,2); P_Class2 = zeros(SIZEX+2,SIZEY+2); Cr2 = zeros(SIZEX+2,SIZEY+2); % for i = 1:SIZEX % for j = 1:SIZEY % if P_Class(i,j) == 0; % P_Class(i,j) = 1; % end % end % end %求图像I的R方向和L方向上的梯度 for time = 1:n str = '已经完成'; str = strcat(str,num2str(time/n*100.0),'%'); disp(str); q0_t(time) = q0*exp(-1*rou*time*0.05); % Delta_T = 5^0.5*q0_t(time)^2; %设定图像边界条件 P_Class2(2:SIZEX+1,2:SIZEY+1) = P_Class; P_Class2(2:SIZEX+1,1) = P_Class(:,1);P_Class2(1,2:SIZEY+1) = P_Class(1,:); P_Class2(2:SIZEX+1,SIZEY+2) = P_Class(:,SIZEY);P_Class2(SIZEX+2,2:SIZEY+1) = P_Class(SIZEX,:); for i = 1:SIZEX for j = 1:SIZEY I_gradR(i,j,:) = [P_Class2(i+1+1,j+1)-P_Class2(i+1,j+1),P_Class2(i+1,j+1+1)-P_Class2(i+1,j+1)]; I_gradL(i,j,:) = [P_Class2(i+1,j+1)-P_Class2(i+1-1,j+1),P_Class2(i+1,j+1)-P_Class2(i+1,j-1+1)]; module_grad = I_gradR(i,j,1)^2+I_gradR(i,j,2)^2+I_gradL(i,j,1)^2+I_gradL(i,j,2)^2; %计算Laplacian算子 I_Laplace(i,j) = P_Class2(i+1+1,j+1)+P_Class2(i+1,j+1+1)+P_Class2(i-1+1,j+1)+P_Class2(i+1,j-1+1)-4*P_Class2(i+1,j+1); %计算扩散系数 q(i,j) = ((1/4*module_grad/(P_Class2(i+1,j+1)+0.0001)^2-1/16*(I_Laplace(i,j)/(P_Class2(i+1,j+1)+0.0001))^2)/... (1+1/4*I_Laplace(i,j)/(P_Class2(i+1,j+1)+0.0001)+0.0001)^2)^0.5; Cr(i,j) = p*1/(2*(1+(q(i,j)-q0_t(time))^2)^0.5);%2.53 % Cr(i,j) = 0.36*exp(log2(1+(q(i,j)-q0_t(time))^2))/((1+(q(i,j)-q0_t(time))^2)*0.3010); % Cr(i,j) = 2.36*1/((2*3.14)^0.5*0.7)*exp(-1*(q(i,j)-q0_t(time)).^2/(2*0.7^2));% 1.88 % Cr(i,j) = 1.8*1/((2*3.14)^0.66*0.4)*exp(-1*10*(q(i,j)-q0_t(time)).^2/(2*0.66^2)); % Cr(i,j) = 2*1/((2*3.14)^1.66*0.4)*exp(-1*(q(i,j)-q0_t(time)).^2/(2*1.66^2)); end end %设定Cr的边界条件 Cr2(2:SIZEX+1,2:SIZEY+1) = Cr; Cr2(2:SIZEX+1,1) = Cr(:,1);Cr2(1,2:SIZEY+1) = Cr(1,:); Cr2(2:SIZEX+1,SIZEY+2) = Cr(:,SIZEY);Cr2(SIZEX+2,2:SIZEY+1) = Cr(SIZEX,:); %计算散度 for i = 1:SIZEX for j = 1:SIZEY DIV(i,j) = Cr2(i+1+1,j+1)*(P_Class2(i+1+1,j+1)-P_Class2(i+1,j+1))+Cr2(i+1,j+1)*(P_Class2(i-1+1,j+1)-... P_Class2(i+1,j+1))+Cr2(i+1,j+1+1)*(P_Class2(i+1,j+1+1)-P_Class2(i+1,j+1))+Cr2(i+1,j+1)*(P_Class2(i+1,j-1+1)-P_Class2(i+1,j+1)); % end % end %开始迭代 % for i = 1:SIZEX % for j = 1:SIZEY Img_filtered(i,j) = P_Class(i,j)+0.05*landa/4*DIV(i,j)+0.05*(P_Class0(i,j)-... P_Class(i,j)); end end P_Class = Img_filtered; end Flag = P_Class;