北大《空间探测信息处理技术(IDL)》第2章 数据、变量、表达式.ppt

表达式 例：创建一个取值范围在-50到50之问的数组MyArray(100, 100)，其各元素的值为随机值，按照下列要求完成相应的运算。 把所有值大于等于30的元素的值改为30；并统计其个数，存入变量sMax。 把所有值小于等于-30的元素的值改为-30；并统计其个数，存人变量sMin。 计算所有取值在-20到20之间的数的个数日SMid，然后对其绝对值分别开平方。 统计数组中所有剩余元素的个数，存入变量SRest，并把其值清O。 输出数组中所有元素的值以及SRest，sMax和SMin的值。 表达式 ; 创建一个取值范围在-50到50之问的数组MyArray(100,100) ; 其各元素的值为随机值 a=-50.0 b=50.0 MyArray = (b-a)*RANDOMU(1, 100, 100)+a ; 1.1 查找大于等于30的元素及个数 IdxGE30 = where(MyArray GE 30, SMax) ; 2.1 查找小于等于-30的元素及个数 IdxLEn30 = where(MyArray LE -30, SMin) ; 3.1 查找-20至20之间的元素及个数 Idx20 = where(abs(MyArray) LE 20, SMid) ; 4.1 查找其余元素及个数 Idx2030 = where((abs(MyArray) GT 20) and (abs(MyArray) LT 30), Srest) ; 1.2 大于等于30的赋值30 MyArray[IdxGE30] = 30 ; 2.2 小于等于-30的赋值-30 MyArray[IdxLEn30] = -30 ; 3.2 -20至20之间的元素被开平方 MyArray[Idx20] = sqrt(abs(MyArray[Idx20])) ; 4.2 -20至-30, 20至30之间的元素赋值0 MyArray[Idx2030] = 0 ; 输出结果 print, SMin, SMid, SMax, SRest, SMin+SMid+SMax+SRest end 内存变量 整型数使用注意 整型数除法为整型数 IDLresult=12/5 IDL print, result 2 例：假如想知道IDL图形窗口的比率，而窗口的大小被储存在两个系统变量 !D.X_Size, !D.Y_Size （整型）中，也许会写出如下的IDL代码： aspect=!D.X_Size / !D.Y_Size 正确的写法应该是： aspect=float(!D.X_Size) / !D.Y_Size IDL的整型数为短整数（2字节）而其他语言的整型数是4字节。 内存变量 矢量和数组 创建常数矢量 矢量就是一维的数组；创建常数矢量用方括号，即：“[ … ]”，… 为元素值，各元素用 “,” 分开，例如： IDLV=[1,2,3] ; 整型数矢量 IDLA=[1.0,2.0,3.0] ; 浮点数矢量 获取变量属性 IDLhelp, 变量名 例如： IDL help,A A FLOAT = Array[3] 内存变量 矢量和数组 下标引用 传统的IDL实用圆括号引用矢量元素，即：“(…)”，… 为下标，例如： IDLI= V(2) ; 引用元素2 IDLI= V(0:2) ; 引用元素0, 1, 2 注意：1）下标排列从0开始， 2）从IDL5以后可以使用方括号。 DLprint, V(1:2) 2 3 DLprint, V[1:2] 2 3 内存变量 矢量和数组 矢量的扩展 IDLV=[V[0:1],5,V[1:2]] IDLPrint, V 1 2 5 2 3 用函数创建矢量 例如：建立一个值在0到50之间的6个元素浮点矢量，可键入： IDLvector=FIndGen(6)*10 IDLPrint,vector 0.000000 10.0000 20.0000 30.0000 40.0000 50.0000 这里：FIndGen(n)函数生成从 0.0 到 n-1 的5个浮点数. 内存变量 用函数创建矢量 array = INTARR( D1[, ..., D8] [, /NOZERO]) 数据类型 初始化函数 产生索引值的函数 字节型 BytArr(...) BIndGen(...) 16位有符号整型 IntArr(...) IndGen(...) 32位有符号长整型 LonArr(...) LIndGen(...) 64位有符号整型 Lon64Arr(...) L64IndGen(...) 16位无符号整型