Another fella asking for Projection Assistance

Justin C

unread,

Oct 27, 2020, 4:37:59 PM10/27/20

to LAStools - efficient tools for LiDAR processing

Hello,

I am aware that what I am about to ask has been asked a lot already, such as:

https://groups.google.com/g/lastools/c/1A36bxhcZUM/m/VkdgBMPVAgAJ

However, I've tried for the last few hours to solve on my own and now have exhausted my efforts. If someone could dumb down the specific steps I need to take for my problem it would be immensely appreciated. Please forgive me in advance if I mix up nomenclature.

I have a publically available LiDAR tile in LAZ format. It's simplified information is as follows:

Projection = NAD_1983_CSRS_UTM_Zone_20N
Horizontal Datum = NAD83_CSRS_2010
Vertical Datum = CGVD2013
Year = 2019
LAS Version = 1.4
Source = LIDAR

Normally I just extract into LAS using the standalone laszip.exe, keeping class 2 and loading into Recap (.rcp) for use in Civil3D (.tin). However, we have recently obtained a project from another firm in a different co-ordinate system (Nova Scotia NAD83 CSRS2010 MTM Z5) and will have to shift our surface model to match their grading plans etc.

I would like to transform the NAD83 LAS from UTM20 into MTM, and maintain the vertical datum of CGVD2013 (If it matters, the full MTM is referenced as NS-NAD83-2010-5).

My understanding is I need to reproject from EPSG 2961 (The usual) into EPSG 8083 (Similar to the old standard EPSG 2295). The laszip standalone GUI has an option for EPSG codes in the projection dropdown but crashes when I select it- so I'm running the LAStools addon through QGIS 2.18 Las Palmas and can't seem to find where to project by EPSG codes through the QGIS GUI.

What steps do I need to take to transform my LAS into the localized co-ordinate grid?

For reference, I'm using the 2019 (not 2013) data file on tile 476_5026 found here:

https://nsgi.novascotia.ca/datalocator/elevation/

Martin Isenburg

unread,

Oct 30, 2020, 1:27:10 PM10/30/20

to LAStools - efficient command line tools for LIDAR processing

Hello,

They have a very nice download portal in Nova Scotia. I got the file in question with ease (see screen shot). Kudos to those in charge for open LiDAR data in Nova Scotia.

Usually I would then run

E:\software\LAStools\bin>las2las -i c:\Users\Martin\Downloads\476_5026_201901.laz -epsg 2961 -target_epsg 8083 -odix _8083 -olaz
EPSG code 8083 not found in 'pcs.csv' file
ERROR: unknown EPSG code in '-target_epsg 8083'.

E:\software\LAStools\bin>las2las -i c:\Users\Martin\Downloads\476_5026_201901.laz -target_epsg 8083 -odix _8083 -olaz
EPSG code 8083 not found in 'pcs.csv' file
ERROR: unknown EPSG code in '-target_epsg 8083'.

but this EPSG code is not something I am familiar with. Are there any Canadian CRS nerds out there? What projection is that exactly?

https://epsg.io/8083

==================================

