Extracting the daily time series from a text file to a single column

Question

Shahram Sahraei 2019-7-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/473118-extracting-the-daily-time-series-from-a-text-file-to-a-single-column

评论： Shahram Sahraei 2019-7-24

D59 - Copy.txt

Dear Matlab users,

I have a text file (attached here) containing daily time series from 2013 to 2019 and I want to take out the daily data points and save it in a column vector in chronological order. I would be really grateful of your help in this regard.

Best,

Shahram

2 个评论
显示无隐藏无

Guillaume 2019-7-23

It would be a lot more useful if you attached a representative text file.

Shahram Sahraei 2019-7-23

Dear Guillaume,

Thanks for you suggestion. I just modified my question and attached the text file to it.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

David Wilson 2019-7-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/473118-extracting-the-daily-time-series-from-a-text-file-to-a-single-column#answer_384481

A quick hack is:

fname = 'D59.txt' 
Flow = []; 
fid=fopen(fname);
Yr = 2012; 
R = []; 
while 1
    tline = fgetl(fid);
    if ~ischar(tline), break, end
       
    s = strtrim(tline); 
    disp(s)
    if strncmp(s,'Day',3)
        disp(s)
        Yr = Yr+1; 
        X = []; 
        i=1; 
        while i < 32
            s = deblank(fgetl(fid));
            if isempty(s) 
                s = deblank(fgetl(fid));
            end 
            %disp(s);
            x = str2num(s); disp(x) 
            n = length(x); 
            if n == 8
                x = [x(1:2), NaN, x(3), NaN, x(4), NaN, x(5:6), NaN, x(7), NaN x(8)]; 
            elseif n == 12
                x = [x(1:2), NaN, x(3:end)];
            end 
            X = [X;x]; 
            i=i+1; 
        end % for 
        Rdat = X(:,[2:end]); 
        r = Rdat(:); 
        r(find(isnan(r)))=[]; % drop bad days (leap year?) 
        R = [R; r]; 
    end 
    
end
fclose(fid);
plot(R) 