IDL渐变散点图

秋之期

查到说8.2.2以后版本的IDL中有引入scatterplot函数，可以完成有颜色变化的散点图，例图如下，但是看help没有特别明白，按照我理解的去做了，程序一直报错，所以想问问有没有人用过这个函数，或者有简便的方法能够实现渐变散点图的。
我的程序段：
set_plot,'ps'
device,xoffset=.2,yoffset=10.8,xsize=11.,ysize=8.,/inches,/landscape,$
    filename=figname,/times,/bold,/color
  levels=[-9999,0,50,100,150,200,250,300,500,700,900,1000]
  contour,ele_CN,lon_CN,lat_CN,pos=pos1,levels=levels,c_colors=colors,/fill,xrange=[112,122],yrange=[28,36],$
    title='Ave Annual Occurrence of Short Heavy Rain',xtickname=['112','114','116','118','120','122E'],$
    ytickname=['28','30','32','34','36N'],xminor=4,yminor=4,xcharsize=1.1,ycharsize=1.1
  s=scatterplot(lon,lat,SYMBOL='o',/sym_filled,rgb_table=32,MAGNITUDE=ave_n_15mm)
  china_map
  plot_map  
device,/close_file
  close,/all
报错：
% Internal error: Unable to set pixel format.
% IDLITWINDOW::ONEXPOSE: Unable to acquire device context.
% Unable to acquire device context.
数据说明：
lon,lat分别是469个站点的经纬度，float；ave_n_15mm是469个站点的降水量，也是float

秋之期

好吧，我倒腾几天以后，用循环的方法解决了，不算复杂，也许scatterplot这个函数的本质也是这样的，贴上我自己写的程序段，希望大家别走我的老路了，另外还是很期待有人能够解答一下scatterplot这个函数，大家也可以学习下嘛。
我的程序段：
  minlev=0. & maxlev=20. ;colorbar分级起始值
  levsize=2.
  nlev=10 ;colorbar分级数量
  levels=minlev+FindGen(nlev)*levsize
  clev=[11,13,15,17,19,20,21,22,23,25]
  ;title='Ave Annual Occurrence of Short Heavy Rain',
  ;-------------- 15mm ----------------------------------------
  loadct,0
  contour,ele_CN,lon_CN,lat_CN,pos=pos1,levels=tlevels,c_colors=tcolors,/fill,xrange=[112,122],yrange=[28,36],$
    title='15mm Short Heavy Rain',xtickname=['112','114','116','118','120','122E'],$
    ytickname=['28','30','32','34','36N'],xminor=4,yminor=4,xcharsize=1.4,ycharsize=1.4

  CS, SCALE=26
  for ilev=0,nlev-1 do begin
    ind15=where(ave_n_15mm ge minlev+ilev*levsize and ave_n_15mm lt minlev+(ilev+1)*levsize,n15)
    if(n15 gt 0)then oplot,lon[ind15],lat[ind15],psym=8,color=clev[ilev],symsize=1.5,thick=1.4
  endfor

秋之期

补充下成果图，没细调，凑合看吧

木木喵

楼主 请教一下色标棒是怎么画的啊

秋之期

木木喵 发表于 2017-4-8 10:41
楼主 请教一下色标棒是怎么画的啊

如果你问的是除了使用函数以外的方法的话……其实是笨方法，完全一点点画出来的。
    ;------------ color bar ----------
    CS, COLS=[224, 231, 331, 419,423], NCOLS=15
    ;CS, SCALE=38
    nlev=n_elements(levels)-1
    nybar = 2 & nxbar = nlev+1
    arr = FltArr(nxbar,nybar)
    yarrb = [0.,1.]
    xarrb = FindGen(nxbar)

    FOR ix = 0, nlev-1 DO BEGIN
      arr[ix,*] = levels[ix]
    ENDFOR
    dlev = levels[1] - levels[0]
    arr[ix,*] = levels[nlev-1]+dlev

    xtickn = StrArr(nxbar)
    xtickn=' '
    FOR i=0,nlev-1,2 DO BEGIN
      xtickn = STRING(levels,format='(f4.1)')
    ENDFOR
    xtickn[nlev] = ' '
    ;  xtickn[0]=' '
    blnk= StrArr(nxbar) & blnk = '  '
    contour,arr,xarrb,yarrb,pos=pos7,yrange=[0.,1.],$
      xrange=[min(xarrb),max(xarrb)],LEVEL=levels,/FILL, $
      yticks=1,xticks=nxbar-1,xminor=0,yminor=0,$
      xtitle='',ytitle='',ytickname=[' ',' '],$
      yticklen=0.001,xticklen=0.001,charsize=1.4,$
      xtickname=blnk,C_COLORS=clev,XSTYLE=9,/NOERASE
    AXIS,XAXIS=0,xrange=[min(xarrb),max(xarrb)], $
      xticks=nxbar-1,xminor=0, xtitle='',$
      xticklen=0.001,charsize=1.7,$
      xtickname=xtickn,XSTYLE=1,color=1
    AXIS,XAXIS=1,xrange=[min(xarrb),max(xarrb)], $
      xticks=nxbar-1,xminor=0, xtitle='',$
      xticklen=0.001,charsize=1.,$
      xtickname=blnk,XSTYLE=1,color=1
    AXIS,YAXIS=1,yticks=1,yminor=0, ytitle='',yticklen=0.001,charsize=1.,$
      ytickname=[' ',' '],YSTYLE=1,color=1
    AXIS,YAXIS=0,yticks=1,yminor=0, ytitle='',yticklen=0.001,charsize=1.,$
      ytickname=[' ',' '],YSTYLE=1,color=1