E:\software\LAStools\bin>lasinfo -i c:\Users\Martin\Downloads\476_5026_201901.laz
lasinfo (201028) report for 'c:\Users\Martin\Downloads\476_5026_201901.laz'
reporting all LAS header entries:
file signature: 'LASF'
file source ID: 0
global_encoding: 17
project ID GUID data 1-4: 00000000-0000-0000-0000-000000000000
version major.minor: 1.4
system identifier: 'LAStools (c) by rapidlasso GmbH'
generating software: 'las2las (version 200131)'
file creation day/year: 15/2019
header size: 375
offset to point data: 1430
number var. length records: 2
point data format: 6
point data record length: 30
number of point records: 0
number of points by return: 0 0 0 0 0
scale factor x y z: 0.01 0.01 0.01
offset x y z: 400000 5000000 0
min x y z: 476000.00 5026000.00 -125.22
max x y z: 476999.99 5026999.99 113.36
start of waveform data packet record: 0
start of first extended variable length record: 0
number of extended_variable length records: 0
extended number of point records: 17514732
extended number of points by return: 14366216 2429785 563140 127123 23629 4178 661 0 0 0 0 0 0 0 0
variable length header record 1 of 2:
reserved 0
user ID 'LASF_Projection'
record ID 34735
length after header 48
description 'by LAStools of rapidlasso GmbH'
GeoKeyDirectoryTag version 1.1.0 number of keys 5
key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected
key 3072 tiff_tag_location 0 count 1 value_offset 2961 - ProjectedCSTypeGeoKey: NAD83(CSRS) / UTM zone 20N
key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter
key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter
key 4096 tiff_tag_location 0 count 1 value_offset 6647 - VerticalCSTypeGeoKey: Canadian Geodetic Vertical Datum of 2013
variable length header record 2 of 2:
reserved 0
user ID 'LASF_Projection'
record ID 2112
length after header 899
description 'by LAStools of rapidlasso GmbH'
WKT OGC COORDINATE SYSTEM:
COMPD_CS["NAD83(CSRS) / UTM zone 20N + CGVD2013",PROJCS["NAD83(CSRS) / UTM zone 20N",GEOGCS["NAD83(CSRS)",DATUM["NAD83_Canadian_Spatial_Reference_System",SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6140"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4617"]],PROJECTION["Tra
nsverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-63],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,A
UTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","2961"]],VERT_CS["CGVD2013",VERT_DATUM["Canadian Geodetic Vertical Datum of 2013",2005,AUTHORITY["EPSG","1127"]],U
NIT["metre",1.0,AUTHORITY["EPSG","9001"]],AXIS["Gravity-related height",UP],AUTHORITY["EPSG","6647"]]]
LASzip compression (version 3.4r3 c3 50000): POINT14 3
LAStiling (idx 205963, lvl 9, sub 0, bbox 187000 4695000 699000 5207000) (size 1000 x 1000, buffer 0)
reporting minimum and maximum for all LAS point record entries ...
X 7600000 7699999
Y 2600000 2699999
Z -12522 11336
intensity 6 21432
return_number 1 7
number_of_returns 1 7
edge_of_flight_line 0 1
scan_direction_flag 0 1
classification 2 18
scan_angle_rank -33 32
user_data 0 255
point_source_ID 10011 30145
gps_time 246362877.637832 248873384.067001
extended_return_number 1 7
extended_number_of_returns 1 7
extended_classification 2 18
extended_scan_angle -5500 5333
extended_scanner_channel 1 3
number of first returns: 14366216
number of intermediate returns: 718552
number of last returns: 14366903
number of single returns: 11936939
overview over extended number of returns of given pulse: 11936939 3733305 1307530 413793 97960 20415 4790 0 0 0 0 0 0 0 0
histogram of classification of points:
2080670 ground (2)
6408505 low vegetation (3)
926788 medium vegetation (4)
8093469 high vegetation (5)
2712 noise (7)
2525 bridge deck (17)
63 Reserved for ASPRS Definition (18)
+-> flagged as withheld: 1186536
+---> 461025 of those are low vegetation (3)
+---> 57107 of those are medium vegetation (4)
+---> 665609 of those are high vegetation (5)
+---> 2712 of those are noise (7)
+---> 20 of those are bridge deck (17)
+---> 63 of those are Reserved for ASPRS Definition (18)
+-> flagged as extended overlap: 5071114
+---> 446228 of those are ground (2)
+---> 1995469 of those are low vegetation (3)
+---> 288884 of those are medium vegetation (4)
+---> 2339986 of those are high vegetation (5)
+---> 542 of those are noise (7)
+---> 5 of those are Reserved for ASPRS Definition (18)

--
Download LAStools at
http://lastools.org
http://rapidlasso.com
Be social with LAStools at
http://facebook.com/LAStools
http://twitter.com/LAStools
http://linkedin.com/groups/LAStools-4408378
Manage your settings at
http://groups.google.com/group/lastools/subscribe
---
You received this message because you are subscribed to the Google Groups "LAStools - efficient tools for LiDAR processing" group.
To unsubscribe from this group and stop receiving emails from it, send an email to lastools+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/49ac3e4e-45ef-4a30-967a-d135b62bce1en%40googlegroups.com.

201030_opendata_download_portal_nova_scotia.jpg

Kirk Waters - NOAA Federal

unread,

Oct 30, 2020, 1:39:19 PM10/30/20

to LAStools - efficient command line tools for LIDAR processing

From www.epsg.org:

PROJCRS["NAD83(CSRS)v6 / MTM Nova Scotia zone 5",
  BASEGEOGCRS["NAD83(CSRS)v6",
    DATUM["North American Datum of 1983 (CSRS) version 6",
      ELLIPSOID["GRS 1980",6378137,298.257222101,
        LENGTHUNIT["metre",1,ID["EPSG",9001]],
        ID["EPSG",7019]],
      ID["EPSG",1197]],
    PRIMEM["Greenwich",0,
      ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9102]],
      ID["EPSG",8901]],
    ID["EPSG",8252]],
  CONVERSION["MTM Nova Scotia zone 5 v2",
    METHOD["Transverse Mercator",
      ID["EPSG",9807]],
    PARAMETER["Latitude of natural origin",0,
      ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9102]]],
    PARAMETER["Longitude of natural origin",-64.5,
      ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9102]]],
    PARAMETER["Scale factor at natural origin",0.9999,
      SCALEUNIT["unity",1,ID["EPSG",9201]]],
    PARAMETER["False easting",25500000,
      LENGTHUNIT["metre",1,ID["EPSG",9001]]],
    PARAMETER["False northing",0,
      LENGTHUNIT["metre",1,ID["EPSG",9001]]],
    ID["EPSG",8081]],
  CS[Cartesian,2,
    ID["EPSG",4400]],
  AXIS["Easting (E)",east,
    ORDER[1]],
  AXIS["Northing (N)",north,
    ORDER[2]],
  LENGTHUNIT["metre",1,ID["EPSG",9001]],
  USAGE[SCOPE["Engineering survey, topographic mapping."],
    AREA["Canada - Nova Scotia - west of 63°W."],
    BBOX[43.41,-66.28,46.02,-63.00]],
  ID["EPSG",8083]]

Kirk

To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/CABSWR-F%2B7FMk281FTGM3Hpqz_inRtCw9D3_-YCBkM%3DXFg2GvRA%40mail.gmail.com.

Martin Isenburg

unread,

Nov 1, 2020, 8:02:45 AM11/1/20

