orientcollect.m

利用ＭＡＴＬＡＢ实现了对测量中的平面网

function [p1,p2,s,dir12,dir21]=orientcollect(e,d,sid,g,f,dir,si,ni,ed,dd,pn)
%实现平面网数据的整合，函数入口为(e,d,sid,g,f,dir,si,ni,ed,dd,pn)，分别表示：
%边长起点点号，终点点号，边长观测值、方向起点点号、终点点号、方向观测值、
%每个测站的第一个方向在总方向观测值中的位置、每个测站的方向数、已知点个数、未知点个数、点号
%函数出口为[p1,p2,s,dir12,dir21]，表示：整合后边长观测和往返方向观测的起点、终点、边长观测值、
%往方向观测值、返方向观测值
sd=ed+dd;
num=0;
for i=1:sd-1
    for j=i+1:sd
        [dist,kdis]=findsid(i,j,e,d,sid);       %判断该边长观测值是否存在并返回之
        [dir1,kdir1]=finddir(i,j,g,f,dir,si,ni);%判断该往测方向观测值是否存在并返回之
        [dir2,kdir2]=finddir(j,i,g,f,dir,si,ni);%判断该反测方向观测值是否存在并返回之
        if kdis+kdir1+kdir2==3
            num=num+1;
            p1(num)=i;
            p2(num)=j;
            s(num)=dist;
            dir12(num)=dir1;
            dir21(num)=dir2;
        end
    end
end
return