| Horizon plot - Problem using vis_horizon.pro function [message #89974] |
Mon, 12 January 2015 05:58  |
atmospheric physics
Messages: 121
Registered: June 2010
|
Senior Member |
|
|
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
|
|
|
|
| Re: Horizon plot - Problem using vis_horizon.pro function [message #89980 is a reply to message #89974] |
Mon, 12 January 2015 12:27  |
Michael Galloy
Messages: 1114
Registered: April 2006
|
Senior Member |
|
|
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 suggestions.
>
> Thanking you in advance,
> Madhavan
>
I've update MG_HORIZON to handle NaN values like you have in your data.
Get updates from the GitHub repo:
http://github.com/mgalloy/mglib
Mike
--
Michael Galloy
www.michaelgalloy.com
Modern IDL: A Guide to IDL Programming (http://modernidl.idldev.com)
Research Mathematician
Tech-X Corporation
|
|
|
|
| Re: Horizon plot - Problem using vis_horizon.pro function [message #89999 is a reply to message #89980] |
Wed, 14 January 2015 02:42 |
atmospheric physics
Messages: 121
Registered: June 2010
|
Senior Member |
|
|
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=ytitles, $
>> xstyle=9, ystyle=8, min=minval, max=maxval, colors=bytarr(12)
>>
>> END
>> ---------------------------------------------------
>>
>> Look forward for your suggestions.
>>
>> Thanking you in advance,
>> Madhavan
>>
>
> I've update MG_HORIZON to handle NaN values like you have in your data.
> Get updates from the GitHub repo:
>
> http://github.com/mgalloy/mglib
>
> Mike
> --
> Michael Galloy
> www.michaelgalloy.com
> Modern IDL: A Guide to IDL Programming (http://modernidl.idldev.com)
> Research Mathematician
> Tech-X Corporation
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