to LAStools - efficient command line tools for LIDAR processing

Hi,

anyone a specialist with Proj4 or the maintenance of the magic 'pcs.csv' file? Is my 'pcs.csv' out of date or why is this Transverse Mercator projection not listed in mine?

EPSG code 8083 not found in 'pcs.csv' file

Martin

To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/CADm%3DQrROd41TU3qCcQRmSYmLoSTK-ADAtDStbT%3D8Rn-V4ngTsg%40mail.gmail.com.

Jeff McKenna

unread,

Nov 3, 2020, 10:42:05 AM11/3/20

to last...@googlegroups.com

Hi Martin, Hi Justin,

As your resident Nova Scotia PROJ (formerly referred to as "PROJ.4")
person, here are my thoughts:

- Nova Scotia's MTM ("Modified Transverse Mercator") has been in use by
the surveying community here since around the 1970's

- it was finally included in the EPSG geodetic database in 2015, as EPSG
8082 (zone 4) and 8083 (zone 5) (details:
https://epsg.org/crs_8083/NAD83-CSRS-v6-MTM-Nova-Scotia-zone-5.html)

- GeoNOVA (our provincial government entity who manages geodata) now
lists both the UTM zone 20 (EPSG: 2961) projection and these 2 MTM
projections as part of their official supported projections in the "Nova
Scotia Coordinate Referencing System" (GeoNOVA's white paper on these
supported projections, that includes the details of MTM:
https://geonova.novascotia.ca/sites/default/files/resource-library/NSCRS%20Technical%20Support%200004%20Projections%2029Feb2016.pdf
)

- it would indeed be good for LAStools to support EPSG 8082 and 8083

- GDAL's 'pcs.csv' file includes both 8082 and 8083, in GDAL verison 2.4.4

- in the GDAL 3.x series, GDAL uses PROJ's SQLite database 'proj.db',
where you can see that EPSG: 8082 and 8083 exist in the 'projected_crs'
table, such as:

sqlite> select * from projected_crs where code = '8083';

EPSG|8083|NAD83(CSRS)v6 / MTM Nova Scotia zone
5||EPSG|4400|EPSG|8252|EPSG|8081||0

Other thoughts of mine:

- as someone pushing the use of OGC standards and FOSS4G for now 20
years, I can tell you that I unfortunately don't see these 2 MTM
projections 'out in the wild' much. For example, note that GeoNOVA's
own OGC services do not support EPSG 8082 or 8083, see an example
GetCapabilities request:
https://nsgiwa.novascotia.ca/arcgis/services/BASE/BASE_NSTDB_10k_Roads_UT83/MapServer/WMSServer?request=GetCapabilities&service=WMS

- these 2 MTM projections are more relied upon by the local surveying
community

In summary, yes I believe these 2 MTM projections should be supported by
LAStools.

Thanks, and hello from a cloudy and cool Blue Rocks, Nova Scotia.

-jeff

--
Jeff McKenna
MapServer Consulting and Training Services
co-founder of FOSS4G
http://gatewaygeo.com/

On 2020-11-01 9:01 a.m., Martin Isenburg wrote:
> Hi,
>
> anyone a specialist with Proj4 or the maintenance of the magic 'pcs.csv'
> file? Is my 'pcs.csv' out of date or why is this Transverse Mercator
> projection not listed in mine?
>
> EPSG code 8083 not found in 'pcs.csv' file
>
> Martin
>
> On Fri, Oct 30, 2020 at 6:39 PM 'Kirk Waters - NOAA Federal' via
> LAStools - efficient tools for LiDAR processing

> <last...@googlegroups.com <mailto:last...@googlegroups.com>> wrote:
>
> From www.epsg.org <http://www.epsg.org>:

> https://epsg.io/8083 <https://epsg.io/8083>

> tile *476_5026* found here:
> https://nsgi.novascotia.ca/datalocator/elevation/

> <https://nsgi.novascotia.ca/datalocator/elevation/>
>
>
> --
> Download LAStools at

> http://lastools.org <http://lastools.org>
> http://rapidlasso.com <http://rapidlasso.com>

> Be social with LAStools at

> http://facebook.com/LAStools <http://facebook.com/LAStools>
> http://twitter.com/LAStools <http://twitter.com/LAStools>
> http://linkedin.com/groups/LAStools-4408378

> <http://linkedin.com/groups/LAStools-4408378>
> Manage your settings at
> http://groups.google.com/group/lastools/subscribe
> <http://groups.google.com/group/lastools/subscribe>
> ---
> You received this message because you are subscribed to the
> Google Groups "LAStools - efficient tools for LiDAR
> processing" group.
> To unsubscribe from this group and stop receiving emails
> from it, send an email to
> lastools+u...@googlegroups.com

> <mailto:lastools+u...@googlegroups.com>.

> To view this discussion on the web visit
> https://groups.google.com/d/msgid/lastools/49ac3e4e-45ef-4a30-967a-d135b62bce1en%40googlegroups.com

> <https://groups.google.com/d/msgid/lastools/49ac3e4e-45ef-4a30-967a-d135b62bce1en%40googlegroups.com?utm_medium=email&utm_source=footer>.
>
> --
> Download LAStools at
> http://lastools.org <http://lastools.org>
> http://rapidlasso.com <http://rapidlasso.com>

> Be social with LAStools at

