IDL实用函数

常用函数

注：
IDL命令行输入：
？函数名
直接查询与该函数有关的所有信息
借鉴博文：https://www.cnblogs.com/alex-bn-lee/p/5782132.html

1.make_array（n,m)

创建一个n列m行浮点型的数组`

IDL> a=make_array(3,2)
IDL> a
      0.00000000      0.00000000      0.00000000
      0.00000000      0.00000000      0.00000000

2.size( )

返回值是维度，列，行，数据类型，数据个数

IDL> size_a=size(a)
IDL> size_a
           2           3           2           4           6

3.typename()

返回变量的类型名

IDL> typename(a)
FLOAT
IDL> typename(size_a)
LONG

4. where

IDL> a=[4,3,6]
IDL> where(a gt 4,count)
           2
IDL> count
           1

（类似于matlab的find函数）

ENVI> a=[[1,2,3],[4,5,6]]
ENVI> a
       1       2       3
       4       5       6
ENVI> index = WHERE(a gt 3, count)
ENVI> index
           3           4           5
ENVI> value=a(WHERE(a gt 3, count))
ENVI> value
       4       5       6
ENVI> 

5.reverse

1-行反转/水平，2-列反转/垂直

IDL> a=[4,3,6]
IDL> where(a gt 4,count)
           2
IDL> count
           1

6.rotate

数组旋转，1-顺时针90°，2-180°中心对称，3-逆时针90°

IDL> a
       1       2       3
       4       5       6
IDL> rotate(a,2)
       6       5       4
       3       2       1      

7.rot

任意角度旋转,第二个参数是角度

IDL> rot(a,90)
% Compiled module: ROT.
       2       2       5
       1       1       4

8.sort

排序，从小到大，返回***索引***数组
上升：a[sort(a)]
下降：a[reverse(sort(a))]

IDL> a=[1,5,3,8,4]
IDL> sort(a)
           0           2           4           1           3
IDL> a(sort(a))
       1       3       4       5       8
IDL> a(reverse(sort(a)))
       8       5       4       3       1

9.UNIQ

相邻唯一值，返回索引数组
唯一上升数组：array[UNIQ(array, SORT(array))]

IDL> a=[1,1,2,5]
IDL> uniq(a)
           1           2           3
IDL> a(uniq(a))
       1       2       5

非相邻 唯一值：

IDL> a=[1,2,1,5,2,3]
IDL> uniq(a)
           0           1           2           3           4           5
IDL> a(uniq(a))
       1       2       1       5       2       3

10.max,min

      最大值，可选参数为最大值的索引
    MAX( Array [, Max_Subscript])

11. total

数据元素总和

IDL> a=[1,2,3]
IDL> total(a)
       6.0000000

12.product

数组元素乘积

IDL> a=[[1,2],[3,4]]
IDL> product(a)
       24.000000000000000

13.mean

均值

IDL> a
       1       2
       3       4
IDL> mean(a)
       2.5000000

14.alog10()

以10为底的对数

IDL> a=[10,100,1000]
IDL> alog10(a)
       1.0000000       2.0000000       3.0000000

15.exp

自然指数

IDL> exp(2)
       7.3890562
IDL> a=[1,2,3]
IDL> exp(a)
       2.7182817       7.3890562       20.085537

16.alog

自然对数
（此处存疑，alog(exp(1))=0.99999999）

IDL> a=[1,2,3]
IDL> alog(exp(a))
      0.99999994       2.0000000       3.0000000
IDL> alog(exp(1))
      0.99999994

17.sqrt

平方根

IDL> a
       4      16       9
IDL> sqrt(a)
       2.0000000       4.0000000       3.0000000

18.指数

IDL> a=[1,2,3]
IDL> a^2
       1       4       9

19.idl_validname

判断变量名是否有效，无效返回值为空或者自动修改

IDL> idl_validname('a')
a
IDL> idl_validname('a_b')
a_b
IDL> idl_validname('a_b c')

IDL> 

20.reform

调整数组的维度，元素个数不变

IDL> a=intarr(10,10,10)
IDL> b=reform(a,20,5)
% REFORM: New subscripts must not change the number elements in A.
% Execution halted at: $MAIN$          
IDL> b=reform(a,200,5)

21.rebin 分辨率变化

类似于分辨率变化，可增大分辨率，可减小分辨率

创建图像D，可见其行列数

IDL> D = SIN(DIST(50)/4)
IDL> 
IDL> im = IMAGE(D)
IDL> size(d)
           2          50          50           4        2500

增大其分辨率，（当然，也可减小分辨率）

IDL> d_rebin=rebin(d,250,250)
IDL> image(d_rebin)
IMAGE <66982>
IDL> size(d_rebin)
           2         250         250           4       62500

22.congrid

重采样，类似REBIN

可视化控件

1.dialog_message

默认的dialog_message的dialog type 是warning

	dialog_message("Right?", /question, title="Look: ", /cancel)

IDL> dialog_message('warning')

IDL> dialog_message('error',/error)

IDL> dialog_message('info',/information)

2.dialog_pickfile

◈ 文件选择，返回值为文件路径
◈ 指定 DIRECTORY 关键字，文件夹选择，返回值为文件夹路径
DEFAULT_EXTENSION 自动增加扩展名；
DIRECTORY 选择文件夹；
FILE 可以使用通配符进行筛选（FILE = ‘sa*’）；
FILTER 扩展名设置（FILTER = [’.jpg’, '.tif’, ‘*.png’]）；
MULTIPLE_FILES 可以选择多个文件；
PATH 设置初始进去的路径，不指定的话就进到当前的工作空间；
READ 标题为 “Select File to Read”。
WRITE 标题为 “Select File to Write”。

IDL>	dialog_pickfile(title="Select", FILTER = ['*.jpg', '*.tif', '*.png']) 
IDL> 	dialog_pickfile(title="Select", FILTER = ['*.jpg', '*.tif', '*.png'],path='D:\',/multiple_files,/read) 

更多信息，请使用:

IDL>  ？dialog_pickfile