フォルダ内のファイル​を「1日分」読み込み​、グラフを作るには？

以下のように一つのフォルダに1分ごとに作成されたファイルがあります。①

ファイル内のデータはテキストとバイナリデータが混在しています。② 簡単なデータ構造は下記に示します。③ データフォーマットはこちらです。

このフォルダ内のファイルを「1日分」読み込み、グラフを作成したいです。

下記のプログラムはエラーが出ます。③ また、下記プログラムはフォルダ内の全てのファイル読み込みになっているが、「1日分」指定して読み込む方法にしたいです。

以前の質問に補足しました。 一人で詰まってしまい大変困っています。よろしくお願い致します。

----①フォルダ内

　　ファイル名　　　　　　作成日時

e16A0319.272413 2016/10/03 19:26

e16A0319.272413 2016/10/03 19:27

・

・

----②一つのファイル内構造

e16A0319.272413

03/10/2016 19:26:23 03/10/2016 19:27:24 0000 0000.0 0000.0 00

0000000 0010 0018952 0010 04

1 0 2 00512 1 0000 0015 00829.o 0 0 00 000 12 150024 0.020 BT1

1 1 2 00512 1 0000 0015 00829.o 0 0 00 000 00 150024 3.1746 BC1

1 0 2 00512 1 0000 0015 00829.o 0 0 00 000 12 150023 0.020 BT1

1 1 2 00512 1 0000 0015 00829.o 0 0 00 000 00 150023 3.1746 BC1

