如何制作拉斯夹子(激光雷达数据)

时间:2016-07-06 18:05:31

标签: lidar

我正在尝试将激光雷达数据剪切成0.5英里的圆形多边形,但它无法正常工作。我使用了las工具" las clip"但它并没有给我一个结果文件。有谁知道如何制作拉斯夹?

2 个答案:

答案 0 :(得分:1)

您可以在matlab中尝试inpolygon功能。 以下是带有随机采样点云和圆的Matlab代码示例。

close all; clear all;clc;
%Creating sample circle.
L = linspace(0,2.*pi,6);
xv = cos(L)';
yv = sin(L)';
%You can change here center of x,y and diameter.
r=1;x=0;y=0;
th = 0:pi/50:2*pi;
xv = r * cos(th) + x;
yv = r * sin(th) + y;
%% random Point cloud sample 3d
rng default
xq = randn(250,1);
yq = randn(250,1);
zq =randn(250,1);
% Main function Finding inside circle
[in,on] = inpolygon(xq,yq,xv,yv);
numel(xq(in))
% numel(xq(on)) % if its on the circle
%% Plot raw point cloud and circle plan view
figure
plot(xv,yv) % circle
axis equal
hold on
plot(xq(in),yq(in),'r+') % points inside
plot(xq(~in),yq(~in),'bo') % points outside
hold off
%%
k = find(in);
exportx =xq(k);
exporty = yq(k);
exportz = zq(k);
l = find(~in);
outx = xq(l);outy = yq(l);outz = zq(l);
%% Topview
figure
scatter3(exportx,exporty,exportz,'g','+')
hold on
scatter3(outx,outy,outz,'r')
title('Top View')
view(2)
hold off
% Perspective view
figure
scatter3(exportx,exporty,exportz,'g','+')
hold on
scatter3(outx,outy,outz,'r')
title('Perspektif')
view(3)
hold off

images from code

答案 1 :(得分:0)

有点晚了,但是在我使用LAStools的过程中,生成的剪辑文件可以写入一些不寻常的地方 - 当然不是你打算放置它们的地方。

尝试查看C驱动器和LAStools BIN文件夹 - 尽管给出了存储的完整文件路径,但这是我剪辑文件的版本已经结束的地方。