ALOS-2/PAL​SAR-2のデータの​読み込み方

ALOS-3,/PALSAR-3 での観測および高分解能 DSM を利用した干渉SAR処理を模擬し，

MATLABでALOS-2/PALSAR-2のSARデータを可視化することを試みております．

初めに，PALSAR-2のL1.1データを読み込みたいと考えています．

読み込みたいCEOSファイルは，VOLデータ（2KB），LEDデータ（1,572KB），IMGデータ（5,939,982KB），TRLデータ（766KB）の順に入っています．

読み込みたいファイルデータは，IMGのみになります．

この際，記述するコードが知りたいです．

お力添えお願いいたします．

Answer 1

[Shunichi Kusano]

こんにちは。

フォーマットがわからないのでコードまではわかりませんが、手順としては

1. fopenでファイルを開く
2. fseekで読みたいところまで移動（←フォーマット説明書で確認）
3. freadで指定したデータ型で指定した分読み込む（←フォーマット説明書で確認）
4. fcloseでファイルを閉じる

でワンセットです。freadのドキュメントを見れば典型的なコードは書いてますので、参考にしてください。

ちなみになのですが、干渉SARをしたい場合、基線長も計算する必要がありますので、画像データだけじゃなく衛星の位置と速度の情報も必要になります。それはLEDに記録されていたような…ちょっと不確かなのでフォーマット説明書を読んでみてください。読み込み方は画像と同じ要領です。

[Ryosuke KIJIMA]

>> fileID = fopen('0000312085_001001_ALOS2096552950-160307.zip');

>> tline = fgetl(fileID)

エラー: fgets

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

エラー: fgetl (line 32)

[tline,lt] = fgets(fid);

>> fileID = fopen('IMG-HH-ALOS2096552950-160307-UBSR1.1__D');

>> tline = fgetl(fileID)

エラー: fgets

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

エラー: fgetl (line 32)

[tline,lt] = fgets(fid);

>> A = importdata('IMG-HH-ALOS2096552950-160307-UBSR1.1__D')

エラー: importdata (line 139)

ファイルを開けません。

うまくいきません。何が原因なのでしょうか。

[Shunichi Kusano]

fgetl, fgets, importdata,どれもバイナリデータの読み出しには使うことはできません。freadを使えば読めます。外部のページですが、バイナリデータについては次のページの「エンディアン」くらいまで読めば一通り勉強できそうです。

https://programming-place.net/ppp/contents/c/042.html