Научная статья на тему 'ТЕОРЕМЫ СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ РЕШЕНИЯ ОБРАТНЫХ ЗАДАЧ ПРОЕКТИВНОЙ ГЕОМЕТРИИ ДЛЯ 3D РЕКОНСТРУКЦИИ ПО ФОТОСНИМКАМ'

ТЕОРЕМЫ СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ РЕШЕНИЯ ОБРАТНЫХ ЗАДАЧ ПРОЕКТИВНОЙ ГЕОМЕТРИИ ДЛЯ 3D РЕКОНСТРУКЦИИ ПО ФОТОСНИМКАМ Текст научной статьи по специальности «Математика»

CC BY
90
13
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Чебышевский сборник
Scopus
ВАК
RSCI
Область наук
Ключевые слова
ЦЕНТРАЛЬНАЯ ПРОЕКЦИЯ / 3D РЕКОНСТРУКЦИЯ / ГЕОМЕТРИЯ ТРЕУГОЛЬНИКА / CENTRAL PROJECTION / 3D RECONSTRUCTION / TRIANGLE GEOMETRY

Аннотация научной статьи по математике, автор научной работы — Клячин Алексей Александрович, Клячин Владимир Александрович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

EXISTENCE AND UNIQUENESS THEOREMS FOR SOLUTIONS OF INVERSE PROBLEMS OF PROJECTIVE GEOMETRY FOR 3D RECONSTRUCTION FROM PHOTOGRAPHS

The paper considers the problem of calculating the parameters of the plane of a spatial triangle from its central projection. Under certain conditions, the existence theorem for a solution to this problem and its uniqueness are proved. Examples of conditions under which a solution does not exist or is not unique are given. An algorithm for the approximate search of all possible solutions to the problem under certain conditions is also proposed. The problem considered in the article arises when constructing three-dimensional models of objects from their photograph.

Текст научной работы на тему «ТЕОРЕМЫ СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ РЕШЕНИЯ ОБРАТНЫХ ЗАДАЧ ПРОЕКТИВНОЙ ГЕОМЕТРИИ ДЛЯ 3D РЕКОНСТРУКЦИИ ПО ФОТОСНИМКАМ»

ЧЕБЫШЕВСКИЙ СБОРНИК

Том 21. Выпуск 4.

УДК 514.142.2+514.174.6 DOI 10.22405/2226-8383-2020-21-4-117-128

Теоремы существования и единственности решения обратных задач проективной геометрии для 3D реконструкции по фотоснимкам1

А. А. Клячин, В. А. Клячин

Алексей Александрович Клячин — доктор физико-математических наук, доцент, Волгоградский государственный университет (г. Волгоград). e-mail: aleksey.klyachin@volsu.ru

Владимир Александрович Клячин — доктор физико-математических наук, доцент, Волгоградский государственный университет (г. Волгоград). e-mail: klchnv@mail.ru

Аннотация

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

Ключевые слова: центральная проекция, 3D реконструкция, геометрия треугольника.

Библиография: 22 названия.

Для цитирования:

А. А.Клячин, В. А. Клячин. Теоремы существования и единственности решения обратных задач проективной геометрии для 3D реконструкции по фотоснимкам // Чебышевский сборник, 2020, т. 21, вып. 4, с. 117-128.

1 Исследование выполнено при финансовой поддержке РФФИ и Администрации Волгоградской области в рамках научного проекта № 19-47-340015, а также при финансовой поддержке Минобрнауки России в рамках научного проекта № 0633-2020-0004 «Развитие методики виртуальной ЗБ реконструкции исторических объектов»

CHEBYSHEVSKII SBORNIK Vol. 21. No. 4.

UDC 514.142.2+514.174.6 DOI 10.22405/2226-8383-2020-21-4-117-128

Existence and uniqueness theorems for solutions of inverse

problems of projective geometry for 3D reconstruction

from photographs

A. A. Klvachin, V. A. Klvachin

Alexey Alexandrovich Klyachin — Doctor of Physical and Mathematical Sciences, Associate Professor, Volgograd State University (Volgograd). e-mail: aleksey.klyachin@volsu.ru

Vladimir Aleksandrovich Klyachin^ Doctor of Physical and Mathematical Sciences, Associate Professor, Volgograd State University (Volgograd). e-mail: klchnv@mail.ru

Abstract