邪ﾗ｡ﾗｼｩﾗOｫﾗそﾗa泰ﾗL釐!TﾘﾗﾖZﾗｮｪﾗ㌘ﾗtKﾗTxﾗ%rﾗ晗o蹟( ﾗ・ﾗｭﾗｿCﾗ 碩>槞遞ﾗM泰%ｬﾗｷﾗCｺﾗdｺﾗ・ﾗﾏﾗﾗﾗﾗ=・"ｮfﾘﾂﾘﾀﾘﾘVﾘN/ﾘｺ4ﾘ6DﾘDHﾘVﾘﾘDﾘﾋ_ﾘﾗaﾘjﾘYnﾘｦyﾘ|ﾘ>ﾘｦ筋X筋ﾄ娩壓s｣ﾘ｡｢ﾘgｩﾘｯﾘyｯﾘ┝ﾘｸﾁﾘｿﾘ捨ﾘ橆ﾘH籀 ﾟﾘﾁ袙､毖・ﾘﾊ."ﾙYﾙ-ﾙﾅﾙ ﾙｱ'ﾙI2ﾙﾄ6ﾙ・ﾙ・ﾙｶ@ﾙ'Jﾙ祟ﾙxVﾙ輯ﾙ(cﾙjdﾙﾅiﾙHlﾙjﾙoﾙｮtﾙﾓgﾙﾐpﾙnﾙ+uﾙ`pﾙ綠ﾙ2・ﾂEj\MｿMS・果_ｯ・Xｵ・T-適・A｣・s=FC>輅筋ﾌfﾍﾅｼk・u%dﾛ0ﾖа4ﾋｾ・ｦ ・・・・・・・・・・・・・・・・

-----③簡単なデータ構造

固定ヘッダ行1

固定ヘッダ行2

固定ヘッダ行3

データセットヘッダ行1

データセットヘッダ行2

データセットヘッダ行3

データセットヘッダ行4

データセット1 (512x4バイト)

CRLF

データセット2 (512x4バイト)

CRLF

データセット3 (512x4バイト)

CRLF

データセット4 (512x4バイト)

------④プログラム

DebugFiles = 'C:\Users\ken\Documents\MATLAB\161004';
%ファイル名をフルパスごと取得 　「161004」というフォルダ
D = dir(DebugFiles);
%フォルダ内のファイルの名前を全取得

for k = 1:length(D)
    % ファイルオープン
    fid = fopen('D(k).name');
    % ヘッダー読み込み
    for i = 1:3
        headers{i} = fgetl(fid);
    end

      % データセットの数を取り出し
      third_header = sscanf(headers{3},'%f');
      num_datasets = third_header(end);              %データセット数

      % データセットのヘッダーを読み込み
      datasetheader = {};

      for i = 1:num_datasets
          datasetheader{i} = fgetl(fid);
      end

      %fread(fid,2,'uint8') % Read CRLF=13d 10d

      % データセットのヘッダーに記載されたデータ分バイナリデータを読み込む
      data = {};
      for i = 1:num_datasets
          dataheader_parsed = sscanf(datasetheader{i},'%d');
          num_read = dataheader_parsed(4);                    %１つのデータセットのデータ数
          [data{i},cout] = fread(fid,num_read,'long');
          fread(fid,2,'uint8') % Read CRLF=13d 10d
      end
      % データファイルのクローズ
      fclose(fid);

      % 行列に整形(各データの数が同じ場合は成功)
      DataMat = [data{:}];
      plot(DataMat);
      hold on
  end

[michio]

関連質問： テキストとバイナリデータ混在のファイルを読み込むには？

[Jiro Doke]

具体的なエラー情報（エラーメッセージとどの行でエラーが起きるか、等）を追加してください。

参考： MATLAB Answersで早く的確な回答を得るためのポイント

Answer 1

[michio]

「1日分」をどこまで厳密に指定されるか次第ですが・・・少なくとも計測開始日はファイル名だけで絞り込めそうです。計測終了までとなると・・ファイル中身（２行目）を見て行く必要がありますね。

 03/10/2016 19:26:23 03/10/2016 19:27:24 0000 0000.0 0000.0 00

ファイル名で絞り込む場合は、いま

 D = dir(DebugFiles);

と実行している所を

 D = dir(DebugFiles, '*16A03*');

と変更すると、2016年10月3日に計測開始されたファイルだけが選ばれますので、その日のデータだけを読み取ることになる・・かと思いますが、いかがでしょう？

ファイルフォーマット を確認するとファイル名の命名規則が冒頭で説明されています。例えば e16A0319.272413 を例に取ると

(略) the filename, which begins with

• a letter -> e (意味はわかりませんが・・）
• two digits showing the years in the century -> 16 (2016年)
• the month (hexadecimal, one digit) -> A (16進数なので10月）
• the day (decimal, two digits) -> 03 (３日)
• the hour (decimal, 24 hours per day, two digits) -> 19 (19時）

といった感じです。

[nknknknk]

ご回答ありがとうございます。

一つのファイルに対する処理 （for k=3:length(D)～end） はお蔭様で出来ました。

161004フォルダのファイルを16A03指定で取得し、全てのファイルに処理をしようとすると以下のエラーが出ます。

『エラー fgets ファイルの識別子が無効です。有効なファイルの識別子を生成するには fopen を使用してください。

エラーの発生場所 fgetl (line 34) [tline,lt] = fgets(fid);

エラーの発生場所 totoko (line 22) headers{k} = fgetl(fid);』

問題は何が考えられるでしょうか。R2012bです。

 %ファイル名をフルパスごと取得
  DebugFiles = 'C:\Users\ken\Documents\MATLAB\161004';                          
  %指定フォルダ内のファイルをリスト化 日にち指定
  D = dir(DebugFiles,'*16A03*');                                                             
  %最初の2つは'.'と'..'だから、k=3からスタート
  for k=3:length(D)
      % ファイルオープン
      fid=fopen(D(k).name,'r');
      % ヘッダー読み込み
      for k = 1:3
          headers{k} = fgetl(fid);
      end
      % データセットの数を取り出し
      third_header = sscanf(headers{3},'%f');
      num_datasets = third_header(end);
      % データセットのヘッダーを読み込み
      datasetheader = {};
      for k = 1:num_datasets
          datasetheader{k} = fgetl(fid);
      end
      fread(fid,2,'uint8') % Read CRLF=13d 10d
      % データセットのヘッダーに記載されたデータ分バイナリデータを読み込む
      data = {};
      for k = 1:num_datasets 
          dataheader_parsed = sscanf(datasetheader{k},'%d');
          num_read = dataheader_parsed(4); % 
          [data{k},cunt] = fread(fid,num_read,'long');
          fread(fid,2,'uint8'); % Read CRLF=13d 10d
      end
      % データファイルのクローズ
      fclose(fid);
      % 行列に整形(各データの数が同じ場合は成功)
      DataMat = [data{:}];
      plot(DataMat)
      hold on 
      %
      D = dir(DebugFiles); 
  end

[michio]

forループ内最後の

 D = dir(DebugFiles);

が疑わしいですね。Dが上書きされているようですが、意図されたものですか？

[nknknknk]

script次第でclear all;してしまうから、 毎回登録するってことで書きました。 ただ、これを無くしてもエラーが出ます、、。

[nknknknk]

下記の変更を加えたところ、プログラムは通りました。

for k=3:length(D)　→　for i=3:length(D)

fid=fopen(D(k).name,'r'); →　fid=fopen(D(k).name);

しかし、

一つのファイル読み込み時は ans= 13 10 となり、図もプロット有りで表示されます(図onefile)が、

今回のプログラム（フォルダ内ファイル一括処理）だと

ans=[],ans=[],ans=[].... となり、図もプロットも無しになりました(図files)。

これはデータが正常に読み込まれてないのでしょうか。 それとも、

% 行列に整形(各データの数が同じ場合は成功) DataMat = [data{:}]; plot(DataMat) hold on の部分に問題があるのでしょうか。 分かりづらい質問で恐縮ですが、よろしくお願い致します。

[nknknknk]

解決しました。丁寧な対応ありがとうございました。

[michio]

解決されましたか。よかったです。タイムリーに対応できず失礼しました。

参考までに教えて頂ければ・・どの辺が問題だったのでしょうか。

[nknknknk]

for k=3:length(D)　→　for i=3:length(D) fid=fopen(D(k).name,'r'); →　fid=fopen(D(k).name); 上記変更で解決しました。

ありがとうございました。グラフについてまた質問させていただきます。 よろしくお願い致します。。