Skip to content

Aligning polygons and point clouds

Ivan Pađen edited this page Oct 7, 2022 · 5 revisions

The most important requirement for the polygon-point cloud reconstruction to work is to have the two datasets aligned, i.e. have them in the same coordinate reference system (CRS), typically defined with a certain EPSG code.

This section of the documentation will lead you through the steps necessary to verify, and eventually fix, the alignment of your datasets.

In which CRS are my datasets? Use QGIS

Start by dragging your datasets into the software QGIS (it is free, (download link)[https://qgis.org]), both polygons and point cloud in LAS/LAZ format. To check the CRS of the layer, double-click on it (which opens up "Layer Properties") and look at the information.

layercrs

Pay attention to units. There is a good chance you will need them in meters. To check how your layers look in their own projection, make sure your layer CRS and project CRS are the same.

projectcrs

The project CRS is written at the bottom-right of the window. To change it, you can click on it and type in the new CRS, or you can right-click on the layer and choose Layer CRS -> Set Project CRS from Layer

If your polygons and point cloud are visually not aligned in QGIS, it means one (or even both) of them are not imported into the program correctly, i.e. in a wrong CRS. Notice that in the first figure below both polygon and point cloud layers are set as visible. However, we can only see polygons. The second figure below shows properly aligned datasets.

not_aligned

aligned

To fix this, you might have to do a bit of guesswork by changing the layer's source CRS:

changesrc

Once you have the alignment right, you can use the reprojection tool to fix the data. Find it in Vector -> Data Management Tools -> Reproject Layer...

Note: It is good practice to double-check the CRS when importing layers. The website/entity that provides the dataset often gives out the information.

(back to top)

Changing CRS'

By changing the project CRS, you can check how your layers look in different projections -- this is just visualization, layer CRS' are the ones that matter for reconstruction. In case you have to change the CRS of a polygon, check the last figure in the Polygon section of the documentation.

To change the CRS of the point cloud, we recommend you use the las2las utility of LAStools. It is compiled as a separate executable in this project under the name city4cfd_las2las. To transform the point cloud with the compiled code, type:

./build/tools/LAStools/city4cfd_las2las -i input_file.laz -target_epsg xxxx -o output_file.laz

If the point cloud is missing CRS information (it will say 'unknown' in QGIS), you can set it with -epsg flag, e.g.

./build/tools/LAStools/city4cfd_las2las -i input_file.laz -epsg xxxx -target_epsg xxxx -o output_file.laz

For macOS and Docker you can directly call the executable, i.e. city4cfd_las2las or ./city4cfd_docker.sh city4cfd_las2las, respectively.

(back to top)

Which CRS should I choose?

Our suggestion would be to first look at the CRS of the point cloud. They are typically localized and therefore stored in CRS that most faithfully describes relations in the scanned region. Otherwise, you can do a location search (typically a country) at https://epsg.io.

Your datasets should ideally be in a local CRS (with units in meters), and you should avoid having them in latitude-longitude (EPSG:4326 is to be avoided!).

(back to top)