Не указана система координат для файлов las LiDAR

[rhot]

Скачал данные LiDAR с сайта http://lidar.cr.usgs.gov/ в формате las.

Во время проверки системы координат с помощью утилиты lasinfo, получаю:

Спойлер

---------------------------------------------------------
Header Summary
---------------------------------------------------------
File Name: salis000001.las
Version: 1.0
Source ID: 0
Reserved: 0
Project ID/GUID: '00000000-0000-0000-0000-000000000000'
System Identifier: ''
Generating Software: 'TerraScan'
File Creation Day/Year: 0/0
Header Size 227
Offset to Point Data 229
Number Var. Length Records 0
Point Data Format 0
Point Data Record Length 20
Number of Point Records 68692084
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.000000 -0.000000 -0.000000
Min X Y Z 335992.940000 3950699.710000 133.300000
Max X Y Z 341772.300000 3957675.350000 587.620000
Spatial Reference
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
End_Of_Tags.
Keyed_Information:
End_Of_Keys.
End_Of_Geotiff.

Видно, что система координат не указана.
В чём может быть проблема? Кто-нибудь пользовался данными с этого ресурса?

[Александр Мурый]

<lasinfo> откуда --- из состава LASTools или libLAS? Они могут работать немного по-разному..

Данными именно с этого ресурса не пользовался. Но вот недавно по работе пришлось иметь дело с LAS-файлами. Использовал libLAS, систему координат оно реально не показывает. В моём случае просто было известно, что это UTM определённой зоны.

Хотя здесь говорится, что надо исп-ть команду:

Код: Выделить всё

lasinfo --no-check srs.las

--> и увидим СК и проекцию

[rhot]

<lasinfo> откуда --- из состава LASTools или libLAS? Они могут работать немного по-разному..

Устанавливал пакет liblas-bin. Проверить сейчас не могу, т.к. данные на рабочей машине.

А как опция --no-check влияет на это?

версия lasinfo (version 1.2.1)
Проверял на файлах вот отсюда - СК определяет нормально. У моих же данных СК неизвестна из LAS файла. В архиве есть также xml файл, не знаю для чего он.

[Александр Мурый]

А как опция --no-check влияет на это?

С ней <lasinfo> выводит только содержимое заголовка файла, без анализа собственно точек в файле.

В архиве есть также xml файл, не знаю для чего он.

КО: может, он и содержит описание СК/проекции?

[rhot]

Проверил. Точно, этот xml-файл является сопутствующими метаданными в стандарте FGDC.
В итоге система координат оказалась следующая:

Спойлер

Spatial_Reference_Information:
Horizontal_Coordinate_System_Definition:
Planar:
Grid_Coordinate_System:
Grid_Coordinate_System_Name: Universal Transverse Mercator
Universal_Transverse_Mercator:
UTM_Zone_Number: 15
Transverse_Mercator:
Scale_Factor_at_Central_Meridian: 0.999600
Longitude_of_Central_Meridian: -93.000000
Latitude_of_Projection_Origin: 0.000000
False_Easting: 500000.000000
False_Northing: 0.000000
Planar_Coordinate_Information:
Planar_Coordinate_Encoding_Method: coordinate pair
Coordinate_Representation:
Abscissa_Resolution: .002
Ordinate_Resolution: .002
Planar_Distance_Units: North American Datum of 1983
Geodetic_Model:
Horizontal_Datum_Name: North American Datum of 1983
Ellipsoid_Name: Geodetic Reference System 80
Semi-major_Axis: 6378137.000000
Denominator_of_Flattening_Ratio: 298.257222
Vertical_Coordinate_System_Definition:
Altitude_System_Definition:
Altitude_Datum_Name: North American Vertical Datum of 1988
Altitude_Resolution: .01
Altitude_Distance_Units: feet
Altitude_Encoding_Method: Implicit coordinate

Вопрос решён, тема закрыта.