テキストファイルへの変換

Question

観測で得られたファイルをMATLABで読み込めるようにテキストファイルに変換したいです。

まず、1分毎に記録されたファイルが、

---

e16A0319.272413 2016/10/03 19:26
e16A0319.272413 2016/10/03 19:27
・・・

---

という風に並んでいます。
その中から、
「e16A0319.272413」
の中身を表示すると、

---

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ﾋｾ・ｦ
・・・・・・・・・・・・・・・・

---

という表記があります。

このファイルをMATLABの関数で読み込めるように、テキストファイルに変換したいのですが、下記のプログラムでエラーが出てしまいます。
下記プログラムで足りない部分、根本的に間違っている部分を教えていただきたいです。

---

fid =fopen('e16A0319.282564','r');

%ヘッダー7行飛ばす
for i = 1:8
header = fgetl(fid);
end
% 各データセットの読み込み
data = cell(1, 4);
for i = 1:4
data{i} = fread(fid, 512, 'long');

% <CRLF> の読み飛ばし
if i ~= 4
    fread(fid, 1, 'short');
end

end

fclose(fid);

---

データフォーマットは「http://licel.com/raw_data_format.html」です。

以前にここで似たような質問をさせていただき、丁寧なご回答をいただいたのですが、解決できずにまた質問させていただきました。
お詳しい方よろしくお願いします。

Answer 1

どのようなエラーかわかりませんが
data={:,i}
で大丈夫そうな予想

Comments

[Tannnnnnnnnnn]

ご回答ありがとうございます。
data ={:,i}
で実行しましたが、下記のエラーが出ました。

エラー fgets
ファイルの識別子が無効です。有効なファイルの識別子を生成するには fopen を使用してください。
エラーの発生場所 fgetl (line 34)
[tline,lt] = fgets(fid);
エラーの発生場所 bynari (line 5)
header = fgetl(fid);

どのような対策が考えられるでしょうか。

[iwanote]

プログラムはこれですべてですか？

[Tannnnnnnnnnn]

はい、すべてです。
全体の構造は、
ーーーーーーーーーー
固定ヘッダ行1
固定ヘッダ行2
固定ヘッダ行3
データセットヘッダ行1
データセットヘッダ行2
データセットヘッダ行3
データセットヘッダ行4
データセット1 (512x4バイト)
<CRLF>
データセット2 (512x4バイト)
<CRLF>
データセット3 (512x4バイト)
<CRLF>
データセット4 (512x4バイト)
ーーーーーーーーーーーーーーーーーーー
になっています。
全てのデータが、データセット4、データ512×４バイト、
32ビットデータをリトルエンディアンの整数と仮定したコードです。

[Tannnnnnnnnnn]

連続ですみません。
fid =fopen('e16A0319.282564','r');⇒fid =fopen('e16A0319.282564','ｗ');
と変更したところ、エラーが消え、
e16A0319.282564　ファイルがMATLABに読み込めたのですが、
中身が空です。
これは正確に読み込めていないということでしょうか。

[iwanote]

w指定だと中身は破棄されてしまいますよ
r+だとどうです？

[Tannnnnnnnnnn]

w指定だから空になってしまうのですね。
r+だと以下のエラーが出てしまいます。

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

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

エラーの発生場所 bynari (line 5)
header = fgetl(fid);

[iwanote]

ファイル名とそのディレクトリは正しいですか？
fidの数値を調べてみてください

[Tannnnnnnnnnn]

ファイル名は正しいのですが、
fidの数値は－１でした。

[Tannnnnnnnnnn]

-1は何を表しているのでしょうか。
無知で申し訳ありません。

[iwanote]

ではファイルの位置がおかしいかもしれません
もしくは拡張子が書いていないとか…？

[iwanote]

-1は読み取れていないものです

[Tannnnnnnnnnn]

ファイルの位置を間違えていました、再度fidの数値を確認したところ　3　でした。
観測データのファイル名が「e16A0319.282564」なので、拡張子がもともとありません。。

[iwanote]

3なら正常に読み取れていると思いますがエラーは出ていますか？

[Tannnnnnnnnnn]

お早い返信ありがとうございます。助かります。

プログラムは通りましたが、そのファイルを開くと下記の表示になり、下部分が変換されていないのですが、、プログラムに問題があるのでしょうか。。
ーーーーーーーーーーーーー
e16A0319.282564
TMU 03/10/2016 19:27:24 03/10/2016 19:28:25 0000 0000.0 0000.0 00
0000000 0010 0018953 0010 04
1 0 2 00512 1 0000 0015 00829.o 0 0 00 000 12 150025 0.020 BT1
1 1 2 00512 1 0000 0015 00829.o 0 0 00 000 00 150025 3.1746 BC1
1 0 2 00512 1 0000 0015 00829.o 0 0 00 000 12 150025 0.020 BT1
1 1 2 00512 1 0000 0015 00829.o 0 0 00 000 00 150025 3.1746 BC1