> http://facebook.com/LAStools <http://facebook.com/LAStools>
> http://twitter.com/LAStools <http://twitter.com/LAStools>
> http://linkedin.com/groups/LAStools-4408378

> <http://linkedin.com/groups/LAStools-4408378>
> Manage your settings at
> http://groups.google.com/group/lastools/subscribe
> <http://groups.google.com/group/lastools/subscribe>
> ---
> You received this message because you are subscribed to the
> Google Groups "LAStools - efficient tools for LiDAR processing"
> group.
> To unsubscribe from this group and stop receiving emails from
> it, send an email to lastools+u...@googlegroups.com

> <mailto:lastools+u...@googlegroups.com>.

> To view this discussion on the web visit

> https://groups.google.com/d/msgid/lastools/CABSWR-F%2B7FMk281FTGM3Hpqz_inRtCw9D3_-YCBkM%3DXFg2GvRA%40mail.gmail.com
> <https://groups.google.com/d/msgid/lastools/CABSWR-F%2B7FMk281FTGM3Hpqz_inRtCw9D3_-YCBkM%3DXFg2GvRA%40mail.gmail.com?utm_medium=email&utm_source=footer>.
>
> --
> Download LAStools at
> http://lastools.org <http://lastools.org>
> http://rapidlasso.com <http://rapidlasso.com>

> Be social with LAStools at

> http://facebook.com/LAStools <http://facebook.com/LAStools>
> http://twitter.com/LAStools <http://twitter.com/LAStools>
> http://linkedin.com/groups/LAStools-4408378

> <http://linkedin.com/groups/LAStools-4408378>
> Manage your settings at
> http://groups.google.com/group/lastools/subscribe
> <http://groups.google.com/group/lastools/subscribe>
> ---
> You received this message because you are subscribed to the Google
> Groups "LAStools - efficient tools for LiDAR processing" group.
> To unsubscribe from this group and stop receiving emails from it,
> send an email to lastools+u...@googlegroups.com

> <mailto:lastools+u...@googlegroups.com>.

> To view this discussion on the web visit

> https://groups.google.com/d/msgid/lastools/CADm%3DQrROd41TU3qCcQRmSYmLoSTK-ADAtDStbT%3D8Rn-V4ngTsg%40mail.gmail.com
> <https://groups.google.com/d/msgid/lastools/CADm%3DQrROd41TU3qCcQRmSYmLoSTK-ADAtDStbT%3D8Rn-V4ngTsg%40mail.gmail.com?utm_medium=email&utm_source=footer>.
>
> --
> Download LAStools at
> http://lastools.org <http://lastools.org>
> http://rapidlasso.com <http://rapidlasso.com>

> Be social with LAStools at

> http://facebook.com/LAStools <http://facebook.com/LAStools>
> http://twitter.com/LAStools <http://twitter.com/LAStools>
> http://linkedin.com/groups/LAStools-4408378

> <http://linkedin.com/groups/LAStools-4408378>
> Manage your settings at
> http://groups.google.com/group/lastools/subscribe
> <http://groups.google.com/group/lastools/subscribe>
> ---
> You received this message because you are subscribed to the Google
> Groups "LAStools - efficient tools for LiDAR processing" group.
> To unsubscribe from this group and stop receiving emails from it, send
> an email to lastools+u...@googlegroups.com

> <mailto:lastools+u...@googlegroups.com>.

> To view this discussion on the web visit

> https://groups.google.com/d/msgid/lastools/CABSWR-HjLUPLS5eguN_ZtzSfAJx-4pwyM6kTEtFTZQKiGvqT9A%40mail.gmail.com
> <https://groups.google.com/d/msgid/lastools/CABSWR-HjLUPLS5eguN_ZtzSfAJx-4pwyM6kTEtFTZQKiGvqT9A%40mail.gmail.com?utm_medium=email&utm_source=footer>.

Justin C

unread,

Nov 9, 2020, 10:59:11 AM11/9/20

to LAStools - efficient tools for LiDAR processing

Hi Martin, Hi Jeff,

Thank you both for your time in providing detailed and informational responses. Although I will admit, a lot of it is over my head!

I will take my brief opportunity to thank both of you for your tireless efforts over your careers in progressing LiDAR. As someone in the younger generation of Civil Industry, LiDAR is mind-blowing in it's nearly limitless potential. Martin I've followed several of your tutorials over the years, and Jeff you've responded to my personal e-mails diligently every time!

Jeff nailed it pretty well. These old MTM projections are heavily used by our land surveying and Civil Engineering communities in NS mostly due to heavy precedent. I'm working in the Civil Engineering industry with a local consulting firm with decades and decades of experience, and we have only seen 3 different projections in total, but the NS-MTM NAD83 and ATS77 systems are the massively majority uses.

For future users experiencing this problem:

I found a "workaround" using another proprietary software, Recap Pro. They currently have Nova Scotia's specific MTM projections loaded for use.

Use LAStools to extract Class 2 data from .LAZ without any addition manipulation
(optional) Use LAStools to stitch the class 2 LAS's into one LAS cloud
Use Recap Pro to reproject into NS-MTM (WARNING: Output .rcp format so limited use)
External Reference .rcp into a Civil3D drawing with the correct target datum defined in Drawing Settings
Manipulate further directly in C3D or export to new point-cloud for further LAStools use

Thanks for all the help, incredibly impressed with the effort you've both given!

Justin C

joual.f...@gmail.com

unread,

Nov 10, 2020, 12:08:49 PM11/10/20

to LAStools - efficient tools for LiDAR processing

It would be great if Lastools can read/transform MTM CSRS projection Datum. Here in Quebec, such as Nova Scotia, we use it.

For example, EPSG 2950, using CGVD2013 vertical datum: PROJCS["NAD83(CSRS) / MTM zone 8",GEOGCS["NAD83(CSRS)",DATUM["NAD83_Canadian_Spatial_Reference_System",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6140"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4617"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-73.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",304800],PARAMETER["false_northing",0],AUTHORITY["EPSG","2950"],AXIS["E(X)",EAST],AXIS["N(Y)",NORTH]]

Softwares like Global Mapper and FME use EPSG CSRS codes. Ex: EPSG 2950.

Thanks

Frank

Terje Mathisen

unread,

Nov 10, 2020, 3:06:03 PM11/10/20

to last...@googlegroups.com, joual.f...@gmail.com

If LAStools is using PROJ then supporting that MTM CSRS would seem automatic as soon as an updated DB was included? I know Martin have added a bunch of these EPSG codes over the last few years. :-)

Terje

joual.f...@gmail.com wrote:

It would be great if Lastools can read/transform MTM CSRS projection Datum. Here in Quebec, such as Nova Scotia, we use it.
For example, EPSG 2950, using CGVD2013 vertical datum: PROJCS["NAD83(CSRS) / MTM zone 8",GEOGCS["NAD83(CSRS)",DATUM["NAD83_Canadian_Spatial_Reference_System",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6140"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4617"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-73.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",304800],PARAMETER["false_northing",0],AUTHORITY["EPSG","2950"],AXIS["E(X)",EAST],AXIS["N(Y)",NORTH]]

Softwares like Global Mapper and FME use EPSG CSRS codes. Ex: EPSG 2950.

Thanks

Frank

Le lundi 9 novembre 2020 Ã 10 h 59 min 11 s UTC-5, Justin C a Ã©critÂ :

> anyone a specialist with Proj4 or the maintenanceÂ of the magic 'pcs.csv'
> file? Is my 'pcs.csv' out of date or why is thisÂ Transverse Mercator

> AREA["Canada - Nova Scotia - west of 63Â°W."],

> Â file signature: Â Â Â Â Â Â 'LASF'
> Â file source ID: Â Â Â Â Â Â 0
> Â global_encoding: Â Â Â Â Â Â 17
> Â project ID GUID data 1-4: Â 00000000-0000-0000-0000-000000000000
> Â version major.minor: Â Â Â Â 1.4
> Â system identifier: Â Â Â Â Â 'LAStools (c) by rapidlasso GmbH'
> Â generating software: Â Â Â Â 'las2las (version 200131)'
> Â file creation day/year: Â Â 15/2019
> Â header size: Â Â Â Â Â Â Â Â 375
> Â offset to point data: Â Â Â 1430
> Â number var. length records: 2
> Â point data format: Â Â Â Â Â 6
> Â point data record length: Â 30
> Â number of point records: Â Â 0
> Â number of points by return: 0 0 0 0 0
> Â scale factor x y z: Â Â Â Â 0.01 0.01 0.01
> Â offset x y z: Â Â Â Â Â Â Â 400000 5000000 0
> Â min x y z: Â Â Â Â Â Â Â Â Â 476000.00 5026000.00 -125.22
> Â max x y z: Â Â Â Â Â Â Â Â Â 476999.99 5026999.99 113.36
> Â start of waveform data packet record: 0
> Â start of first extended variable length record: 0
> Â number of extended_variable length records: 0
> Â extended number of point records: 17514732
> Â extended number of points by return: 14366216 2429785 563140

> 127123 23629 4178 661 0 0 0 0 0 0 0 0
> variable length header record 1 of 2:

> Â reserved Â Â Â Â Â Â 0
> Â user ID Â Â Â Â Â Â Â 'LASF_Projection'
> Â record ID Â Â Â Â Â Â 34735
> Â length after header Â 48
> Â description Â Â Â Â Â 'by LAStools of rapidlasso GmbH'
> Â Â GeoKeyDirectoryTag version 1.1.0 number of keys 5
> Â Â Â key 1024 tiff_tag_location 0 count 1 value_offset 1 -
> GTModelTypeGeoKey: ModelTypeProjected
> Â Â Â key 3072 tiff_tag_location 0 count 1 value_offset 2961 -

> ProjectedCSTypeGeoKey: NAD83(CSRS) / UTM zone 20N

> Â Â Â key 3076 tiff_tag_location 0 count 1 value_offset 9001 -
> ProjLinearUnitsGeoKey: Linear_Meter
> Â Â Â key 4099 tiff_tag_location 0 count 1 value_offset 9001 -
> VerticalUnitsGeoKey: Linear_Meter
> Â Â Â key 4096 tiff_tag_location 0 count 1 value_offset 6647 -

> VerticalCSTypeGeoKey: Canadian Geodetic Vertical Datum of 2013
> variable length header record 2 of 2:

> Â reserved Â Â Â Â Â Â 0
> Â user ID Â Â Â Â Â Â Â 'LASF_Projection'
> Â record ID Â Â Â Â Â Â 2112
> Â length after header Â 899
> Â description Â Â Â Â Â 'by LAStools of rapidlasso GmbH'
> Â Â WKT OGC COORDINATE SYSTEM:
> Â Â COMPD_CS["NAD83(CSRS) / UTM zone 20N +

> CGVD2013",PROJCS["NAD83(CSRS) / UTM zone
> 20N",GEOGCS["NAD83(CSRS)",DATUM["NAD83_Canadian_Spatial_Reference_System",SPHEROID["GRS
> 1980",6378137,298.257222101,
> AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6140"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4617"]],PROJECTION["Tra
> nsverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-63],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,A
> UTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","2961"]],VERT_CS["CGVD2013",VERT_DATUM["Canadian
> Geodetic Vertical Datum of 2013",2005,AUTHORITY["EPSG","1127"]],U
> NIT["metre",1.0,AUTHORITY["EPSG","9001"]],AXIS["Gravity-related
> height",UP],AUTHORITY["EPSG","6647"]]]
> LASzip compression (version 3.4r3 c3 50000): POINT14 3
> LAStiling (idx 205963, lvl 9, sub 0, bbox 187000 4695000 699000
> 5207000) (size 1000 x 1000, buffer 0)
> reporting minimum and maximum for all LAS point record entries ...

> Â X Â Â Â Â Â Â 7600000 Â Â 7699999
> Â Y Â Â Â Â Â Â 2600000 Â Â 2699999
> Â Z Â Â Â Â Â Â Â -12522 Â Â Â 11336
> Â intensity Â Â Â Â Â 6 Â Â Â 21432
> Â return_number Â Â Â 1 Â Â Â Â Â 7
> Â number_of_returns Â 1 Â Â Â Â Â 7
> Â edge_of_flight_line 0 Â Â Â Â Â 1
> Â scan_direction_flag 0 Â Â Â Â Â 1
> Â classification Â Â Â 2 Â Â Â Â 18
> Â scan_angle_rank Â -33 Â Â Â Â 32
> Â user_data Â Â Â Â Â 0 Â Â Â Â 255
> Â point_source_ID 10011 Â Â Â 30145
> Â gps_time 246362877.637832 248873384.067001
> Â extended_return_number Â Â Â Â Â 1 Â Â Â 7
> Â extended_number_of_returns Â Â Â 1 Â Â Â 7
> Â extended_classification Â Â Â Â 2 Â Â 18
> Â extended_scan_angle Â Â Â Â -5500 Â 5333
> Â extended_scanner_channel Â Â Â Â 1 Â Â Â 3
> number of first returns: Â Â Â Â 14366216

> number of intermediate returns: 718552

> number of last returns: Â Â Â Â 14366903
> number of single returns: Â Â Â 11936939

> overview over extended number of returns of given pulse:
> 11936939 3733305 1307530 413793 97960 20415 4790 0 0 0 0 0 0 0 0
> histogram of classification of points:

> Â Â Â Â Â 2080670 Â ground (2)
> Â Â Â Â Â 6408505 Â low vegetation (3)
> Â Â Â Â Â 926788 Â medium vegetation (4)
> Â Â Â Â Â 8093469 Â high vegetation (5)
> Â Â Â Â Â Â 2712 Â noise (7)
> Â Â Â Â Â Â 2525 Â bridge deck (17)
> Â Â Â Â Â Â Â 63 Â Reserved for ASPRS Definition (18)
> Â +-> flagged as withheld: Â 1186536
> Â +---> Â Â Â Â Â 461025 of those are low vegetation (3)
> Â +---> Â Â Â Â Â 57107 of those are medium vegetation (4)
> Â +---> Â Â Â Â Â 665609 of those are high vegetation (5)
> Â +---> Â Â Â Â Â Â 2712 of those are noise (7)
> Â +---> Â Â Â Â Â Â Â 20 of those are bridge deck (17)
> Â +---> Â Â Â Â Â Â Â 63 of those are Reserved for ASPRS
> Definition (18)
> Â +-> flagged as extended overlap: 5071114
> Â +---> Â Â Â Â Â 446228 of those are ground (2)
> Â +---> Â Â Â Â 1995469 of those are low vegetation (3)
> Â +---> Â Â Â Â Â 288884 of those are medium vegetation (4)
> Â +---> Â Â Â Â 2339986 of those are high vegetation (5)
> Â +---> Â Â Â Â Â Â 542 of those are noise (7)
> Â +---> Â Â Â Â Â Â Â 5 of those are Reserved for ASPRS

http://twitter.com/LAStools

http://linkedin.com/groups/LAStools-4408378
Manage your settings at
http://groups.google.com/group/lastools/subscribe

---
You received this message because you are subscribed to the Google Groups "LAStools - efficient tools for LiDAR processing" group.

To unsubscribe from this group and stop receiving emails from it, send an email to lastools+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/5f7eb025-67c9-43f7-b267-34f9df148271n%40googlegroups.com.

-- 
- <Terje.M...@tmsw.no>
"almost all programming can be viewed as an exercise in caching"

Kirk Waters - NOAA Federal

unread,

Nov 10, 2020, 4:24:31 PM11/10/20

to LAStools - efficient command line tools for LIDAR processing

I don't think PROJ is in the back-end. As I recall, the *.csv files from GDAL (pre-GDAL3?) were incorporated.

Kirk

To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/901ef5d2-2c31-a896-d8f2-b73d7823132e%40tmsw.no.

Martin Isenburg

unread,

Nov 11, 2020, 7:32:57 AM11/11/20

to LAStools - efficient command line tools for LIDAR processing

Hello,

indeed. Thanks to Jeff for connecting me to the latest *.csv files. They will be part of the next release. Here they are until then

https://www.dropbox.com/s/magapow18hakmlr/LAStools_bin_serf_geo.zip

Simply replace the three files in the ./LAStools/bin/serf/geo folder with these ones. Then the EPSG codes are recognized,

E:\LAStools\bin>las2las -i 476_5026_201901.laz -epsg 2961 -target_epsg 8083 -odix _8083 -olaz
SERIOUS WARNING: horizontal datum of source and target incompatible.
WARNING: reprojecting file with lastiling VLR invalidates this VLR

But there are two issues.

(1) The horizontal datum of source and target are incompatible. Do you need a more complex transformation? Do you need to use a grid of shift vectors? Your use case and your national geospatial agency (or other professionals) should have an answer for you here.

(2) The pre-existing LAStiling VLR in this LAZ file (added by lastile) no longer makes sense. It should be removed with option '-remove_tiling_vlr'.

Below the resulting lasinfo reports with and without option '-remove_tiling_vlr'.

E:\LAStools\bin>lasinfo -i 476_5026_201901_8083.laz
lasinfo (201111) report for '476_5026_201901_8083.laz'

reporting all LAS header entries:

file signature: 'LASF'
file source ID: 0
global_encoding: 17
project ID GUID data 1-4: 00000000-0000-0000-0000-000000000000
version major.minor: 1.4

system identifier: 'LAStools (c) by rapidlasso GmbH'
generating software: 'las2las (version 201030)'

file creation day/year: 15/2019
header size: 375

offset to point data: 1205

number var. length records: 2

point data format: 6
point data record length: 30
number of point records: 0

number of points by return: 0 0 0 0 0

scale factor x y z: 0.01 0.01 0.01

offset x y z: 25500000 5000000 0
min x y z: 25593433.53 5028155.58 -125.22
max x y z: 25594452.33 5029174.38 113.36

start of waveform data packet record: 0

start of first extended variable length record: 0

number of extended_variable length records: 0

extended number of point records: 17514732

extended number of points by return: 14366216 2429785 563140 127123 23629 4178 661 0 0 0 0 0 0 0 0
variable length header record 1 of 2:

reserved 0
user ID 'LASF_Projection'
record ID 34735

length after header 40
description 'by LAStools of rapidlasso GmbH'
GeoKeyDirectoryTag version 1.1.0 number of keys 4

key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected

key 3072 tiff_tag_location 0 count 1 value_offset 8083 - ProjectedCSTypeGeoKey: NAD83(CSRS)v6 / MTM Nova Scotia zone 5

key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter

key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter

variable length header record 2 of 2:

reserved 0
user ID 'LASF_Projection'
record ID 2112

length after header 682
description 'by LAStools of rapidlasso GmbH'
WKT OGC COORDINATE SYSTEM:
PROJCS["NAD83(CSRS)v6 / NAD83(CSRS)v6 / MTM Nova Scotia zone 5",GEOGCS["NAD83(CSRS)v6",DATUM["North American Datum of 1983 (CSRS) version 6",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORIT
SG","7019"]],AUTHORITY["EPSG","1197"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","8252"]],PROJECTION["Transverse_
tor"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-64.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",25500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHO
"EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","8083"]]

LASzip compression (version 3.4r3 c3 50000): POINT14 3
LAStiling (idx 205963, lvl 9, sub 0, bbox 187000 4695000 699000 5207000) (size 1000 x 1000, buffer 0)
reporting minimum and maximum for all LAS point record entries ...

X 9343353 9445233
Y 2815558 2917438

Z -12522 11336
intensity 6 21432
return_number 1 7
number_of_returns 1 7
edge_of_flight_line 0 1
scan_direction_flag 0 1
classification 2 18
scan_angle_rank -33 32
user_data 0 255
point_source_ID 10011 30145
gps_time 246362877.637832 248873384.067001
extended_return_number 1 7
extended_number_of_returns 1 7
extended_classification 2 18
extended_scan_angle -5500 5333

extended_scanner_channel 1 3
number of first returns: 14366216

number of intermediate returns: 718552

number of last returns: 14366903
number of single returns: 11936939

overview over extended number of returns of given pulse: 11936939 3733305 1307530 413793 97960 20415 4790 0 0 0 0 0 0 0 0
histogram of classification of points:

2080670 ground (2)
6408505 low vegetation (3)
926788 medium vegetation (4)
8093469 high vegetation (5)
2712 noise (7)
2525 bridge deck (17)

63 Reserved for ASPRS Definition (18)

+-> flagged as withheld: 1186536

+---> 461025 of those are low vegetation (3)
+---> 57107 of those are medium vegetation (4)
+---> 665609 of those are high vegetation (5)
+---> 2712 of those are noise (7)
+---> 20 of those are bridge deck (17)
+---> 63 of those are Reserved for ASPRS Definition (18)

+-> flagged as extended overlap: 5071114

+---> 446228 of those are ground (2)
+---> 1995469 of those are low vegetation (3)
+---> 288884 of those are medium vegetation (4)
+---> 2339986 of those are high vegetation (5)
+---> 542 of those are noise (7)
+---> 5 of those are Reserved for ASPRS Definition (18)

E:\LAStools\bin>las2las -i 476_5026_201901.laz -epsg 2961 -remove_tiling_vlr -target_epsg 8083 -odix _8083 -olaz
SERIOUS WARNING: horizontal datum of source and target incompatible.

E:\LAStools\bin>lasinfo -i 476_5026_201901_8083.laz
lasinfo (201111) report for '476_5026_201901_8083.laz'

reporting all LAS header entries:

file signature: 'LASF'
file source ID: 0
global_encoding: 17
project ID GUID data 1-4: 00000000-0000-0000-0000-000000000000
version major.minor: 1.4

system identifier: 'LAStools (c) by rapidlasso GmbH'
generating software: 'las2las (version 201030)'

file creation day/year: 15/2019
header size: 375

offset to point data: 1205

number var. length records: 2

point data format: 6
point data record length: 30
number of point records: 0

number of points by return: 0 0 0 0 0

scale factor x y z: 0.01 0.01 0.01

offset x y z: 25500000 5000000 0
min x y z: 25593433.53 5028155.58 -125.22
max x y z: 25594452.33 5029174.38 113.36

start of waveform data packet record: 0

start of first extended variable length record: 0

number of extended_variable length records: 0

extended number of point records: 17514732

extended number of points by return: 14366216 2429785 563140 127123 23629 4178 661 0 0 0 0 0 0 0 0
variable length header record 1 of 2:

reserved 0
user ID 'LASF_Projection'
record ID 34735

length after header 40
description 'by LAStools of rapidlasso GmbH'
GeoKeyDirectoryTag version 1.1.0 number of keys 4

key 1024 tiff_tag_location 0 count 1 value_offset 1 - GTModelTypeGeoKey: ModelTypeProjected

key 3072 tiff_tag_location 0 count 1 value_offset 8083 - ProjectedCSTypeGeoKey: NAD83(CSRS)v6 / MTM Nova Scotia zone 5

key 3076 tiff_tag_location 0 count 1 value_offset 9001 - ProjLinearUnitsGeoKey: Linear_Meter

key 4099 tiff_tag_location 0 count 1 value_offset 9001 - VerticalUnitsGeoKey: Linear_Meter

variable length header record 2 of 2:

reserved 0
user ID 'LASF_Projection'
record ID 2112

length after header 682
description 'by LAStools of rapidlasso GmbH'
WKT OGC COORDINATE SYSTEM:
PROJCS["NAD83(CSRS)v6 / NAD83(CSRS)v6 / MTM Nova Scotia zone 5",GEOGCS["NAD83(CSRS)v6",DATUM["North American Datum of 1983 (CSRS) version 6",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY
SG","7019"]],AUTHORITY["EPSG","1197"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","8252"]],PROJECTION["Transverse_M
tor"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-64.5],PARAMETER["scale_factor",0.9999],PARAMETER["false_easting",25500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHOR
"EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","8083"]]

LASzip compression (version 3.4r3 c3 50000): POINT14 3

reporting minimum and maximum for all LAS point record entries ...

X 9343353 9445233
Y 2815558 2917438

Z -12522 11336
intensity 6 21432
return_number 1 7
number_of_returns 1 7
edge_of_flight_line 0 1
scan_direction_flag 0 1
classification 2 18
scan_angle_rank -33 32
user_data 0 255
point_source_ID 10011 30145
gps_time 246362877.637832 248873384.067001
extended_return_number 1 7
extended_number_of_returns 1 7
extended_classification 2 18
extended_scan_angle -5500 5333

extended_scanner_channel 1 3
number of first returns: 14366216

number of intermediate returns: 718552

number of last returns: 14366903
number of single returns: 11936939

overview over extended number of returns of given pulse: 11936939 3733305 1307530 413793 97960 20415 4790 0 0 0 0 0 0 0 0
histogram of classification of points:

2080670 ground (2)
6408505 low vegetation (3)
926788 medium vegetation (4)
8093469 high vegetation (5)
2712 noise (7)
2525 bridge deck (17)

63 Reserved for ASPRS Definition (18)

+-> flagged as withheld: 1186536

+---> 461025 of those are low vegetation (3)
+---> 57107 of those are medium vegetation (4)
+---> 665609 of those are high vegetation (5)
+---> 2712 of those are noise (7)
+---> 20 of those are bridge deck (17)
+---> 63 of those are Reserved for ASPRS Definition (18)

+-> flagged as extended overlap: 5071114

+---> 446228 of those are ground (2)
+---> 1995469 of those are low vegetation (3)
+---> 288884 of those are medium vegetation (4)
+---> 2339986 of those are high vegetation (5)
+---> 542 of those are noise (7)
+---> 5 of those are Reserved for ASPRS Definition (18)

To view this discussion on the web visit https://groups.google.com/d/msgid/lastools/CADm%3DQrQPGBUfaWN2eeJOJKwNgtKAYAo7KXi8nZpewVL%3DZdO9XA%40mail.gmail.com.

Reply all

Reply to author

Forward