Функция surfnorm более эффективный способ Matlab

Question

После построения облака точек я хочу получить нормальное значение для каждой точки, и я использовал встроенную функцию surflorm для функции matlab, но для этого требуется много времени обработки. Поэтому, если бы кто-нибудь мог помочь мне сделать это лучше и эффективнее.

Answer 1

Интересно, поможет ли вам следующий код. Здесь есть три шага.

1. Создать 500 случайно расположенные точки (х, у), и вычислить соответствующее значение г (высоту поверхности), для которого я выбрал sinc как функции
2. ресэмплировать случайные точки с использованием функции TriScatteredInterp - это позволяет мне получать точки на сетке с равномерным отображением, которые «грубо соответствуют» исходной поверхности.
3. Вычислить нормаль к «некоторым точкам» на этой сетке (так как есть 480x640 точек, вычисление нормали в каждой точке будет просто создавать невероятно плотный «лес векторов», путем выборки «каждый десятый пункт», вы действительно можете видеть, что делаете

код я был следующим:

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]); 

И полученный участок: