WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Сикорский Д.А. Метод оценки уровня турбулентности атмосферы

Научная статья

 

Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         824                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

Метод оценки уровня турбулентности атмосферы

Сикорский Д. А. fold@tut.by)

НПО "Лептон" Введение

В настоящее время широкое распространение получили космические оптико-электронные системы наблюдения (ОЭСН), решающие различные задачи обзора земной поверхности и околоземного пространства. К основным из них можно отнести системы космического картографирования, исследования природных ресурсов земли, экологического мониторинга, специального наблюдения и т. д.

Все ОЭСН требуют методического аппарата, позволяющего проводить априорные и апостериорные оценки их целевых характеристик, т.е. оценки качественных показателей формируемой информации. Оптико-электронный тракт (ОЭТ) таких систем включает в себя не только собственно оптико-электронную аппаратуру (ОЭА), но и турбулентную среду, в которой распространяется электромагнитное излучение. На рис.1 представлена обобщенная структурная схема ОЭТ.


 


Объект наблюдения


Среда распространения


Зрительный анализатор


Рис. 1. Обобщенная структурная схема ОЭТ

В качестве турбулентной среды выступает атмосфера. Ее оптические турбулентные свойства обусловлены случайными изменениями коэффициента преломления по трассе наблюдения, что вызывает искривление волнового фронта и, следовательно, приводит к размытию изображения в картинной плоскости ОЭА. Как правило, аппаратура наблюдения вышеотмеченных систем работает с временами накопления Тн=1-5 мс. Эти времена значительно меньше временного радиуса корреляции случайных турбулентных возмущений. В этом случае можно считать, что съемка объекта наблюдения осуществляется через так называемую "замороженную атмосферу". На рис.2 в одномерной геометрической интерпретации иллюстрируются варианты формирования суммарной функции рассеяния точки оптической системы и "замороженной атмосферы".

Р. Хафнагель [1] оценил функцию распределения интенсивности кружка рассеяния турбулентной атмосферы при наблюдении в полдень через всю её толщу для видимого диапазона. В дальнейшем кру-


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         825                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

жок рассеяния будем также называть функцией рассеяния точки (ФРТ) или апертурной функцией (АФ). Результат, полученный Р. Хафнаге-лем, представлен на рис.3, откуда видно, что максимальный угловой раствор усредненной апертурной функции турбулентной среды составляет около 16 угл. сек.

Полученные Р. Хафнагелем данные свидетельствуют о возможности значительного влияния оптической неоднородности среды распространения на качество формируемых изображений. Так, например, при наблюдении точечных источников с поверхности Земли по вертикальной трассе с использованием оптических систем (ОС)


Турбулентная среда


Плоскийфронт, формируемыйточечнымисточником

Неискаженныйфронтпри отсутствиитурбулентности


Оптическаясистема


ФРТоптическойсистемы



¦

/:-".u- Л.ур^улентн-ая1 среда *&*#

¦¦ ¦ ¦__________________

liii___ !_


V5^3

к'¦¦'.'-•' "Турбулентная'

_"J____ I___ "   I |Г ' '   '_______ ~^1        ¦    '¦_________ ^    "¦   ¦ 'ч-


Фронты, искаженные "замороженной" турбулентнойатмосферой


/V   /\


aJ±


А


Суммарныефункциирассеянияточкиоптическойсистемыитурбулентнойсреды

Рис. 2. Геометрическая одномерная интерпретация процесса формирования функции рассеяния точки при различных турбулентных искажениях, вносимых "замороженной атмосферой"


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         826                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

с фокусным расстоянием более 0,5 м, увеличение раствора кружка рассеяния будет составлять более 35 мкм, что значительно превышает уровень дифракционных и аберрационных искажений собственно оптической системы с D/F= 1/5-1/10.


Относительная

интенсивностьв

среднемкружке

рассеяния

Угловойрастворкружкарассеяния

Угловоерасстояниеотцентра усредненногокружкарассеяния

угл. сек

-4"    '      0      '   +4"        +8"

Рис. 3. Распределение интенсивности в кружке рассеяния турбулентной среды при наблюдении в полдень через всю её толщу

Автором с целью юстировки зеркального телескопа с фокусным расстоянием 1 м осуществлялось наблюдение удаленного (L«2,5 км) квазиточечного объекта на поверхности Земли по слабонаклонной (ф=5°) открытой трассе. Результаты измерений ФРТ телескопа показали, что раствор этой функции значительно превышает расчетное значение. По мнению автора это связано с негативным влиянием турбулентности среды распространения.

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

Для проведения этих оценок был использован известный инструментарий [2], описание которого приводится ниже.

1. Метод оценки пространственных свойств функциональной пары «оптическая система+фотоприемник» по обобщенной пограничной кривой

Для пояснения сути метода определим распределение выходного сигнала фотоприёмного устройства (ФПУ) с пространственно-регулярной дискретной структурой его элементов.

Пусть Н(х,у) - двумерное распределение яркости объекта наблюдения. Тогда распределение сигнала в непрерывном (гипотетическом) фотоприемнике


Q*(x,y)=к Н(х,у)**0(х,у)**Р(х,у),


(1.1)


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         827                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

где 0(х,у), Р(х,у) - апертурные функции ОС и элемента ФПУ, соответственно;

к - коэффицент преобразования светоэнергетических величин в электрические.

Реально Q*(x,y) дискретизируется в соответствии с пространственной "решеткой" используемого фотоприемника.

Выходной сигнал ФПУ

Q(x,y)= Q*(x,Y) S5(x,y),                                 (1.2)

+00       +оо

где   S5(x,y)= ?    Z6(x-Ax-m + AFx)-5(y-^Y-n + AFY)AxAY        (1.3)

П = -ооГП = -оо


здесь Ax, AY - периоды решетки ФПУ по соответствующим осям;

Afx, Afy - случайные фазы выборки.

Функция S5(x,y) называется дискретизирующей функцией или функцией выборки.

Суть методики сводится к оценке суммарной апертурной функции Т(х,у), равной свертке апертурных функций ОС и элемента ФПУ. Будем считать, что все приведенные апертурные функции пространственно инвариантны.

Учитывая (1.1), а также то обстоятельство, что преобразование Фурье свертки нескольких функций равно произведению преобразований Фурье свертываемых функций, можно записать:

T(v,m)=0(v,M) P(v,M).                                        (1.4)

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

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

Математическая запись нормированной гауссоиды вращения представляется в виде:

X

~~2

,2 \(     „2 Л

ехр

(1.5)

2

J

Т(х,у)=ехр

j


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»


828


http://zhurnal.ape.relarn.ru/articles/2002/075.pdf


Сечение гауссоиды вращения по оси х представлено на рис.1.1.

^-

1,0   м


Рис. 1.1. Сечение гауссоиды вращения. Здесь о - параметр гауссоиды

Наличие Т(х,у) изображающей системы с о>0 приводит к тому, что резконтрастное распределение яркости объекта Н(х,у) преобразуется в размытое изображение. Степень размытия тем больше, чем больше величина о.

Пусть объект Н(х,у) имеет двумерное распределение яркости резкого края:

0 при х<х0,

1 при хп >х,

НР.к.(х,у) =

(1-6) 1/2 при х0=0,

для любых у.

Функция распределения сигнала резкого края в непрерывном фотоприемнике Qp.K.(x) п0 оси х за счет наличия апертурной функции Т(х,у) однозначно связана с ней и определяется её видом.

Определим распределение Qp.K.(y) по оси У и Qp.k.(x) п0 оси х-Функции Qp.K.(x) и Qp.K.(y) называются пограничными кривыми резкого края (по соответствующим направлениям).

Значение двумерной функции распределения сигнала изображения в точке с координатами х, у в непрерывном фотоприемнике находится в соответствии с выражением (1.1).

Найдем производную по оси х, т.е. в направлении, перпендикулярном границе резкого края.


5Qp.K.(x>y)

дх


_д_ дх


+00+00

{  jHp.K.(xi,yi)-exp


(x-Xl)2

2&


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»


829


http://zhurnal.ape.relarn.ru/articles/2002/075.pdf



ехр


(y-yi)2

2а2


chqdy]


(1.7)


Так как НРК (х,у) не зависит от у и интегрирование осуществляется по переменным Xi yi, а также учитывая (1.5), последнее выражение запишется в виде:


Эх         I


ехр


(x-Xl)2


/


+оо

6х1 J ехр


(У-Yi)2 2о2


dy1 =



= V2noexp


х

2о'


(1.8)


где {QJx-обозначает операцию дифференцирования функции ( ) по

переменной х.

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

( х ^ л/2о.

X         (      t2   N\

ехр

-оо           V

(1.9)

dt = erf

<&c(x)=4uJ

2а'

2по

Найдем производную функции распределения <Зр.к.(х'У) по оси у,

совпадающей                с                границей               резкого               края


dQp.K.(*,yL д

дуду


+ GO+GO

| |Hp.K.(xi,yi)exp


2а2


ехр


(У-Yi)2 2а2


dXidy!


(1.10) Учитывая те же обстоятельства, отмеченные при нахождении производной <Зр.к.(х/У) по оси х, запишем:


+оо

dXlJ

+оо

{ехр

о |п

SQ?.K.(x,y)   +

ах

-о • ехр


(x-xt)2 2а2

" (у-уд)2

2а2


ехр

О


(y-yi)2

2о2


JJ Y


dy1 =

(1.11)


Поскольку производная функции <Зр.к.(х/У) резкого края по направлению у равна нулю, то сама функция Qp.K.(y) есть величина, постоянная для любых х.

Если в поле зрения изображающей системы присутствует объект, описанный выражением (1.5), то, измерив одномерную функцию распределения сигнала от этого объекта, формируемого в непрерывном фотоприемнике в направлении, перпендикулярном границе резкого края,  и  взяв от нее производную,  можно определить апертурную


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         830                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

функцию, то есть при приведенном условии (1.5) отыскать параметр гауссоиды вращения.

Нахождение параметра гауссоиды по распределению сигнала резкого края Q*(x,y) путем непосредственного определения пограничной кривой по ее выборкам, взятым с периодом дискретизации приемника, связано в определенных ситуациях со значительной неоднозначностью.

Так как в реальной ОЭА практически всегда присутствуют в той или иной степени искажения первого рода, то восстановление Q*(x,y) по ее отсчетам Q(x,y) никогда не бывает абсолютно точным. Указанное обстоятельство и приводит к неоднозначности определения QpK(x,y)  и, следовательно, определяет наличие ошибки в оценке

апертурной функции Т(х,у).

Эта неоднозначность еще более усугубляется при анализе реальной ОЭА, когда учитывается шумовая компонента, складывающаяся с отсчетами Q(x,y) и вызывающая изменение их значений.

Рис. 1.2 поясняет суть проявления указанной неоднозначности при определении Qp.k.(x/Y) п0 отсчетам пограничной кривой.

Н(х)

Х0

%¦¦

q;j

X

х)

......... 4

*.........

......... 4

Е^-"*""

""*"'...... i

i

--------------- <

Р"

->>

,...               QpAx)

т...

------------- >


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         831                 http://zhumal.ape.relarn.ru/articles/2002/075.pdf

Рис. 1.2. Два графика отмеченные пунктиром - результат восстановления функции Qp.K.(x) п0 ее отсчетам в условиях действия шума

Указанную неоднозначностьвозможновзначительнойстепени уменьшитьзасчетполучениябольшогоколичествакоррелированных отсчетовпограничной кривой на участке ее изменения. Искажения первого рода при этом существенно уменьшаются. Более того, усреднение некоррелированной шумовой составляющей по полученным отсчетам в еще большей степени устраняет отмеченную неоднозначность. Рис. 1.3 графически иллюстрирует данное положение.

Суть предлагаемого способа сводится к следующему.

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


QpjJx)

<Ж.М


Рис.1.3. Сглаживание функции  QpЈ(x) (пунктирная линия) по большому количеству ее отсчетов в условиях действия шума

Из выделенных объектов отбираются такие, у которых их прямолинейная граница расположена под малым углом ср по отношению к одному из направлений пространственной структуры приемника (рис.1.4).

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

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

Ниже рассматривается случай расположения границы смежных полей изображения резкого края под малым углом к направлению


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         832                http://zhumal.ape.relarn.ru/articles/2002/075.pdf


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

Из всей совокупности столбцов фотоприемника выбирают такой j-й столбец, элементы которого полностью пересекают границу изображения выбранного объекта (рис.1.4). Выборкиэтогостолбцаобразуют, такназываемую, обобщеннуюпограничнуюкривую (ОПК).

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

Поиск параметра о гауссоиды вращения осуществляется в следующем порядке.

Определяется область последовательных выборок j-ro столбца, которая несет информацию о степени размытия резкого края.

В соответствии с (1.9) огибающая сигналов элементов выбранного столбца аппроксимируется (например, по минимуму СКО) функцией

вида erf —-   путем подбора параметра о' (x'=x/sin^).

] -йстолбец

Одиночныйэлементфотоприемника

Размытаяобласть' резкогокрая

v2o')

Рис.1.4. Структура приемника определяет пространственно регулярный характер выборки распределения сигнала, формируемого в непрерывном приемнике. Из всей совокупности столбцов приемника


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         833                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

выбирается такой столбец, элементы которого пересекают размытую границу резкого края

Искомое значение параметра о, качественно характеризующее Т(х,у), определяется в направлении, перпендикулярном границе смежных полей объекта, т.е.

о = o'-sin (1-12)

На рис. 1.5 для пояснения процесса нахождения параметра о апертурной функции приведена его графическая интерпретация.

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

Аппроксимирующая ОПКфункция

ПериоданализирующейдискретизацииА

Огибающая ОПК

Размытаяобласть резкогокрая

]-йстолбец


Рис. 1.5. Геометрическая интерпретация процесса определения параметра апертурной функции


ОПК


2. Метод оценки уровня турбулентности атмосферы по


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         834                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

Суть метода поясняет рис. 2.1, на котором представлена блок-схема, отражающая этапы получения и обработки экспериментальных данных.

Л^

*-#*

*-'

Турбулентная среда отсутствует

. Турбулентная <Ереда \- ., ¦. присутствует. '


 


Алгоритм

обработки

отсчетов ОПК

ЭВМ

Результаты измерений


Рис. 2.1. Блок-схема эксперимента

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

Оператор проводит анализ полученных изображений с целью определения отсчетов ОПК. Полученные данные (отсчеты ОПК) вводятся в ЭВМ, где обрабатываются по алгоритму, представленному на рис.2.2.


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         835                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

Обработка осуществляется с целью определения раствора апер-турной функции Т(х,у). Разница между растворами апертурных функций, оцененных по ОПК для турбулентной среды и без нее, позволяет выделить вклад турбулентности среды в формирование Т(х,у).

Ввод отсчетов ОПК

*

Фильтрация ОПК медианой с ядром из пяти элементов

i

НЧ фильтрация ОПК линейным

звеном с апертурной функцией

вида ехр(-х2/2сГф)

Дифференцирование ОПК

Определение раствора апертурной функции по уровню 0.01

а)