The paper considers the problem of calculating the parameters of the plane of a spatial triangle from its central projection. Under certain conditions, the existence theorem for a solution to this problem and its uniqueness are proved. Examples of conditions under which a solution does not exist or is not unique are given. An algorithm for the approximate search of all possible solutions to the problem under certain conditions is also proposed. The problem considered in the article arises when constructing three-dimensional models of objects from their photograph.

Keywords: central projection, 3D reconstruction, triangle geometry. Bibliography: 22 titles.

For citation:

A. A. Klyachin, V. A. Klyachin, 2020, "Existence and uniqueness theorems for solutions of inverse problems of projective geometry for 3D reconstruction from photographs", Chebyshevskii sbornik, vol. 21, no. 4, pp. 117-128.

1. Введение

Одна из сложных задач 3D моделирования - моделирование по плоскому изображению, различного рода фото и видео материалам реальных объектов. Это особенно актуально в задачах 3D реконструкции исторических архитектурных сооружений и целых комплексов, имеющих как историческую, так и культурную ценность. Подобные задачи реконструкции предполагают, в частности, решение ряда задач нацеленных на развитие методов и технологий 3D моделирования по плоскому изображению объекта при отсутствии какой-либо возможности доступа к реальному объекту. Здесь возникает ряд геометрических задач, связанных с особенностями и условиями подобных 3D реконструкций. Примерами таких задач могут быть: задачи определения характерных размеров реальных объектов и их составных частей, задачи определения их местоположения в пространстве (например, в терминах географических координат), задачи определения настроек камеры (ориентация камеры в пространстве), с помощью которой было получено изображение. Фундаментальная научно-техническая дисциплина, занимающаяся определением формы, размеров, положения и иных характеристик

объектов по их фотоизображениям носит название фотограмметрия. Однако существенным отличием классических задач фотограмметрии (см., например [1],[2]) от задач рассматриваемых нами является практически неизвестные характеристики элементов внутреннего и внешнего ориентирования снимков (настроек камеры). Частичное нахождение этих характеристик и составляет определенную задачу настоящей статьи. В основе предлагаемой методики лежит поиск характерных наборов точек снимка, которые определяются существенно исходя из информации геометрического строения объекта съемки. Например, для архитектурных сооружений характерными точками могут быть угловые точки зданий, оконных проемов, точки, лежащие на параллельных прямых, пространственных окружностях и т.п.

В настоящей статье предпринята попытка решения задачи вычисления параметров плоскости пространственного треугольника по его центральной проекции. Надо отметить, что несколько подобных задач ранее были решены в статье [3].

Предполагается, что найденные решения могут быть использованы в задачах ЗГ) реконструкции архитектурных комплексов довоенного Сталинграда по имеющимся фотоснимкам. Описание задачи этой реконструкции было дано в работе [4].

В плане задач определения тех или иных характеристик объектов по их снимкам имеется ряд публикаций, в которых решаются прикладные задачи в самых разнообразных отраслях науки и техники. Так в работах [5] -[7] исследуется совокупность характерных точек на изображении, образующая размытость. В работах [8] - [11] решается задача вычисления геометрических характеристик объекта по нескольким специально подготовленным его снимкам. В работах [12], [13] анализ характерных точек применяется в задачах ЗБ реконструкции в медицине. Отметим, что похожие задачи возникают в прикладной области компьютерного зрения. В частности, можно указать интересную работу [14], в которой решается задача 3Г) реконструкции лица по фотографии. Методика реконструкции основана на использовании обученной сверточной нейронной сети. В дополнении к этому отметим еще ряд публикаций, посвященных восстановлению трехмерных поверхностей по их плоскому изображению [15] -[18] на основе методов машинного обучения и накопленной базы данных изображений с данными буфера глубины.

В [19] предложен метод восстановления дефектных областей карт глубины, полученных 3 Г) сканерами при анализе реальных объектов сцены для подготовки к реконструкции трёхмерных моделей поверхностей объектов.

В [20] рассмотрен класс задач о реконструкции внутренней структуры плоских и объемных объектов по внешним данным, находящимся в определенной структуре и имеющим возможность оцифровки. Так же отметим работы [21], [22], в которых решаются задачи восстановления пространственных множеств точек и поверхностей в несколько иной постановке задачи, чем в настоящей работе.

Всюду в статье используется модель камеры, состоящая из трех уравнений центральной проекции

(х,у,г) ^ (Х,¥) = (^Л),

