如何制作 las 剪辑（激光雷达数据）答案

【问题标题】：How to make a las clip (lidar data)如何制作 las 剪辑（激光雷达数据）
【发布时间】：2016-11-08 21:46:50
【问题描述】：

我正在尝试将激光雷达数据裁剪为 0.5 英里的圆形多边形，但它不起作用。我使用了 las 工具“las clip”，但它没有给我一个结果文件。有谁知道如何制作las剪辑？

【问题讨论】：

gis.stackexchange.com/questions/86075/…

标签： lidar

【解决方案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

【讨论】：

为了将 las 数据导入 matlab，你可以找到它 FEX Matlab。 [link] (mathworks.com/matlabcentral/fileexchange/48073-lasdata) 将las dat 导入Matlab 真是很棒的功能

【解决方案2】：

这有点晚了，但根据我使用 LAStools 的经验，生成的剪辑文件可以写入一些不寻常的地方——当然不是你打算放置它们的地方。

尝试查看 C 驱动器和 LAStools BIN 文件夹 - 尽管已提供完整的存储文件路径，但我的剪辑文件的版本已结束。

【讨论】：