鱼眼镜头校正matlab程序

% 功能：标定椭圆参数 % clc;close all;clear all; a = 400; %b = 160; Hnine = 1; data = imread('C:\Documents and Settings\Administrator\桌面\a\90.bmp'); [M N ] = size(data); %data(287:289,359:361) = 255; data(288,:) = 255; data(:,360) = 255; figure;imshow(data); test1 = 0; test2 = 0; test3 = 0; test4 = 0; rec1 = 0; rec2 = 0; data_then = zeros(M,N); datatest = zeros(M,N); ftest = zeros(M,N); ftest = data; x = 0; y = 0; %figure; for i = 72 for j = 1:N x_pie = j-N/2; % 矩阵坐标转论文（图像中点即原点）坐标 y_pie = M/2-i; %datatest(i,j) = (x_pie^2 + y_pie^2)/(210*210+145*145)*255; b = a*(y_pie/x_pie); if(x_pie<0) a = -abs(a); end if(y_pie<0) b = -abs(b); end end % x = a*(1 - y_pie*y_pie/b/b).^0.5; % y = b*(1 - x_pie*x_pie/a/a).^0.5; bb = 10; aa = 308; tuoyuan = @(xx) (bb*(1-xx.*xx/aa/aa).^0.5); x = x_pie; if(x_pie>aa | x_pie<-aa) x_pie = aa; end y = tuoyuan(x_pie); % x = x_pie/; % y = y_pie/2; if(x<-360 | x>360) rec1 = rec1+1; end if(y<-288 | y>288) rec2 = rec2+1; end i_index = M/2 - y; % 论文坐标转矩阵坐标 j_index = x + N/2 ; i_index = uint32(i_index); j_index = uint32(j_index); if(i_index>M) i_index = M; test1 = test1 + 1; elseif(i_index<1) i_index = 1; test2 = test2 + 1; end if(j_index>N) j_index = N; test3 = test3 + 1; else if(j_index<1) j_index = 1; test4 = test4 +1; end ftest(i_index,j_index) = 0; ftest(i,j) = 0; imshow(uint8(ftest)); data_then(i,j) = data(i_index,j_index); end end %figure;imshow(uint8(datatest),[]); figure;imshow(uint8(data_then)); % % %% -----------------透镜形变校正------------------------ % A = [ % 13,59,-1,0,0,0,585,2655 % 0,0,0,13,59,-1,338,1534 % -68,58,-1,0,0,0,7344,-6264 % 0,0,0,68,-58,1,1496,-1276 % 51,124,-1,0,0,0,2601,6324 % 0,0,0,51,124,-1,6324,15376 % -113,122,-1,0,0,0,12769,-13786 % 0,0,0,113,-122,1,13786,-14884 % %-30,106,-1,0,0,0,900,-3180,30 % % ]; % B = [45 26 -108 -22 51 -124 -113 -122 ]'; % % value = det(A); % Ali = inv(A); % Htmp = Ali*B; % % H = zeros(3,3); % for i = 1:8 % H(i) = Htmp(i)*Hnine; % end % H(9) = Hnine; % % H = H'; % % h = inv(H); % % % data_then = zeros(M,N,NN); % test1 = 0; % test2 = 0; % test3 = 0; % test4 = 0; % for i = 1:M % for j = 1:N % X = j-N/2; % 矩阵坐标转论文（图像中点即原点）坐标 % Y = M/2-i; % % x_pie = (h(1,1)*X + h(1,2)*Y + h(1,3))/(h(3,1)*X + h(3,2)*Y + h(3,3)); % y_pie = (h(2,1)*X + h(2,2)*Y + h(2,3))/(h(3,1)*X + h(3,2)*Y + h(3,3)); % % i_index = M/2 - y_pie; % 论文坐标转矩阵坐标 % j_index = x_pie + N/2; % % i_index = uint32(i_index); % j_index = uint32(j_index); % if(i_index>M) % i_index = M; % test1 = test1 + 1 % elseif(i_index<1) % i_index = 1; % test2 = test2 + 1 % end % % if(j_index>N) % j_index = N; % test3 = test3 + 1 % elseif(j_index<1) % j_index = 1; % test4 = test4 +1 % end % % % data_then(i,j,:) = data(i_index,j_index,:); % end % end % % figure;imshow(uint8(data_then));