Визуализация плотности облака точек с помощью PDAL

Question

У меня есть наземное лазерное облако точек сканирования, собранное с помощью системы координат сканера. Я хотел бы создать поверхность плотности, используя Pdal и запустив следующий код в оболочке OSGeo4W. Но получаю ошибку.

C:\>pdal density ^
More? /9A-1B_subset15m.las ^
More? -o /9A-1B_sub15m_den.sqlite ^
More? -f SQLite
(pdal density Error) GDAL failure (6) No translation for an empty SRS to PROJ.4 
format is known.

Можно ли вообще использовать PDAL для обработки наземного лидара с однократным возвратом? Я хотел бы создать равномерно распределенное облако точек (истончение), извлекать оголенную поверхность, устраняя шум.

Answer 1

На вопрос, может ли PDAL, вообще говоря, обрабатывать наземный лидар с однократным возвратом, ответ - да. PDAL не делает никаких предположений относительно того, доступны ли множественные возвраты. (Некоторые фильтры могут использовать возвращаемую информацию, но должны либо использовать по умолчанию какое-либо другое поведение, либо иным образом жаловаться, если возвращаемая информация недоступна.)

Что касается вашей ошибки, я бы предположил, что входное облако точек LAS не имеет назначенного SRS, и что это требуется для создания выходных данных SQLite. Если вы знаете SRS, вы можете назначить его, используя pdal translate (или, возможно, установив --readers.las.a_srs=<your SRS> в призыве к pdal density).