Find the last non-nan in large array
조회 수: 14 (최근 30일)
이전 댓글 표시
I have very large 4-dimensional arrays with dimensions being time,depth, lat, lon. I want to find the linear index in the depth dimension to the last non-nan data value to I can evaluate the values at the "bottom". I am looking for the output to be 3-D (time, lat, lon) with linear indexes to the last depth with a value (non-nan).
Any idea how to do this efficiently? I have a bunch of very large files to run this on.
댓글 수: 0
채택된 답변
Guillaume
2019년 5월 30일
One way:
%m: a 4D matrix (time x depth x lat x lon
[~, idx] = max(cumsum(~isnan(m), 2), [], 2)
This relies on the fact that max returns the index of the first max value if several are identical. The first max value of cumsum(~isnan) is the last number before zero or more nan.
댓글 수: 3
Pepe Grillo
2021년 2월 24일
Hi, I was looking exactly the same, but when I run this code I get a 4D array and when plt is the number of values before the last NaN, but I want to plot the value,
what I am doing wrong?
[~, idx] = max(cumsum(~isnan(m), 3), [], 3) %??
Thanks
Walter Roberson
2021년 2월 24일
I am not clear as to what size your m is ? You should expect the output, idx, to be the same size as m, except with the third dimension now scalar.
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Matrix Indexing에 대해 자세히 알아보기
제품
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!