Journal of Siberian Federal University. Engineering & Technologies 2020 13(3): 328-337
DOI: 10.17516/1999-494X-0224 УДК 519.816
A Mathematical Model of the Detection Probability of the Point Target by the Operator of the Payload Electro-Optical System for Unmanned Aerial Vehicle
Aleksey I. Tishchenko and Stepan V. Artyshchenko*
Military Education and Research Centre of Military-Air Forces
«Military-Air Academy Named After Professor N.E. Zhukovsky and Yu.A. Gagarin»
Voronezh, Russian Federation
Received 27.09.2019, received in revised form 04.12.2019, accepted 21.01.2020
Abstract. The paper deals with a mathematical model of the probability of detection of a point target by the operator of the payload of an optoelectronic device installed on an unmanned aerial vehicle, based on the use of the Johnson criterion and the binomial law of repeated tests.
Keywords: a mathematical model, a probability of detection, optoelectronic systems, unmanned aerial systems, a payload operator, rescue operations.
Citation: Tishchenko A.I., Artyshchenko S.V. A mathematical model of the detection probability of the point target by the operator of the payload electro-optical system for unmanned aerial vehicle, J. Sib. Fed. Univ. Eng. & Technol., 2020, 13(3), 328-337. DOI: 10.17516/1999-494X-0224
© Siberian Federal University. All rights reserved
This work is licensed under a Creative Commons Attribution-Non Commercial 4.0 International License (CC BY-NC 4.0). Corresponding author E-mail address: aleksei.tishenko@yandex.ru
*
Математическая модель вероятности обнаружения точечной цели оператором полезной нагрузки оптико-электронной системы беспилотного летательного аппарата
А.И. Тищенко, С.В. Артыщенко
Военный учебно-научный центр Военно-воздушных сил
«Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» Российская Федерация, Воронеж
Аннотация. В работе рассматривается математическая модель вероятности обнаружения точечной цели оператором полезной нагрузки оптико-электронного средства, установленного на беспилотном летательном аппарате, на основе использования критерия Джонсона и биноминального закона повторных испытаний.
Ключевые слова: математическая модель, вероятность обнаружения, оптико-электронные системы, беспилотные летательные системы, оператор полезной нагрузки, аварийно-спасательные работы.
Цитирование: Тищенко, А.И. Математическая модель вероятности обнаружения точечной цели оператором полезной нагрузки оптико-электронной системы беспилотного летательного аппарата / А.И. Тищенко, С.В. Артыщенко // Журн. Сиб. федер. ун-та. Техника и технологии, 2020. 13(3). С. 328-337. DOI: 10.17516/1999-494Х-0224
Введение
Подразделения поисково-спасательной службы используют беспилотные летательные аппараты (БпЛА) при проведении аварийно-спасательных работ в тех случаях, когда необходимо оперативно, достоверно без риска для жизнедеятельности личного состава спасательного подразделения, экипажей самолетов и вертолетов оценить ущерб от чрезвычайных ситуаций природного или техногенного характера, а также при поиске и спасении пропавших или терпящих бедствие людей в труднодоступных районах.
В зависимости от решаемой задачи на БпЛА устанавливается соответствующая целевая нагрузка. Для поиска людей используют видео и инфракрасные (ИК) камеры, которые в комплексе с оборудованием беспилотного летательного аппарата можно рассматривать как оптико-электронные системы (ОЭС).
Для планирования мероприятий, связанных с применением БпЛА в ходе выполнения аварийно-спасательных работ, штабу поисково-спасательных работ целесообразно иметь расчетные программы, позволяющие с высокой степенью достоверности, исходя из наличия и возможностей ОЭС, парка беспилотных летательных аппаратов прогнозировать сроки выполнения поставленной задачи.
Реализация расчетных программ возможна при наличии математических моделей, соответствующих задачам, стоящим перед расчетом БпЛА в ходе проведения аварийно-спасательных работ.
Целью данной работы является разработка математической модели расчета вероятностей обнаружения и распознавания наземной цели оператором полезной нагрузки оптико-электронной системы, установленной на БпЛА, осуществляющем в заданном районе случайный поиск. Обнаружение и распознавание оператором осуществляется по монитору наземного пункта управления БпЛА.
Для решения поставленной цели определимся с некоторыми понятиями, которые будем использовать в дальнейшем.
Под наземным объектом поиска будем понимать простые, малоразмерные, стационарные или малоподвижные объекты, находящиеся на поверхности Земли. Площадь этих объектов, как правило, в сотни тысяч раз меньше площади районов поиска. Поэтому в дальнейшем объекты поиска будем называть точечными.
Для подразделения поисково-спасательных работ, осуществляющего поиск людей (объектов поиска) с помощью ОЭС БпЛА, основными визуальными задачами являются обнаружение и распознавание объекта поиска. Под обнаружением будем понимать уровень восприятия, при котором осуществляется выделение на экране монитора пятна, соответствующего объекту поиска, на фоне помех. Под распознаванием будем понимать уровень восприятия, при котором объект поиска различается уровнем детализации, позволяющим определить его класс - человек, автомобиль и т.д.
Тактическими характеристиками задачи являются высота, маршрут полета БпЛА в районе поиска, а также время обнаружения и распознавания объекта поиска.
На решение задачи по обнаружению и распознаванию объектов поиска могут оказать влияние следующие факторы: качество изображения объекта поиска на экране ОЭС; объем априорной информации, на основе которой осуществляется выделение объекта поиска на фоне окружающей среды; степень подготовки оператора полезной нагрузки. Первые три фактора определяются параметрами ОЭС и БпЛА, привлекаемых к поиску. К числу технических характеристик ОЭС, оказывающих влияние на вероятности обнаружения и распознавания наземных объектов, на выбор маршрута и профиль полета БпЛА, относятся: фокусное расстояние, углы поля зрения по горизонтали и вертикали и разрешающая способность объектива; размер пикселя и матрицы приемника излучения.
Математическая модель вероятности обнаружения цели оператором полезной нагрузки оптико-электронной системы беспилотного летательного аппарата
Решение поставленной цели выполним применительно к следующим условиям. Расчет беспилотного летательного аппарата осуществляет поиск точечного наземного объекта площадью Sт в заданном районе поиска площадью Sр с помощью ОЭС видимого диапазона, установленной на БпЛА. ОЭС включает в себя видеокамеру и монитор. Монитор находится на наземном пункте управления беспилотного летательного аппарата. Формат изображения монитора полностью соответствует формату изображения приемника излучения. Видеокамера установлена непосредственно на летательном аппарате под углом визирования а. Технические характеристики камеры с установленным объективом: вертикальный угол поля зрения ав; горизонтальный угол поля зрения аг; фокусное расстояние объектива /; скорость съемки т (для
фотосъемки серийность съемки, для видеосъемки частота кадра); линейное разрешение ОЭС на местности Л; размер пикселя матрицы д; размер матрицы приемника излучения по горизонтали а и по вертикали Ь.
Проекция поля зрения видеокамеры на контролируемую в процессе полета БпЛА земную поверхность при угле визирования 0° имеет форму прямоугольника, а при угле визирования, отличном от 0°, - трапецию (рис. 1) с параметрами и площадью, определяемыми следующими выражениями:
1, = АВ ■■
2KJg\
а.
L = DC =
cos\ а + -
-<90°;
,а--- < -90°;
2
2 cos а,
а II а, cos\ а +--I cos\ а--
2 J I 2
,а ±—а(+90°; 2 4
14 = PG =
[2_А.
(1)
(2)
(3)
(4)
So63 = l4l3,
(5)
где 11 - основание трапеции (сторона прямоугольника), ближнее относительно БпЛА; 12 - основание трапеции (сторона прямо}'гольника), дальнее относительно БпЛА; 13 - высота трапеции (прямоугольника); 14 - средняя линяя трапеции; Sобз - площадь зоны обзора ОЭС земной поверхности; кпол - высота полета беспилотного летательного аппарата.
1
l3 = МК = ипол
+
2
Рис. 1. Проекция поля зрения видеокамеры на земную поверхность Fig. 1. Projection of the field ofview of the camera on the earth's surface
- 331 -
Связь между размерам и матрицы приемника излучения и углами поля зрения ОЭС такова:
аг = аГС^ 27; (6)
ь
а = аг^ —. (7)
* 2 7
Решение задачи рассмотрим применительно к случаю, когда количество БпЛА, привлекаемых к поиску, ограничено, район поиска обширен, данные о месте и маршруте движения терпящих бедствие отсутствует. В подобных условиях для детального осмотра участков местности в пределах района поиска используют прямолинейные взаимно параллельные маршруты - параллельное галсирование. Маршруты прокладывают таким образом, чтобы края поля зрения ОЭС БпЛА двух соседних маршрутов при установленной высоте полета hпол летательного аппарата, как минимум, соприкасались. Это обеспечит 100%-ный просмотр ОЭС БпЛА района поиска. Следовательно, расстояние между маршрутами должно быть равным ¡¡. Облет района поиска летательным аппаратом осуществляется с крейсерской скоростью Укр.
В процессе полета БпЛА полачаемое с помощью видеокамеры видеоизображение участка местности передается по радиоканалу на наземный пункт управления, где отображается на мониторе оператора целевой нагрузки. Оператор в ходе анализа полученной информации обнаруживает цель.
В системах видеорегистрации для определения вероятностей обнаружения и распознавания объектов используют критерий Джонсона. Его суть применительно к решаемой задаче заключается в том, что существует соответствие между числом разрешаемых периодов эквивалентной миры, укладывающихся на минимальном размере наблюдаемого объекта поиска, и уровнем его восприятия. В качестве эквивалентной миры взята штриховая мира прямоугольной формы, ширина которой равна минимааьному размеру объекта поиска, а длина соответствует размеру, перпендикулярному минимальному размеру. Этот критерий т.кже применяется для оптико-электронных сист ем с выходным видеосигналом видеокамер с тем отличием, что вместо пар линий штриховой миры значения указывают в телевизионных линиях или, что то же самое, в элементах изображения (пикселях). Графическое выражение критерия Джонсона применительно к решаемой задаче обнаружения и распознавания изображений объектов поиска представлено на рис. 2 [1].
Используя метод наименьших квадратов и аппроксимируя кривые вероятностей обнаружения и распознавания объектов от количества элементов разрешения занимаемыми ими на видеоматериалах экспоненциальными функциями, получим аналитические зависимости:
Роб = 1 и ти0'2(-1)2, (8)
Ррас = 1 и т и°,™2, (9)
где Роб, Ррас - вероятности обнаружения и распознавания и идентификации объектов; N - количество элементов разрешения занимаемых объектом поиска на мониторе наземного пункта управления.
Ь 10 15 20 25
N
Рис. 2. Зависимости вероятностей о бнаружения (1) и распознавания (2) объектов поиска от количества элементов разрешения, занимаемых ими на мониторе
Fig. 2. Dependences of probabilities of detection (1) and recognition (2) of search objects the number of resolution elements they occupy on the montio
Зависимость между размерами объекта поиска на местности и количеством N определяется выражением [2]
т /
N = i
RK
(10)
где 1т - линейные размеры объекта поиска по направлению полета БпЛА; R - линейное разрешение ОЭС на местности, связанное с размерами д соотношением вида
К8
R =-
f
(11)
При попадании объекта поиска в поле зрения ОЭС он появится на экране монитора наземного пункта управления и оператор полезной нагрузки обнаружит и распознает объект поиска с вероятность р, равной вероятности распознавания Ррас.
Район поиска применительно к ОЭС БпЛА можно рассматривать как серию видеокадров п=п1п2 - независимых испытаний, определяемых произведением количества элементов разрешения п1, приходящимся на полный маршрут полета БпЛА, и количеством кадров п2, приходящихся на элемент разрешения, определяемых выражениями (12), (135). Элемент разрешения представляет собой поле зрения ОЭС на земную поверхность.
1' 2 2
3
V..
(12)
(13)
(14)
где 1пол - полная протяженность полета БпЛА в районе поиска; |~х] - наибольшее целое, большее или равное х.
Количество кадров т, на которых имеется изображение объекта поиска, определяется произведением количества элементов разрешения п3, приходящимся на размер объекта поиска, и
количеством кадров п2, приходящихся на элемент разрешения, определяемых выражениями (13), (14).
(15)
(16)
Если обозначить вероятность обнаружения и распознавания объекта поиска в i-м снимке черезр , равное р, то для расчета вероятности того, что в п снимках района поиска обнаружение и распознавание объекта поиска произойдут не менее чем в т снимках можно использовать биноминальный закон повторных испытаний [2]:
Рп (=) = С=рт (1 - р)п
(17)
где Рп(т) - вероятность обнаружения объекта поиска оператором полезной нагрузки т раз в серии из п видеокадров; С% - число сочетаний из п по т, равное
Ст -
n\
т\(п —т )|
(18)
Используя Рп(т) можно вычислить вероятность обнаружения и распознавания объекта поиска по видеоматериалам хотя бы 1 раз):
Run= i-I-^1 -п qj,
no 1=1
где q=1 - Pi, Pi=0 для (=1 - С p= Ррар для i= z+1 - z+m; p¿=0 для i= z+m - n. Отсюда
R — 1 — (P P • • • P )
1 ,n n pacz+l pacz+2 pacz+m /*
Так как РРаС+= РПа—+2= ■■■ = Ррас+m, то выражение (199) примет вид
R =1- рт
1 M,n * 1 рас •
(19)
(20)
Подставляя в (20) выражения (3), (6), (7), (9), (12) - (16), получим выражение, позволяющее оценить вероятность обнаружения и распознавания оператором полезной нагрузки цели с помощью установленных на БпЛА ОЭС в первом зале те.
( г 1
о= -1 — ехр
I f
— 0,015(——--1)
Rh
2hno. c°s(E)sm| Е- I
(21)
кр
sini a--isini «н—
I 2 I I 2
Время, необходимое для выполнения первого залета БпЛА в район поиска, определяется выражением
(
t = \L +1 + л±
l, 1 l,.
2 IIV
¿■J lV кр
2hnJg \ ^
Л
f
sin
hp +1+ Я-
a
\
a + -
l,
2KJg \ a2r
(22)
f
sin
a+
2
V
N кр
где 1д, 1ш - размеры района поиска по длине и ширине соответствен но.
Выражения (21) и (22) показывают, что вероятность обнаружения и распознавания объектов поиска в заданном районе и время, необходимое для выполнения этого поиска, можно рассчитать по известным характеристикам ОЭС, полета БпЛА, размерам района поиска и объекта поиска.
Эти зависимости позволяют решить обратную задачу, когда по заданному значению вероятности обнаружения и распознавания объекта поиска, известным характеристикам ОЭС, размерам района поиска и объекта поиска можно рассчитать параметры полета БпЛА и количество вылетов в район поиска, обеспечивающих минимальное время поиска.
С этой целью устанавливается пороговое значение Rl^. Из условия Rl,n*< Rl,n применительно к размерам района и объекта поиска определяют характеристики ОЭС, полета БпЛА, при которых время поиска объекта будет минимальным.
Если в результате залета условие Rl|n*< Rl|n не выполнимо, то в этом случае необходимо выполнить последующие залеты БпЛА в район поиска до тех пор, пока неравенство не будет выполнено.
Расчет Rl,k в последующих залетах осуществляется с использован ием выражения (23), а время / определяется как произведение времени выполнения одного залета на количества залетов:
*и = 1 -П (1 - Rin )
(23)
где k - количество залетов БпЛА в район поиска.
Численное исследование математической модели
Штат подразделения поисково-спасательных работ может содержать расчет БпЛА Phanton 3 Professional [3]. Поэтому численное исследование математической модели проведем применительно к данному беспилотному летательному аппарату, состоящему из непосредственно летательного аппарата, видеокамеры и пульта дистанционного управления.
Основные технические характеристики рассматриваемой модели:
беспилотный летательный аппарат: максимальная скорость полета 16 м/с (~60 км/ч); максимальное время полета ~23 мин; максимальная высота полета 6000 м; дальность передачи видеосигнала до 2000 м; подвес обеспечивает охват с углом наклона от -90 до 30°;
2
2
1=1
видеокамеры: оптика EXMOR 1/2,3"; эффективные пиксели 12,4 млн (всего пикселей 12,76 млн); объектив с углом обзора 94°; фокусное расстояние 20 мм; эквивалентный формат 35 мм f/2,8; максимальный размер изображения 4000x3000; режим видеосъемки FHD: 1920х1020р 24/25/30/48/50/60 fps [4].
В качестве объекта поиска выберем человека, площадь поперечного сечения которого описывается овалом, длина и ширина которого составляют соответственно 0,75 и 0,5 м. Район поиска представим в виде квадрата со сторонами 2000 м, исходя из дальности передачи видеосигнала. Скорость полета БпЛА 60 км/ч. Пороговое значение вероятности обнаружения и распознавания объекта поиска по видеоматериалам хотя бы 1 раз определим равным 0,99.
Результаты расчета времени, необходимого для обнаружения и распознавания расчетом БпЛА потерявшегося в районе поиска человека в зависимости от высоты полета летательного аппарата (50 - 1000 м), углов визирования видеокамеры 0, 10, 20, 30 градусов и технических характеристик, приведенных выше, представлены на рис. 3.
Н. м 1000 900 800 700 600 500 400 300 200 100 о
5 10 15 20 25 20 35
t. мин
Рис. 3. Расчетные значения времени, необходимого для обнаружения и распознавания расчетом БпЛА потерявшегося в районе поиска человека при различных углах визирования 1 - 0°, 2 - 10°, 3 - 20°, 4 - 30°
Fig. 3. Estimated time required for detection and calculation of the UAV lost in the search area of the person at different viewing angles 1 - 0°, 2 - 10°, 3 - 20°, 4 - 30°
Анализ полученных результатов показывает, что время поиска потерявшихся людей в целом зависит от высоты полета летательного аппарата в ходе поиска и угла визирования видеокамеры. Наиболее существенна эта зависимость в диапазоне высот от 50 до 200 м и от 600 до 1000 м. В диапазоне высот от 200 до 600 м зависимость практически несущественно зависит от высоты полета летательного аппарата и угла визирования видеокамеры. Разность результатов не превышает 3 - 5 мин. Следовательно, можно сделать вывод, что поиск целесообразно вести на указанных высотах при угле визирования 0 градусов, так как при данном угле погрешности определения координат объекта по данным воздушной разведки минимальны. Полученные результаты согласуются с практическими рекомендациями, изложенными в [3], и указывают на адекватность представленной математической модели.
Заключение
Таким образом, получена математическая модель вероятностей обнаружения и распознавания точечной цели оператором полезной нагрузки оптико-электронной системы беспилотного летательного аппарата. Она может быть реализована с применением программ высокого уровня. Результаты моделирования целесообразно использовать при планировании поисково-спасательных работ, связанных с поиском пропавших людей или объектов поиска в труднодоступных районах.
Список литературы / References
[1] Тищенко А.И., Ищук И.Н., Громов Ю.Ю. Определение высоты полета беспилотного летательного аппарата для выполнения видеонаблюдения с требуемой степенью детализации информации при поиске тепловых объектов, Информационно-сенсорные системы в теплофи-зических исследованиях: сб. науч. ст. (Тамбов, 6 - 9 ноября 2018 г.). Тамбов: Издательский центр ФГБОУ ВО «ТГТУ», 2018, 2, 190 - 194. [Tischenko A.I., Ischuk, I.N., Gromov Yu.Yu. Determination of the height of the flight of unmanned aircraft for video monitoring with the required degree of detail of information when searching for thermal objects, Information and sensor systems in thermophysical studies: collection of scientific works. article (Tambov, 6 - 9 November 2018). Tambov: Publishing house FGBOU VO «TSTU», 2018, 2, 190 - 194 (In Russian)]
[2] Корн Г., Корн Т. Справочник по математике для научных работников и инженеров, М.: Главная редакция физико-математической литературы, 1984. 832 с. [Korn G., Korn T. Handbook of mathematics for scientists and engineers, Moscow.: Home edition of physical and mathematical literature, 1984. 832 p. (In Russian)]
[3] Поисково-спасательные работы при помощи беспилотных летательных аппаратов [Электронный ресурс] - Режим доступа: https://fireman.club/statyi-polzovateley/poiskovo-spasatelnyih-rabotyi-pri-pomoshhi-bespilotnyih-letatelnyih-apparatov - Заглавие с экрана. [Search and rescue operations using unmanned aerial vehicles [Electronic resource] - Access: https://fireman. club/statyi-polzovateley/poiskovo-spasatelnyih-rabotyi-pri-pomoshhi-bespilotnyih-letatelnyih-apparatov]
[4] Обзор квадрокоптера DJI Phanton 3 Professional. [Электронный ресурс] - Режим доступа: https://mirquadrocopterov.ru/populyarnye-modeli/universalnye-drony/obzor-dji-phantom-3-professional.html - Заглавие с экрана. [DJI Phantom 3 Professional. [Electronic resource] - Access: https://mirquadrocopterov.ru/populyarnye-modeli/universalnye-drony/obzor-dji-phantom-3 -professional.html]