n 차원 배열에서 선행 및 후행 NaN 값을 잘라내는 방법은 무엇입니까?

Question

이것은 예를 들어, 2 차원 쉽다 : NaN의 모든 NaN의 유지된다없는 행과 열로 값 것을

>> A = NaN(5,4) 
>> A(2:4,2:3) = [1 2; 3 4; 5 6] 
>> A(2,2) = NaN     
>> A(4,3) = NaN 

A = 

    NaN NaN NaN NaN 
    NaN NaN  2 NaN 
    NaN  3  4 NaN 
    NaN  5 NaN NaN 
    NaN NaN NaN NaN 

>> A(~all(isnan(A),2),~all(isnan(A),1)) 

ans = 

    NaN  2 
    3  4 
    5 NaN 

참고.

여러 치수로 확장하는 방법은 무엇입니까? 그때

ans = 

    NaN  2 
    3  4 
    5 NaN 

내가 네 가지 차원에서이 작업을 수행하고 싶습니다

, 여기에 예 행렬 A보다 훨씬 더 매트릭스로 얻을 어떻게

>> A = NaN(5,4,3)   
>> A(2:4,2:3,2) = [1 2; 3 4; 5 6] 
>> A(2,2,2) = NaN     
>> A(4,3,2) = NaN     

A(:,:,1) = 

    NaN NaN NaN NaN 
    NaN NaN NaN NaN 
    NaN NaN NaN NaN 
    NaN NaN NaN NaN 
    NaN NaN NaN NaN 

A(:,:,2) = 

    NaN NaN NaN NaN 
    NaN NaN  2 NaN 
    NaN  3  4 NaN 
    NaN  5 NaN NaN 
    NaN NaN NaN NaN 


A(:,:,3) = 

    NaN NaN NaN NaN 
    NaN NaN NaN NaN 
    NaN NaN NaN NaN 
    NaN NaN NaN NaN 
    NaN NaN NaN NaN 

예를 들어 A는 세 가지 차원이있는 경우.

Answer 1

입력 A에 따라 문제에 대한 내 솔루션 영업 이익에 의해 게시 된 : 당신은 당신이 내가의 성능에 대한 확실하지 않다 훨씬 더 큰 행렬이 작업을 수행 할 수 있음을 언급 한 이후

>> [i,j,k] = ind2sub(size(A),find(~isnan(A))); 
>> l = min([i j k]); 
>> u = max([i j k]); 
>> B=A(l(1):u(1),l(2):u(2),l(3):u(3)) 
B = 
    NaN  2 
    3  4 
    5 NaN 
>> size(B) 
ans = 
    3  2 

@ ronalchn의 솔루션 -이 모든 것이 all- 호출입니다. 그러나 나는 그 문제를 확장하는 것에 대해 아무 것도 모른다. 어쩌면 누군가가 주석을 달 수있다.

Answer 2

이 시도 :

2 차원

A(~all(isnan(A),2),~all(isnan(A),1)) 

3 차원

A(~all(all(isnan(A),2),3),... 
    ~all(all(isnan(A),1),3),... 
    ~all(all(isnan(A),1),2)) 

4 차원

01 2,311,

기본적으로 규칙은 N 차원입니다 : 당신이 때는 isNaN() 일을 모든 N 크기에

• 합니다. ,

• 는 그 후, 전체() 함수는 N-1 회

• 그것을에 감아 제 2 인수는 i 번째 차원의 전체() 함수의 각각은 임의의 순서로 N에 숫자 1이되어야 그러나 i 제외.

  ---

  테오 드 로스 테 Zelleke가 그 방법을 빠른 (그는 그의 방법은 너무 빨리하다고 생각 말하는 좋은 방법)입니다보고 싶어하기 때문에

, 여기에 벤치 마크입니다.행렬 A를 같이 정의 : 결과와

tic; 
for i=1:100 
isnanA=isnan(A); 
ai34=all(all(isnanA,3),4); 
ai12=all(all(isnanA,1),2); 
B=A(~all(ai34,2),~all(ai34,1),~all(ai12,4),~all(ai12,3)); 
end 
toc 

: 정의

Elapsed time is 0.533932 seconds. 
Elapsed time is 0.519216 seconds. 
Elapsed time is 0.575037 seconds. 
Elapsed time is 0.525000 seconds. 

내 테스트 : 결과와

tic; 
for i=1:100 
[i,j,k,z] = ind2sub(size(A),find(~isnan(A))); 
l = min([i j k z]); 
u = max([i j k z]); 
B=A(l(1):u(1),l(2):u(2),l(3):u(3),l(4):u(4)); 
end 
toc 

:

A = NaN*ones(100,400,3,3); 
A(2:4,2:3,2,2) = [1 2; 3 4; 5 6]; 
A(2,2,2,2) = NaN;A(4,3,2,2) = NaN; 
A(5:80,4:200,2,2)=ones(76,197); 

그의 테스트로 정의

Elapsed time is 0.224869 seconds. 
Elapsed time is 0.225132 seconds. 
Elapsed time is 0.246762 seconds. 
Elapsed time is 0.236989 seconds.