Исходный массив

'          к

Результат медианной фильтрации

уРаствор апертуры- 27

Апертура фильтра

Результат низкочастотной фильтрации

^

Результат дифференцирования ОПК

^

W


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»                       836                 http://zhumal.ape.relarn.ru/articles/2002/075.pdf

б)

Рис.2.2. Блок-схема обработки экспериментально полученных отсчетов ОПК (а); пример результата обработки (б)

В соответствии с описанным методом автором был проведен натурный эксперимент, позволивший количественно оценить уровень турбулентности "замороженной атмосферы". В эксперименте использовалась промышленная оптическая система с фокусным расстоянием f=75 мм. Шаг решетки ФПУ составлял 4,65 мкм.

Однократная съёмка резкого края в помещении проводилась на расстоянии 1=20 м (l/f=266).

Относительные значения

сигнала в элементах

анализируемого столбца                                       Обобщенная ПОГраНИЧНаЯ КрИВаЯ

Угол- 8.333083 град

Оптическая система - ОКС ( фокус 75 мм)

Шаг решетки - 4.65 мкм

Диафрагма - 5.6

в)

Рис. 2.3. Изображение резкого края, полученное в помещении:

Съёмки через турбулентную атмосферу осуществлялись по слабонаклонной трассе. Объект наблюдения находился на расстоянии 1_=2,5 км. На рис. 2.3, 2.4 приведены примеры полученных экспериментальных данных.


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»        837          http://zhumal.ape.relarn.ru/articles/2002/075.pdf

