RTKPOST导出的pos文件转tum格式

clear;
clc;
data = importdata('street_02.pos');
outfile = 'street_02.tum';
data = data(20:end,:);
addpath('FUNCTIONS');

datamatrix = zeros(size(data,1),9);
for i = 1:size(data,1)
    b = data(i,1);
    c = strsplit(b{1,1},{'/',' ',':'});
    datamatrix(i,1)=str2double(c{1,1});
    datamatrix(i,2)=str2double(c{1,2});
    datamatrix(i,3)=str2double(c{1,3});
    datamatrix(i,4)=str2double(c{1,4});
    datamatrix(i,5)=str2double(c{1,5});
    datamatrix(i,6)=str2double(c{1,6});
    datamatrix(i,7)=str2double(c{1,7});
    datamatrix(i,8)=str2double(c{1,8});
    datamatrix(i,9)=str2double(c{1,9});
end


row = size(datamatrix,1);
tum = zeros(row,8);
for i=1:row
    CT = datamatrix(i:i,1:6);
    Mjd = CommonTimeToMjdtime(CT);
    GT = MjdTimeToGPSTime(Mjd);
    GT2second = GT(1)*7*24*60*60+GT(2);
    tum(i,1) = (GT2second  +  315964785-2.92);
    tum(i,2) = datamatrix(i,7);
    tum(i,3) = datamatrix(i,8);
    tum(i,4) = datamatrix(i,9);
    tum(i,8) = 1;
end
matTotum(tum, outfile);
clear;