查看“IDL8”的源代码


'''相对于IDL7.1的IDL8系列的新功能'''
==功能==
*python bridge
:关键是要设置IDL的库文件在[[动态链接库文件]]的目录中
*jupyter(pyhton), IDL kernel, 参见[http://www.harrisgeospatial.com/docs/idl_kernel.html]
#将复制到<IDL_DIR>/lib/bridges/kernels/IDL/kernel.json 复制到~/share/jupyter/kernels/IDL
#设置，PATH包含IDL的可执行命令文件夹
#设置环境变量PYTHONPATH包含IDL和PYTHON的可执行命令文件夹

==函数==
*idl_convert （单位转换，很实用），参见IDLUNIT() 

==变量==
*变量的运算 [http://www.harrisgeospatial.com/docs/whats_new_8_4.html]
 var = RANDOMU(seed, 200, 100)
 PRINT, var.length
 PRINT, var.ndim
 PRINT, var.dim
 PRINT, var.typecode
 PRINT, var.typename
 PRINT, var1.Mean()
 PRINT, var1.Total()
 var2 = var1.Sort()

*biginteger
:任意位数的整数
===系统变量===
*!const
:常用物理常数


==函数==
*FOLDERWATCH 
:监控文件夹
* CODE_COVERAGE
:监控那些CODE被执行 （选择判断）

*FFT Power Spectrum




==画图==

===IDL graphics===
新的图形系统相比于之前的Direct Graphics系统的差别是可以动态的调整图形元素

*P1.window.setcurrent ; 将某个图形窗口设为当前窗口

====plot函数====
[[IDL8中的plot函数]]

====legend====
  a = findgen(11)
  b = 20./(a+1)^2
  p1 = PLOT(a, b, '4', /HISTOGRAM, NAME='Theoretical',  /YLOG, YRANGE=[0.01,100])
  p2 = PLOT(RANDOMU(-1, 11), '4r', NAME='Experimental', TRANSPARENCY=50, /HISTOGRAM, /OVERPLOT) 
  h = LEGEND(POSITION=[9,50],label=['haha','papa'] /DATA)
 P3=plot(/test,/overplot)
 h1.add,p3   ;自动增加一个标题，不能是数组

====text函数 ====
*和原来的xyouts对应, 直接支持latex符号
* 注意，在多个panel的图上，使用/data定位的时候，要用target关键词制定图形对象
  fn = exp(SIN(FINDGEN(100))^FINDGEN(100))
  p1 = PLOT(fn, LAYOUT = [1, 2, 1], /CURRENT, TITLE = '1')
  p2 = PLOT(fn, LAYOUT = [2, 2, 3], xrange=[25, 42], /CURRENT, TITLE = '2')
  t=text(30,2.5,'test1',/data,target=P1) 
  t=text(30,2.5,'text2',/data,target=P2)

====errorplot====
完整的误差棒程序

====barplot====
:棒状图

====boxplot====

====polarplot====
:极坐标

====vector====
:画矢量图（箭头）

====symbol====
:画符号，支持latex


====axes====
*AXES可以直接在image上画坐标轴

===image===
*函数
:RGB_VALUE : 选择color bar
:max：可以用最大值归一化
:axis: 图像是否显示坐标轴
:p.axes 可定制
       p=image(Agez,position=[0.15,0.1,0.9,0.5],max=max(Agez),RGB='rainbow',axis=1,$
                xtitle='Age (Gy)',ytitle='$Z/Z_\odot$',/current)
        ax=p.axes
        ax[0].tickvalues=[0,2,4,7,9,11,13,15,17,19]
        ax[0].tickname=['0.001','0.01','0.1','1','3','5','7','9','11','13']
        ax[1].tickname=['0.005','0.002','0.02','0.04','1','2.5']


===地图相关===
*map (可以投影天球)
 mymap = MAP('hammer', FILL_COLOR="light_blue")
 ; Convert from longitude/latitude to meters
 longitudes = [-83.02, 88.50]
 latitudes = [42.42, 47.17]
 xy = mymap.MapForward(longitudes, latitudes)
 ; Convert back to longitude/latitude and plot
 EXPAND, xy, 2, 20, xycoordinates
 lonlat = mymap.MapInverse(xycoordinates)
 p = PLOT(lonlat[0,*], lonlat[1,*], 'rS',   /OVERPLOT, /SYM_FILLED)

==数组运算==
*array_equal