где (х, у, г) - точка в пространстве, величины X, У моделируют пиксельные координаты на изображении, которое располагается в плоскости г = 5,5 > 0. Предельным случаем центральной проекции является ортогональная проекция, заданная формулами

(х,у, г) ^ (X, V) = (х, у)

2. Математическая модель

В начале мы рассмотрим случай, когда на входном изображении предполагается наличие изображения пространственного треугольника с заданными углами = 0,1, 2. Задача состоит в том, чтобы определить параметры плоскости этого треугольника.

Пусть в пространстве заданы три точки р^ = (хг, х^, г^ > 0, г = 0,1, 2 не лежащие на одной прямой и являющиеся вершинами треугольника. При этом, мы предполагаем, что нумерация соответствует некоторому обходу вдоль его сторон. Обозначим через П - плоскость этого треугольника. Рассмотрим так же плоскость ф, заданную уравнением г = 1. Построим точки д^, г = 0,1, 2, 3 как проекции точек р^ на плоскость ф при центральной проекции с центром в начале координат. Положим qi = (Хг, Уг, Е^). Несложно видеть, что

X, = ^

¿г

У = ж

г -г,-

гг = i

Поскольку треугольник qoqiq2 считается известным, то известными будут величины углов ao = Z qoOqi, ai = Z qiOq2, a2 = Z q2Oqo. Для углов в треугольнике Р0Р1Р2 введем обозначения

Po = ZP2P0P1, Pi = ZP0P1P2, P2 = ZP1P2P0.

Величины этих углов так же считаются известными. Ясно, что ориентация плоскости этого треугольника может быть определена с точностью до гомотетичного преобразования относительно начала координат: при этом центральная проекция треугольника не измениться. Так, что достаточно определить радиус-векторы точек Po,Pi,P2 также с точностью до преобразования гомотетии. Для определенности будем считать, что углы в треугольнике пронумерованы так, что р0 > Pi > P2.

Согласно теореме синусов имеют место равенства

| PoPil = | P0P2I = | PiP2l sin P2 sin Pi sin P0 '

Откуда

l PoPil = Sin P2 = x цч

| PiP2l sin Po ,

lPoP2l = sin Pi = \PVP2I sin Po .

В силу выбора нумерации углов треугольника Apopip2 эти величины удовлетворяют неравенствам \,р, < 1.

По теореме косинусов можно записать равенства

2 2 2

+ |pi|2 - 2|poIIpiI cosao = ipopil Ipil2 + |p212 - 2|p 11|p21 cos ai = |pml2, |p212 + |po|2 - 2|p2||po| cos a2 = ip2poi2.

Учитывая обозначения для величин Л, y, получим

Í I p o|2 + |pi|2 - 2| po||pi| cosao = Л2(| p 112 + |p212 - 2| p 1 ||p21 cos ai), 1 | p 212 + | p o|2 - 2| p 21| po | cos (12 = y2(|p i|2 + | p 2|2 - 2| p 11| p21 cos ai).

Полагая u = |pi\/\ро\, v = |Р2\/\Ро\, приходим к следующей системе квадратичных уравнений

J 1 + и2 — 2и cos a0 = X2 (и2 + v2 — 2uv cos a\), \ V2 + 1 — 2v cos а2 = ß2(u,2 + v2 — 2uv cos a\).

И, окончательно,

(1 — X2)u2 — X2v2 + 2X2uv cos a\ — 2u cos a0 + 1 = 0, —¡j?u2 + (1 — ß2)v2 + 2/j?uv cos a\ — 2v cos a2 + 1 = 0.

В следующем разделе будет приведен алгоритм поиска всех приближенных решений этой системы уравнений. Перед этим мы докажем теорему, в которой будут представлены достаточные условия существования и единственности рассматриваемой задачи. С этой целью мы ее немного переформулируем. Именно, требуется для треугольника Арор\р2 с углами @о,@1, @2 и заданного набора величин углов ао, а\, а2 найти точку р е М3, такую, что

¿Р0РР1 = ао, ¿р\рр 2 = а1, Ар2рро = а2. (4)

Имеет место

Теорема 1. Если для, треугольника Ар0р1р2 с углами г = 0,1,2 выполнено Д < к/2, а, заданные величины а удовлетворяют условиям а^ > к/2 и а0 + а1 + а2 < 2тг, то найдется единственная точка р такая, что выполнены равенства (4).

Доказательство. Для каждой стороны рiрj,г,] = 0,1,2,г < ] заданного треугольника построим множество точек

Бц = [де М3 : Ар^хр^ = аг},

