Igor Pro 日变化图 Diurnal Plot

点为中位数，柱为1/4和3/4分位数

 

Function MakeDiurnalCycle(TimeWave,DataWave)
wave TimeWave, DataWave
variable i
Make/O/N=25 HourofDay
HourofDay=p
duplicate/O DataWave temp
duplicate/O HourofDay Data_Q1
duplicate/O HourofDay Data_Q2
duplicate/O HourofDay Data_Q3

for (i=0;i<24;i+=1)
	temp = (floor(mod(TimeWave,86400)/3600)==i ? DataWave : NaN)
	StatsQuantiles/iNaN temp
	Data_Q1[i]=V_Q25
	Data_Q2[i]=V_Median
	Data_Q3[i]=V_Q75
endfor
Data_Q1[24]=Data_Q1[0]
Data_Q2[24]=Data_Q2[0]
Data_Q3[24]=Data_Q3[0]
KillWaves Temp
//Make graph
Display Data_Q3 vs HourofDay
AppendtoGraph Data_Q1 vs HourofDay
ModifyGraph mode=5,useBarStrokeRGB=1,hbFill(Data_Q3)=4,hbFill(Data_Q1)=1,rgb(Data_Q1)=(65535,65535,65535),barStrokeRGB=(65535,65535,65535)
AppendtoGraph Data_Q2 vs HourofDay
ModifyGraph mode(Data_Q2)=4,lsize(Data_Q2)=2,marker(Data_Q2)=19,msize(Data_Q2)=4,offset(Data_Q2)={0.5,0}
ModifyGraph mirror=2,fSize=16
ModifyGraph manTick(bottom)={0,4,0,0},manMinor(bottom)={3,50}
Label bottom "Hour of Day"
End