Problems with readPointCloud function

I have used the following functions to read data from a las file

lasReader = lasFileReader('DataCloud.las');

ptCloud = readPointCloud(lasReader);

The obtained pointCloud object has ptCloud.Location stored in single precision even if it also accepts values in double precision. Since for my data it is quite important to have double precision locations (and I am sure that in my las file locations are in double precision), it is possible to suggest to readPointCloud to save locations in double precision?

Thanks a lot for your help

[Helge Nysæter]

I have the same problem. The decimals in the UTM-north-values are lost, and only one remains for the eastings. I then tried to move the origin with pctransform() before retrieving the coordinates from the Location-property. This resulted in more decimals for the eastings, but only one decimal for the northings.

Answer 1

[Helge Nysæter]

Found a workaround with use of lasdata - File Exchange - MATLAB Central (mathworks.com):

Read utm-coordinates from one .las-file, shift origin, and then write to a new .las-file

pkt = lasdata('kd_1.las') 
origoX = floor(min(pkt.x)) 
origoY = floor(min(pkt.y)) 
pkt.x = pkt.x-origoX 
pkt.y = pkt.y-origoY
pkt.z = pkt.z
write_las(pkt,'kd_2.las')