Восстановление эллипса по распределениям хорд по направлениям Текст научной статьи по специальности «Естественные и точные науки»

Б01 10.24412/1829-0450-£т-2024-1-17-24 УДК 519.218.5

Поступила: 20.02.2024г. Сдана на рецензию: 26.02.2024г. Подписана к печати: 19.05.2024г.

ВОССТАНОВЛЕНИЕ ЭЛЛИПСА ПО РАСПРЕДЕЛЕНИЯМ

ХОРД ПО НАПРАВЛЕНИЯМ

Р.Г. Арамян, Л.А. Апинян

Российско-Армянский (Славянский) университет apmlev00@gmailcom, [email protected]

АННОТАЦИЯ

Задача нахождения информации или полного восстановления выпуклых тел по распределениям характеристик к-мерных сечений -одна из основных задач стохастической томографии. Данная работа посвящена восстановлению эллипса на плоскости через распределения случайных хорд по направлениям. Предлагается алгоритм восстановления.

Ключевые слова: эллипс, выпуклость, функция распределения, опорная функция, восстановление эллипса.

1. Введение

Геометрические характеристики выпуклого тела К, инвариантные относительно Евклидовых движений, могут быть вычислены статистическими методами. Решение таких задач выходит далеко за рамки теоретического интереса.

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

Восстановление выпуклых тел (компактных множеств) является главной задачей геометрической томографии. Геометрическая томография - это математическая область, которая фокусируется на проблемах восстановления однородных (часто выпуклых) объектов из томографических данных (это могут быть рентгеновские лучи, проекции, сечения, функции яркости или кова-риограммы). По словам Ричарда Гарднера, «Геометрическая томография посвящена поиску информации о геометрическом объекте из данных о его проекциях (тенях) на плоскостях или сечениях плоскостями». Геометрическая томография в рамках математики в основном связана с выпуклой геометрией, интегральной геометрией, геометрией Минковского, функциональным ана-

лизом и комбинаторикой, а также вне математики, с компьютерной томографией, медицинской томографией, дискретной томографией, электронной микроскопией, компьютерным зрением и распознаванием образов.

Стохастическая томография посвящена поиску информации о геометрическом объекте из случайных данных о его проекциях (тенях) на плоскостях или распределения случайных сечений к-мерными плоскостями.

Постановка задачи

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

Проблему восстановления выпуклых тел по распределению длин случайных хорд впервые рассмотрел Вильгельм Бляшке [1]. Пусть задана выпуклая фигура. Случайная прямая (с инвариантной мерой), пересекающая фигуру, задает случайную хорду х(^). В. Бляшке поставил задачу: определяется ли выпуклое тело по распределению длины случайных хорд с точностью до параллельных переносов и отражений? Позже Колин Л. Маллоус и Джон Морис Кларк опровергали вышесказанное [2], построив два неконгруэнтных 12-угольника с одинаковыми распределениями длин хорд.

После была поставлена вторая гипотеза: восстанавливается ли выпуклая фигура по распределению длин случайных хорд по направлению с точностью до параллельных переносов и отражений? Также эту проблему рассматривал Жорж Матерон. В 1975г. он ввел понятие ковариограммы, показал связь между ней и функцией распределения длины случайной хорды по направлению [6], а в 1986г. предложил гипотезу [4]: ковариограмма однозначно определяет выпуклое тело в Мп, с точностью до параллельных переносов и отражений. В 1993г. Вернер Нагель доказал эту гипотезу для всех выпуклых многоугольников [3]. Матерон также предложил положительный ответ при п = 2, который был доказан в 2007 г. Аверковым и Г. Бьянки [4]. Габриэль Бьянки построил контрпримеры в Мп для каждого п> 4 и показал, что задача ковариограммы в общей постановке имеет отрицательный ответ [5]. Для п = 3 проблема пока не решена. Но несмотря на то, что задача имеет положительный ответ при п = 2, алгоритм восстановления выпуклых фигур с помощью распределения длин хорд по направлению до сих пор остается открытым.

