EPSG Help: The Highlands at Dove Mountain Golf Course

[rocket2]

Apparently I'm probably not using the lastools properly.  I'm using the app window with the point and click options and not the command line approach.  I tried many searches to find your approach to using lastools to get the right EPSG.  The only thing I found was your 24 step process for course development that didn't contain details for finding the right EPSG.  Can you point me in the right direction to get your detail way of using LAStools?  Thanks for all your help.

[mal]

In the lastools package there are readme files for each program. You want to use lasinfo.exe to create a text file of the laz file header, it is a command prompt program, and there is a corresponding readme that instructs you how to use it. Once you have the text file, you open it with notepad.

[BaconJunkie1]

You are killing me Smalls This should help

tgctours.proboards.com/thread/29827/use-lastools-lasinfo-epsg-code

[mal]

Nov 8, 2021 11:02:36 GMT -5 BaconJunkie1 said:

You are killing me Smalls This should help

tgctours.proboards.com/thread/29827/use-lastools-lasinfo-epsg-code

or you can read the readme file and make a batch file that automatically does this for you.
put lasinfo.exe in a folder
make a txt file named lidar.bat in this folder and then edit it
paste the following four lines and then save the .bat file

ren *.laz *.las
lasinfo -i *.las -o lidar_info.txt
START "" notepad.exe lidar_info.txt
ren *.las *.old 

(or you can put "del *.las" in place of the "ren *.las *.old" if you would rather delete the .las file after you have created the lidar_info.txt file)

From then on, you just copy/paste a las or laz file into this folder and double click the lidar.bat



You can also add commands to automagically find and highlight the correct EPSG code if you wanted to... but why?

[rocket2]

I know it all sounds easy for you guys; but for me it'll be a satisfying accomplishment. So I'll dive into it and make it work. Thanks!

[rocket2]

OK I went with mal's approach and got results.  But these results were the same that I got by running the lasinfo windows gui screen.   It's obvious with the results on the AZ LAZ info the X & Y EPSG is 6405 (see AZ txt file below).  But when I run it with the WI LAZ file I can't find the X & Y EPSG to save my life (see WI txt file below).  Am I blind, nuts or just plain can't read???  Or all three!?  I've also included the bat file in case something in there is not right.  Please take a quick look and get a good laugh ( at my expense) for the obvious mistake I'm making.       Again, thanks for all your help and patience.

Well I guess I can't even post a simple text file so here the full text from each:

AZ info txt:

lasinfo (210720) report for 'USGS_LPC_AZ_Eastern_PimaCO_2015_955523_LAS_2017.las'

reporting all LAS header entries:

file signature: 'LASF'

file source ID: 0

global_encoding: 17

project ID GUID data 1-4: 1DC295D3-7619-4035-F78C-BF1C04E7360F

version major.minor: 1.4

system identifier: 'LAStools (c) by rapidlasso GmbH'

generating software: 'las2las (version 150622)'

file creation day/year: 307/2015

header size: 375

offset to point data: 2764

number var. length records: 3

point data format: 6

point data record length: 30

number of point records: 11115069

number of points by return: 10791453 323230 386 0 0

scale factor x y z: 0.01 0.01 0.01

offset x y z: 0 0 0

min x y z: 955112.76 523427.81 2548.69

max x y z: 960112.75 528427.80 2812.13

start of waveform data packet record: 0

start of first extended variable length record: 48632242

number of extended_variable length records: 1

extended number of point records: 11115069

extended number of points by return: 10791453 323230 386 0 0 0 0 0 0 0 0 0 0 0 0

variable length header record 1 of 3:

reserved 43707

user ID 'LASF_Spec'

record ID 7

length after header 459

description 'OGC WKT Coordinate System'

variable length header record 2 of 3:

reserved 43707

user ID 'NIIRS10'

record ID 4

length after header 10

description 'NIIRS10 Timestamp'

variable length header record 3 of 3:

reserved 43707

user ID 'NIIRS10'

record ID 1

length after header 26

description 'NIIRS10 Tile Index'

extended variable length header record 1 of 1:

reserved 0

user ID 'LASF_Projection'

record ID 2112

length after header 832

description 'WKT Projection'

OGC COORDINATE SYSTEM WKT:



PROJCS["NAD83(2011) / Arizona Central (ft)",

GEOGCS["NAD83(2011)",

DATUM["NAD83 (National Spatial Reference System 2011)",

SPHEROID["GRS 1980",6378137,298.257222101,

AUTHORITY["EPSG","7019"]],

AUTHORITY["EPSG","1116"]],

PRIMEM["Greenwich",0,

AUTHORITY["EPSG","8901"]],

UNIT["degree",0.0174532925199433,

AUTHORITY["EPSG","9122"]],

AUTHORITY["EPSG","6318"]],

PROJECTION["Transverse_Mercator"],

PARAMETER["latitude_of_origin",31],

PARAMETER["central_meridian",-111.9166666666667],

PARAMETER["scale_factor",0.9999],

PARAMETER["false_easting",700000],

PARAMETER["false_northing",0],

UNIT["foot",0.3048,

AUTHORITY["EPSG","9002"]],

AXIS["X",EAST],

AXIS["Y",NORTH],

AUTHORITY["EPSG","6405"]]

the header is followed by 1732 user-defined bytes

LASzip compression (version 3.1r0 c3 50000): POINT14 3

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

X 95511276 96011275

Y 52342781 52842780

Z 254869 281213

intensity 256 65280

return_number 1 3

number_of_returns 1 3

edge_of_flight_line 0 1

scan_direction_flag 0 1

classification 1 17

scan_angle_rank -27 19

user_data 0 1

point_source_ID 208 210

gps_time 108560714.319444 108561929.861713

extended_return_number 1 3

extended_number_of_returns 1 3

extended_classification 1 17

extended_scan_angle -4500 3167

extended_scanner_channel 0 0

number of first returns: 10791453

number of intermediate returns: 386

number of last returns: 10791442

number of single returns: 10468212

WARNING: point type is 6 but (legacy) number of point records in header is 11115069 instead zero.

WARNING: point type is 6 but (legacy) number of points by return [1] in header is 10791453 instead zero.

WARNING: point type is 6 but (legacy) number of points by return [2] in header is 323230 instead zero.

WARNING: point type is 6 but (legacy) number of points by return [3] in header is 386 instead zero.

overview over extended number of returns of given pulse: 10468212 645699 1158 0 0 0 0 0 0 0 0 0 0 0 0

histogram of classification of points:

6296097 unclassified (1)

4801681 ground (2)

17158 noise (7)

9 rail (10)

124 bridge deck (17)

+-> flagged as withheld: 17158

+---> 17158 of those are noise (7)

+-> flagged as extended overlap: 3852687

+---> 3852687 of those are unclassified (1)


WI info txt:

lasinfo (210720) report for 'USGS_LPC_WI_Racine_2017_LAS_11752100_LAS_2019.las'

reporting all LAS header entries:

file signature: 'LASF'

file source ID: 0

global_encoding: 17

project ID GUID data 1-4: A3D933DA-63E9-4972-0000-000000000000

version major.minor: 1.4

system identifier: 'NIIRS10'

generating software: 'GeoCue LAS Updater'

file creation day/year: 78/2019

header size: 375

offset to point data: 1471

number var. length records: 1

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: 0 0 0

min x y z: 1117500.00 2210000.00 646.99

max x y z: 1120000.00 2212500.00 792.60

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: 15929734

extended number of points by return: 15929734 0 0 0 0 0 0 0 0 0 0 0 0 0 0

variable length header record 1 of 1:

reserved 0

user ID 'LASF_Projection'

record ID 2112

length after header 1040

description 'OGC WKT Coordinate System'

WKT OGC COORDINATE SYSTEM:

COMPD_CS["NAD83(2011) / Illinois East (ftUS) + NAVD88 height - Geoid12B (ftUS)",PROJCS["NAD83(2011) / Illinois East (ftUS)",GEOGCS["NAD83(2011)",DATUM["NAD83_National_Spatial_Reference_System_2011",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","1116"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","6318"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",36.66666666666666],PARAMETER["central_meridian",-88.33333333333333],PARAMETER["scale_factor",0.999975],PARAMETER["false_easting",984250.0000000002],PARAMETER["false_northing",0],UNIT["US survey foot",0.3048006096012192,AUTHORITY["EPSG","9003"]],AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","6455"]],VERT_CS["NAVD88 height - Geoid12B (ftUS)",VERT_DATUM["North American Vertical Datum 1988",2005,AUTHORITY["EPSG","5103"]],UNIT["US survey foot",0.3048006096012192,AUTHORITY["EPSG","9003"]],AXIS["Gravity-related height",UP],AUTHORITY["EPSG","6360"]]]

the header is followed by 2 user-defined bytes

LASzip compression (version 3.1r0 c3 50000): POINT14 3

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

X 111750000 112000000

Y 221000000 221250000

Z 64699 79260

intensity 0 56855

return_number 1 1

number_of_returns 1 1

edge_of_flight_line 0 0

scan_direction_flag 0 0

classification 1 18

scan_angle_rank 15 15

user_data 0 0

point_source_ID 0 0

gps_time 176455728.000000 176456896.000000

extended_return_number 1 1

extended_number_of_returns 1 1

extended_classification 1 18

extended_scan_angle 2500 2500

extended_scanner_channel 0 0

number of first returns: 15929734

number of intermediate returns: 0

number of last returns: 15929734

number of single returns: 15929734

overview over extended number of returns of given pulse: 15929734 0 0 0 0 0 0 0 0 0 0 0 0 0 0

histogram of classification of points:

11525337 unclassified (1)

4398558 ground (2)

4273 noise (7)

1566 Reserved for ASPRS Definition (18)



Bat file:

ren *.laz *.las

lasinfo -i *.las -o lidar_info.txt

START "" notepad.exe lidar_info.txt

ren *.las *.old

[mal]

You are blind. do a find (control f) and search for AXIS["Y",NORTH]. Should then be the next EPSG code.

[rocket2]

Operator error on my part. I assumed the line with all the EPSG info ended at about a full screens length; didn't realize the line was so long. I'm officially embarrassed! Thanks again. I'll try to stay out of your hair for awhile.