双像空间前方交会（点投影系数法

clcclear;%清理命令窗口和工作区间
format long g;%规定输出数据为长整型
x=[17.383;-4.706]/1000;y=[8.396;6.444]/1000;%代入左右像片的同名像点坐标
f=0.1014;%主距f
x0=-0.180/1000;y0=0;%内方位元素
Xs=[50766.27048297550;51406.58175310970];
Ys=[85092.60204012240;85098.60787072870];
Zs=[2764.51598579742;2763.02759736577];%代入左右两张像片外方位元素中的线元素值
phi=[0.00749335445825581;-0.01158295898298480];
omi=[-0.00699425949527119;0.01003605897533120];
kap=[-0.00213195906904753;0.00849024086417501];%代入左右两张像片外方位元素中的角元素值
for i=1:2
a1（i）=cos（phi（i））*cos（kap（i））-sin（phi（i））*sin（omi（i））*sin（kap（i））;
a2（i）=-cos（phi（i））*sin（kap（i））-sin（phi（i））*sin（omi（i））*cos（kap（i））;
a3（i）=-sin（phi（i））*cos（omi（i））;
b1（i）=cos（omi（i））*sin（kap（i））;b2（i）=cos（omi（i））*cos（kap（i））;
b3（i）=-sin（o