Given a small sample input file like:

0,0,10
10,10,10

and running p2g with a very large search radius and resolution, one would expect a single output point with a min, mean, and max of 10. This works. But what is the definition of north, south, east, and west in this case?

The output file says they are all 10, but I would expect them to be the cell edges [0,10] or cell center 5.

P2G by itself presently only support LAS v 1.0-1.2 files. With USGS requiring LAS 1.4 deliverables, would be nice to support later LAS versions (1.3 and 1.4) without having to link against PDAL.

I would like to use points2grid to make a DEM from a pointcloud output from openMVS.

Is this possible or does anyone have experience with doing this?

Hey! I'm having installing with cmake. Here is the error log associated to it. Thanks for any help!

my_error.zip

Should be relatively straightforward. Just need to sit down and do it.

Points2Grid should support metric units.

There is a decent argument that this should be the default, but supporting legacy code probably trumps that argument. Perhaps if a ~/.p2g config file existed it for defaults it could be set there. Either way, a command line argument should enable metric units.

I'm not expecting y'all to work on this, but file this anyway as a wishlist and something I may implement at some point myself, and/or others if more people start using this toolchain.

Hi,

I sent this in to the [email protected] list but there was no reply. Perhaps this is the correct place to submit it?

I seem to be running into a user error (or maybe bug?) with points2grid and am seeking some help.

I have a high density pointcloud where the scene FOV is about 1m^2, and the points are about 0.5 to 1 mm apart. The data output (mean,min,max) look good, but the values reported in the density DEM seem a bit high.

As shown below, grid size is 1 mm, search radius is the default, in this case, sqrt(2 mm), full DEM is about 1 sq. m.. Shouldn't GRID_X and GRID_Y be about 1000 (1 m), and shouldn't the point density be around 1. I'm getting GRID_X and GRID_Y of 303 and 236, and point densities of ~ 50.

Thanks for any help you can provide

points2grid --resolution 0.0032808401 -i test.xyz -o test.xyz

stdout:

Parameters ************************
inputName: 'test.xyz'
input_format: 0
outputName: 'test.xyz'
GRID_DIST_X: 0.00328084
GRID_DIST_Y: 0.00328084
searchRadius: 0.00463981
output_format: 0 type: 69905
fill window size: 0 ************************************
inputName: 'test.xyz'
Min/Max searching time: 0.15
min_x: -0.454966, max_x: 0.533365, min_y: -0.357392, max_y: 0.411059
GRID_SIZE_X 303
GRID_SIZE_Y 236
Using incore interp code
InCoreInterp created successfully Interpolation uses in-core algorithm
InCoreInterp::init() done
Interpolation::init() done successfully
Init + Min/Max time: 0.15 Interpolation
Starts Output Execution time: 0.32
Interpolation::interpolation() done successfully
DEM generation + Output time: 0.67

of data: 296104

dimension: 303 x 236

stderr: nothing.

Given a small sample input file like:

0,0,10
10,10,10
10,0,10
0,10,10

and running p2g with a very large search radius and resolution, one would expect a single output point with a min, mean, and max of 10.

I run p2g like this:

points2grid -i test.xyz --output_format grid -r 1E4 --mean --input_format ascii --resolution 1E4 --output_file_name test.xyz --min --max

The output for min, mean, and max are all:

north: 10.000000
south: 0.000000
east: 10.000000
west: 0.000000
rows: 2
cols: 2
10.000000 -9999
10.000000 10.000000

Note that with an even simpler input file of just two points, the two corner points, of [0,0,10], and [10,10,10], then the result is more correct. There output is a 1x1 grid with min, mean, and max of 10.

Hi,

I am using points2grid to generate point count raster as follows:
points2grid -i IN.las -o OUT -r 100 --den --input_format "las" --output_format "arc" --resolution 200

Each grid cell in this output shows the number of points found within the search radius at that node. However, I would like the output to contain number of points per grid cell (i.e. in this example number of points in each 200m cell). Does anybody know if it's possible to do this using Points2Grid?

Thanks for any advice,
Yuri

Given:

0,0,10
10,10,0

And processing with:

points2grid -i test.xyz --output_format grid -r 1E4 --input_format ascii --resolution 1E4 --output_file_name test.xyz --min --max --mean

The mean should be 5, the min 0, and the max 10. But all values are just 0. If the order is switched, all values are 10. North, East, South, and West are also all reported as just the last line.

Pace #17, we have been writing asc and grid file headers incorrectly. @atdyer indicates that the problem exists in geotiff files as well.

I merged #17 to get the fix in, so this ticket is for the geotiff errors only.