单向空间后方交会matlab编程(共3页).doc

精选优质文档-倾情为你奉上%读取包含点号的像点坐标、控制点坐标分别到矩阵IP、CP%IP=load(f:imgpoint.txt);CP=load(f:ctlpoint.txt);%删除像点坐标、控制点坐标矩阵中的点号列%IP(:,1)=;CP(:,1)=;%像片大小%H=4008;W=5344;IP1=IP(:,1)-W/2;IP2=H/2-IP(:,2);IP=IP1,IP2;%焦距%fx=4547.99665;fy=4547.87373;%内方位元素%x0=47.48571;y0=12.02756;%畸变参数%k1=-5.00793e-009;k2=1.83462e-016;p1=-2.24419e-008;p2=1.76820e-008;%外方位元素初值%Xs=3000;Ys=-100;Zs=100;Phi=0;Omega=0;Kappa=0;m=size(IP,1);AIP=zeros(m,2);L=zeros(2*m