近期，需要用到modis数据批量预处理，于是重新学习idl,感谢郭师兄推荐，以及张洋老师的详细教导。特以此为学习笔记，望学有所成。

IDL注释为；

数据类型

数据类型

创建数组

Result = 函数（[D1[,…,D8]）至多八个维度
常见函数：
BYTARR: 字节型数组
DBLARR:双精度浮点型
FLTARR：浮点型数组
FINDGEN:浮点型索引数组
indgen:整型索引数组
STRARR:字符串
…
函数不区分大小写

BYTARR,DBLARR,FLTARR,INTARR, LON64ARR,LONARR,UINTARR,ULON64ARR,ULONARR,STRARR
BINDGEN,DCINDGEN, DINDGEN，FINDGEN,CINDGEN,LINDGEN,-SINDGEN,UINDGEN, UL64INDGEN, ULINDGEN

类型转换函数

BYTE, FIX, COMPLEX, DCOMPLEX, DOUBLE, FLOAT, LONG, LONG64, STRING，UINT, ULONG, ULONG64

整型 转换为 浮点型
a = indgen(5,5)
浮点型 转换为 整型
b = float(a)

FIX：转换为整型

IDL 中第一维是列，第二维是行，下标从0开始

代码输出

print

往文件中输出

openw,1,'D:/out.txt'

这里的1是硬盘的一个硬件存储位置，用1代替，运行后会创建一个txt

b=intarr(3,4)
printf,1,b
printf,1,'hello'

这个意思是往1当中，把b写进去。这里的1代指硬件地址。此时打开文件是看不到东西的，因为文件被idl占用了，所以需要归还文件,代码为：

free_lun,1

此时可以看到内容了。

注意:
此时再次，使用

openw,1,'D:/out.txt'
printf,1,'hello'
free_lun,1

此时文件中只剩下hello了，这是因为openw这个操作，不加其他关键字的时候，代表新建一个文件，若存在，则清空后新建，写入新的东西，原本东西不复存在。

那么如何保留之前的，继写呢？加关键字！

IDL> openw,1,'D:/out.txt',/append

这样就可以续写了！

help查询命令

HELP[,Expression1,…,ExpressionN]
查询对象相关信息
查询未知对象a的信息

help,a

可以输出数据类型 、行数、列数。

基础运算

创建数组的时候，一定要记住他的下标（几列几行）。且下标从0开始的。

a = findgen(3,4)
print,a[3,4]

上述代码是错误的，a没有第三列第四行，因为从0开始，下标最大只有[2,3]。

索引数组是从0开始，按行蛇形排列的。

IDL> a = findgen(3,4)
IDL> print,a 

结果为：
0.00000 1.00000 2.00000
3.00000 4.00000 5.00000
6.00000 7.00000 8.00000
9.00000 10.0000 11.0000

关系运算

a = findgen(3,4)
print,a gt 5

结果为：
0 0 0
0 0 0
1 1 1
1 1 1
返回的是0或者1

(a gt 5) * a
print，c

这样就可以把大于5的保存，小于5的赋值为0，结果为
0.0000000 0.0000000 0.0000000
0.0000000 0.0000000 0.0000000
6.0000000 7.0000000 8.0000000
9.0000000 10.000000 11.000000

控制语句

第一次程序编写

pro  program_name_firstClassprint,'helloa = 'cccc'print,aend

最基本的是:

pro 程序名程序题end

并且保存的文件名必须要 和 程序名 保持一致。

函数编写

如果需要函数，可以再pro前面先定义好function，比如

function fun_nameendpro  program_name_firstClassprint,'helloa = fun_name
end

第一次作业

pro  program_name_firstClassa = findgen(4,6)b = 3c = [3]d = [9,3,1]print,'a[3,4]'print,a[3,4]print,'a[15]'print,a[15]print,'a+b'print,a+bprint,'a[1,1]+b'print,a[1,1]+bprint,'a+c'print,a+cprint,'a+d'print,a+dend

输出：
a[3,4]
19.0000
a[15]
15.0000
a+b
3.00000 4.00000 5.00000
6.00000
7.00000 8.00000 9.00000
10.0000
11.0000 12.0000 13.0000
14.0000
15.0000 16.0000 17.0000
18.0000
19.0000 20.0000 21.0000
22.0000
23.0000 24.0000 25.0000
26.0000
a[1,1]+b
8.00000
a+c
3.00000
a+d
9.00000 4.00000 3.00000

解析：
① a[3,4]指的是第四列的第五行，行列不要搞混
② a+c等于一个单数3。是因为数组第二个元素维度不相同，则只有和【较小数组】一样大的元素进行相加后，得到结果。即c[0]+3 = 3.0000 ( 因为c是浮点数，所以默认是float )

pro  program_name_firstClassa = [[3,9,10],[2,7,5],[4,1,6]]print,a
end

结果为：

       3       9      102       7       54       1       6

pro  program_name_firstClassa = [[3,9,10],[2,7,5],[4,1,6]]print,'a'print,ab = [[7,10,2],[5,8,9],[3,1,6]]print,'b'print,bc = a+bprint,'c'print,cd = a*bprint,'d'print,d
end

A*B是对应位置元素相乘。

结果为：

a3       9      102       7       54       1       6
b7      10       25       8       93       1       6
c10      19      127      15      147       2      12
d21      90      2010      56      4512       1      36

除法的时候，处于整数要小心，是否结果需要的是float

pro  program_name_firstClassa = [[0,5,3],[4,0,2],[0,7,8]]b = [[0,0,1],[9,7,4],[1,0,2]]c = (a gt 3)*aprint,'a'print,aprint,'c'print,cprint,'b'print,bd = (b le 4)*b +(b gt 4)*9print,'d'print,de = (a+b)/2.0f = (a+b)/2print,'e'print,eprint,'f'print,fg = float(a+b)h = (a gt 0)+ (b gt 0)j = g/hprint,'j'print,j

结果为：

a0       5       34       0       20       7       8
c0       5       04       0       00       7       8
b0       0       19       7       41       0       2
d0       0       19       9       41       0       2
e0.00000      2.50000      2.000006.50000      3.50000      3.000000.500000      3.50000      5.00000
f0       2       26       3       30       3       5
j-NaN      5.00000      2.000006.50000      7.00000      3.000001.00000      7.00000      5.00000

解析：

1. d = (b le 4)*b +(b gt 4)*9
2. g = float(a+b)
   h = (a gt 0)+ (b gt 0)
   j = g/h

这两个逻辑关系要复习，在纸上画一下就可以知道为什么了

---