数组统计

wzf171

idl中有没有提供对数组中各个数据出现的频次进行统计的函数？

chenhongfengye

我建议你可以使用where函数

925450204

我有个问题为什么我的循环进行不下去？

;data=make_array(144,73,17,718,/float)
cdfid=ncdf_open('atmprfzhuan.nc')
varid2=ncdf_varid(cdfid,'refrac')
varid1=ncdf_varid(cdfid,'alt_refrac')
ncdf_varget,cdfid,varid2,atm1
ncdf_varget,cdfid,varid1,atm2
;ncdf_attget,cdfid,varid,'scale_factor',a   ;系数
;ncdf_attget,cdfid,varid,'add_offset',b      ;偏值
ncdf_close,cdfid

help,atm1,atm2

;openw,lun,'d:\1.dat',/get_lun
;writeu,lun,data1,data2
;free_lun,lun

;print,data[*,*,1,1]
;print,atm1[0],atm2[0]

;print,atm2,ecm2

;a=replicate(0.000,50)
;for i=1L,50L,1L do begin
;a(i)=i*1000
;endfor

;print,a
cdfid=ncdf_open('ucar1.nc')
varid2=ncdf_varid(cdfid,'refrac')
varid1=ncdf_varid(cdfid,'alt_refrac')
ncdf_varget,cdfid,varid2,ropp1
ncdf_varget,cdfid,varid1,ropp2
ncdf_close,cdfid

help,atm1,atm2
a=[0,2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36,38,40,42,44,46,48,50]
atmrec=interpol(atm1,atm2/1000,a)
ropprec=interpol(ropp1,ropp2/1000,a)
;print,atmrec,ecmrec
bias=atmrec-ropprec
print,bias
;print,b
plot,bias,a,xtitle="bias",ytitle="alt(km)",xrange=[-10,10],yrange=[0,50],$
title='refrac COSMIC-ROPP',ycharsize=0.85,xcharsize=0.85
plotS,[0,0],[0,50],LineStyle=2,color='0000FF'XL
Window, size=400,xsize=400
LoadCT, 3
image = LoadData(7)
TVImage, image, Position=[0.1, 0.1, 0.75, 0.9]
Colorbar, Position=[0.8, 0.1, 0.86, 0.9], /Right, $
/Vertical, Division=5, Format='(F5.1)'
image = TVRD()
image = BytArr(3, 400, 400)
TVLCT, r, g, b, /Get
image3D[0, *,*]=r[image]
image3D[1, *,*] =g [image]
image3D[2, *, *] = b[image]
image3D=TVRD(true=1)

    Device, Decomposed=1
    !P.Background='FFFFFF'XL
    !P.Color='000000'XL

Write_JPEG,' COSMIC-ROPP.jpg', image3D, true=1, quality=75
end