Эффективное понижение разрешения файлов LAS/LAZ с помощью Laspy: пошаговое руководство

Понимание процесса субдискретизации данных LAS с помощью Laspy

При работе с большими файлами LAS или LAZ в Python понижающая дискретизация необходима для эффективной обработки и анализа. Ласпи, пакет Python для чтения, записи и изменения данных LAS, предлагает множество способов манипулирования данными облака точек, таких как создание и редактирование заголовков LAS.

В этом примере показано, как уменьшить выборку набора данных, извлекая каждую десятую точку из лаз файл и повторно использовать существующий Ласхедер. Это требует понимания того, как заголовки взаимодействуют с данными, особенно при работе с различным количеством точек.

При создании нового ЛасДата объект из существующего заголовка, пользователи часто сталкиваются с несоответствием размера массива. Это несоответствие возникает потому, что заголовок point_count может не соответствовать автоматически новым данным.

Задача состоит в том, чтобы оценить, требуется ли вручную изменять атрибуты заголовка, такие как смещения, масштабы и point_count, или существует более автоматическое решение. В этом посте объясняется, как правильно обновить эти значения при понижении разрешения с помощью Ласпи, что приводит к эффективному процессу.

Оптимизация понижения разрешения облака точек с помощью Laspy

Первый скрипт в этом посте фокусируется на понижении разрешения ЛАЗ файл, который необходим для управления большими наборами данных облаков точек. Импортировав исходный файл с помощью laspy.read() мы можем получить доступ к данным точек и заголовку, который содержит метаданные об облаке точек. Метод пониженной выборки предполагает выбор каждой десятой точки, что минимизирует размер набора данных, сохраняя при этом важные географические свойства. Это делается с помощью np.arange() построить массив индексов. После выбора точек скопируйте заголовок из исходного файла, чтобы обеспечить совместимость метаданных, например point_format и версия.

Однако распространенная проблема возникает, когда количество точек в исходном заголовке не соответствует данным с пониженной дискретизацией. Чтобы исправить это, мы используем копировать() функция, позволяющая сделать неглубокую копию исходного заголовка и вручную изменить point_count поле, чтобы отразить количество точек пониженной дискретизации. После создания нового заголовка точки пониженной дискретизации распределяются по новому ЛасДата объект, содержащий реальные координаты x, y и z. Наконец, ЛасДата сохраняется как новый файл LAZ с помощью писать() метод. Этот сценарий эффективен для пользователей, которым необходимо извлечь меньшие наборы данных из более крупных облаков точек.

Второй скрипт расширяет первый, автоматически пересчитывая смещения и масштабы для данных с пониженной дискретизацией. При работе с облаками точек наличие точных смещений имеет решающее значение, поскольку они указывают происхождение данных в трехмерном пространстве. заголовок.смещения Атрибут обновляется с использованием минимальных координат x, y и z из точек с пониженной дискретизацией. Аналогично, масштабные коэффициенты, влияющие на точность точечных данных, устанавливаются с помощью заголовок.масштабы атрибут. Этот сценарий не только минимизирует размер облака точек, но также гарантирует точность и выравнивание данных, что делает их более подходящими для практического использования.

Наконец, окончательный сценарий демонстрирует модульное тестирование с использованием Python. юниттест рамки. В этом сценарии тестовый пример определяет, соответствует ли количество точек пониженной дискретизации прогнозируемому значению. Это имеет решающее значение для обеспечения того, чтобы процедура понижающей выборки выполнялась последовательно в разных контекстах и ​​наборах данных. Тестовый пример определяется с помощью Тесткейс класс, и сравнение проводится с использованием self.assertEqual() метод. Включив тестирование в рабочий процесс, мы можем гарантировать, что процедура понижения разрешения работает правильно, прежде чем развертывать ее в более крупных проектах или конвейерах. Этот сценарий помогает пользователям избежать проблем и несоответствий при работе с несколькими файлами облаков точек.

import laspy
import numpy as np
from copy import copy
# Load the existing LAZ file
las = laspy.read("input_file.laz")
# Downsample by taking every 10th point
indices = np.arange(0, len(las.points), 10)
downsampled_points = las.points[indices]
# Copy the header and adjust the point count
header = copy(las.header)
header.point_count = len(downsampled_points)
# Create new LasData with downsampled points
d_las = laspy.LasData(header)
d_las.points = downsampled_points
# Write to a new LAZ file
d_las.write("downsampled_output.laz")

import laspy
import numpy as np
# Load the original LAZ file
las = laspy.read("input_file.laz")
# Downsample by taking every 10th point
indices = np.arange(0, len(las.points), 10)
downsampled_points = las.points[indices]
# Create new header and adjust offsets/scales
header = laspy.LasHeader(point_format=las.header.point_format, version=las.header.version)
header.offsets = np.min([las.x[indices], las.y[indices], las.z[indices]], axis=1)
header.scales = np.array([0.01, 0.01, 0.01])  # Set new scales
# Create new LasData and write to file
downsampled_las = laspy.LasData(header)
downsampled_las.points = downsampled_points
downsampled_las.write("downsampled_with_scales.laz")

import unittest
import laspy
import numpy as np
class TestDownsampling(unittest.TestCase):
    def test_downsample_point_count(self):
        las = laspy.read("input_file.laz")
        indices = np.arange(0, len(las.points), 10)
        downsampled_points = las.points[indices]
        self.assertEqual(len(downsampled_points), len(indices))
if __name__ == "__main__":
    unittest.main()

Обработка метаданных файла LAS и расширенные методы понижения разрешения

При работе с огромными наборами данных с лесбиянкаУправление метаданными столь же важно, как и управление фактическими данными облака точек. Поддержание точности Ласхедер значений после понижения дискретизации представляет собой значительную трудность. Поскольку координаты данных облака точек (x, y и z) изменяются, заголовок должен отражать эти изменения. Перерасчет смещения требует пересчета минимальных значений для каждого измерения, тогда как весы определить точность точечных данных, особенно для хранения.

Еще одним фактором, который необходимо оценить, является целостность дополнительных измерений в файле LAS. Дополнительные байты обычно используются для хранения информации, отличной от обычных координат x, y и z, например интенсивности или времени GPS. Если набор данных содержит эти дополнительные измерения, их необходимо обработать при понижении разрешения. Вы должны гарантировать, что количество точек в дополнительных измерениях соответствует уменьшенному количеству точек в первичных данных. add_extra_dim функциональность в лесбиянка позволяет добавлять пользовательские размеры в заголовок LAS.

Наконец, оптимизация скорости является важным фактором, который следует учитывать при субдискретизации облаков точек. Хотя обычно требуется внесение изменений в заголовок человеком, автоматизация процесса осуществляется за счет эффективного индексирования и применения операций с массивами с помощью бестолковый может значительно ускорить этот процесс. Используя силу бестолковый, вы можете быстро управлять огромными наборами данных, не жертвуя при этом производительностью. Это позволяет расширять решения для более крупных проектов или даже автоматизировать конвейеры для обработки нескольких файлов LAZ.