[SciPy-user] converting hourly series to annual unneccessaryly masks data

Timmie timmichelsen at gmx-topmail.de
Thu Jan 8 13:49:15 EST 2009


Hello Pierre,

> The documentation is still a bit lacking, sorry. Still, in the  
> docstring of convert, you can see that if you don't precise a func  
> input parameter, the series is converted to 2D, as stated:
> `
>          If ``func`` is not given, the output series group the points  
> of the
>          initial series that share the same new date. For example, if  
> the
>          initial series has a daily frequency and is 1D, the output  
> series is
>          2D.
No problem here.
We discussedit already here:
aggregation of long-term time series
http://article.gmane.org/gmane.comp.python.scientific.user/15584

 
> 1. You're not starting at 01/01:00-00, but 8 days later
Yes, I am aware of it.

> 2. We are using this 366d year: as there are no leap year in your  
> range of years, the last 24 data of each line will be masked.
This explains what I was looking for.
Because it affects how I handle the data later.

I need averages for all hours over the years:

atest.mean(0)
=> this the data array for the new one-year hourly time series (8760 h).
And since the data is masked at the and, I am lacking a day when I build the
timeseries.
Is there a way to handle this generically?

I mean if my long-term years contain a leap year I neeed the masked points but
normally not.

How would you suggest to build the one-year hourly average time series in a
flexible way?

A example case what I am aiming at:
Averge hourly temperatures over 20 years of data.

> 3. You don't finish at 12/31-23:00, but (365-8) days earlier.
I also know this here.

 
> So all is well and works as expected (developer-wise), no need for a  
> ticket (good reflex, though).
Sorry, too fast.
 
> Now, of course, you need to tell us what you were expecting, and what  
> kind of average you wanted to calculate.
See above.




More information about the SciPy-User mailing list