Цель данной статьи - найти какие-то характеристики эллипса от его распределений хорд и выяснить способ восстановления эллипса через его распределения длин хорд по направлениям.

2. Предварительные понятия

2.1. Координаты прямой на плоскости

Пространство прямых на плоскости R2 будем обозначать через G. Прямую g Е G параметризуем парой координат (^ ф), p > 0, 0 < ф < 2п, то есть полярными координатами основания перпендикуляра опущенного на прямую g из начала координат O. Введем в G инвариантную относительно евклидовых движений, меру dg = dpdф.

2.2. Распределение длины хорды по направлению

Пусть D - компактное множество в Евклидовом пространстве М2. С пространство прямых в М2, д Е С, (р, —полярные координаты прямой д, где р > 0, <р Е 51, 51 — единичная окружность на плоскости с центром в начале координат. Прямая д пересекает множество И и образует хорду, которую обозначим х(в) (Рис. 1).

( * ; Г

Рисунок 1.

Зафиксируем направление Рассмотрим все прямые в этом направлении. Каждая такая прямая однозначно определяется ее пересечением с осью 0x.

Функция распределения длины хорды по направлению - это

^(0 = Р{1х(д<р)1 < ь: х(д<р) = 9<р^о} ,

где д<р имеет направление (р. 3. Основные результаты

Здесь рассмотрен эллипс на R2 и выявлена его функция распределения длин хорд по направлению, описана опорная функция и представлен способ восстановления.

3.1. Опорная функция эллипса

Опорная функция тела K - это 5ир(<р, х) = к(ф), где (р есть единичный

* Е К

вектор в направлении Известно, что опорная функция выпуклого тела K однозначно задает К. Пусть K эллипс с полуосями a и Ь:

2 2 а2 + Ь2 1

Известно, что опорная функция эллипса имеет следующий вид:

к(ф) = ^а2соз2<р + Ь2зт2ф

3.2. Функция распределения длин хорд для эллипса по направлению

Для эллипса К обозначим через 1 длину хорды (А, В) отсекаемой произвольной прямой g брошенной по направлению ф, обозначим через р длину перпендикуляра, проведенного от центра эллипса до хорды 1. Найдем функцию распределения длины хорды К по направлению ф, то есть ^ (Ь).

Из определения функции распределения имеем:

РФ(1) =

КФ)-Р Чф)

Чтобы найти координаты А = (х1,у1) и В = (х2,у2), рассмотрим следующую систему:

( у2 _

{ а2 + Ь2 1

[хсояф + узтф = р

Заметим, что

I = \АВ\ = ^(Х2-Х1)2+ (У2-У1)2 Подставляя выражение для х из второго уравнения системы в первое, по-

лучим

р- уБтф \

2

СОБф ) У_ _ ^

Ь2

Откуда получим

, ^2 _ 4р2Ь4зт2ф 4Ь2р2-4Ь2а2соз2ф

(У2 - У1) = ^(ф) Ь2(ф)

Аналогичными образом получим

а

(Х2 - Х1)

2 _ 4р2а4соз2р 4а2р2-4Ъ2а2зт2р

Л4(р) л^СФ)

Так как

,.2 _ („ л2 _ 4Р2 (а4с°з2р-1}4зт2р , иА

I -(Х2-Х!) +(У2-У1) -ЛЙРР)--а -Й)+ЛЧР)

Окончательно получается

■2 4а2Ъ2\ _ 4р2а2Ъ2 - К2(р)) - - Л4(р)

Отсюда и получается выражение для p

р — к(ф)7 ' 1 1

Л(ф)2 4 а2Ъ2

Подставив его в Fр ( £) и сделав соответствующие упрощения, был получен окончательный вид для функции распределения длин хорд по направлению:

ко —1-./1-^

Отсюда можно выразить Н(р) через Рр(£)

К<р)— ^ ^р(О)2

Получается следующая лемма:

Лемма. Для любого эллипса и для любого направления р выражение

^ ./1-(1-*р(0)2

не зависит от t.

Обратное утверждение мы представляем как гипотезу:

Гипотеза. Если для некой фигуры и любого направления р выражение

^ — (1- ^р(О)2 не зависит от t, то эта фигура является эллипсом.

3.3. Восстановление неповернутого эллипса

Мы предлагаем следующий алгоритм для восстановления неповернутого эллипса

Для восстановления необходимы 2 самые длинные хорды в двух не кол-линеарных направлениях. Самую длинную хорду в направлении ф обозначим Найти ее можно из данной системы:

— + = 1

а2+ ь2 1

xcos<p + ysm<p = 0 1<р = V(х2 - Х\)2 + (У2 - У1)2

аЬсозф

У\,2 = ±

2,1 ± Vа2соз2ф+ Ь2зт2ф

2=

аЬзтф

Откуда получаем, что

2аЬ

Ф ^а2соз2ф+Ь2зШ2ф

Остаётся решить следующую систему выразив а и Ь через эти хорды и их направления

^ _ 2аЬ

Va2cos2Фl+ Ь2зт2ф1 ^ _ 2аЬ

Ф2 V а2 саз2 Ф2+ Ъ2з1п2ф2

Окончательное выражение для полуосей выглядят так:

_ 1^11^2созф1 \ зт2ф2- 1д2ф1

сиь2ф1- 1ф2

2 \11юл2соз2ф1- и22соз2ф2

£ _ ¡^¡^зтф! соз2ф2-ад2ф!

2 ^ (^12^1П2ф1- 1р22з(п2ф2

3.4. Восстановление повернутого эллипса

Мы предлагаем следующий алгоритм для восстановления эллипса повернутого на угол Аналогичным образом составляется следующая система:

Ipi

=

L =

2 ab

*Ja2cos2(p1- ф) + b2sin2(pp1- Ф)

2 ab

Ja2cos2(p2- ф)+ b2sin2(pp2- Ф)

2 ab

< ^а2соз2(ф3-1р)+ Ь2зт2(ф3- ф)

Решив систему аналогичным образом, мы получим a, Ь и ^ в терминах Фъ (2, (3, , , ^з. Получаются следующие выражения для параметров эллипса:

ф =

а =

Ъ =

CI+C3-I(CI-C3)2+C22

Ci+C3+I(CI-C3)2+C22

где

Cl = Тй2'

С3 = Т2' 1 0

8

С2 = С1 + С3

4

Заключение

Таким образом мы показали, способ восстановления произвольного эллипса (повернутого) по его распределению случайных хорд по трем направлениям. Вопрос восстановления произвольной выпуклой фигуры по распределению случайных хорд по направлениям открыт

ЛИТЕРАТУРА

1. Luis A. Santalo. Integral Geometry and Geometric Probability, Addision-Wesley, Reading, MA (2004).

2. Mallows C. L., Clark J. M.C. Linear-intercept distributions do not characterize plane sets. J. Appl. Probab. 7, 240-244 (1970).

3. Nagel W. Orientation-dependent chord length distributions characterize convex polygons. J. Appl. Probab. 30, 730-736 (1993).

2

2

2

2

24

P.r. ApaMMH, flA. Anmm

4. Averkov G., Bianchi G. Confirmation of Matheron's conjecture on the covariogram of a planar convex body, J. Eur. Math. Soc., in press, arXiv:0711.0572 [math.MG].

5. Bianchi G. The covariogram determines three-dimensional convex polytopes. Advances in Mathematics 220 (2009), 1771-1808; arXiv:0805.1605v1 [math.MG].

6. Matheron G. Random Sets and Integral Geometry. Willey, New York (1975).

RECONSTRUCTION OF AN ELLIPSE FROM THE DISTRIBUTIONS OF CHORDS IN DIRECTIONS

R. Aramyan, L. Apinyan

Russian-Armenian(Siavonic) University

ABSTRACT

The problem of finding information or completely reconstructing convex bodies from the distributions of characteristics of k-dimensional sections is one of the main problems of stochastic tomography. This work is devoted to the reconstruction of an ellipse on a plane through the distribution of random chords in directions. A reconstruction algorithm is proposed. Keywords: ellipse, convexity, distribution function, support function, ellipse restoration.