mfbonfigli / gocesiumtiler

A Cesium.js point cloud 3D tiles generator from LAS files written in Golang
Mozilla Public License 2.0
195 stars 40 forks source link

The projection caused the height offset ?(My input is EPSG:2326) #30

Closed LeoLiJingWen closed 5 months ago

LeoLiJingWen commented 1 year ago

image According to the above image:

The point cloud at the bottom was transformed using gocesiumtiler

The mesh at the top is a gltf file which was generated from the same point cloud file I use pyproj(proj==9.0.0) to transform the coordinates from EPSG:2326 to EPSG:4978

The final display on Cesium shows that the gltf is correct in position and height, while the point cloud is correct in position but incorrect in height. What causes it?

My conversion command is as follows:

gocesiumtiler -i ./ hk.las -o ./ cesium -e 2326 -n 1.5

mfbonfigli commented 1 year ago

Hello,

I'm not quite sure what could be causing it.

Did you try to add the -g parameter to the converter? Z coordinates could be expressed as elevation above the ellipsoid or above the geoid. gocesiumtiler assumes by default that the elevation is defined as above the ellipsoid as that's the Cesium expectation. If you pass the -g flag it will interpret the Z as above the geoid and attempt to convert it to ellipsoid elevation.

I'm not however sure if this is the issue at play here.

LeoLiJingWen commented 1 year ago

Thank you very much for your reply.

I have tried to add the -g parameter to convert the point cloud, but I waited for a long time and the command was still executing (as shown below). Eventually I gave up, I am not sure if running for a sufficient amount of time would allow the command to complete executing. Without the -g parameter, the command completed executing very quickly.

image image

Here are some basic information about the point cloud(open in CloudCompare), the file size is 113MB. image

LeoLiJingWen commented 1 year ago

image the command executed for over 16 hours and the status was still "exporting data", so I eventually had to terminate it.

mfbonfigli commented 1 year ago

Hi, 16 hours for a 100mb file is definitely too much, a file of that size on a modern computer should normally take minutes to complete. Without data to debug the issue is however hard for me to figure out what is happening here. Please consider if you could send me a sample to troubleshoot

mfbonfigli commented 5 months ago

Closing this one. I followed up offline with LeoLiJingWen and the issue was solved with the version v1.2.3 which improved LAS 1.4 compatibility.