자료의 구간이 등간격이 아니어서 hist_2d()를 사용할 수 있는 상황은 아닌 것 같습니다. 2차원 구조의 각 계급구간에 대해 COUNT를 수행하는 과정이 필요할 것 같습니다.
bd=[0, 1, 3, 9, 20]
bv=[0, 0.4, 0.9, 2]
nbd=n_elements(bd)
nbv=n_elements(bv)
count=lonarr(nbd-1, nbv-1)
;make sample data
;
d=randomn(seed, 1000)*20 ;0~20
v=randomn(seed, 1000)*2 ;0~2
for i=0, nbd-2 do begin
for j=0, nbv-2 do begin
ok=where(d ge bd[i] and d lt bd[i+1] and v ge bv[j] and v lt [j+1], ct)
count[i, j] = ct
endfor
endfor
end
이후에 같은 계급 구간으로 box를 그려가면서 box 내부를 도수에 따라 색으로 채워 넣으면 될 것 같습니다.