GPS RINEX格式读取MATLAB版

%========================================
%     Cycle-slip
%     Zhen Dai
%========================================
% Functions:
%      Read  the RINEX observation file and save the measurement
% Input parameters:
%      filename --> RINEX observation file name
%      sampling_interval --> Sampling interval in Seconds （0.1 s == 10 Hz ......）
%      totalepoch -> total epochs to be processed. -999 implies processing
%               all the data in this file
% Output:
%      mSatID -> Satellite visibility of each epoch
%      vTime-> Time of each epoch scaled to “second of week“
%      vTimeStr -> Time expressed in STRING
%      interval -> Sampling interval of the data
%      mCode -> C/A code on  L1 signal
%      mPhase -> Carrier phase data on L1 signal
%  Remarks:
%      1. Only read the data from the GPS satellites
%      2. Only consider the L1 C/A and L1 phase data
%      3. Outputs are saved in a data file
%      4. It may provide some unexpected result when processing some
%      specially formatted RINEX observation file.

function [datafilenametotalepoch_real]=GUIMain_LoadRinexOBS（filenamesampling_intervaltotalepoch）
% function [datafilenametotalepoch_real]=GUIMain_LoadRinexOBS（）
%% GUI Initialization
f1=figure;
close（f1）;
figure;
title=sprintf（‘Reading RINEX observation file‘）;
set（gcf‘MenuBar‘‘none‘‘Numbertitle‘‘off‘‘Name‘title...
    ‘Position‘[400 500 450 150]‘Visible‘‘off‘）;
movegui（gcf‘center‘）;
hlistbox = uicontrol（gcf‘style‘ ‘listbox‘ ‘String‘ ‘Clear‘...
    ‘Position‘ [0 0 450 150]‘String‘‘‘‘Foregroundcolor‘‘k‘‘Backgroundcolor‘‘w‘）;
set（gcf‘Visible‘‘on‘）;

datafilename=sprintf（‘RawDataTemp‘）;
 

%% Start reading the RINEX file
totalGNSSsatellite=32;
vSatID=zeros（totalGNSSsatellite1）;
interval=0;
%% Initialize the data matrix.
if totalepoch==-999 %% When the total number of epochs are unknown
    vTime=[];  mCode=[]; mPhase=[]; vTimeStr=[];mDoppler1=[];mDoppler2=[];
elseif   totalepoch>0%% With the known number of epochs
    vTime=zeros（totalepoch1）;
    vTimeStr=cell（totalepoch1）;
    mP1=zeros（totalGNSSsatellitetotalepoch）;
    mC1=zeros（totalGNSSsatellitetotalepoch）;
    mL1=zeros（totalGNSSsatellitetotalepoch）;
    mP2=zeros（totalGNSSsatellitetotalepoch）;
    mL2=zeros（totalGNSSsatellitetotalepoch）;
    mD1=zeros（totalGNSSsatellitetotalepoch）;
    mD2=zeros（totalGNSSsatellitetotalepoch）;
else
    errordlg（‘Invalid total number of epochs‘）
end

%% Open the file
fid = fopen（filename）;
if fid==-1
    errmsg=sprintf（‘Can not open this file%s‘filename）;
    errordlg（errmsg）;
    return;
end

line80num=0;

obsC1pos=0; %% Indicates which observation is C/A code
obsL1pos=0; %% Indicates which observation is phase on L1
obsP1pos=0; %% Indicates which observation is P1 code
obsL2pos=0; %% Indicates which observation is phase on L2
obsP2pos=0; %% Indicates which observation is P2 code
obsD1pos=0; %% Indicates which obs