а) изображение, имеющее резкую границу смежных полей; б) анализируемый фрагмент изображения; в) отсчеты ОПК




а)


в)


Рис. 2.4. Одно из серии изображений резкого края, полученное через турбулентную среду: а)изображение, имеющее резкую границу смежных полей; б) анализируемый франмент изображения; в) отсчеты ОПК

В процессе эксперимента было получено двадцать аналогичных рис. 2.4.а изображений. Съёмка объекта осуществлялась через равные промежутки времени в течение 10 минут. Результаты, количественно отражающие уровень турбулентности атмосферы на момент съемки, сведены в гистограмму, представленную на рис.2.5.


Частота появления

^


i------------- 1----------- 1---------- 1----------- 1---------- 1----------- 1у гл. сек

8.0                                          11.0                                          14.0

Угловой размер раствора апертурной функции турбулентной "замороженной атмосферы"


Электронный журнал «ИССЛЕДОВАНО В РОССИИ»         838                http://zhumal.ape.relarn.ru/articles/2002/075.pdf

Рис. 2.5. Оценка дифференциальной плотности распределения уровня турбулентности "замороженной атмосферы"

Полученные результаты могут быть трансформированы на ОЭСН, включающие ОС с любыми фокусными расстояниями, при условии равенства углов полей зрения исследуемой и измерительной ОС.

Выводы:

  1. Предложенный метод позволяет на количественном уровне оценить дифференциальную плотность распределения углового раствора апертурной функции турбулентной "замороженной атмосферы" для конкретных условий наблюдения.
  2. Полученные экспериментальные результаты показывают, что для горизонтальных приземных оптических трасс уровень турбулентности значителен. Это в сильной степени ограничивает качество формируемых ОЭСН изображений, особенно если фокусные расстояния оптических систем более 0.3-0.5 метра.

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

  1. Фриден Б. Компьютеры в оптических исследованиях/ Пер. с англ. под ред. С. А. Ахманова. М.: Мир, 1983. С. 190-198.
  2. Бобылев В. И. Методика апостериорной оценки ЧКХ бортовой части ОЭТ по характеристике резкого края/ Депонированная статья, вып. 5(32), сер. А, 1991.
  3. Прэтт У. Цифровая обработка изображений. Пер с англ. под ред. канд. техн. наук Д. С. Лебедева. М.: Мир, 1982.
  4. Исимару А. Распространение и рассеяние волн в случайно-неоднородных средах/ Пер. с англ. канд. физ.-мат. наук Л. А. Апресяна, А. Г. Виноградова, 3. И. Фейзулина. М.: Мир, 1981
  5. Госсорг Ж. Инфракрасная термография/ Пер. с франц. канд. техн. наук Н. В. Васильченко под ред. чл.-корр. АН СССР Л. Н. Курбатова. М.: Мир, 1988.
  6. Лансберг Г. С. Оптика. М.: Наука, 1976.
Мирошников М. М. Теоретические основы оптико-электронных
 



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.