We are trying to write a workflow to generate voxel metrics for imputation over large areas, the voxelisation is at 5x5x1 but the summary metric (as a toy example let us say filled voxels/total voxels) is to be computed at a 25 x 25 m resolution (so 5x5 voxels). We usually achieve summary metrics at 25 m resolution using lascanopy. Do you have a suggested workflow for using the outputs of LASvoxel with LAScanopy to get the summary metrics at this resolution? We can do this in intermediate steps with other tools but we would like to start using lasvoxel + lascanopy to generate more metrics to better define our forest structure.
Example of the normalised LiDAR pixel (attached)
We have tried two approaches:
lasvoxel.exe -i %NORMALISE% -step_xy 5 -step_z 1 -odir %VOXEL% -odix _xy5_z1_Count_Int -olaz
Here is what we obtain, coloured by intensity, I think in the intensity field we have the number of points (count) for each particular voxel
From this file we could compute our voxel metrics, the voxels are 5*5*1 but they will need to be aggregated in 25 by 25m, at this stage we can run lascanopy on the intensity field, so this is not a problem.
We are trying as well, that I think is where we actually want to go to extract the voxel metrics.
lasvoxel.exe -i %NORMALISE% -step_xy 5 -step_z 1 -compute_IDs_and_voxel_table -store_IDs_in_intensity -odir %VOXEL% -olaz -odix _xy5_z1_vID_Int
where the voxel id is saved in the intensity field.
Coloured by Voxel_ID (intensity)
It would be interesting to be available to compute the standard metrics by voxel id and then aggregate them in 25 by 25m, ready for the LiDAR imputation.
Tell me your thoughts.
Susana Gonzalez - Forest Engineer, LiDAR Science
Interpine Group Ltd
+64 7 350 3209
Australia 0280113645 ext 722
Address : 99 Sala Street, PO Box 1209, Rotorua 3010, New Zealand
Website : www.interpine.co.nz
Interpine Innovation is Shaping Today’s Forests with the Technology of Tomorrow
He rangahau tenei ra he hangarau apopo