Как сделать клип las (данные лидара)

Question

Я пытаюсь закрепить лидарные данные в круге полигона 0,5 мили, но он не работает. Я использовал инструмент las las clip, но он не дает мне результирующего файла. Кто-нибудь знает, как сделать клип las?

Answer 1

Вы можете попробовать inpolygon функцию в Matlab. Вот пример кода 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

Answer 2

Немного поздно, но в моем expereince использования LAStools, результирующие файлы клипов могут быть записаны в некоторых необычных местах - конечно, не там, где вы собираетесь разместить их.

Попробуйте взглянуть на диск C и папку LAStools BIN - это то, где версии моих обрезанных файлов закончились, несмотря на то, что им был предоставлен полный путь к файлу.