IDL之构造垂直风速色条

15195775117

垂直风速，下沉为蓝，上升为红，
如果红蓝函数线性配色，红蓝过渡的黑色段很长，
于是我想着把过渡区的蓝红曲线做陡一些，这些过渡区就短些

pro vertical_wind
a=findgen(1000)
vc=vcolor(a)
help,vc
fig=plot(a,a,vert_colors=vc,thick=10)
end
;构造垂直风色条
;根据向量x从纯蓝到纯红进行颜色匹配
;返回值是3列匹配好的RGB颜色
function vcolor,x
n=n_elements(x)
n2=fix(n/2)
s=min(x)
t=max(x)
m=(s+t)/2
a=-255.0/((m-s)^2)
b=-2.0*s*a
c=0.0-a*m*m-b*m
bb=[a*x[0:n2-1]^2+b*x[0:n2-1]+c,fltarr(n-n2)]
a=-255.0/((m-t)^2)
b=-2.0*a*t
c=0.0-a*m*m-b*m
rr=[fltarr(n2),a*(x[n2:-1])^2+b*x[n2:-1]+c]
return,transpose([[rr],[fltarr(n)],[bb]])
end