IDL (The Interactive Data Languege) 覚書き

初期設定

方法

IDL の設定ファイルを読み込むため idl_setup.bash を読み込むよう .bashrc に書き込んでおく。
```
. /usr/local/rsi/idl_5.5/bin/idl_setup.bash
```

（最終更新日：2003/12/09）

起動と終了

方法

IDL Development Environment （以下 DE ）を起動。
```
       $ idlde
       
```
IDL Development Environment を終了。 IDL プロンプト上で
```
       IDL> exit
       
```
demo を実行。IDL プロンプト上で
```
       IDL> demo
       
```

（最終更新日：2003/12/09）

ASCII データ読み込み

方法

ascii_template 関数を使ってデータの型を定義。以下を実行するとファイル選択ダイアログが表示される。
```
       IDL> datatemplate=ascii_template()
       
```
テンプレートダイアログが表示される。
データファイルの区切りがわかっていれば Delimited を選択。
Comment String to Ignore: でコメント記号を指定すればそれは無視される。
Data Starts at Line: でデータの開始行を指定。
Next >
Delimiter Between Data Elements: で区切り文字を指定。
Next >
各列をクリックし（field1など) Name: データ名を入力する。
Finish >
テンプレートが valuechec window に表示される。
ファイル選択ダイアログを呼び出し、変数にファイルのパスを入れる。
```
      
      IDL> file=dialog_pickfile()
      
      
```

read_ascii 関数でファイルを読み込む。

      
      IDL> data=read_ascii(file, template=datatemplate)

logmessege window に出力。

      
      IDL> print, data

特定の列だけを表示するには以下のようにする。
```
      
      IDL> print, data.field1
      
      
```
構造体の下層データのアクセスは . で変数名をつなぐ。 field1はテンプレートダイアログで決めた列の名前。

dialog_pickfile() を以下のようにするとファイル選択ダイアログで選択候補にフィルタをかけることができる。
```
       
       IDL> file=dialog_pickfile(filter='*.dat')
       
       
```

参考

無いからつくった IDL まにゅある
Getting Started with IDL

（最終更新日：2003/12/09）

読み込んだデータから 2D グラフを表示

方法

上記のようにファイルからデータを読み込んでおく。

データをプロットする。

       
       IDL> plot, data.field1, data.field2

グラフにタイトルや軸のタイトルを入れるには以下のようにする。

       
       IDL> plot, data.field1, data.field2, $
       title='title', xtitle='xtitle', ytitle='ytitle'

参考

IDL Basics
IDL Reference Guide
無いからつくった IDL まにゅある

（最終更新日：2003/12/09）

テキストデータの書き込み

方法

ファイル選択ダイアログを呼び出し、変数 file にファイルのパスを入れる。

       
       IDL> file=dialog_pickfile(/write, filter='*.dat')

ファイル file を開く。
```
       
       IDL> openw, 1, file
       
       
```
'1' は適当に割り当てたは理論ユニット番号で、以降このファイルはこの理論番号で参照される。

変数 data を file （1）に書き込む。

       
       IDL> printf, 1, data

ファイルを閉じる。

       
       IDL> close, 1

データのフォーマットを指定するには printf で以下のようにする。
```
       
       IDL> printf, 1, 'datafile.dat', format='(5f8.3)'
       
       
```
5 は値を横に並べる数。f は実数を表す。 8 は全桁数（小数点も含む）で 3 は小数点以下の桁。ここでは直接ファイルを指定している（ファイル名を "'" で囲む）。
以下 lun はユニット番号。/get_lun でユニット番号を自動的に割り当てる。 /get_lun でユニット番号を割り当てたときは free_run で閉じる。(lun も開放される。)
```
       
       openw, lun, mod_file, /get_lun
       printf, lun, format='(i4, f15.2)', j, k
       free_lun, lun
       
       
```

参考

IDL Basics
無いからつくった IDL まにゅある

（最終更新日：2003/12/10）

プログラミング

書式

IDL には 4 つのタイプのコードがある。

プロシージャ
関数
メインレベルプログラム
バッチファイル

以下、プロシージャのについて。

pro で始まり、end で終わる。

       
       pro proname [,argument ,...]
       ...
       ...
       end

proname はプロシージャ名。 argument は引数。

