surfnorm 더 효율적인 방법 Matlab

Question

포인트 클라우드를 구성한 후에 각 포인트의 법선을 얻고 자 할 때 내장 된 matlab 함수 surfnorm을 사용했지만 처리 시간이 많이 걸립니다. 그래서 누군가가 나를 도울 수 있다면 더 좋고 더 효율적인 방법입니다.

Answer 1

다음 코드가 도움이 될지 궁금합니다. 여기에는 세 단계가 있습니다.

1. 500 개 무작위 이격 된 지점 (X, Y)를 생성하고 I 기능
2. 리샘플링 랜덤 포인트 TriScatteredInterp 함수를 사용하여 같은 sinc을 선택하는 대응하는 z 값 (표면의 높이)를 계산 - 이렇게하면 초기 표면에 "대략 일치"하는 균등하게 샘플링 된 격자의 점을 얻을 수 있습니다.
3. 그리드의 "일부 점"에 대한 법선을 계산합니다 (480x640 점이 있기 때문에 모든 점에서 법선을 계산하면 생성됩니다). 엄청나게 조밀 한 "벡터 숲", "매 10 번째 포인트"를 샘플링하여 실제로 무엇을하는지 볼 수 있습니다.

다음과 같이 내가 사용하는 코드였다

:

randomX = rand(1,500); 
randomY = rand(1,500); 
r = 5*sqrt(randomX.^2 + randomY.^2); 
randomZ = sin(r) ./ r; 

% resample the data: 
[xx yy] = meshgrid(linspace(0,1,640), linspace(0,1,480)); 
F = TriScatteredInterp(randomX(:), randomY(:), randomZ(:)); 
zz = F(xx, yy); 

%% at each point, the normal is cross product of vectors to neighbors 
xyz=reshape([xx yy zz],[size(xx) 3]); 
xv = 10:30:479; yv = 10:30:639; % points at which to compute normals 
dx = xyz(xv, yv+1, :) - xyz(xv, yv, :); 
dy = xyz(xv+1, yv, :) - xyz(xv, yv, :); 

normVecs = cross(dx, dy); % here we compute the normals. 
normVecs = normVecs ./ repmat(sqrt(sum(normVecs.^2, 3)), [1 1 3]); 

figure; 
quiver3(xx(xv, yv), yy(xv, yv), zz(xv, yv), ... 
    normVecs(:,:,1), normVecs(:,:,2), normVecs(:,:,3)); 
axis equal 
view([56 22]); 

을 그리고 결과 플롯 :