[pdal] Riegl extrabytes stripped

[Benjamin Kabelik]

Dear mailing list members,

I read the thread about the extra bytes and wrote a pdal translate command for a ntv2 projection from UTM33n ETRS89 to Austria MGI M34:


For %i in (*.las) do (docker run -v c:/temp/trafo_wien2/:/data pdal/pdal:1.3 pdal translate reprojection --filters.reprojection.in_srs="EPSG:25833" --filters.reprojection.out_srs="+proj=tmerc +lat_0=0 +lon_0=16.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs +nadgrids=/data/AT_GIS_GRID.gsb" --writers.las.forward=vlr --writers.las.extra_dims="Amplitude=float,Pulse width=float" --writers.las.minor_version=4 --writers.las.dataformat_id=6  --writers.las.offset_x=auto --writers.las.offset_y=auto --writers.las.offset_z=auto -i /data/%i -o /data/%i_gk_gisgrid.las)


There are two things that came up:

1.       I had to switch the data type of amplitude and pulse width from uint16 to float because with uint16 I get integer values instead of float.

2.       Is there any option to forward the description and scale of the extra bytes?



Source file UTM33n ETRS89:
data type: 3 (unsigned short), name "Amplitude", description: "Echo signal amplitude [dB]", min: 0, max: 10000, scale: 0.01, offset: 0 (not set)
                data type: 3 (unsigned short), name "Pulse width", description: "Full width at half maximum [ns]", min: 1, max: 10000, scale: 0.1, offset: 0 (not set)

                Target file Austria MGI M34:
                data type: 9 (float), name "Amplitude", description: "", scale: 1 (not set), offset: 0 (not set)
                data type: 9 (float), name "Pulse width", description: "", scale: 1 (not set), offset: 0 (not set)

Thanks in advance!

Mit freundlichen Grüßen / Kind regards

Benjamin Kabelik
Head of Data Solutions | CIO

[cid:image003.png at 01D22E1F.743CA100]

Airborne Technologies GmbH
Viktor Lang Straße 8
2700 Wiener Neustadt | Austria
P   +43 2622 347 18 206
F   +43 2622 347 18 300
M  +43 664 8845 3023
b.kabelik at airbornetechnologies.at<mailto:b.kabelik at airbornetechnologies.at>
www.airbornetechnologies.at<http://www.airbornetechnologies.at/>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20161024/19e09805/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image003.png
Type: image/png
Size: 4211 bytes
Desc: image003.png
URL: <http://lists.osgeo.org/pipermail/pdal/attachments/20161024/19e09805/attachment.png>