[GRASS-dev] Segmentation fault using v.in.lidar
Doug_Newcomb at fws.gov
Doug_Newcomb at fws.gov
Wed Jul 13 08:59:24 EDT 2011
Pierre,
I think the main issue with doing *.las would be that although LAS
data is coming in tiled but is also available as overlapping swaths. I
imagine to impliment the option for multiple las files you would need to
build an internal index of the "footprint" of each input las file that
tells r.in.lidar and v.in.lidar which file(s) to draw the input las data
from to create the output. In the case of r.in.lidar, you would need to
specify which files to simultaneously open to gather the point values to
calculate the statistics that are the value of the output cell. You would
need something like the way gdaltindex builds a polygon shape file from a
list of input georeferenced images. ( Now if I only had some skill with C
to throw together some code :-( )
In the short term, use the las2las command from liblas.org to
merge your .las files into a single las file and then process the one
file. You will have to watch out for the 4 billion point limit that is
currently in the LAS file format, but for most folks that's not an issue.
Doug
Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 14 doug_newcomb at fws.gov
---------------------------------------------------------------------------------------------------------
The opinions I express are my own and are not representative of the
official policy of the U.S.Fish and Wildlife Service or Dept. of the
Interior. Life is too short for undocumented, proprietary data formats.
Pierre Roudier <pierre.roudier at gmail.com>
Sent by: grass-dev-bounces at lists.osgeo.org
07/13/2011 03:05 AM
To
Markus Metz <markus.metz.giswork at googlemail.com>
cc
grass-dev <grass-dev at lists.osgeo.org>
Subject
Re: [GRASS-dev] Segmentation fault using v.in.lidar
Thanks for the quick answer Markus,
That would be a nice feature to add though. A lot of the LAS files are
coming "tiled", and it'd be nice to be able to do something like:
v.in.lidar in=*.las out=test_input_lidar -trb
or
v.in.lidar in=zone_32_*.las out=test_input_lidar -trb
to import a special subset of LAS files.
I've very few coding abilities, so this is just meant as another line
on the wishlist ;)
Thanks heaps for your work on *.in.lidar, it is working well otherwise,
Pierre
2011/7/13 Markus Metz <markus.metz.giswork at googlemail.com>:
> On Wed, Jul 13, 2011 at 2:08 AM, Pierre Roudier
> <pierre.roudier at gmail.com> wrote:
>> Hi,
>>
>> I've been trying to use v.in.lidar. It yields good results on one LAS
>> file, but I get a segfault when trying it on several files:
>>
>>> v.in.lidar in=BD32_1610.las,BD32_1611.las out=test_input_lidar -trb
>> Segmentation fault
>>> v.in.lidar in=BD32_161*.las out=test_input_lidar -trb
>> Segmentation fault
>>
>> Am I missing something, or is it a bug?
>
> v.in.lidar takes only one input file at a time.
>
> Markus M
>
--
Scientist
Landcare Research, New Zealand
_______________________________________________
grass-dev mailing list
grass-dev at lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-dev/attachments/20110713/4aa27635/attachment-0001.html
More information about the grass-dev
mailing list