относительно плоскости треугольника Арор1р2, что и начадо координат. Множество представляет собой половину поверхности вращения дуги окружности, опирающейся на соответствующую сторону этого треугольника вокруг этой стороны. Угловая величина этой дуги в точности равна соответствующей величине а^. Поскольку аг > к/2, то эти поверхности являются выпуклыми и их можно рассматривать как графики функций Д,- (д), 0 < г < ] < 1, где точка д принадлежит плоскости треугольника Арор1р2■ У каждой такой функции своя область определения, которая ограничена двумя дугами окружностей, опирающихся на соответствующую сторону р^ треугольника Арорда. Заметим, что искомая точка р является точкой пересечения трех поверхностей Этой точке соответствует решение системы уравнений

Г ¡12(д) = fоl(g), г.,

I Ыд) = Ыд).

Заметим, что пересечение областей определения функций Д,- не пусто только если ао+а1+а2 < < 2к. Причем, если сумма этих углов в точности равна 2к, то пересечением является единственная точка. Она же в этом случае будет искомой. Так, что в дальнейшем считаем, что выполнено строгое неравенство а0 + а1 + а2 < 2тг. Тогда пересечением областей определения функций Д,- является выпуклая область О С П ограниченная тремя дугами окружно-

А о 1 2

ствующее угловое значение а^ и расположенных в той же полуплоскости относительно сторон треугольника Аро'р1'р2, что и сам треугольник. Пусть до, д1, д2 - точки, в которых указанные дуги окружностей соединяются. Нумерация точек выполнена так, что бы Д,-(х) = 0, если х принадлежит дуге ад.,. Покажем, что внутри О имеется по крайней мере одно решение системы (5). Действительно, заметим, что на дуге д1д2 найдется точка д', в которой !01(д') = 1'02(д') > 0. Рассмотрим непрерывную кривую ^(Ъ) вдоль которой ¡'01 = /02 и такую,

что 7(0) = ро,7(1) = ч'■ Положим <p(t) = /oi(t(î)) = f02(l(t)) ш ф(£) = f12(j(t)). Поскольку Ф(0) > 0,ф(1) = 0 и ^>(0) = 0,^(1) > 0, то найдется такое to, что будет выполнено равенство Ф(£0) = p(t0) Не трудно видеть, что точка 7(t0) является решением системы (5). Положим

Dc = {х : foi(x) > с, fo2(x) > с, /12(ж) > с}.

Эти множества являются выпуклыми, причем Dcu с Dc' при с' < с". Ясно, что Do = ^ содержит все решения системы (5). Найдется такое Со > 0 чт0 ПРИ всяком 0 < с < Со множество Dc будет содержать все решения системы (5). В дальнейшем считаем, что 0 < с < со-

Рассмотрим множества

Kl = {х : foi > с, /12 > с},

е<02 = {х : foi > с, fo2 > с},

Ес!2 = {х : ¡02 > С, f12 > с}.

Каждое из этих множеств выпукло и имеет границу состоящую из двух выпуклых дуг линий уровня соответствующих функций. Эти линии уровня явяются выпуклыми кривыми, опирающимися на соответствующие стороны треугольника ApoPiP2- Поэтому они пересекаются в двух точках. В результате мы получаем шесть точек пересечения линий уровня функций Д,-, причем три из них ограничивают дуги образующих криволинейный треугольник, являющийся границей выпуклого множества Dc. Рассмотрим выпуклое множество

D0 = П0<c<c0Dc.

Это множество замкнуто и выпукло. Предположим, что его внутренность не пуста. Тогда Do = DC0 и по выбору значения со на границе этого множества имеется решение q системы (5). Получаем, что с одной стороны Do = DC0 и граница этого множества состоит из трех выпуклых дуг линий уровня функций fij (х), а сдругой стороны эти дуги перескаются в точке q1 G dDo на его границе. Получим противоречие с предположением, что внутренность множества Do не пуста. Таким образом Do состоит из единственной точки q' или является отрезком. Последнее невозможно, так как функции fij строго выпуклы вверх. Таким образом, q' - единственное решение системы (5). Теорема доказана.

