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

 

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

СОЮЗ СОВЕТСКИХ

СОЦИАЛИСТИЧЕСКИХ

РЕСПУБЛИК

ÄÄSUÄÄ 1597816,А 1 (51)5 G 01 W 1/00

ОПИСАНИЕ ИЗОБРЕТЕНИЯ

H АBTGPCHOMY СВИДЕТЕЛЬСТВУ

ГОСУДАРСТВЕННЫЙ КОМИТЕТ

ПО ИЗОБРЕТЕНИЯМ И ОТКРЫТИЯМ

ПРИ ГКНТ СССР (21) 4449896/24-10 (22) 05.07.88 (46) 07.10.90. Бюл. № 37 (71) Научно-исследовательский институт прикладных физических проблем нм. А. И. Севченко (72) 1О. A. Полканов (53) 551,501.8:551.521.3(088.8) (56) Полканов Ю. А. О регулярности структуры неоднородностей рассеяния оптического излучения и их связи с метеопараметрами. — Известия АН СССР. Физика атмосферы и океана, т. 21, № 7, 1985, с. 720—

726. (54) СПОСОБ ОБНАРУЖЕНИЯ СТРУКТУРЫ ОПТИЧЕСКИХ НЕОДНОРОДНОСТЕЙ АТМОСФЕРЫ (57) Изобретение относится к измерению о««тичсских характеристик неоднородных сред и позволяет повысить надежность обнаружен и я элементов периодичности в структу ре неоднородностей. В атмосферу посылают имИзобретение относится к измерению оптических характеристик неоднородных сред и может применяться для целей метеорологии, физики плазмы, оптики атмосферы и океана, лазерной локации в условиях замутненной среды.

Цель изобретения — повышение надежности выявления периодичности оптических

««еоднородностей среды flo временной структуре сигнала обратного рассеяния в условиях воздействия шумов путем нелинейной фильтрации и анализа характера выделен-«oA структуры неоднородностей в зависимости от растующего нелинейного порога амплитудной дискриминации.

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

«/ У

ЛС (С+ С-)(а (С +2СО) где С вЂ” значение сигнала в локальном минимуме между максимумами С+« и Сч, а — коэффициент обеспечивающий вероятность (P««) соответствия рассчитанного перепада (Лс) перепаду обусловленному флуктуациями самого сигнала;

C0 — средний уровень шумов, от««eve««ный к интервалу временной д««скретизации сигнала, 1597816 находят в последовательности локальных минимумов пары, соседствующие с отброшенными максимумами и отбрасывают меньший минимум пары, формируют новые последовательности оставшихся локальных максимумов, минимумов, определяют наличие максимального и минимального максимумов (С .«., С. ), находят их положение в послед= довательности, отбрасывают все экстремумы

+ перед С - и после С, задают начальный уровень дискриминации локальных максимумов

d=(C + )

С формируют новую модельную последовательность отсчетов (C+,= 1, С, = — 1), на ее основе формируют распределение полученного числа интервалов между максимумами, минимумами определенной протяженности в зависимости от величины интервала, определяют величину, положение и число относительных максимумов (N„) в полученном распределении, увеличивают уровень дискриминации на заданный шаг (Ado), повторяют действия по определению величи ны, положения и числа относительных максимумов, запоминают эти числа, определяют перепад уровней дискриминации, соответствующий переходу к числу относительных максимумов распределения менше на одно, по их максимальному числу (N) судят о числе характерных элементов структуры неоднородностей по отношению выявленных перепадов уровней дискриминации к уровню фильтрации (Со+АС) судят о контрасте соответствующего элемента структуры относительно фоновых неоднородностей, по отношению соседних перепадов уровней судят об относительном контрасте элементов структуры неоднородностей.

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

Схема включает пульт 1 управления, импульсный источник света (ОКГ) 2, фотоприемник 3, усилитель 4, накопитель 5, амплитудно-цифровой преобразователь (АЦП)

6, блоки памяти 7, 8, 17, блок 9 вычитания фона, блок 10 выделенных локальных экстремумов, блок 11 задания уровня дискриминации, блок 12 нелинейной фильтрации и дискриминации локальных максимумов, блок

13 формирования распределения числа интервалов от протяженности интервала, блок 14 выделения локальных максимумов распределения, блок 15 определения числа локальных максимумов распределения (типов интервалов), блок 16 соотнесения типа интервала уровню дискриминации, блок 18 определения максимального и минимального числа типов интервалов, вычислитель 19 параметров обнаруженной структуры. С пульта 1 в накопитель 5 заносят интервал времени измерений, 5

55 шаг дискретизации и число циклов, по которому проводят накопление приходягцего сигнала. После этого с пульта 1 подают команду на излучение импульса света источником 2. Одновременно подают команду на усилитель 4 для усиления приходящего сигнала пропорционально квадрату текущего времени, начиная с момента излучения импульса. Кроме того, с пульта 1 задают число шагов дискретизации, запоминаемое в блоке 7, и последовательность предварительно заданных команд обработки сигнала.

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

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

В блоке 9 производится вычитание уровня фона соотнесенного к шагу дискретизации сигнала (по данным из блока 7) из полученных значений сигнала. Полученные рабочие значения сигнала обратного рассеяния рабочая последовательность вновь запоминаются в блоке 8.

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

В зависимости от выставленного уровня фильтрации (С ), выдаваемого пультом 1, здесь происходит фильтрация локальных максимумов, когда значения большие ЛС сохраняются. Если в результате фильтрации отфильтровывается локальный максимум, то из двух соседствующих минимумов за действительный принимается наибольший, т. е. в блоке 11 производится дополнительная операция (схемой сравнения и выбора) .

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

1597816

20

5 новании данных об уровне шумов (Сю), поступающих из блока 8 памяти. На выходе блока 12 фильтрации имеем набор локальных максимумов, минимумов, разделенных определенными интервалами, число которых меньше их исходного числа (без дискриминации) . В блоке 13 происходит сортировка этих максимумов и минимумов в две независимые последовательности, генерация нормированной последовательности, отбор интервалов между импульсами нормированной амплитуды (+1 или — 1) их сортировка по величине и суммирование числа интервалов определенной протяженности. В результате по команде с пульта 1 получаем на выходе блока 13 распределение числа интервалов в зависимости от протяженности интервала («спектр» интервалов). Следующий этап анализа — анализ зависимости спектра интервалов от уровня дискретизации. Исходный спектр и отфильтрованный сигнал — в блоке 12. С пульта 1 и блока 11 на входы блока 12 поступают команды на выставление уровня дискриминации, который растет от цикла к циклу анализа с заданным шагом (Лйю).

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

+ у глобального максимума (С- ), затем такой же следующий и т. д. пока не исчезнет экстремум у глобального максимума (Сч") ..

Происходит прореживание локальных экстремумов сигнала, граница которого перемещается в направлении С- с — С" . Причем исчезновение локальных экстремумов идет с оп редел ен ным интер валом (последова тел ьно) и уменьшение их числа на одно соответствует определенному приращению уровня дискриминации Ad. После фильтрации определенного типа локальных экстремумов, соответствую щи х существованию на чального элемента периодической структуры, существенный рост уровня дискриминации (Лд))

))Адю) не приводит к изменению общей картины в последовательности следования локальных экстремумов. С некоторого значения d опять начинается селекция экстремумов и их прореживание. Это вызывает появление нового типа интервалов между оставшимися экстремумами, число которых растет с ростом уровня дискриминации.

