时间:2023-10-09 09:14:31 下载该word文档
相对定向—绝对定向解法实验报告1、实验代码1.1根据所给同名像点的像平面坐标进行相对定向,求解相对的相对定向元素functionxP=xiangduidingxiang(Lxy,Rxy,f%UNTITLEDSummaryofthisfunctiongoeshere%Detailedexplanationgoeshere%设置相对定向元素是初始值u=0;v=0;w=0;q=0;k=0;bu=Rxy(1,1-Lxy(1,1;while(1%求解余弦元素a1=cos(q*cos(k-sin(q*sin(w*sin(k;a2=-cos(q*sin(k-sin(q*sin(w*cos(k;a3=-sin(q*cos(w;b1=cos(w*sin(k;b2=cos(w*cos(k;b3=-sin(w;c1=sin(q*cos(k+cos(q*sin(w*sin(k;c2=-sin(q*sin(k+cos(q*sin(w*cos(k;c3=cos(q*cos(w;R=[a1,a2,a3;b1,b2,b3;c1,c2,c3];[n,m]=size(Lxy;u2=[];v2=[];w2=[];fori=1:nu2(i=a1*Rxy(i,1+a2*Rxy(i,2-a3*f;v2(i=b1*Rxy(i,1+b2*Rxy(i,2-b3*f;w2(i=c1*Rxy(i,1+c2*Rxy(i,2-c3*f;endfori=1:nu1(i=Lxy(i,1;v1(i=Lxy(i,2;
w1(i=-f;endbv=bu*u;bw=bu*v;fori=1:nN1(i=(bu*w2(i-bw*u2(i/(u1(i*w2(i-u2(i*w1(i;N2(i=(bu*w1(i-bw*u1(i/(u1(i*w2(i-u2(i*w1(i;endfori=1:n