2016-07-06 2 views
0

ライダーデータを0.5マイルの円ポリゴンでクリップしようとしていますが、動作しません。私はlasツール "las clip"を使用しましたが、それは私に結果ファイルを与えません。誰もがラスクリップを作る方法を知っていますか?lasクリップ(ライダーデータ)の作成方法

+0

http://gis.stackexchange.com/questions/86075/how-to-clip-las-data-using-shapefile-polygons-and-open-source-software/86160#86160 –

答えて

1

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

+0

インポート用のlasデータあなたはそれがFEX Matlabを見つけることができるmatlab。 [リンク](https://www.mathworks.com/matlabcentral/fileexchange/48073-lasdata)Matlabにlas datをインポートするのは本当に素晴らしい機能です – Muzo583

0

が、LAStoolsを使用しての私のexpereinceで、得られるクリップファイルは、いくつかの珍しい場所に書き込むことができます - 確かにあなたがそれらを配置することを予定している場所。

CドライブとLAStools BINフォルダを見てみてください。これは、ストレージへの完全なファイルパスが与えられているにもかかわらず、クリップされたファイルのバージョンが終了した場所です。

関連する問題