Реферат
Введение 4
1 Метод восстановления Пуассона 8
1.1 Основные сложности 10
2 Описание метода 14
2.1 Определение пространства функций 14
2.1.1 Регулярная сетка 15
2.1.2 Октодерево 16
2.2 Выбор сглаживающей функции 18
2.3 Определение векторного поля 19
2.4 Решение задачи Пуассона 23
2.4.1 Метод сопряженных градиентов 27
2.5 Получение изозначения 29
2.6 Последующее возможное построение поверхности 30
2.6.1 Метод марширующих кубов 30
3 Реализация 32
3.1 Входные данные 32
3.2 Выходные данные 33
3.3 Принцип работы программы 33
3.4 Примеры визуализованных входных и выходных данных 34
4 Результаты 36
4.1 Вариативность реализации, возможные улучшения 36
Заключение 37
Список использованных источников 38
Приложение A 40
На сегодняшний день лазерное сканирование уже широко применяется для изысканий на местности, в промышленных зонах, архитектуре. Характерной особенностью данного вида изысканий является тип генерируемых данных - гигантский массив геометрических точек, с известными координатами в пространстве и некоторой атрибутивной информацией (интенсивность сигнала, номер отражения, иногда цвет), также называемое облаком точек. Преимуществом такой съёмки является крайне высокая степень детализации, позволяющая получить представление об объекте удалённо.
В результате лазерного сканирования модели (архитектуры, рельефа земной поверхности с прилегающими объектами: машинами, светофорами, деревьями, ...) объектов трехмерного пространства вычисляются данные об этой модели, представленные облаком точек. Облака точек могут обрабатываться в различных целях, например, для выделения отдельных объектов, классификации объектов, нахождения плоских горизонтальных или вертикальных поверхностей (рельефа, стен зданий), и т.п. И как одну из основных целей в обработке облаков точек можно выделить восстановление поверхности.
Задача восстановления поверхности делится на основные этапы, которые представлены на рисунке 1. Каждый из этапов обрабатывает данные полученные в результате предыдущего этапа. На рисунке этапы идут от начального до заключительного в порядке слева направо.
Этапы восстановления поверхности:
1. На вход (анализ, удаление шума, классификация) подается облако точек, снятое с помощью лазерных сканеров, которое проходит алгоритмы предварительного анализа и обработки данных. Так, например, для совмещения облаков применяется алгоритм ICP, для сегментации облаков точек, например, с целью выделения плоскостей рельефа, и т. д.
2. На вход подается очищенное облако точек, для которого в дальнейшем предстоит определить поверхность. И на выходе должно выйти то же облако, но только практически для каждой из точек должны иметься рассчитанные направления нормалей. Вычисление направления нормалей можно выполнить, например, с помощью алгоритма principal component analysis (PCA) .
3. На вход подается очищенное от посторонних объектов облако точек с ориентированными нормалями. В результате обработки этого облака вычисляется функция-индикатор, которая показывает нахождение любой точки трехмерного пространства относительно поверхности модели (внутри или снаружи поверхности). Так же могут передаваться данные, параметры, необходимые на этапе построения контура, подробнее описание таких данных приводится в п.2.6.1.
4. При имеющейся функции-индикаторе и параметрах, ей сопутствующих, выполняется построение контура поверхности, например, с помощью метода марширующих кубов, подробнее в п.2.6.1.
В данном отчете будет рассматриваться этап 3: «Восстановление поверхности» на примере современного метода восстановления поверхности Пуассона . В результате чего, для очищенного от посторонних объектов облака точек с ориентированными нормалями будет вычислена функция-индикатор, определяющая расположение любой точки относительно поверхности модели на определенном изоуровне.
Были достигнута, поставленная в начале работы цель, и выполнены соответствующие задачи. А именно, написана программа, реализующая основные шаги метода восстановления поверхности Пуассона, по статье его автора . В п.4.1 описаны возможности для оптимзации метода.
1. Велижев А. В. Автоматическая сегментация облаков точек на основе элементов поверх
ности [Электронный ресурс] / Велижев А. В., Шаповалов Р. В., Потапов Д., Третьяк Л., Конушин А. - 01.2009. - URL:
http://www.graphicon.ru/html/2009/conference/se10/114/114_Paper.pdf (дата обращения: 24.05.2018).
2. Медведев В.И. Программы для обработки данных лазерного сканирования местности [Текст] / В.И. Медведев, Л.С. Райкова // САПР и ГИС автомобильных дорог. — 2017. — № 2(9). — С. 10-31. — DOI: 10.17273/CADGIS.2017.2.2.
3. Метод сопряжённых градиентов (для решения СЛАУ) [Электронный ресурс]: свободная энциклопедия / Wikipedia. - 20.03.2017. - URL: 1Ш|^:/т1.1к1рес11а.огдт'1к1А1е- тод_сопряжённых_градиентов_(для_решения_СЛАУ) (дата обращения: 24.05.2018).
4. Наземное лазерное сканирование [Электронный ресурс]: свободная энциклопедия / Wikipedia. - 24.05.2018. - URL: 1йр8://ги.1к1ред1а.огд/1к1/Наземное_лазерное_скани- рование (дата обращения: 27.05.2018).
5. Некипелов Н. Метод сопряженных градиентов — математический аппарат. - 09.08.2001. - URL: https://basegroup.ru/community/articles/conjugate (дата обращения: 24.05.2018).
6. Нормальное распределение [Электронный ресурс]: свободная энциклопедия / Wikipedia. - 11.01.2018. - URL: 1йр8://ги.1к1ред1а.огд/1к1/Нормальное_распределение (дата обращения: 25.05.2018).
7. Понятие обусловленности [Электронный ресурс]. -
URL: http://orloff.am.tpu.ru/chisl_metod/Lab3/cond.htm (дата обращения: 24.05.2018).
8. Сравнительный анализ прямых и итерационных методов решения систем линейных алгебраических уравнений [Электронный ресурс]. - 16.01.2016. - URL: http://entropiya- blog.ru/sravnitelnyj-analiz-pryamyx-i-iteracionnyx-metodov-resheniya-sistem-linejnyx- algebraicheskix-uravnenij.html (дата обращения: 23.05.2018).
9. 3D reconstruction [Electronic resource]: the free encyclopedia / Wikipedia. - 27.04.2018. - URL: https://en.wikipedia.org/wiki/3D_reconstruction (access date: 27.05.2018).
10. Alqudah A. Survey of Surface Reconstruction Algorithms [Electronic resource] // Journal of
Signal and Information Processing - 2014. - №5. - P. 63-79. -
URL: http://dx.doi.org/10.4236/jsip.2014.53009 (access date: 22.05.2018).
11. Berger M. State of the art in surface reconstruction from point clouds /Berger M., Tagliasacchi A. // Proc. Eurographics 2014 - State of the Art Reports. - 04.2014. - Vol. 1, № 1. - P. 161185.
12. Besl P. J. A Method for Registration of 3-D Shapes / Besl P. J., McKay N. D. // IEEE Transactions on Pattern Analysis and Machine Intelligence - Special issue on interpretation of 3-D scenes—part II. - 01.02.1992. - Vol 5, № 2. - P. 239-256. - URL: http://www-evasion.in- rialpes.fr/people/Franck.Hetroy/Teaching/ProjetsImage/2007/Bib/besl_mckay-pami1992.pdf (access date: 24.05.2018).
13. Bourk P. Polygonising a scalar field [Electronic resource]. - 05.1994. - URL: http://paulb- ourke.net/geometry/polygonise/ (access date: 22.05.2018).
14. Condition number [Electronic resource]: the free encyclopedia / Wikipedia. - 28.11.2017. - URL: https://en.wikipedia.org/wiki/Condition_number (access date: 24.05.2018).
15. Fattal R. Gradient domain high dynamic range compression / Fattal R., Lischinksi D., Wer- man M. // SIGGRAPH. - 2002. - P. 249-256.
... всего 23 источников