comp.lang.idl-pvwave archive - RDF feed
type: channel
description: Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado
items: Horizon plot - Problem using vis_horizon.pro function
items: Re: Horizon plot - Problem using vis_horizon.pro function

Horizon plot - Problem using vis_horizon.pro function
type: item
description: 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
creator: atmospheric physics
date: 2015-01-12T13:58:09-00:00

Re: Horizon plot - Problem using vis_horizon.pro function
type: item
description: On 1/12/15, 6:58 AM, Madhavan Bomidi wrote: > 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 su…
creator: Michael Galloy
date: 2015-01-12T20:27:11-00:00

Re: Horizon plot - Problem using vis_horizon.pro function
type: item
description: Hello Mike, Thanks and it works fine. Regards, Madhavan On Monday, January 12, 2015 at 9:27:13 PM UTC+1, Mike Galloy wrote: > On 1/12/15, 6:58 AM, Madhavan Bomidi wrote: >> 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=yt…
creator: atmospheric physics
date: 2015-01-14T10:42:46-00:00