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