【心碎】有谁知道这种支离破碎的图问题在哪？？

15195775117

前提：con是个2*350000的数组，里面有2列，第一列是高度，第二列是雷达dbz的值
我想统计下不同高度和dbz值对应的频数
即，途中横坐标是dbz,纵坐标是高度ht,中间等值线是频数
我取的ht的范围是0~15km，dbz的值范围-30~20
我设定connum数组来放频数，起始时是全零阵
connum的横坐标的序号对应fix((dbz+30)*2)
connum的纵坐标的序号对应fix(ht*10)
比如，（7.97km，-18.94）对应着connum中（79,22）的位置
算法中，对应（79,22）位置的情况每出现一次，该处的connum值就+1
connum=fltarr(102,152)
      for j=0,s-1 do begin
        if(con[0,j] gt 0)then begin
          connum[fix((con[1,j]+30)*2),fix(con[0,j]*10)]++
        endif
      endfor
;以下这段将ht和dbz的真实值还原，不会影响图的形态
      x=fltarr(102,152)
      y=fltarr(102,152)
      z=fltarr(102,152)
      for i=0,100 do begin
        for j=0,150 do begin
          x[i,j]=i/2-30
          y[i,j]=j/10
          z[i,j]=connum[i,j]
        endfor
      endfor
画图：
      contour,z,x,y,ystyle=1,xstyle=1,xtitle='dbz',xrange=[-30,20],ytitle='ht/km',nlevels=40,/fill,title=figurename[a]
      colorbar,ncolors=250,division=10,range=[0,10],/vertical,/right,pos=[0.95,0.2,0.97,0.8]
但是如果我把connum设为52*15，即1km平均以下，图像将是这样：

我的要去就是图像平滑一点，其他的我可以自己该，没见过这么破碎的图啊，救命！

longlivehj

还原的时候，i/2，j/10出来都是整数，改成i/2.和j/10.试试！

15195775117

longlivehj 发表于 2015-4-6 08:08
还原的时候，i/2，j/10出来都是整数，改成i/2.和j/10.试试！

对哦，大神果然是火眼金睛

舒服多了，不过还是有横向拉丝的样子，大神以为，这是何故？？