Замечание 1. Из этой теоремы следует однозначное восстановление плоскости пространственного треугольника с точностью до гомотетии по его проекции, при условии, что углы ai удовлетворяют указанным в теореме требованиям. Действительно, предположим, что можно найти два равных треугольника, но расположенных в различных плоскостях, не сводящихся одна к другой преобразованием гомотетии относительно начала координат и имеющих совпадающие проекции в плоскости z = 1. Таким образом мы имеем два тетраэдра в основаниях которых лежат выбранные треугольники, а вершинами является начало координат. Движением совместим основания этих тетраэдров, так, что бы вершины оказались по одно сторону от общей плоскости их оснований. Эти вершины не будут совпадать поскольку треугольники лежали в разных плоскостях. Но, с другой стороны, это противоречит утверждению теоремы 1. Таким образом предположении о существовании двух указанных выше треугольников ошибочно.

Замечание 2. В общем случае теорема существования не справедлива. Действительно, рассмотрим три произвольных, взаимно ортогональных вектора ei, в2, &г> направленных вверх по отношению оси Oz. Построим три луча в направлениях ei,i = 1, 2, 3. Плоскость пересекающая эти три луча пересекает их в вершинах всегда остроугольного треугольника. Таким образом, не существует пространственного треугольника с углами Pi, г = 0,1, 2 с хотя бы одним углом не меньше чем ж/2 и проектирующегося в треугольник на плоскости z = 1 так, что ai = ж/2, i = 0,1, 2.

