编写高效率的IDL程序

1、程序的效率问题，在大数据或复杂运算的时候是不能忽略的。但在IDL程序的编写方式上，不能按照常规的循环 for依次处理方式写，简单归纳下，提高效率的运行的写法注意下 面两种方式。1、尽量避免或少用循环2、多用 where 和 Histogram；说起来很容易，但实际写的时候一定要多斟酌斟酌。 举例1:对2000*2000的数组中大于100的值进行累加。 PRO TEST_TIME ;a = DIST ( 2000 , 2000 ) sum =0.sum1 =0.start = SYSTIME ( 1) FOR i= 0L , N_ELEMENTS (a)- 1L DO BEGIN IF (ai

2、GT 100.0 ) THEN BEGIN sum = sum +aiENDIF ENDFOR fortime = SYSTIME ( 1 )-start PRINT , for time: fortime i= 0L start = SYSTIME ( 1) WHILE i LT N_ELEMENTS (a)- 1L DO BEGINIF (ai GT 100.0 ) THEN sum = sum +ai i+ ENDWHILE whiletime = SYSTIME ( 1 )-start PRINT , while time: ,whiletime start = SYSTIME ( 1

3、) sum = TOTAL (a * (a GT 100.0 ) funtime = SYSTIME ( 1 )-start PRINT , function time:,funtime;倍数 print , result for VS function:, fortime/funtimeprint , result while VS function:,whiletime/funtimeEND 运行后的输出 IDL test_time % Compiled module: TEST_TIME. for time: 1.4510000 while time: 2.0430000 functio

4、n time: 0.0 result for VS function: 26.870301 result while VS function: 37.833235看到差别了吧，循环比函数直接运算慢至少一个数量级！举例2：对一个图像中的特定值，若存在，则以该像素为中心，特定半径内的元素统一修改 为某值。以IDL自带的一个图像为例，将数据值等于142的赋为0.源码如下：:；;计算两个点的距离；function CalDistance , pointl, point2compile_optidl2Return , SQRT(point1end；0-point20)A 2+(point1 1 -poi

5、nt21)A 2)搜索当前坐标周围 Distance 内的下标，注意输入 x和y方向的坐标范围xRange和yRangefunctioncalldxInDistance,curLoc, distance;,xRange,yRange;初始化临时下标suitLoc = 0, 0；循环一次，计算矩形范围内的符合要求下标for xLoc = curLoc 0-distance, curLoc0+distance do beginfor yLoc = curLoc1 -distance, curLoc1 +distancedo beginif calDistance (curLoc, xLoc,yLo

6、c)LE distance then suitLoc =suitLoc,xLoc,yLoc endfor end ； return , suitLoc*, 1:( N_Elements (suitLoc)/ 2-1) end ; ； ;测试即调用主函数 pro test_process ;原数据 file = FILEPATH ( rbcells.jpg , $ SUBDIRECTORY = examples , data) READ_JPEG , file, data ； if size (data,/n_dimensions)ne 2 then returndims =size (data

7、,/dimension)startTime =systime(1)；当值等于142时，半径5内的元素赋值为0eqValue =142repValue = i;符合要求的坐标suitIdx =calldxInDistance(0, 0, 5)；新数据nData = dataidxs =where (data eq eqValue,count)for curldx =0 ,count- 1 do beginnSuit = suitIdxsuitLoc =ARRAY_INDICES (data, idxscurIdx)nSuit0,*= nSuit0 ,*+suitLoc0nSuit1 ,*= nS

8、uit;下标要在数组自身范围内1 ,*+suitLoc1nSuit0 ,* = dims0 nSuit1 ,* = dims;符合要求的位置赋值1 000,*,nSuit1 ,* = repValue;转换为二维坐标nDatanSuit endfor；输出花费时间 print , systime;常规的循环写法(1) - startTimens = dims nl = dims01new = datanew1 = datastartTime =for i= 0 ,ns-systime ( 1)1 do beginfor j= 0 ,nl- if newi,j1 doeqfor m= 0,ns- 1begineqValue then begin do beginfor n= 0,nl-if (m-i*endfor1 do begin1L2+(n-j*1L)A 2 le25then new1m,n= 0endforendifendforendfor1print , common time: window , 1 ,xSize = dims tvscl ,Data, 0,systime (1) - startTime 0* 2,ySize = dims0* 2,ySize = dims1tvscl ,ndata, 1window , 2 ,xSi