how to take out specific month from data

Question

ravi 2014-4-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/125177-how-to-take-out-specific-month-from-data

评论： dpb 2019-6-24

sir i have 95 year daily model output rainfall data form 2006-1-1 to 2100-12-31.i want to take out nov to feb which include leap years. i have program for this.this is working perfectly for who does not include leap year but i want to include leap years then please suggest me what will change in this program for this. in between 2006 to 2100 have 24 leap years but my data include only 23 leap years,this is not not include end(last) leap year. this is daily data.sir my precipitaion data size is ( 20x22x34698) (lat,lon.time) for india region.sir please help me. this is include leap years.please help me.thank you

dten = datenum([2006 1 1 ;2100 12 31])

[~,M] = datevec(dten(1):dten(2));

xx = find(ismember(M,[11 12 1 2]));

[m,n,k] = size(prec_all_years1);

ii = cumsum(eomday(2014,1:12));

x0 = rem(0:k-1,365)+1;

xx_logic = x0 <= ii(2) | x0 > ii(10);

p = prec_all_years1(:,:,xx_logic);

max_rainfall = max(reshape(p,m*n,[]));

4 个评论
显示 2更早的评论隐藏 2更早的评论

dpb 2014-4-10

在 MATLAB Online 中打开

What do you mean by "take out", precisely? You mean you want to remove (or selectively keep) those months' data and ignore the rest? Your M variable and ismember will find those. If you're looking to find the leap years and either keep or discard only them, then do the same kind of thing keeping Y as well and build the logical vector of which year is a leap year for the addressing. That's easily done by

isLpYr=(eomday(Y,2)==29); % logical vector T for Y==leap year

ravi 2014-4-10

编辑：ravi 2014-4-10

it's mean i want to make another data set for nov to feb from 12 month data (12x16x34675 (lat,lon.time) for india region who include leap years).i want to cut only nov to feb my data set but this program is right only for those who does not include leap year then suggest me what will change in this program for my requirement.please help me,thank you

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

dpb 2014-4-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/125177-how-to-take-out-specific-month-from-data#answer_132934

编辑：dpb 2014-4-11

在 MATLAB Online 中打开

As noted above, you were already there--

dten = datenum([2006 1 1 ;2100 12 31]);   
[~,M] = datevec(dten(1):dten(2));
p = prec_all_years1(:,:,ismember(M,[11 12 1 2]));

IOW, just throw away all the leap year stuff if don't want it. But if want to go back, the preceding logical vector is much simpler way to select/not select those years than the initial logic.

[Y,M] = datevec(dten(1):dten(2));
p = prec_all_years1(:,:,ismember(M,[11 12 1 2]) & ~(eomday(Y,2)==29));

will select Nov-Feb for the non-leap years in the range of years Y

Look up Using Logicals in Array Indexing in the documentation for more details.

14 个评论
显示 12更早的评论隐藏 12更早的评论

dpb 2014-4-10

编辑：dpb 2014-4-10

在 MATLAB Online 中打开

A) It is, after all, your program/work/job/research so it shouldn't be too much to expect you put in some effort to debug/fix any donated suggestions, and

B) What does "is not working" mean, specifically? I can't see your terminal from here.

One observation -- 95*365 --> 34675

whereas

length(datevec(dten(1):dten(2))) --> 34698

It thus appears that your model doesn't know anything about leap years? Is this the problem that your resulting vectors aren't the commensurate length? If so, either fix the calendar in the model to generate full time histories or remove the leap-year day from the time vector you've generated.

ADDENDUM:

Again on the assumption it's the above length mismatch that's the problem and that you'd rather use the data from the model as is rather than fix its broken time stamp, one way to proceed would be to eliminate the leap days from the M vector. The simplest way to do this means returning the day as well as month but don't really need the year altho I returned it anyway on general principles...

>> [Y,M,D] = datevec(dten(1):dten(2));  % create the date vector over period
>> ix29= (M==2) & (D==29);              % find Feb 29 locations
>> sum(ix29)    % just observe how many there are
ans =
  23