Замечание 3. Приведем пример не единственности решения поставленной задачи. Рассмотрим точки р0 = (0, 0, (), р1 = (к, 0, (), р2 = (0, Н, (), лежащие на плоскости г = (1> 1. Соответствующие проекции этих точек на плоскость г = 1 будут до = (0, 0,1) д1 = (^, 0,1), 42 = (0, , 1). Будем предполагать, что угол Ар1Ор2 < к/4. В треугольнике АроОр2 проведем высоту роЬ и рассмотрим точку р'2, лежащую на отрезке ОЬ и такую, что |рор2| = |рор21-Очевидно, что треугольники Ар1рор2 и Ар^ор^ равны так как оба прямоугольные и катеты у них равны. Таким образом, мы получили два треугольника с равными углами, лежащими в не параллельных плоскостях и центральные проекции которых равны одному и тому же треугольнику А д 1дод2-

3. Алгоритм поиска приближенных решений задачи

Для приближенного поиска всех решений системы (3) мы будем использовать двумерный аналог метода деления отрезка. Поэтому, в первую очередь, нужно определить прямоугольник на плоскости (и, и), в пределах которого находятся все решения этой системы.

Ясно, что и > 0, V > 0. Следовательно, достаточно получить верхнюю оценку границы и

венства

| рор1| > || ро| - |р 1||, | р ор21 > || ро| - | р21|, | р 1'р21 < | р 11 + | р21,

получаем

||ро| - < х Ы<ц

|р11 + |р21 ' |р11 + |р21 и

Л—^ < А, И—

и + V и + У

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

( и, )

1 — А + ц 1 + А — /1 — А — ц 1 + А + ц

1 + А + ц' 1 + А + ц/' \ 1 + А — ц' 1 + А — ц 1 + А + ц 1 — А — /1 + А — ц 1 — А + ц

1 — А + ц' 1 — А + ц/' \1 — А — ц' 1 — А — ц)

Отметим, что из неравенства треугольника и равенств (1) и (2) имеем

х . _ |рор11 + [ротЫ ^ 1

А + ц = -;-;- > 1,

| 1 2|

о, 1 , 2

и > 0 > 0

прямоугольнике [0, а] х [0, Ь], где

а = тах

Ь = тах

1 — А + ц 1 + А + ц 1 + А + ц, 1 — А + ц/'

1+А—ц 1+А+ ц

1+А+ ц' 1+А—ц

Теперь перейдем к описанию алгоритма. Обозначим через Гд и Г^ кривые, задаваемые уравнениями системы (3). Будем предполагать, что Гд и Г^ не являются эллипсами.

1. Разделим прямоугольник [0, а] х [0, Ь] на четыре прямоугольника прямыми и = | и V =

2. Выберем из получившихся прямоугольников те, которые пересекают обе кривые Гд и

Г-

3. Для каждого из найденных прямоугольников повторяем те же действия: разбиваем на четыре прямоугольника прямыми, проходящими через середины сторон и вбираем из получившихся прямоугольников те, которые пересекают обе кривые Гд и Г^.

4. Далее шаг 3 повторяется до тех пор, пока размеры получающихся прямоугольников не будут отвечать требуемой погрешности вычисления корней системы.

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

Теперь поясним некоторые моменты алгоритма. Во-первых, отбрасывая прямоугольники, которые не пересекаются с обеими кривыми Гд и Г^, мы на каждом шаге будем иметь не более 16 прямоугольников. Поэтому количество рассматриваемых на каждом шаге прямоугольников не будет неограниченно возрастать. Во-вторых, предположение, что кривые Гд и Г^ не эллипсы дает возможность свести проверку пересечения этих кривых с прямоугольником к проверке их пересечения со сторонами прямоугольника. Данная проверка осуществляется следующим образом.

Пусть, для определенности, нужно проверить пересекает ли кривая Гд отрезок (как одну из сторон некоторого прямоугольника) и = Uo, V\ <v < Это равносильно тому, что нужно

f(v) = (1 - \2)и0 - \2v2 + 2X2u0v cos а1 - 2и0 cos а0 + 1

хотя бы один корень на отрезке [v i, г^]. Обозначим через v* = Uo cos а\ точку, соответствующую вершине параболы. Тогда f(v) имеет хотя бы один корень на отрезке [vi, V2] только в том случае, если выполняется одно из условий

1. f(v 1) f(v2) < 0;

2. v* £ [vi, v2] и либо f(vi) f(v*) < 0 либо f(v*) f(v2) < 0.

Аналогично осуществляется проверка для второй кривой и других сторон прямоугольника. Отметим, что описанный алгоритм за п шагов дает погрешность вычисления корней системы

/0\ maxi а,Ь\

(3) не превосходящую величины —2п ■

Замечание 4. В силу того, что измерения углов треугольника qoQiQ2 для последующего вычисления углов ао,а1,а2 несут в себе определенную погрешность, возникает вопрос об устойчивости найденного приближенного решения относительно этих углов. Отметим, что углы ßo, ßi и ß2 не требуют измерений, так как они считаются известными. Для исследования этого вопроса мы предлагаем следующий подход.

Пусть (и*, V*) - точное решение системы (3), а (ип, vn) - приближенное решение, сходяще-

( и* , * ) п

F(u,v, cos а0, cos а1, cosа2) = 0, G(u,v, cos а0, cos а1, cosа2) = 0.

Для проверки устойчивости можно воспользоваться теоремой об обратном отображении. Вычислим якобиан этой системы

т/ \ 0(Г, G) . 2 2,\2 2,

J(и, v) = —;-- = 4(ц cos а\и + Л cos aiv +

о(и, v)

+(1 — Л2 — ц2)иу — (Л2 cos ai cos a2 + (1 — ц2) cos a0)v—

—((1 — Л2) cos a2 + cos aoai/j2)u + cos a0 cos a2).

Если J(un, vn) не сходится к 0 при п ^ œ, то точка (и*, v*) будет непрерывно зависеть от a0,ai и a2. Рассмотрим, например, частный случай, когда a0 = ai = a2 = к/2. Тогда из

и* = 0 * = 0

J (и*, и*) = (1 — Л2 — ц2 )u*v* = 0

только если Л2 + ц =1, что соответствует прямоугольному треугольнику. Таким образом, может наблюдаться неустойчивость решения, если сам треугольник p0'pip2 прямоугольный и углы a0 = ai = a2 = к/2.

СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ

1. Михайлов А. П., Чибуничев А. Г. Курс лекций по фотограмметрии. МИИГАиК Ракурс, 2013.

2. Лобанов А. Н., Куприна Н. Т., Чумаченко 3. Н. Фотограмметрия. М.: «Недра», 1984. - 552 с.

3. Клячин В.А., Григорьева Е.Г. Алгоритм автоматического определения параметров ориентации камеры в пространстве на основе характерных элементов фотоснимка// Тенденции развития науки и образования. 2018. № 45-6. С. 10-20.

4. Grigorieva Е. G., Klvachin V. A. Mathematical model of 3D maps and design of information system for its control//Journal of Computational and Engineering Mathematics, 2016., № 4 p. 51-58.

5. Локтев Д.A., Быков Ю.А., Коваленко H. Использование метода анализа размытия изображения для определения внешних дефектов железнодорожного пути // Наука и техника транспорта. 2016. № 1. С. 69-75.

6. Локтев А.А., Локтев Д.А. Метод определения расстояния до объекта путем анализа размытия его изображения // Вестник МГСУ. 2015. № 6. С. 140-151.

7. Локтев А.А., Бахтин В.Ф., Черников И.Ю., Локтев Д.А. Методика определения внешних дефектов сооружений путем анализа серии его изображений в системе мониторинга // Вестник МГСУ. 2015. № 3. С. 7-16.

8. Локтев Д.А., Локтев А.А. Определение расстояния до объекта по серии его изображений //В сборнике: Математика, информатика, естествознание в экономике и обществе (МИ-ЕСЭКО - 2015) труды Всероссийской научной конференции: в 2-х томах. 2015. С. 52-57.

9. Локтев Д.А. Определение геометрических параметров объекта с помощью анализа серии его изображений // T-Comm: Телекоммуникации и транспорт. 2015. Т. 9. № 5. С. 47-53.

10. Локтев А.А., Локтев Д.А. Оценка измерений расстояния до объекта при исследовании его графического образа // Вестник МГСУ. 2015. № 10. С. 54-65.

11. Локтев Д.А. Определение параметров объекта по серии его изображений в комплексной системе мониторинга //Путь и путевое хозяйство. 2015. № 2. С. 31-33.

12. Недзьведь О.В., Абламейко С.В., Белоцерковский A.M. Определение объемных характеристик динамических медицинских объектов // «Искусственный интеллект» 2010. № 4. С. 262-270.

13. Миров С., Иванов А., Огурцова Т., Дюкенджиев Е. Применение данных дистанционного зондирования в подометрии // 4-я Международная конференция пользователей ЦФС PHOTOMOD : Сборник тезисов докладов. - 2004. - С. 25-28.

14. Jackson A. S. et al. Large pose 3D face reconstruction from a single image via direct volumetric CNN regression //2017 IEEE International Conference on Computer Vision (ICCV). - IEEE, 2017. - p. 1031-1039.

15. Ferkova, Z., Urbanova, P., Cernv, D., Zuzi, M., k, Matula, P. (2018). Age and gender-based human face reconstruction from single frontal image. // Multimedia Tools and Applications, 1-26.

16. Pang, G., Qiu, R., Huang, J., You, S., Neumann, U. (2015, May). Automatic 3D industrial point cloud modeling and recognition. In 2015 14th IAPR international conference on machine vision applications (MVA) (pp. 22-25). IEEE.

17. Kamvab, S., Ghodsi, A., Azimifar, Z. Deep Structure for end-to-end inverse rendering. // ArXiv, abs/1708.08998 (2017).

18. Kamvab S. and Zohreh A. End-to-end 3D shape inverse rendering of different classes of objects from a single input image. // ArXiv, abs/1711.05858 (2017).

19. Левина О.С., Воронин В.В., Письменскова М.М., Гапон Н.В., Куркина А.В. Анализ основных этапов метода реконструкции трехмерных моделей поверхностей объектов // Информационные технологии. Радиоэлектроника. Телекоммуникации. 2016. № 6-2. С. 25-32.

20. Тагиров Т.С. Алгоритмические методы решения задач реконструкции объектов в 2D и 3D областях // «Современные проблемы науки и образования» № 2, 2014.

21. Penczek P. A. Fundamentals of three-dimensional reconstruction from projections. // Methods Enzvmol. 2010;482:1?33. doi:10.1016/S0076-6879(10)82001-4

22. Molnar, J., Frohlich, R., Chetverikov, D., Kato, Z. 3D reconstruction of planar patches seen by omnidirectional cameras. //In 2014 International Conference on Digital Image Computing: Techniques and Applications (DICTA) . IEEE. 2014. p. 1-8.

REFERENCES

1. Mikhailov A. P. & Chibunichev A. G. 2013. Kurs lekcij po fotogrammetrii [Course of lectures on photogrammetry], Moscow State University of Geodesy and Cartography, Racurs, Moscow.

2. Lobanov A. N., Kuprina N. T. & Chumachenko 3. N. 1984. Fotogrammetriya /Photogrammetryj, «Nedra», Moscow, pp. 552.

3. Klvachin V.A. k, Grigorieva E.G. 2018. "Algorithm for automatic determination of camera orientation parameters in space based on characteristic elements of a photograph", Tendencii razvitiya nauki i obrazovaniya., no. 45-6, pp. 10 - 20.

4. Grigorieva E. G.& Klvachin V. A. 2016. "Mathematical model of 3D maps and design of information system for its control", Journal of Computational and Engineering Mathematics, no. 4, pp. 51 - 58.

5. Loktev D.A., Bvkov Yu.A. k Kovalenko N. 2016. "Using the method of image blur analysis to determine external defects of the railway", Nauka i texnika transporta, no. 1, pp. 69 - 75

6. Loktev A.A. k Loktev D.A. 2015. "Method of determining the distance to the object by analyzing the blurring of its image", Vestnik MGSU, no. 6, pp. 140 - 151.

7. Loktev A.A., Bakhtin V.F., Chernikov I.Y. k Loktev D.A. 2015. "Methods for determining the external defects of structures by analyzing a series of its images in the monitoring system", Vestnik MGSU, no. 3, pp. 7 - 16.

8. Loktev D.A. k Loktev A.A. 2015. "Determining the distance to an object from a series of its images", V shornike: Matem,atika, inform,atika, estestvoznanie v ekonomike i, obshhestve trudy Vserossijskoj nauchnoj konferencii: v 2-x torn,ax ( In the collection: Mathematics, Informatics, Science in Economics and Society, works All-Russian scientific conference: in 2 volumes). Moscow. 2015. pp. 52 - 57)

9. Loktev D.A. 2015. "Determination of geometric parameters of the object by analyzing a series of its images", T-Comm: Telekommunikacii i, transport, vol. 9, no. 5, pp. 47 - 53

10. Loktev A.A. k Loktev D.A. 2015. "Estimation of measurements of distance to object at research of its graphic image", Vestnik MGSU, no. 10, pp. 54 - 65.

11. Loktev D.A. 2015. "Determination of object parameters by a series of its images in a complex monitoring system", Put' i, putevoe xozyajstvo, no. 2, pp. 31 - 33.

12. Nedzved O.V., Ablameiko S.V. k Belotserkovskv A.M. 2010. "Determination of volumetric characteristics dynamic medical objects", Iskusstvenny'j intellekt, no. 4, pp. 262 - 270.

13. Mirov S., Ivanov A., Ogurtsova T. k Dvukendzhiev E. 2004. "Application of remote sensing data in podometrv", 4-ya Mezhdunarodnaya konferenciya pol'zovatelej CzFS PHOTOMOD : Sbornik tezisov dokladov (4th International Conference of users of CFS PHOTOMOD: Collection of abstracts.) pp. 25 - 28.

14. Jackson A. S. et al. 2017. "Large pose 3D face reconstruction from a single image via direct volumetric CNN regression" IEEE International Conference on Computer Vision (ICCV). IEEE, pp. 1031 - 1039.

15. Ferkova, Z., Urbanova, P., Cernv, D., Zuzi, M., k Matula, P. 2018. "Age and gender-based human face reconstruction from single frontal image", Multimedia Tools and Applications, pp. 1 - 26.

16. Pang, G., Qiu, R., Huang, J., You, S., k Neumann, U. 2015, "Automatic 3D industrial point cloud modeling and recognition", In 2015 14-th IAPR international conference on machine vision applications (MVA). IEEE. pp. 22 - 25.

17. Kamvab S., Ghodsi A. k Azimifar, Z. 2017, "Deep Structure for end-to-end inverse rendering", ArXiv, abs/1708.08998.

18. Kamvab S. and Zohreh A. 2017, "End-to-end 3D shape inverse rendering of different classes of objects from a single input image", ArXiv, abs/1711.05858.

128

A. A. K. imhiiii. B. A. K. imhiiii

19. Levina O.S., Voronin V.V., Pis'menkova M.M., Gapon N.V.,& Kurkina A.V. 2016, "Analysis of the main stages of the method of reconstruction of three-dimensional models of object surfaces", Informacionnye tehnologii. Radiojeiektronika. Telekommunikacii, no. 62. pp. 25 - 32.

20. Tagirov T.S. 2014, "Algorithmic methods for solving problems of reconstruction of objects in 2D and 3D regions", Sovremennye problemy nauki i obrazovanija, no. 2.

21. Penczek P. A. 2010, "Fundamentals of three-dimensional reconstruction from projections", Methods Enzymol., no. 482, pp. 1 - 33. doi:10.1016/S0076-6879(10)82001-4

22. Molnar, J., Frohlich, R., Chetverikov, D. k, Kato, Z. 2014, "3D reconstruction of planar patches seen by omnidirectional cameras", In 2014 International Conference on Digital Image Commuting: Techniques and Applications (DICTA). IEEE. pp. 1-8.

nojivHeno 11.06.2020 r. IlpHHaTO B nenaib 22.10.2020 r.

i Надоели баннеры? Вы всегда можете отключить рекламу.