编写高效率IDL程序

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

2、EN BEGINsum = sum +aiENDIFENDFORfortime =SYSTIME ( 1 )-startPRINT , for time:fortimei= 0Lstart =SYSTIME ( 1)WHILE iLTN_ELEMENTS (a)- 1L DO BEGINIF (aiGT 100.0 ) THEN sum = sum +aii+ENDWHILEwhiletime =SYSTIME ( 1 )-startPRINT ,while time:,whiletimestart =SYSTIME ( 1)sum =TOTAL (a * (a GT 100.0 )fun t

3、ime =SYSTIME ( 1 )-startPRINT ,fun cti on time:,fun time;倍数print,result for VS function:,fortime/fu ntimeprint,result while VS function:,whiletime/fu ntimeEND运行后的输出IDL test time % Com piled module: TEST_TIME.for time:1.4510000while time:2.0430000 fun ctio n time: 0.054000139result for VS fun ctio n:

4、26.870301result while VS fun ctio n:37.833235看到差别了吧，循环比函数直接运算慢至少一个数量级!举例2:对一个图像中的特定值，若存在，则以该像素为中心，特定半径内的元素统一修改 为某值。以IDL自带的一个图像为例，将数据值等于142的赋为0.源码如下:SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS;计算两个点的距离fun cti onCalDista nee, poin t1, point2comp ile_ opt idl2Return , SQRT(point10-point20)人 2+(point11-point21

5、)人 2) endSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS;搜索当前坐标周围Distanee内的下标，注意输入x和y方向的坐标范围 xRange 和yRangefun cti oncalldxl nDista nee,curLoc, dista nee;,xRa ngeyRa nge;初始化临时下标suitLoc =0, 0；循环一次，计算矩形范围内的符合要求下标for xLoc = curLoc0-dista nee, curLoc0+dista needo beginfor yLoc = curLoc1 -dista nee, curLoc1 +dista

6、needo beginif calDista nee(curLoc, xLocyLoc)LE distaneethen suitLoc =suitLoc,xLoc,yLocen dforendretur n, suitLoc*,1:( N_Elements(suitLoc)/2-1)endssssssssssssssssssssssssssssssssssss;测试即调用主函数 pro test_processfile =;原数据FILE PATH ( rbcells.j pgSUBDIRECTORY =exa mp les,data)READ_JPEG , file, dataif size

7、 (data,/n_dimensions)ne 2 thenreturndims = size (data,/dimension)startTime =systime ( 1);当值等于142时，半径5内的元素赋值为 0eqValue =142rep Value =;符合要求的坐标suitIdx =calIdx In Dista nee(0, 0，5)；新数据n Data = dataidxs = where (dataeq eqValue,count)for curIdx =0 ,cou nt-1 do beginn Suit = suitIdxsuitLoc =;转换为二维坐标ARRAY_

8、INDICES (data, idxscurIdx)nSuit0 ,*= n Suit0,*+suitLoc0n Suit1 ,*= n Suit1 ,*+suitLoc1;下标要在数组自身范围内nSuit0,* = dims0 n Suit1 ,* = dims1 ;符合要求的位置赋值n Data nSuit0,* ,n Suit1 ,* = rep Valueen dfor;输岀花费时间 print , systime ( 1) - startTime;常规的循环写法ns = dims0nl = dims1new = datan ew1 = datastartTime =systime (

9、1)for i= 0,ns-do beg infor j= 0,nl-1 do beg inifn ewi,jeq eqValuethe n beg infor m=0,ns- 1do beginfor n=0,nl-1 do beg inif (m-i*1L)A 2+(n-j*1L)A 2 le25then n ew1 m,n=0en dforen dforen difen dforen dforprint , common time:,systime ( 1) - startTimewindow , 1 ,xSize = dims0* 2,ySize = dims1tvscl ,Data, 0tvscl,n data,1window , 2,xSize = dims0* 2,ySize = dims1tvscl ,n ew, 0tvscl,n ew1, 1I en