Загрузите набор 3-D координатных точек в рабочую область.

load('xyzPoints.mat');

Создайте объект облака точек.

ptCloud = pointCloud(xyzPoints);

Укажите точку запроса и радиус, в пределах которого должны быть идентифицированы соседи.

point = [0,0,3];
radius = 0.5;

Получите индексы и расстояния точек, которые находятся в заданном радиусе.

[indices,dists] = findNeighborsInRadius(ptCloud,point,radius);

Получите данные облака точек радиальных соседей.

ptCloudB = select(ptCloud,indices);

Отобразите облако точек. Постройте график точки запроса и соответствующих радиальных соседей.

figure
pcshow(ptCloud)
hold on
plot3(point(1),point(2),point(3),'*')
pcshow(ptCloudB.Location,'r')
legend('Point Cloud','Query Point','Radial Neighbors','Location','southoutside','Color',[1 1 1])
hold off

Найдите радиальные соседи точки запроса в данных организованного облака точек с помощью матрицы проекции камеры. Вычислите матрицу проекции камеры из выборочных точек данных облака точек и их соответствующих координат точки изображения.

Загрузите организованные данные облака точек в рабочую область. Облако точек генерируется с помощью датчика глубины Kinect.

ld = load('object3d.mat');
ptCloud = ld.ptCloud;

Укажите размер шага для выборки данных облака точек.

stepSize = 100;

Дискретизируйте входное облако точек и сохраните выбранные координаты точек 3-D как объект облака точек.

indices = 1:stepSize:ptCloud.Count;
tempPtCloud = select(ptCloud,indices);

Удалите недопустимые точки из выборочного облака точек.

[tempPtCloud,validIndices] = removeInvalidPoints(tempPtCloud);

Задайте 3-D мировых координат точки входного облака точек.

worldPoints = tempPtCloud.Location;

Найдите 2-D координаты изображения, соответствующие 3-D координатам точки входного облака точек.

[Y,X] = ind2sub([size(ptCloud.Location,1),size(ptCloud.Location,2)],indices);
imagePoints = [X(validIndices)' Y(validIndices)'];

Оцените матрицу проекции камеры из изображения и мировых координат точки.

camMatrix = estimateCameraMatrix(imagePoints,worldPoints);

Укажите точку запроса и радиус, в пределах которого должны быть идентифицированы соседи.

point = [0.4 0.3 0.2];
radius = 0.05;

Получите индексы и расстояния радиальных соседей. Используйте метод облака точек select для получения данных облака точек соседних точек.

[indices,dists] = findNeighborsInRadius(ptCloud,point,radius,camMatrix);
ptCloudB = select(ptCloud,indices);

Отобразите облако точек и радиальные соседи, найденные вокруг точки запроса.

figure
pcshow(ptCloud);
hold on;
pcshow(ptCloudB.Location, 'b');
legend('Point Cloud','Radial Neighbors','Location','southoutside','Color',[1 1 1] )
hold off;