Причем, как и в начальной последовательности, уменьшение локального экстремума на один соответствует приращению Ad, т. е. к появлению еще одного интервала нового типа. С ростом d число этих интервалов растет. Они появляются последовательно в

1- + направлении С - -+ С -. Максимальное число этих интервалов равно отношению протяженности всего анализируемого участка на протяженность такого интервала. Если есть еще периодические компоненты структуры (последующие ее элементы), то с ростом картина прореживания экстремумов сигнала и появления интервалов все большей протяженности повторяется. В случае отсутствия локальных экстремумов блок 12 работает как обычный ограничитель по амплитуде.

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

Распределение интервалов («спектр» интервалов) формируется блоком 13 на основе полученных после фильтрации, дискриминации и нормировки данных. В блоке 14 анализируется это распределение на предмет выявления относительных максимумов по методике поиска локальных максимумов блоком 10. Данные о величине положения каждого относительного максимума поступают на блок 15, где число относительных максимумов суммируется и запоминается.

Эти данные поступают на вход блока 16 одновременно с набором взятых уровней дискриминации от начальнс о до текущего (с блока 11). Здесь сравнивается номер уровня и номер интервала и в случае их совпадения в блоке 17 запоминается пара значений: уровень дискриминации н интервал выявленного типа.

Рост уровня дискриминации при периодичности структуры неоднородностей приводит к росту числа типов интервалов и их последующему уменьшению до одиночного.

Текущее число типов интервалов для определенного уровня d с блока 15 и пары значений уровень — дискриминации — интервал выявленного типа с блока 17 поступают на входы блока 18 и запоминаются. Число типов интервалов сравнивается с предшествующим числом для меньшего уровня дискриминации. При его неизменности или росте с блока 18 подается команда на блок 11 о выставлении следующего уровня дискриминации и повторении цикла анализа. При уменьшении типов интервалов до одного в блоке 18 определяется наибольшее число типов интервалов и запоминается. Одновременно анализируются пары уровень дискриминации — + интервал выявленного типа для всего диапазона уровней дискриминации и определяются перепады уровней Л при изменении числа типов интервалов на единицу.

Эти данные поступают на вычислитель 1 одновременно с данными блока 17 и по максимальному числу типов интервалов уста1597816 навливается число элементов периодичности ст ру кту р ы неоднородностей.

Каждому элементу структуры соотносится перепад уровней дискриминации Л„. Из перепадов вычитается предварительно запоминаемый уровень, соответствующий фону, и полученные значения (Л вЂ” Л,) являются характеристикой контрастности элементов структуры относительно фоновых неоднородностей. Здесь же вычисляются отношения полученных перепадов (Л вЂ” Л ) / (Л вЂ” Л ), которые принимаются за характеристики относительного контраста элементов структуры неоднородностей.

Формула изобретения

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

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

20 максимумов меньше на одно, причем порог фильтрации, определяют по формуле

+ (С +2Сф

С где С+ — минимально допустимая величина локального максимума;

С вЂ” значение сигнала в локальном минимуме между максимумами каждой пары;

Со — средний уровень шумов, нормиро30 ванный на интервале дискретизации сигнала; а — коэффициент, обеспечивающий требуемую вероятность соответствия максимума перепаду, обусловленному флуктуациями самого сиг35 нала, а начальный уровень дискриминации определяют по формуле — С - „1/ 2

d = (Cìèí + -), 1597816

Составитель С. Скуратов

Редактор С. Патрушева Техред A. Кравчук Корректор А. Обручар

Заказ 3053 Тираж 34 l Подписное

ВНИИПИ Государственного комитета по изобретениям и открытиям при ГКНТ СССР

I !3035, Москва, Ж вЂ” 35, Раушская наб., д. 4/5

Производственно-издательский комбинат «Патент», r Ужгород, ул. Гагарина, 10!