Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

PointCloud - Converting LAZ to COPC.LAZ : CRS information is lost in COPC.LAZ file #59662

Open
2 tasks done
P-SLISSE opened this issue Nov 29, 2024 · 3 comments
Open
2 tasks done
Labels
Bug Either a bug report, or a bug fix. Let's hope for the latter! Point Clouds Projections/Transformations Related to coordinate reference systems or coordinate transformation Regression Something which used to work, but doesn't anymore

Comments

@P-SLISSE
Copy link

What is the bug or the crash?

Hi,
I have some pointcloud tiles in LAZ format.
LAZ files contain CRS information.
When I drag and drop the LAZ file into QGis, LAZ file is automaticaly transform in COPC.LAZ
Because my LAZ files contain CRS information, QGis don't ask for CRS when importing it (my files are in EPSG:3948 and my QGis project too)
In this case : everything's OK
image

I decide to delete the pointcloud entity in Qgis and I now drag and drop the copc.laz file created by Qgis in the preceding step
This entity as no CRS associated and Qgis ask me for one
image
So, I think that the CRS information is lost during the conversion LAZ to COPC, and I don't expect this.

Here is my original LAZ pointcloud file for testing :
2034000_7275500_LIDAR2021.zip

Steps to reproduce the issue

I investigate a little bit more.
I remove the pointcloud entity from Qgis

I use the processing Tool 'Create COPC'
As entry, I set the original LAZ file (2034000_7275500_LIDAR2021.laz on my local hard drive)
I get the transform file 2034000_7275500_LIDAR2021.copc.laz
image
Then I drag and drop the copc.laz file into QGis with the same result : CRS is lost

I decided to transform my laz file with PDAL (as a command line), completely independant from QGis
The command line is simple : pdal translate D:\TMP\LIDAR2025\DATA_LNC\2034000_7275500_LIDAR2021.laz D:\TMP\LIDAR2025\DATA_LNC\2034000_7275500_LIDAR2021_PDAL.copc.laz
When I drag and drop the copc.laz file obtained with PDAL, there is no loss of CRS

So, I assume that Qgis processing don't read ? don't use ? don't write ? the CRS information in copc.laz files

I hope this will help to improve pointcloud processing in QGis

Versions

Version de QGIS
3.40.1-Bratislava
Révision du code
989853a

Libraries
Version de Qt
5.15.13
Version de Python
3.12.7
Version de GDAL/OGR
3.9.3
Version de Proj
9.5.0
Version de la base de données du registre EPSG
v11.016 (2024-08-31)
Version de GEOS
3.13.0-CAPI-1.19.0
Version de SQLite
3.46.1
Version de PDAL
2.8.1
Version du client PostgreSQL
16.2
Version de SpatiaLite
5.1.0
Version de QWT
6.3.0
Version de QScintilla2
2.14.1
Version de l'OS
Windows 10 Version 2009

Extensions Python actives
db_manager
0.1.20
grassprovider
2.12.99
MetaSearch
0.3.6
processing
2.12.99

Supported QGIS version

  • I'm running a supported QGIS version according to the roadmap.

New profile

Additional context

No response

@P-SLISSE P-SLISSE added the Bug Either a bug report, or a bug fix. Let's hope for the latter! label Nov 29, 2024
@agiudiceandrea
Copy link
Contributor

agiudiceandrea commented Nov 29, 2024

Previously also reported at #57983. Not sure why @PeterDrexel closed such previous issue report: running which QGIS version later than 3.38.0 the issue didn't occur?

@P-SLISSE
Copy link
Author

Could anyone tag the issue as Regression ?

@agiudiceandrea agiudiceandrea added Regression Something which used to work, but doesn't anymore Projections/Transformations Related to coordinate reference systems or coordinate transformation labels Nov 29, 2024
@agiudiceandrea
Copy link
Contributor

agiudiceandrea commented Nov 30, 2024

Hi @wonder-sk, may it be the case that the PR #57717 is the root cause of this issue?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Bug Either a bug report, or a bug fix. Let's hope for the latter! Point Clouds Projections/Transformations Related to coordinate reference systems or coordinate transformation Regression Something which used to work, but doesn't anymore
Projects
None yet
Development

No branches or pull requests

3 participants