Hello !

 

I am working with georeferenced LiDAR data with R and sp. I am trying to
normalize height and remove the slope effect from a Terrestrial Laser
Scanner (TLS) point cloud. 

The available data is :

-          a TLS point cloud in ASCII format converted to sp objet of class
SpatialPoint (variable named "tls.z.sp")

-          a raster DTM in TIF format (1 m resolution) converted to sp
object of class SpatialGridDataFrame (variable named "tls.dtm$z")

 

The method consists in subtracting the DTM z value from each TLS point z
value. I proceeded to some tests with the "over()" function with :

x = "SpatialPoints"

y = "SpatialGridDataFrame"

 

nz <- over(tls.z.sp, tls.dtm$z, tls.z.sp-tls.dtm$z)

 

I am not sure to be on the right track. Any suggestion would be very much
appreciated !

 

Cheers,

Nathalie

 


	[[alternative HTML version deleted]]