h�S2{4�ﾟ糠�ｰ�c�ｸ�4@ﾛ�'�ﾁH熬
ｼ)
ｴﾁ
b�ﾛH/ﾛ
ｾ&ﾄ;A()xF広tX鉄ﾐjw芭礼��ﾎ�7�ｧ�:･ｰ嫉鮃ｾ|ﾀﾘﾈﾎ﨓=ﾙKﾘ8�晨ﾝ�w��ﾙ3�Z
ｴN!�ﾎ(/z::6DsH!LJｪUﾏV_”0e(o nu�ｮ�ｰ�ﾍ�逶棗ｺｵ(ｸﾍｸjﾃﾈ�@ﾎｳﾔSﾜﾕ�杢x�>�T��ﾅ�ﾜ�瑪H�?�S�纏K:"�r g'.M�ｾ�v�ve^ﾟ {ﾞﾋ+今|5
･ﾙｺ�@w�+�ｽ�Gヱｷｱ ﾁlｺ1聾ﾄﾍ?�ｷﾊｽ-ﾏﾔﾌr(ﾗI"mﾅ
ｍ)ﾓﾏ6@ｦ�cｹ�;ﾊtl>ｸｺg
'Yn��WﾁXOﾉ;k8ﾓ�ｼ
e6;
`Jﾋ*得<K
�吋f7A 9U愰�
Oﾐ蕭ｷ�4�t
�
y�7xt�ﾁ�[�佐
~�K�F�ﾉ五 #
h_  ��ｵ,U)ﾄn
ｳz  �u
ﾉｹ ﾄﾃ
!
s｢ lm
2+樟
'oﾒ枚4�m｡ @+
ｙ季
ｵe ;ﾌﾆﾇ�� 赳
ｴax
6�7�7=
J
Wﾂ
･ﾜｴZﾎﾁ-#:(m
�
pｼ
ﾍ�ﾍ]Pｭ�L~
ﾊk
ｭﾋ
oD楓+ｨｶ

4u
ｻ
S}4'Zﾌ
�x
?�ﾌ
�ﾍ [�G�c
U
Tv
ﾙｯ
ｫﾒ
Bﾐ
ｫｨ
{
Q
a
ｷs
ｻ｣
､ｷ
愀
n�ｭq

ﾂ_
ﾘ|
ｮ�擯
ｫ
囿
ﾄ�｢i
考
從
2�n｣
ﾍｭ
?ｮ
?�
�
縒
<�?�箙
德
ﾑｳ


[Tannnnnnnnnnn]

ちなみに最初の質問文と今実行したファイルは別物です。。

[iwanote]

CRLF読み飛ばしのfreadの第二引数は2だと思います
あとは…もしかしたらdata{i, :}かもしれないです…

[Tannnnnnnnnnn]

ご回答ありがとうございます。
freadの第二引数２でも結果は変わらず、
date{i,:}ではエラーになりました...。

[iwanote]

こちらでも確認してみましたが
data = cell(1, 4);
for i = 1:4
data{:,i} = fread(fid, 512, 'long');
こちらでdata配列に数値が入っているのを確認しました。
正常なデータがどのようなものかわからないため何とも言えませんが...
ステップ実行をすれば内部データがわかるので確かめで見てはどうでしょうか？

[Tannnnnnnnnnn]

ご回答ありがとうございます。
下記のプログラムを実行したところ、date配列の中身に数値が入っていました。

ーーーーーーーーーーーーーーーーーー
fid =fopen('e16A0319.292710','r+');

for i = 1:8
header = fgetl(fid);
end
% 各データセットの読み込み
data = cell(1, 4);
for i = 1:4
data{:,i} = fread(fid, 512, 'long');

% <CRLF> の読み飛ばし
if i ~= 4
fread(fid, 2, 'short');
end
end

celldisp(data)
fclose(fid);
ーーーーーーーーーーーーーーーーーーーーーー

いまいち理解仕切れていないのですが、
このdate配列に数値が入っているということは、
元の’e16A0319.292710’ファイルをテキストファイルに変換できたということでしょうか。

[Tannnnnnnnnnn]

このプログラムを実行しても、’e16A0319.292710’ファイルの、’h�S2・・・・’の部分は変わっていません。

[Tannnnnnnnnnn]

このプログラムの正しい実行後、’e16A0319.292710’ファイルの、’h�S2・・・・’の部分が数値に変わる。というそもそもの理解が間違っているのでしょうか。

分かりづらい質問で申し訳ございませんが、よろしくお願い致します。

[iwanote]

ファイルへアクセスしている関数は今のところfgetlとfreadだけです。
この二つはファイルからの読み出しをしているだけです。
ファイルへの書き込みをするばあいはfprintfなどをする必要があります。

[Tannnnnnnnnnn]

ご回答ありがとうございます。
最後に質問なのですが、
「e16A0319.272413」のファイルを、拡張子がtxtに変換して出力したいのですが、
どのような方法があるのでしょうか。
長々と質問した後に恐縮ですが、よろしくお願い致します。

[iwanote]

fopen関数に作りたいファイル名を入力すればそのファイルが作られます。
a = fopen('b.txt','w');
このようにすればb.txtが作られ、a(fidと同じ種類の変数)にそのfileIDが渡されます
詳しくはfopenやfprintfなどのMATLABのドキュメンテーションを調べてみるといいかもです。