ASCII ファイルからデータを読み込み、ファイルに書き込むプロシージャの例。

       
       pro readwrite    ; ';' から文末まではコメント

       datatemplate=ascii_template()
       file=dialog_pickfile(filter='*.DAT')
       data=read_ascii(file, template=datatemplate)

       file=dialog_pickfile(/write, filter='*.txt')
       openw, 1, file
       printf, 1, data.field1, format='(1f9.5)'
       close, 1

       end

保存・コンパイル・実行

保存時のファイル名は プロシージャ名.pro にする。プロシージャ名と同じにしたほうがよい。
コンパイルは DE では Run メニューの Compile 。または Ctrl+F5 。コマンドライン上では
```
       
       IDL> .run proname [,argument, ...]
       
       
```
もしくは
```
       
       IDL> .rnew proname [,argument, ...]
       
       
```
2 つ目の方法は以前にコンパイルしたプログラムの変数がすべて消去される。

実行は Run メニューの Run。または F5 。コマンドライン上では

       
       IDL> proname [,argument, ...]

参考

IDL Basics
無いからつくった IDL まにゅある

（最終更新日：2003/12/10）

Tips （または未整理）

IDL コマンドライン上では「上矢印」キーで履歴を呼べる。
window ：ウィンドウを開く。
window, xsize=800, ysize=600 ：大きさを指定してウィンドウを開く。
wdelete ：ウィンドウを閉じる。
2次元配列要素は
```
      
      arr[列:行]
      
      
```
の順。

複数の配列要素に値を代入

      
      IDL> arr[0:9] = 0

配列の全要素に 1 を足す。
```
      
      arr=arr+1
      
      
```

多次元配列の定義例

      
      IDL> arr = [[1,2,3],[4,5,6]]   ; 2 次元配列

可変長配列

      
      i = 0
      repeat begin
        arr = [arr,i]
        i = i + 1
      endrep until (i=10)

初めに arr を定義しておく必要がある。e.g. arr=0

ベクトルの要素数を調べる

      
      print, N_ELEMENTS(arr)

xyouts, 100, 50, 'hoge' ：ウィンドウの (100,50) に文字列 hoge を表示する。
spawn, 'command' ： UNIX コマンドを実行。

プロンプトから数値を読み込む。


      IDL> var=0
      IDL> read, var

プロンプトの変更は


      IDL> read, var, PROMPT='hoge'

とする。

プロンプトから文字列を読み込む。


      IDL> var=''
      IDL> read, var

ascii file からデータを読み込む。
```
      IDL> data=read_ascii(file, delimiter=',')
      
```
delimiter は指定しなければスペース。
ascii fileの最後の行末で改行されていないとascii fileとして認識してくれない。
1 byte (定数)は "1b" 。
空白を表すバイトコードは "32b" 。

DEVICE coordinates

       +----------------+(Vx-1,Vy-1)

       |                |

       |     region     |

       |                |

  (0,0)+----------------+

!D.X_SIZE, !D.Y_SIZE
比で表す normal

if文を入れ子にする。


      if (x eq 1) then statement1 else $
      if (x eq 2) then statement2 else $
      if (x eq 3) then statement3 else $
      if (x eq 4) then statement4 else $
      if (x eq 5) then statement5 else statementx

Building IDL Applications Chapter12 : Program Control 参照。

case文。


      case var of
        'hoge'    :  statement1
        'hero'    :  statement2
        else      :  begin
                       statementx
	             end
      endcase

Building IDL Applications Chapter12 : Program Control 参照。
begin ... end は命令文のまとまりを作るコマンド。

x 軸の数値を表示しない。
```
      plot xtickformat='(a1)'
      
```
idl では int は 2 byte 整数 (-32,768～32,767) なので、ループ変数などは long 4 byte 整数にしておいた方がよい。
```
      i = 0L
      
```

数値を文字列に変換


      num = 100
      str = string(num,format='(i3)')

文字列から一部を切り出す


      date  = '2005/01/10'
      year  = strmid(date, 0, 4)
      month = strmid(date, 5, 2)
      day   = strmid(date, 7, 2)

      print, 'year:', year
      print, 'month:', month
      print, 'day:', day

      (output)
      year:2005
      month:01
      day:10

（最終更新日：2005/02/23）

IDL (The Interactive Data Languege) 覚書き

目次

初期設定

方法

起動と終了

方法

ASCII データ読み込み

方法

参考

読み込んだデータから 2D グラフを表示

方法

参考

テキストデータの書き込み

方法

参考

プログラミング

書式

保存・コンパイル・実行

参考

Tips （または未整理）