>> Y(ix29)=[];M(ix29)=[];D(ix29)=[];    % remove those from Y,M,D 
>> length(M)          % and observe resulting length is same as your data
ans =
     34675

Now just select for the months desired as before; the length of the logical vector will now match the number of planes in the data array. Since we've previously eliminated the leap days, now no longer any need for it in the selection logic

p = prec_all_years1(:,:,ismember(M,[11 12 1 2]));

dpb 2014-4-10

编辑：dpb 2014-4-11

在 MATLAB Online 中打开

Your precise words have been "sir my precipitaion data size is 12x16x34675". That's not enough to include the leap day for those 24 years (of which as noted before, one is NOT a leap year). That would require as noted above an array size of 34698 planes.

The selection I gave earlier that you claimed "does not work" does either depending on what you do want but the vector sizes are not going to be commensurate unless the data are consistent.

You need to show explicitly what "doesn't work" including the error message in context of the code you tried.

It would be good if you also included the result of

size(prec_all_years1)
size(M)

for absolute confirmation of what we're dealing with. We need hard data/code to debug with; not descriptions as we can't see your terminal or know what you know inherently and the language barrier doesn't help so it's got to be something in a common unambiguous language and that's Matlab code and data.

ADDENDUM:

And, what do you mean by "cut"? Is that what you want to retrieve or not retrieve?

dpb 2014-4-11

编辑：dpb 2014-4-11

在 MATLAB Online 中打开

OK, now we're getting somewhere...so, indeed, your raw data is 34698 in length...in that case the first response above is all you need. What you don't seem to recognize is the vectorization of operations in Matlab. The result of

ismember(M,[11 12 1 2])

is a logical array of the length( M) (commensurate with the number of planes in prec_all_years) which is T for the months in the selection vector. Logical addressing in Matlab will return the entries in the array for the dimension in which that vector is used based on whether the value is T or F; hence there's no need to loop over the year--that's what the vector does automagically.

By putting it in a loop, you've done the same thing 95 times but have also included the code to exclude the leap year which isn't needed since you do want those years as well. All you need is the one line

p = prec_all_years(:,:,ismember(M,[11 12 1 2]));

Note if I just use the time vector M at the command line--

sum(ismember(M,[11 12 1 2]))
ans =
     11423
>> 95*120+23
ans =
     11423
>>

the sum() of the result of ismember is precisely the number of days in 95 years of the months of [N D J F]=[30+31+31+28 -->120]*95+23. Again, remember that 2100 is not a leap year. The sum of the logical vector is the number of T elements in it so you've got exactly the right number selected.

Oh, and while you're not wanting to exclude the leap years after all, the reason you got no results from your selection is that I had a misplaced parens in the exclusion condition on the leap year flag--

(~eomday(Y,2)==29) should have been ~(eomday(Y,2)==29)

to negate the overall test vector, not just the function value. As it was written, it returns no T values at all.

So, the exclusion condition if it were to be used would be

p = prec_all_years(:,:,ismember(M,[11 12 1 2]) & ~(eomday(Y,2)==29))

We can also test this easily enough at the command line --

>> sum(ismember(M,[11 12 1 2]) & ~(eomday(Y,2)==29))
ans =
      8640
>> (95-23)*120
ans =
      8640
>>

Poulomi Ganguli 2019-6-24

编辑：Poulomi Ganguli 2019-6-24

在 MATLAB Online 中打开

Hello:

I came across the same issue. My data are arranged in a matrix containing year, month, day and value. If I run this code, I am encountering foll. error message:

The logical indices in position 3 contain a true value outside of the array bounds

dpb 2019-6-24

Well, the above code is for retrieving from a 3D array, not 2D...you'll need to write the proper expression to match your different storage arrangement. But, can't see your terminal from here so have no definite knowledge to go on.

请先登录，再进行评论。

Answer 2

dpb 2014-4-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/125177-how-to-take-out-specific-month-from-data#answer_132869

2100 is NOT a leap year; it's the "divided by 100" exception of the "divisible by four" rule. Then there's the 400 yr correction as well but we'll likely not be seeing it.

I am wondering, however, about how you've managed to accumulate daily precipitation data for years from 2015-2100 as yet...