[gdal-dev] weeding out x,y,z points

[Jose Gomez-Dans]

Hi,

2009/2/26 Ari Jolma <ari.jolma at tkk.fi>

> nicholas.g.lawrence at mainroads.qld.gov.au wrote:
>
>>
>>> <nicholas.g.lawrence at mainroads.qld.gov.au>
>>>>
>>>> Can GDAL or OGR weed out 3D points based on proximity to each other?
>>>>
>>>>
>>>  I would use plain Perl (or Python or high level language X) to read in
> the points sequentially and discard new ones that are "too close" to
> existing ones. It may take a while & require a bit of memory it is probably
> doable that way (millions of points but not 100s of millions in the end
> set). The program shouldn't need to be more than 20 lines I guess.


scipy (python's numeric/scientific workhorse) now has a Kdtrees
capabilities, which are an efficient way of working out distances between
points:
<http://docs.scipy.org/doc/scipy/reference/spatial.html>

Best regards,
Jose
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20090226/1b14c172/attachment.html