Dear All,
Greetings. I was trying to make horizon plot to represent wavelet multi-resolution analysis (MRA) of irradiance with each row denoting the different detail of the MRA. Following Michael Galloy's vis_horizon.pro function (http://michaelgalloy.com/2009/05/05/horizon-graph-code.html), I tried to make the visualization code for my data. I found that I could not get anything except a black figure window. Will it be possible for anyone to correct me if I am using vis_horizon function wrongly?
My ASCII input file contains the following columns:
[ UTC_time, raw_data(1s), wd_5s, wd_10s, wd_20s, wd_40s, wd_1m20s, $
wd_2m40s, wd_5m20s, wd_10m40s, wd_21m20s, wd_42m40s, $
wd_1h25m20s, wd_2h50m40s ]
The first column represent the UTC time (in hours), the second column represents the raw irradiance data. From third column onwards, the each column represent the wavelet detail for different smoothing scales. My intention was to represent the wavelet details in the form of horizon plot similar to the figure shown in page 6 (http://oa.upm.es/4953/1/Perpinan.Lorenzo2010.pdf).
I have written the following lines in the IDL code:
-----------------------------------------------------
PRO HORIZONPLOT
infile = 'pyr43_rsds_wdj_20130413.txt'
nrows = FILE_LINES(infile)
allData = FLTARR(14,nrows)
OPENR, lun0, wdjfile, /GET_LUN
READF, lun0, allData
CLOSE, lun0 & FREE_LUN, lun0
utcTime = REFORM(allData[0,*]) ; UTC Time (hours)
ws0 = REFORM(allData[1,*]) ; Raw irradiance (@1 sec)
wdj = allData[2:13,*] ; Wavelet Details
ytitles = ['5s','10s','20s','40s','1m20s','2m40s','5m20s','10m40s','21 m20s', $
'42m40s','1h25m20s','2h50m40s']
minval = MIN(wdj,MAX=maxval)
!P.Multi=0
vis_horizon,utcTime, wdj, nbands=12,titles=ytitles, $
xstyle=9, ystyle=8, min=minval, max=maxval, colors=bytarr(12)
END
---------------------------------------------------
Look forward for your suggestions.
Thanking you in advance,
Madhavan
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
