function?[headobs]=ReadObsData
%讀接收機觀測數據文件
%HeadODat?:a?structure?stores?header?information?if?o-file
%???????????.ApproXYZ[3];??????????//approximate?coordinate
%???????????.interval;??????????//intervals?of?two?adjacent?epochs
%???????????.SiteName;??????????//point?name
%???????????.Ant_H;?????????????//antenna?height
%???????????.Ant_E;?????????????//east?offset?of?the?antenna?center
%???????????.Ant_N;?????????????//north?offset?of?then?antenna?center
%???????????.TimeOB;????????????//first?epoch?time?in?modifuied?Julian?time
%???????????.TimeOE;????????????//last?epoch?time?in?modifuied?Julian?time
%???????????.SumOType;??????????//number?of?types?of?observables
%???????????.SumOO[SumOType];??//type?of?observables?0-L11-L22-C13-P14-P25-D16-D2

%ObsODat?:a?structure?stores?observables?by?one?and?one?epoch
%???????????.TimeOEpp[2];???????//reciever?time?of?epoch
%???????????.SatSum;????????????//number?of?satellites
%???????????.SatCode[SatSum];???//satellites‘?PRN
%???????????.Obs_FreL1[SatSum];?
%???????????.Obs_FreL2[SatSum];
%???????????.Obs_RangeC1[SatSum];
%???????????.Obs_RangeP1[SatSum];
%???????????.Obs_RangeP2[SatSum];
global?HeadODat;
global?ObsODat;

????????[fnamefpath]=uigetfile（‘*.*‘‘選擇一個O文件‘）;
????????O_filename=strcat（fpathfname）;
???
????fid1=fopen（O_filename‘rt‘）;
????if?（fid1==-1）
????????msgbox（‘file?invalide‘‘warning‘‘warn‘）;
????????return;
????end
????%將文件頭數據存入結構體HeadODat中
????t=0;
????while（t<100）
????????s=fgets（fid1）;
????????t=t+1;
????????L=size（s2）;
????????if?L<81
????????????s（L+1:81）=‘?‘;
????????end
??????????
????????instrS=s（61:81）;
????????%測站點近似坐標
????????if?strncmp（instrS‘APPROX?POSITION?XYZ‘19）
????????????HeadODat.ApproXYZ=zeros（13）;
????????????HeadODat.ApproXYZ（11）=str2num（s（1:14））;
????????????HeadODat.ApproXYZ（12）=str2num（s（15:28））;
????????????HeadODat.ApproXYZ（13）=str2num（s（29:42））;
????????%歷元間隔;
????????elseif?strncmp（instrS‘INTERVAL‘8）
????????????HeadODat.interval=str2num（s（5:11））;
????????%測站點號????
????????elseif?strncmp（instrS‘MARKER?NAME‘11）
????????????HeadODat.SiteName=s（1:4）
????????%天線中心改正????
????????elseif?strncmp（instrS‘ANTENNA:?DELTA?H/E/N‘20）
????????????HeadODat.Ant_H=str2num（s（1:14））;
????????????HeadODat.Ant_E=str2num（s（15:28））;
????????????HeadODat.Ant_N=str2num（s（29:42））;
????????%第一個歷元時間????
????????elseif?strncmp（instrS‘TIME?OF?FIRST?OBS‘17）
????????????year=str2num（s（1:6））;
????????????month=str2num（s（7:12））;
????????????day=str2num（s（13:18））;
????????????hour=str2num（s（19:24））;
????????????minute=str2num（s（25:30））;
????????????second=str2num（s（31:42））;
????????????HeadODat.TimeOB=TimetoJD（yearmonthdayhourminutesecond）;
????????%最后一個歷元時間????
????????elseif?strncmp（instrS‘TIME?OF?LAST?OBS‘16）
????????????year=str2num（s（1:6））;
????????????month=str2num（s（7:12））;
????????????day=str2num（s（13:18））;
????????????hour=str2num（s（19:24））;
????????????minute=str2num（s（25: