Научная статья на тему 'Расчет параметров тепловизионной съемки объектов с беспилотных авианосителей'

Расчет параметров тепловизионной съемки объектов с беспилотных авианосителей Текст научной статьи по специальности «Физика»

CC BY
213
57
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ТЕПЛОВИЗИОННАЯ СЪЕМКА / БЕСПИЛОТНЫЕ АВИАНОСИТЕЛИ / МИНИМАЛЬНАЯ ДЕТЕКТИРУЕМАЯ И РАЗРЕШАЕМАЯ ТЕМПЕРАТУРЫ / МОДУЛЯЦИОННАЯ ПЕРЕДАТОЧНАЯ ФУНКЦИЯ / THERMAL IMAGING / UNMANNED AERIAL VEHICLES / MINIMUM DETECTABLE AND RESOLVED TEMPERATURE / MODULATION TRANSFER FUNCTION

Аннотация научной статьи по физике, автор научной работы — Катковский Л. В.

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

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

Похожие темы научных работ по физике , автор научной работы — Катковский Л. В.

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

CALCULATION OF OBJECTS THERMAL IMAGING PARAMETERS FROM UNMANNED AERIAL VEHICLES

This article is aimed at studying the influence of the entire complex parameters on the characteristics of the obtained thermal images when shooting the Earth’s surface from unmanned aerial vehicles. The values of the minimum detectable and minimum resolvable temperature differences are calculated depending on the parameters of the thermal imager, the survey and the size of the object (spatial frequency) for three serial miniature thermal cameras used in aerial surveys of the Earth's surface from unmanned aerial vehicles. Analytical formulas for estimating the minimum resolvable temperature difference are obtained on the basis of a mathematical model of the thermal imager as a linear system of individual components based on the technique that differs from the generally accepted one. Estimates were made for two cases: observation of a thermal image by an operator on a display screen and for the case when an electronic image is analyzed by a threshold algorithm with no operator engaged. For the first time, the influence of the carrier velocity on the overall modulation transfer function of the system and, accordingly, the temperature and spatial resolution of thermal imagers was taken into account. The main components that must be considered when calculating the full modulation transfer function of the system are: a thermal imager lens, a radiation detector, carrier movement and the observer's visual system. Moreover, the parameters of the detector and the speed of the carrier have the greatest influence on the temperatures resolvable by the system.

Текст научной работы на тему «Расчет параметров тепловизионной съемки объектов с беспилотных авианосителей»

Щ

http://dx.doi.org/10.35596/1729-7648-2020-18-2-53-61

Оригинальная статья Original paper

УДК 621.384.3

РАСЧЕТ ПАРАМЕТРОВ ТЕПЛОВИЗИОННОЙ СЪЕМКИ ОБЪЕКТОВ С БЕСПИЛОТНЫХ АВИАНОСИТЕЛЕЙ

Научно-исследовательское учреждение «Институт прикладных физических проблем им. А. Н. Севченко» Белорусского государственного университета (г. Минск, Республика Беларусь)

© Белорусский государственный университет информатики и радиоэлектроники, 2020

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

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

Конфликт интересов. Автор заявляет об отсутствии конфликта интересов.

Для цитирования. Катковский Л.В. Расчет параметров тепловизионной съемки объектов с беспилотных авианосителей. Доклады БГУИР. 2020; 18(2): 53-61.

КАТКОВСКИЙ Л.В.

Поступила в редакцию 27 января 2020

CALCULATION OF OBJECTS THERMAL IMAGING PARAMETERS FROM UNMANNED AERIAL VEHICLES

LEONID V. KATKOVSKY

Scientific research institution «Institute of Applied Physics Problems named after A. N. Sevchenko» of Belarusian State University (Minsk, Republic of Belarus)

Submitted 27 January 2020

© Belarusian State University of Informatics and Radioelectronics, 2020

Abstract. This article is aimed at studying the influence of the entire complex parameters on the characteristics of the obtained thermal images when shooting the Earth's surface from unmanned aerial vehicles. The values of the minimum detectable and minimum resolvable temperature differences are calculated depending on the parameters of the thermal imager, the survey and the size of the object (spatial frequency) for three serial miniature thermal cameras used in aerial surveys of the Earth's surface from unmanned aerial vehicles. Analytical formulas for estimating the minimum resolvable temperature difference are obtained on the basis of a mathematical model of the thermal imager as a linear system of individual components based on the technique that differs from the generally accepted one. Estimates were made for two cases: observation of a thermal image by an operator on a display screen and for the case when an electronic image is analyzed by a threshold algorithm with no operator engaged. For the first time, the influence of the carrier velocity on the overall modulation transfer function of the system and, accordingly, the temperature and spatial resolution of thermal imagers was taken into account. The main components that must be considered when calculating the full modulation transfer function of the system are: a thermal imager lens, a radiation detector, carrier movement and the observer's visual system. Moreover, the parameters of the detector and the speed of the carrier have the greatest influence on the temperatures resolvable by the system.

Keywords: thermal imaging, unmanned aerial vehicles, minimum detectable and resolved temperature, modulation transfer function.

Conflict of interests. The author declares no conflict of interests.

For citation. Katkovsky L.V. Calculation of objects thermal imaging parameters from unmanned aerial vehicles. Doklady BGUIR. 2020; 18(2): 53-61.

Введение

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

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

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

В работе проводятся расчеты пространственного и температурного разрешений тепловых изображений серийных тепловизоров при съемке с БПЛА. При этом впервые учитывается влияние на них скорости носителя. Расчеты проведены для двух случаев: наблюдение изображений на дисплее оператором [3-5] и в случае отсутствия наблюдателя и дисплея, когда анализ изображения проводится на борту беспилотного носителя [6].

Уравнения и расчеты температурного разрешения тепловизоров

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

Таблица 1. Характеристики исследуемых тепловизоров Table 1. Characteristics of the studied thermal cameras

Параметры камеры 1 2 3

Camera options ICI8640 P FLIR Vue Pro R thermoMap

Спектральный диапазон, мкм 7-14 7,5-13,5 7,5-13,5

Частота кадров, Гц 30 30 7,5

Погрешность, °C 1 (+/-) 5 (+/-) 5

Количество пикселей 640x512 640x512 640x512

Радиометрическое разрешение, бит 14 14 14

Шаг пикселя, мкм 17 17 17

Эквивалентная шуму разность температур, °C 0,02 0,05 0,1

Фокусное расстояние, мм 13 13 9

Диафрагменное число 1,0 1,25 1,4

При относительно длительном времени наблюдения объекта (неподвижная система) используют обычно две характеристики: минимальную обнаруживаемую разность температур ATd (Detectable) и минимальную разрешаемую разность температур ATr (Resolvable) [4]. Первая - как мера обнаружения объекта определенных размеров, вторая - как мера распознавания объекта на шумовом фоне. Часто ограничиваются расчетами минимальной разрешаемой разности температур [5, 6], которая является более информативной и фактически устанавливает связь между пространственным и температурным разрешениями. Указанные разности температур, в свою очередь, зависят от размеров объекта (пространственной частоты) и эквивалентной шуму разности температур ATnetd (Noise Equivalent Temperature Difference), которая характеризует температурную чувствительность тепловизора и является, как правило, паспортным параметром, определяемым производителем. Как известно, ATnetd может быть рассчитана по следующей формуле [4]:

дт Kr24AD¥ __W_ (i)

D ч Л (х, т) . , ч2 г—h.. ,дИх (х,т) , ()

Af TobAob f D (X) ) dX (Dob/Fob )2 TobJAD f D (X) ^ т ) dX

где Aob, Dob, Fob - площадь, диаметр и фокусное расстояние входного объектива соответственно, R - расстояние до объекта, AD - площадь чувствительной площадки приемника (пикселя), если пиксель квадратный со стороной ав, то AD _ a2D. Вторая часть формулы (1) получается из первой более общей, если изображение протяженного объекта полностью заполняет пиксель, т. е. необходимо положить: площадь наблюдаемого объекта

Af _ Ad (RjFob) и площадь входного объектива Ал _%DlblA . Af - эквивалентная

шумовая полоса частот, D (Я) - удельная обнаружительная способность фотоприемника, тоЬ - пропускание объектива, X - длина волны излучения, T - абсолютная температура. Обычно считают, что объект излучает по закону Ламберта с коэффициентом излучения s . В этом случае спектральная функция светимости М% (Я,T) определяется формулой Планка.

Заданные величины ATnetd позволяют рассчитать по формуле (1) величины D исследуемых камер и сравнить чувствительности их фотоприемников. Обычно ATNETD рассчитывается в предположении, что тоЬ =1, и для эффективной ширины полосы электронного тракта используется соотношение [4] (t0 - время экспозиции одного элемента матрицы, ff-частота кадров):

Af = l/2t0 = ff/ 2. (2)

С учетом данных табл. 1 для значений соответствующих величин D получаем dc 8640 p =1,5109 смТц1/2-Вт-1, Dflr Vue pr0r =1,0109 смТц12^Вт-1, D**herm0Mtap =3,2^108 смТц1/2-Вт-1.

Оценим минимальную обнаруживаемую разность температур ATd для этих камер и неподвижного относительно камеры объекта. Для этого воспользуемся следующей формулой [4] для случая наблюдения изображения оператором на экране дисплея:

ATd (S) = (2v • ATnetd • SNRd)(«DP^ff*e f L r). (3)

Выражение (3) получено в предположении, что тест-объект в пространстве изображений представляет собой квадрат с угловым размером £ на однородном фоне. В этом выражении v = 1/2^ - угловая пространственная частота (мрад-1), соответствующая угловому размеру объекта £=l/H (l - линейный размер объекта, H - расстояние до объекта - высота съемки), отношение сигнал/шум для обнаружения объекта, SNRD, зависит от задаваемой вероятности обнаружения. Для вероятности 0,9 принято полагать SNRD =4,5. aD, PD - угловые размеры пикселя фотоприемника в пространстве изображений, aD = PD = aD/Fob , tE = 0,2 с -временная постоянная (временное разрешение) глаза человека, осуществляющего интегрирование кадров (величина ^fftE - количество кадров, интегрируемых зрительной

системой), Ls r) - среднее нормированное значение яркости изображения объекта на экране

дисплея в зависимости от нормированного углового размера тест-объекта, его находят обычно на основе заранее рассчитанных таблиц значений, r - угловой размер радиуса кружка рассеяния камеры (оптической системы и приемника) - параметр функции рассеяния точки. Обычно модуляционная передаточная функция (МПФ) всей оптико-электронной системы h (v) хорошо аппроксимируется гауссовой функцией (по одной из пространственных координат) h(v) = exp(-2л2г2v2), параметр этой функции r, при условии согласования МПФ оптики

и приемника, определяется размером пикселя приемника и равен [4] r = 0,282 • aD .

Результаты расчетов по формуле (3) минимальной обнаруживаемой разности температур ATd для трех рассматриваемых камер и 4-х размеров объекта (l=5, 10, 20, 50 см) с высоты съемки H=100 м приведены в табл. 2 (во второй колонке приведены соответствующие размерам объекта l пространственные частоты v).

Таблица 2. Минимальные обнаруживаемые разности температур ATD для трех камер с высоты H=100 м Table 2. Minimum detectable temperature differences ATd for three cameras from the height of H = 100 m

Параметр Parameter Частота, v=H/2l, мрад-1 Frequency, v=H/2l, mrad-1 1 ICI8640 P ATd 2 FLIR Vue Pro R ATd 3 thermoMap ATd

atnetd , K 0,02 0,05 0,1

l = 5 см 1 0,83 2,1 18,3

l = 10 см 0,5 0,21 0,52 4,23

l = 20 см 0,25 0,059 0,15 1,12

l = 50 см 0,1 0,015 0,037 0,25

Понятно, что формула (3) не может применяться для любых размеров объекта, поскольку по своему физическому смыслу ATd не может быть меньше ATNlTD . Ограничение применения (3) обусловлено выполнением предположений, при которых она была получена [4]. Убывание ATd в (3) определяется как убыванием v с увеличением размера объекта, так и ростом безразмерной функции Ls (£/r) с ростом £ = 1/2v . Из условия ATD > ATnetd

можно установить, что формула (3) применима для угловых размеров объекта £,=l/H < 3,1aD. Заметим также, что выражение для минимальной обнаруживаемой разности температур не содержит МПФ, поскольку тест-объектом при ее получении является квадрат, а не мира, и отношение сигнал/шум вычисляется через величину яркости объекта в пространственной области (v здесь введена как параметр объекта для единообразия).

Перейдем к расчетам минимальной разрешаемой разности температур ATr, которая представляет собой минимальную разницу температур между штрихами миры определенной пространственной частоты и фоном, при которой возможно различение глазом на экране штрихов миры. Эта функция частоты устанавливает связь между пространственным и температурным разрешениями. При вычислении здесь отношения сигнал/шум используется МПФ, и вкачестве объекта рассматривается тепловая тест-мира (мира Фуко), состоящая из 4-х прямоугольных полос, высота которых в семь раз превышает их ширину [4].

Для ATr при указанных условиях в [3, 4] получено выражение (периодичная структура миры выбирается вдоль оси x и частота v относится к этой же оси)

2aöß

( . о Г'5 ^ ( о V'5

= 0,66

ATNEtd • SNRRYsv

aDßD

Д/• to'./We

(4)

8-h(v) ^7А/Ч-ff -tE J h(v)

В этом выражении в дополнение к определенным ранее величинам: h (v) - МПФ всей системы на частоте v, равная произведению МПФ отдельных звеньев тепловизионной системы наблюдения, Г - угловое увеличение системы камера - дисплей. Для матричного приемника

Af -t0 =0,5. Примем также, что Г =1 (оно обычно невелико), aD =рв, тогда (4) перепишется так (выделен множитель С, независящий от частоты):

ATR (v) = Сv/h(v), С = (0,93-ATNetd-SWRr-aD)/(ff^ f . (5)

В большинстве практических случаев электронный блок и дисплей не искажают изображения, которое создает тепловизор [5]. Поэтому считают, что МПФ электронного блока и дисплея равняются единице в пределах рабочей полосы частот. Таким образом, для МПФ тепловизора h (v) достаточно учесть МПФ оптической системы (объектива) hQ (v),

пространственную hB (v) и временную hDt (v), МПФ матричного приемника излучения, а также МПФ глаза наблюдателя hE (v) . Кроме этих компонентов МПФ, при движении объекта относительно камеры необходимо ввести в МПФ еще множитель, связанный с движением объекта (motion) hm (v), который не рассматривался ранее, в том числе при аэрокосмической

съемке [6], но который оказывает существенное влияние (как будет видно ниже) на МПФ и разрешаемую температуру системы. Таким образом, для полной МПФ имеем выражение h (v) = h0 (v) hD (v) hDt (v) hE (v) hm (v). (6)

Вычислим последовательно все компоненты МПФ. Для МПФ объектива воспользуемся приближением линейной зависимости от частоты [4]:

ho (v) = 1 -(1 -h)(v/vi), (7)

где Yj - пространственная частота, при которой значение МПФ уменьшается до значения hi.

Пространственная часть МПФ матричного приемника hD (v) определяется в основном

его чувствительным элементом - угловым размером пикселя aD , т. е.

h (v) = sin c(aDv) = sin (rcaDv)/(rcaDv). (8)

Эта МПФ зависит, таким образом, от соотношения угловой пространственной частоты наблюдаемой картины (миры) и углового мгновенного поля зрения прибора. Параметры

Доклады БГУИР Doklady BGUIR № 18 (2) (2020)_No. 18 (2) (2020)

h\ и vi функции (7) можно найти из условия согласования пространственной МПФ объектива и приемника по принципу равнозначности, считая что они одинаково ухудшают качество изображения, так что на частоте разрешения v контраст каждой подсистемы равен 0,5:

1 "(1 -h) = sin(™Dvres ) = 0,5 . (9)

Решение второго из уравнений (9) есть aDvres = 0,6 , тогда получаем h0 (v) = 1 -(ад/ 1,2)v. (10)

Для МПФ глаза наблюдателя h (v) используем аппроксимацию Шульца [7], которая, как считается [5], наилучшим образом аппроксимирует экспериментальную МПФ глаза: hE (v) = 2,63 • [exp(-1,92v) - exp (-5,67v)] . (11)

Временная МПФ hDt (v) приемника излучения, учитывающая его инерционность [6] hDt (v) = (1 + 4n%f2 )"0,5 =(1 + 4^2t¿af2v2 )"0,5, (12)

где tD - постоянная времени приемника излучения (для микроболометра обычно tD = 10 мс), и подставлена связь между пространственной угловой частотой v и временной частотой f в соответствии с формулой f = (aD/t0)v = aDfv ; время формирования одного элемента изображения t0 для микроболометрической матрицы определяется частотой кадров t0 = 1jf .

Очевидно, что для МПФ, учитывающей движение объекта (носителя), hm (v), можно

применить аналогичную (12) передаточную функцию эталонного фильтра, в которую войдут скорость движения носителя V и период миры l:

hm (v) = [1 + 4%2tD (V/l)2J^ = [1 + (2vV¡H)2j"'. (13)

Отметим, что обычно при расчетах ATR (v) используется методика [5, 6], когда для

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

учета МПФ электронного тракта и глаза вводится коэффициент ширины полосы, показывающий отличие реальной полосы пропускания от эффективной шумовой полосы эталонного фильтра А/ в (4). В данной работе используется иной подход, при котором МПФ всех звеньев системы учитываются непосредственно в общей МПФ (6), а эффективная ширина полосы соответствует эталонному фильтру (формула (2)).

Наконец, перейдем к случаю оценки ATr, когда в системе отсутствует оператор и анализируется электронное изображение миры Фуко с помощью порогового устройства. При этом исключается не только МПФ глаза из формулы (6), но и из (4) исключаются факторы, учитывающие временное и пространственное интегрирование изображений глазом оператора. В результате вместо (5), (6) получаем следующее выражение для минимальной разрешаемой разности температур [6]:

ATR (v) = (я • ATnetd • SNRr )/(4 • h* (v)), h' (v) = ho (v)hD (v)ha (v)hm (v) (14)

Результаты и их обсуждение

На рис. 1 изображены МПФ отдельных звеньев тепловизора ICI 8640 P с наименьшей эквивалентной шуму разницей температур 0,02 К и общая МПФ для трех камер в зависимости от угловой пространственной частоты наблюдаемой миры Av (мрад-1). Расчеты проведены, как и в табл. 2, для высоты съемки #=100 м и скорости движения БПЛА V = 5 м/с. Заметим, что типичный диапазон скоростей БПЛА составляет 0-15 м/с и диапазон высот 0-500 м. Общие МПФ для камер 1 и 2 совпадают ввиду одинаковых параметров, от которых зависит МПФ.

Для того чтобы оценить влияние МПФ глаза оператора и МПФ движения БПЛА на общую МПФ системы, на рис. 2 показаны также МПФ камеры ICI 8640 P без учета МПФ глаза и МПФ движения. Как видим, влияние МПФ глаза существенно только при малых частотах (< 0,2 мрад-1), тогда как влияние движения носителя значительно во всем диапазоне частот.

Расчет минимальных разрешаемых температур для ИК камер выполняется по формулам (5), (6). На рис. 3 показаны зависимости детектируемых и разрешаемых камерами

разностей температур в зависимости от частоты, причем для камеры 1 приведены еще кривые: без учета МПФ глаза наблюдателя и без МПФ движения носителя. Включение в схему глаза наблюдателя несущественно меняет разрешаемый температурный контраст, тогда как движение носителя, даже с относительно небольшой скоростью (5 м/с), существенно снижает разрешаемую разность температур. Заметим, что на графиках для малых частот, если ДТя становилось меньше АТттв (что физически неверно), полагалось ДТ&= АТттв . Детектируемые разности температур в основном меньше, чем разрешаемые, хотя для малых частот они могут быть и больше, что, возможно, связано с улучшением разрешения при сближении угловых размеров пикселя и объекта (пересечение кривых ДТп и ДТо наблюдается при частоте, соответствующей угловому размеру пикселя). На рис. 4 показаны для сравнения разрешаемые разности температур камер для случая наблюдения оператором (5) и для случая машинной обработки (14). При наблюдении на экране разрешаемые температуры ниже, что объясняется временным и пространственным интегрированием изображений глазом (сравните выражения (4) и (14) без учета И (у) и И* (у)) и что для рассматриваемых параметров заметно улучшает отношение сигнал/шум.

Упповая частота, 1/мрад

Рис. 1. МПФ отдельных звеньев для камеры ICI 8640 P и общие МПФ для трех камер: 1 - МПФ глаза; 2 - временная МПФ пикселя; 3 - МПФ объектива; 4 - пространственная МПФ пикселя; 5 - МПФ движения носителя; 6 - общая МПФ камер 1 и 2 (табл. 1); 7 - общая МПФ камеры 3 Fig. 1. MTF of individual components for the ICI 8640 P camera and the entire MTF for all three cameras: 1 - MTF of eye; 2 - temporary MTF of pixel; 3 - MTF of lens; 4 -spatial MTF of pixel; 5 - MTF of carrier motion; 6 - total MTF of cameras 1 and 2 (Table 1); 7 - total MTF of camera 3

Рис. 2. МПФ камеры ICI 8640 P: 1 - общая без учета МПФ глаза; 2 - общая без учета МПФ движения БПЛА; 3 - общая с учетом

всех звеньев системы Fig. 2. MTF of the ICI 8640 P camera: 1 - MTF without taking into account the MTF of the eye; 2 - MTF without taking into account the MTF of UAV movement; 3 - total MTF taking into account all parts of the system

Рис. 3. Разрешаемые ATR и детектируемые ATD разности температур для камер 1, 2, 3 (из табл. 1): 1 - ATd камеры 3; 2 - ATr камеры 3; 3 - ATr камеры 2; 4 - ATr камеры 1; 5 -ATr камеры 1 без учета МПФ глаза; 6 - ATr камеры 1 без учета МПФ движения; 7 - ATd камеры 2; 8 - ATd камеры 1 Fig. 3. Resolvable ATr and detectable ATd temperature differences for cameras 1, 2, 3 (from Table 1): 1 - ATd of camera 3; 2 - ATr of camera 3; 3 - ATr of camera 2; 4 - ATr of camera 1; 5 - ATr of camera 1 excluding eye MTF; 6 - ATr of camera 1 excluding motion MTF; 7 - ATd of camera 2; 8 - ATd of camera 1

Рис. 4. ATr для камер 1, 2, 3: при компьютерном анализе - кривые 1, 3, 4 для камер 3, 2, 1 соответственно; при наблюдении оператором на экране дисплея - кривые 2, 5, 6 для камер 3, 2, 1 Fig. 4. ATr for cameras 1, 2, 3: by computer analysis - curves 1, 3, 4 for cameras 3, 2, 1, respectively; observed by the operator on the display screen - curves 2, 5, 6 for cameras 3, 2, 1 respectively

Заключение

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

Расчеты проведены как для случая наблюдения изображения оператором на дисплее, так и для случая отсутствия дисплея и наблюдателя (машинная обработка). Показано, что сама по себе МПФ глаза незначительно влияет на разрешаемую температуру, тогда как влияние пространственной и временной интегрирующей способности глаза может значительно уменьшить разрешаемые температуры. Впервые учтено влияние на разрешаемые температуры скорости движения носителя и показано, что оно значительно даже при небольшой скорости движения (5 м/с при высоте съемки 100 м), и эту компоненту МПФ необходимо учитывать при тепловизионной съемке с движущихся платформ.

Список литературы

1 Uto K., Seki H., Saito G., Kosugi Y. Characterization of Rice Paddies by a UAV-Mounted Miniature Hyperspectral Sensor System. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 2013;6(2):851-860. DOI: 10.1109/JSTARS.2013.2250921.

2 Cen H. Y., Wan L., Zhu J. P. Dynamic monitoring of biomass of rice under different nitrogen treatments using a lightweight UAV with dual image-frame snapshot cameras. Plant Methods. 2019;5(1). DOI: 10.1186/s 13007-019-0418-8.

3 Ллойд Дж. Системы тепловидения. Москва: Мир; 1978.

4 Колобродов В.Г., Лихолгг М.1. Проектування тепловЫйних i телевЫйних систем спостереження. Кшв: НТУУ «КП1»; 2007.

5 Аль-Мзирави А., Колобродов В.Г., Микитенко В.И. Расчет минимальной разрешаемой разности температур тепловизионных приборов. Приборы и методы измерений. 2015;6(1):64-69.

6 Колобродов В.Г., Лихолгг М.1., Тягур В.М.. Мшмальна роздшювана рiзниця температур тепловiзора аерокосмiчного базування. Космiчна наука i технологiя. 2014;20(1):23-27. DOI: 10.15407/knit2014.01.023.

7 Карасик В.Е., Орлов В.М. Лазерные системы видения. Москва: Изд. МГТУ им. Н.Э. Баумана; 2001.

References

1 Uto K., Seki H., Saito G., Kosugi Y. Characterization of Rice Paddies by a UAV-Mounted Miniature Hyperspectral Sensor System. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 2013;6(2):851-860. DOI: 10.1109/JSTARS.2013.2250921.

2 Cen H. Y., Wan L., Zhu J. P. Dynamic monitoring of biomass of rice under different nitrogen treatments using a lightweight UAV with dual image-frame snapshot cameras. Plant Methods. 2019;5(1). DOI: 10.1186/s13007-019-0418-8.

3 Lloyd J. [Thermal imaging systems]. Moscow: Mir; 1978. (In Russ.)

4 Колобродов В.Г., Лихотт М.1. [Design of thermal and television surveillance systems]. Кшв: НТУУ «КП1»; 2007. (In Ukranian)

5 Al-Mziravi A., Kolobrodov V.G., Mikitenko V.I. [Calculation of the minimum resolved temperature difference of thermal imaging devices]. Pribory i metody izmereniy=Devices and methods of measurements. 2015;6(1):64-69. (In Russ.)

6 Колобродов В.Г., Лихолгт М.1., Тягур В.М. [Minimum resolvable temperature difference for thermal imager of space basing]. Kosmichna nauka tehnologoii =Space science and technology. 2014;20(1):23-27. DOI: 10.15407/knit2014.01.023. (In Ukranian)

7 Karasik V.E., Orlov V.M. [Laser vision systems]. Moscow: Publishing house MGTU im. N.E. Bauman; 2001. (In Russ.)

Сведения об авторах

Катковский Л.В., д.ф.-м.н., доцент, заведующий лабораторией дистанционной фотометрии НИИПФП им. А. Н. Севченко Белорусского государственного университета.

Адрес для корреспонденции

220045, Республика Беларусь, г. Минск, ул. Курчатова, д. 7, Научно-исследовательское учреждение «Институт прикладных физических проблем им. А. Н. Севченко» Белорусского государственного университета тел. +375-17-396-44-09; [email protected], katkovskyl@gmail. com Катковский Леонид Владимирович

Information about the authors

Katkovsky L.V., D.Sci., Associate Professor, Head of the Laboratory of Remote Photometry NIIPFP named by A. N. Sevchenko of Belarusian State University.

Address for correspondence

220045, Republic of Belarus, Minsk, Kurchatova str., 7, Research Institution "Institute of Applied Physical Problems named by A. N. Sevchenko" of Belarusian State University tel. + 375-17-396-44-09; [email protected], [email protected] Katkovsky Leonid Vladimirovich

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