Average of irregular range of data

2 次查看(过去 30 天)
David
David 2013-8-15
I have a time series of satellite data over a bout 6 years. With each pass of the satellite a number of measurements are taken. The number of measurements taken is non-uniform, so anywhere between 2 and 7 for a given pass of the satellite. This occurs every 10 days but there are occasions when the satellite misses a pass. I want to get the average of each pass without having to go through the data and inputting the correct ranges manually. Is there a way of setting a rule so that when the time between two measurements is greater than say a day a new range of data is selected? Any help would be appreciated.
  2 个评论
Andrei Bobrov
Andrei Bobrov 2013-8-15
编辑:Andrei Bobrov 2013-8-15
short example of your data, please
David
David 2013-8-15
Hi Andrei, here are the first few passes of the satellite, with the figure on the right being the wave height:
09/01/2003 21:05 1.73
09/01/2003 21:05 2.196
09/01/2003 21:05 1.979
19/01/2003 19:03 6.507
19/01/2003 19:03 6.676
29/01/2003 17:02 6.132
29/01/2003 17:02 6.379
29/01/2003 17:02 5.714
08/02/2003 15:00 5.019
08/02/2003 15:00 4.671
08/02/2003 15:00 4.814
08/02/2003 15:00 4.783
08/02/2003 15:00 4.334

请先登录,再进行评论。

回答(2 个)

Azzi Abdelmalek
Azzi Abdelmalek 2013-8-15
编辑:Azzi Abdelmalek 2013-8-15
data={'09/01/2003 21:05 1.73'
'09/01/2003 21:05 2.196'
'09/01/2003 21:05 1.979'
'19/01/2003 19:03 6.507'
'19/01/2003 19:03 6.676'
'29/01/2003 17:02 6.132'
'29/01/2003 17:02 6.379'
'29/01/2003 17:02 5.714'
'08/02/2003 15:00 5.019'
'08/02/2003 15:00 4.671'
'08/02/2003 15:00 4.814'
'08/02/2003 15:00 4.783'
'08/02/2003 15:00 4.334'}
dates= datenum(data,'dd/mm/yyyy');
weights=cellfun(@(x) str2double(x{3}),regexp(data,'\s+','split'));
[a,b,c]=unique(dates);
average_values=accumarray(c,weights,[],@mean);
out=[cellstr(datestr(dates(b))) num2cell(average_values)]
%or
dates=datenum(data,'dd/mm/yyyy');
s=regexp(data,'\s+','split');
s=[s{:}];
weights=str2double(s(3:3:end))';
[a,b,c]=unique(dates);
average_values=accumarray(c,weights,[],@mean);
out=[cellstr(datestr(dates(b))) num2cell(average_values)];

Andrei Bobrov
Andrei Bobrov 2013-8-15
编辑:Andrei Bobrov 2013-8-15
tresh = 2; % Let greate than two days
f = fopen('PathAndNameYourdata.txt');%'E:\wtdy2.txt'
c = textscan(f,'%s %s %f','CollectOutput',true);
fclose(f);
k = strcat(c{1}(:,1),{' '},c{1}(:,2));
n = datenum(k,'dd/mm/yyyy HH:MM');
ii = [true;diff(n)>tresh];
i1 = cumsum(ii);
i2 = find(ii);
out = [{'start','finish','mean of value'};
[k([i2,[i2(2:end)-1;numel(n)]]), accumarray(i1,c{2},[],@(x){mean(x)})]];
  3 个评论
David
David 2013-8-15
Just changed it around a little and it works perfectly. Thanks for the help Andrei, much appreciated.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Simulink 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by