how do I extract time series?

28 ビュー (過去 30 日間)
Lilya
Lilya 2020 年 3 月 31 日
コメント済み: Lilya 2020 年 4 月 1 日
Hi all, I am trying to extract the data from the nc list of files using the following loop.
1- extract the whole parameter
2- set them in a 3d matrix
3- calculate the mean
4- extract the data in a given domain
5- calculate the mean to have the time series
I am not convinced with the final result (useries)
any help would be appreciated.
Thank you
useries = [];
for day=1:length(datenum_array);
filename=strcat('./MERRA2_400.tavg1_2d_ocn_Nx.2019',num2str(datevec_array(day,2),'%02i'),num2str(datevec_array(day,3),'%02i'),'.SUB.nc');
% disp(filename)
time = ncread(filename,'time');
lon = ncread(filename,'lon');
lat = ncread(filename,'lat');
%this box domain
lat_range=find(lat<= 28 & lat >= 26);
lon_range=find(lon<= 36.5 & lon >= 33);
time_range = datenum_array(32:12);
%read U10
ucom = ncread(filename,'U10M');
umean = nanmean(ucom,3);
U10 (:,:,day) = umean;
Ubox = U10(lon_range,lat_range,time_range);
umean2 = nanmean(nanmean(nanmean(Ubox)));
useries = [useries,umean2];
end

採用された回答

dpb
dpb 2020 年 3 月 31 日
>> which time
C:\ML_R2019b\toolbox\matlab\bigdata\@tall\time.m % tall method
>> which day
C:\ML_R2019b\toolbox\matlab\timefun\@datetime\datetime.m % datetime method
>>
Better to not alias builtin/ML-supplied functions...changed names below.
useries = [];
lat_low=26; lat_up=28; % don't bury data in code--use variables so can change easily
lon_low=33; lon_up=36.5; % ditto
t_low=???; t_up=???; % dunno here but datenum_array(32:12) will return empty array reference
No way to tell what you really want for the time ranges but
time_range = datenum_array(32:12);
will return a null array as the indexing colon expression is equivalent to 32:1:12 which is an increasing difference of +1 over a decreasing interval which is empty. That is undoubtedly the biggest cause of your problems.
Set a time value and use logical indexing within in similar as demonstrated for lat, lon except using datetime
for d=1:size(datenum_array,1); % length unreliable as is max(size()), not specific dimension
filename=sprintf('./MERRA2_400.tavg1_2d_ocn_Nx.2019%02d%02d.SUB.nc',datevec_array(d,2:3)); % somewhat simpler
t=ncread(filename,'time'); % above comment re: builtin function names
lon = ncread(filename,'lon');
lat = ncread(filename,'lat');
%this box domain
lat_range=iswithin(lat,lat_low,lat_up);
lon_range=iswithin(lon,lon_low,lon_up);
time_range=iswithin(datenum_array,t_low,t_up); % t_low, t_up datetime() values begin, end desired
%read U10
ucom = ncread(filename,'U10M');
umean = nanmean(ucom,3);
U10(:,:,d) = umean;
Ubox = U10(lon_range,lat_range,time_range);
umean2 = nanmean(nanmean(nanmean(Ubox)));
useries = [useries,umean2];
end
Above presumes all files have same lat, lon and time data so the result of the locations will be same length--otherwise, the concatenation will fail on mismatched sizes.
Would be more efficient to determine that number first and then preallocate and store, but presuming the number of files is relatively small, the extra overhead of dynamic allocation may not be excessive and worth the extra coding.
iswithin is my "syntactic sugar" utility routine:
function flg=iswithin(x,lo,hi)
% returns T for values within range of input
% SYNTAX:
% [log] = iswithin(x,lo,hi)
% returns T for x between lo and hi values, inclusive
flg= (x>=lo) & (x<=hi);
which just moves the logical into a function for somewhat cleaner code at the user level.
  1 件のコメント
Lilya
Lilya 2020 年 4 月 1 日
This is great!!..
The data is huge, so I couldn't attach them. But I really appreciate that you get the idea and help with the problem.
Thank you.

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeCreating and Concatenating Matrices についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by