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

СОЮЗ СОВЕТСКИХСОЦИАЛИСТИЧЕСКИХРЕСПУБЛИК 81 19) (11) 51)5 б 01 Ч 1/00 ОПИСАНИЕ ИЗОБРЕТЕН ский институт прием гулярности струк- сеяния оптическос метеопарамет- ССР. Физика ат- О 7, 1985, с. 720 -ЕНИЯ СТРУК ОД НОР ОДНОС(54) СГ 1 ОСОБ ОБНАРУ ТУРЫ ОПТИЧЕСКИХ Н ТЕЙ АТМОСФЕРЫ (57) Изобретение относитс тичсских характеристик не позволяет повысить наде ния элементов периодичн неоднородностей. В атмос я к измерению оподнородных сред и жность обнаружеости в структуре феру посылают имов о брасыва ютсусловие те, для которых выпол(С +2 СО) 7 Сгде игнала в локальном миежду максимумами С 1 ве-,Ь рас- епа- ями циент обеть (Р)ого перепловленносигнала;уровеньинтервалуции сигна спечивающииоответствияада (ЛС) пе 1му флхктчац коэффи роятно считанду обу самого среди и ный к кретиза шумов, отнвременной ГОСУДАРСТВЕННЫЙ КОМИТЕТПО ИЗОБРЕТЕНИЯМ И ОТКРЫТИПРИ ГКНТ СССР К АВТОРСКОМУ СВИДЕТЕЛЬСТВУ(56) Полканов Ю. А. О ретуры неоднородностей расго излучения и их связирами. - Известия АН Смосферы и океана, т. 21, М726. Изобретение относится к измерению опгических характеристик неоднородных сред и может применяться для целей метеорологии, физики плазмы, оптики атмосферы и океана, лазерной локации в условиях замутненной среды.Цель изобретения - повышение надежности выявления периодичности оптических неоднородностей среды по временной структуре сигнала обратного рассеяния в условиях воздействия шумов путем нелинейной фильтрации и анализа характера выделенной структуры неоднородностей в зависимости от растующего нелинейного порога ампл итуд ной дискрими нации.Поставленная цель достигается тем, что из последовательности локальных максимупульсы света малои длительности, регистрируют сигналы обратного рассеяния, усиливают их пропорционально квадрату текущего времени, компенсируют влияние фоновой засветки. Последовательность отсчетов фильтруют, исключая локальные максимумы, не превышающие порогового значения, которое зависит от значения принятого сигнала и уровня шумов. Исключают также наименьшие локальные минимумы, соседствующие с исключенными максимумами и экстремумами перед наибольшим и после наименьшего максимумов. Задают начальный уровень дискриминации. Последовательно увеличивая уровень, дискриминируют последовательность отсчетов. Для каждого уровня дискриминации формируют распределение числа интервалов между экстремумами от величины а интервалов. Характерные размеры неоднородностей, устойчивость, число и контраст элементов структуры неоднородностей определяют, анализируя зависимость распределений от уровня дискриминации. 1 ил.+С+С)С находят в последовательности локальных минимумов лары, соседствующие с отброшенными максимумами и отбрасывают меньший минимум пары, формируют новые последовательности оставшихся локальных максимумов, минимумов, определяют наличие максимального и минимального максимумов (Смакс, Счиц), находят их положение в после+довательности, отбрасывают все экстремумы+перед Сяке и после С 1 цн, задают начальный уровень дискриминации локальных мак- симумов формируют новую модельную последовательность отсчетов (С = 1, С= - 1), на ее основе формируют распределение полученного числа интервалов между максимумами, минимумами определенной протяженности в зависимости от величины интервала, определяют величину, положение и число относительных максимумов (М) в полученном распределении, увеличивают уровень дискриминации на заданный шаг (ЛЙ), повторяют действия по определению величины, положения и числа относительных максимумов, запоминают эти числа, определяют перепад уровней дискриминации, соответствующий переходу к числу относительных максимумов распределения менше на одно, по их максимальному числу (В) судят о числе характерных элементов структуры неоднородностей по отношению выявленных перепадов уровней дискриминации к уровню фильтрации (Со+АС) судят о контрасте соответствующего элемента структуры относительно фоновых неоднородностей, по отношению соседних перепадов уровней судят об относительном контрасте элементов структуры неоднородностей. На чертеже приведена блок-схема устройства, реализу ющего предлагаемый способ.Схема включает пульт 1 управления, импульсный источник света (ОКГ) 2, фото- приемник 3, усилитель 4, накопитель 5, амплитудно-цифровой преобразователь (АЦП) 6, блоки памяти 7, 8, 17, блок 9 вычитания фона, блок 10 выделенных локальных экстремумов, блок 11 задания уровня дискриминации, блок 12 нелинейной фильтрации и дискриминации локальных максимумов, блок 13 формирования распределения числа интервалов от протяженности интервала, блок 14 выделения локальных максимумов распределения, блок 15 определения числа локальных максимумов распределения (типов интервалов), блок 16 соотнесения типа интервала уровню дискриминации, блок 18 определения максимального и минимального числа типов интервалов, вычислитель 19 параметров обнаруженной структуры. С пульта 1 в накопитель 5 заносят интервал времени измерений,5 10 15 20 25 30 35 40 45 50 55 шаг дискретизации и число циклов, по которому проводят накопление приходягцего сигнала. После этого с пульта 1 подают команду на излучение импульса света источником 2. Одновременно подают команду на усилитель 4 для усиления приходящего сигнала пропорционально квадрату текущего времени, начиная с момента излучения импульса. Кроме того, с пульта 1 задают число шагов дискретизации, запоминаемое в блоке 7, и последовательность предварительно заданных команд обработки сигнала. Приходящий сигнал с фотоприемника 3 поступает на усилитель 4 и после усиления накапливается в блоке 5, за требуемый интервал времени по заданному числу измерений, с заданным шагом дискретизации. В блоке 6 (АЦП) он преобразуется в цифровую форму и запоминается в блоке 8. Причем в блоке 8 запоминаются как фоновые отсчеты, полученные в цикле предварительных измерений (при отключении источника излучения 2), так и сигнальные отсчеты полученные в рабочих циклах измерений. В блоке 9 производится вычитание уровня фона соотнесенного к шагу дискретизации сигнала (по данным из блока 7) из полученных значений сигнала. Полученные рабочие значения сигнала обратного рассеяния рабочая последовательность вновь запоминаются в блоке 8.После завершения измерений по команде с пульта 1 в блоке 1 О производится сравнение отсчетов сигнала по три и в случае, если для них выполняются условия 1(2)3 или 1)2(3, отсчет 2 принимается за локальный максимум (минимум) и его положение (номер) запоминается отдельно для максимумов и минимумов. После этого отсчеты сигнала и номера локальных экстремумов с блоков 8 и 10 поступают на блок 11. В зависимости от выставленного уровня фильтрации (Со), выдаваемого пультом 1, здесь происходит фильтрация локальных максимумов, когда значения большие ЛС сохраняются. Если в результате фильтрации отфильтровывается локальный максимум, то из двух соседствующих минимумов за действительный принимается наибольший, т. е. в блоке 11 производится дополнительная операция (схемой сравнения и выбора) .Кроме того, из выявленных локальных максимумов здесь выделяется наибольший и наименьший (Серакс и См") и определяется их номер в последовательности экстремумов, Эти номера соответственно принимаются за первый и последний, делается вывод о возможности анализа структуры и задается начальный уровень амплитудной дискриминации, который запоминается в блоке 11 и выделяется только по команде с пульта 1. В отсутствие этой команды на блок 12 поступает команда на фильтрацию выбросов сигнала, обусловленных средой от выбросов сигнала, обусловленных его природой, на ос 159781610 20 25 30 35 40 45 50 5 б 5новании данных об уровне шумов (Сп), поступающих из блока 8 памяти. На выходе блока 12 фильтрации имеем набор локальных максимумов, минимумов, разделенных определенными интервалами, число которых меньше их исходного числа (без дискриминации). В блоке 13 происходит сортировка этих максимумов и минимумов в две независимые последовательности, генерация нормированной последовательности, отбор интервалов между импульсами нормированной амплитуды (+1 или - 1) их сортировка по величине и суммирование числа интервалов определенной протяженности. В результате по команде с пульта 1 получаем на выходе блока 13 распределение числа интервалов в зависимости от протяженности интервала (спектр интервалов). Следующий этап анализа - анализ зависимости спектра интервалов от уровня дискретизации. Исходный спектр и отфильтрованный сигнал - в блоке 12. С пульта 1 и блока 11 на входы блока 12 поступают команды на выставление уровня дискриминации, который растет от цикла к циклу анализа с заданным шагом (Лдв).В соответствии с заданным уровнем дискриминации из отфильтрованной последовательности локальных максимумов, поступающей с блока 13, начинают исключаться локальные максимумы и минимумы определенных интервалов, начиная с минимального интервала. Причем с ростом уровня с 1 из последовательности локальных интервалов сначала исчезает экстремум, расположенный+у глобального максимума (С- ), затем такой же следующий и т. д. пока не исчезнет экстремум у глобального максимума (С"). Происходит прореживание локальных экстремумов сигнала, граница которого перемещается В направлении Смаке -Счин. Причем исчезновение локальных экстремумов идет с оп редел ен ным интер валом ( последова тел ьно) и уменьшение их числа на одно соответствует определенному приращению уровня дискриминации Лд. После фильтрации определенного типа локальных экстремумов, соответствую щи х существованию на чального элемента периодической структуры, существенный рост уровня дискриминации (Лд Яо) не приводит к изменению общей картины в последовательности следования локальных экстремумов. С некоторого значения д опять начинается селекция экстремумов и их прореживание. Это вызывает появление нового типа интервалов между оставшимися экстремумами, число которых растет с ростом уровня дискриминации. Причем, как и в начальной последовательности, уменьшение локального экстремума на один соответствует приращению Лд, т. е. к появлению еще одного интервала нового типа. С ростом с( число этих интервалов растет. Они появляются последовательно в1- +направлении Серакс -+ Смин, Максимальное число этих интервалов равно отношению протяженности всего анализируемого участка на протяженность такого интервала. Если есть еще периодические компоненты структуры (последующие ее элементы), то с ростом картина прореживания экстремумов сигнала и появления интервалов все большей протяженности повторяется. В случае отсутствия локальных экстремумов блок 12 работает как обычный ограничитель по амплитуде.Следующий этап - анализ динамики спектра интервалов структуры неоднородностей сигнала в зависимости от выставляемого порога дискриминации. Он предусматривает выделение интервалов определенной протяженности преобладающих по частоте их встречи, т. е. поиск относительных максимумов в распределении числа интервалов в зависимости от его величины.Распределение интервалов (спектр интервалов) формируется блоком 13 на основе полученных после фильтрации, дискриминации и нормировки данных, В блоке 14 анализируется это распределение на предмет выявления относительных максимумов по методике поиска локальных максимумов блоком 10. Данные о величине положения каждого относительного максимума поступают на блок 15, где число относительных максимумов суммируется и запоминается. Эти данные поступают на вход блока 16 одновременно с набором взятых уровней дискриминации от начальнсо до текущего (с блока 11). Здесь сравнивается номер уровня и номер интервала и в случае их совпадения в блоке 17 запоминается пара значений: уровень дискриминации н интервал выявленного типа. Рост уровня дискриминации при периодичности структуры неоднородностей приводит к росту числа типов интервалов и их последующему уменьшению до одиночного. Текущее число типов интервалов для определенного уровня д с блока 15 и пары значений уровень - дискриминации - интервал выявленного типа с блока 17 поступают на входы блока 18 и запоминаются. Число типов интервалов сравнивается с предшествующим числом для меньшего уровня дискриминации. При его неизменности или росте с блока 18 подается команда на блок 11 о выставлении следующего уровня дискриминации и повторении цикла анализа. При уменьшении типов интервалов до одного в блоке 18 определяется наибольшее число типов интервалов и запоминается. Одновременно анализируются пары уровень дискриминации - + интервал выявленного типа для всего диапазона уровней дискриминации и определяются перепады уровней Л при изменении числа типов интервалов на единицу. Эти данные поступают на вычислитель 1 одновременно с данными блока 17 и по максимальному числу типов интервалов уста 1597816навливается число элементов периодичностиструктуры неоднородностей.Каждому элементу структуры соотносится перепад уровней дискриминации Л,. Из перепадов вычитается предварительно запоминаемый уровень, соответствующий фону, и полученные значения (Л - Л,) являются характеристикой контрастности элементов структуры относительно фоновых неоднородностей. Здесь же вычисляются отношения полученных перепадов (Л - Л) / (Л - Л), которые принимаются за характеристики относительного контраста элементов структуры неоднородностей.Формула изобретенияСпособ обнаружения структуры оптических неоднородностей атмосферы путем посылки импульсов света малой длительности по набору трасс зондирования, приема и преобразования фотоприемником рассеянного в обратном направлении света в электрические сигналы, регистрации дискретных значений приходящего сигнала, их накопления и усиления пропорционально квадрату текущего времени, разбиения последовательности отсчетов сигнала на участки двух типов, ограниченные парами локальных максимумов и минимумов соответственно, формирования модельной последовательности отсчетов путем соотнесения нормированной амплитуды, равной единице, значениям локальных максимумов, отрицательного значения нормированной амплитуды, равной единице, значениям локальных минимумов и исключения промежуточных сигналов, определения по модельной последовательности интервалов между локальными максимумами, минимумами, максимумами и минимумами и определения по средним величинам этих интервалов для всего набора трасс зондирования, характерных размеров неоднородностей, а по среднеквадратичным отклонениям устойчивости структуры неоднородностей, отличающийся тем, что, с целью повышения надежности, измеряют величины сигналов фона, нормируют их на интервал времени дискретизации и исключают из дискретныхуровень дискриминации;значение наименьшего из неотфильтрованных локальных максимумов;значение соответствующего локального минимума. отсчетов приходящих сигналов, фильтруют последовательность отсчетов сигналов путем последовательного исключения локальных максимумов, не превышающих порога, меньших из каждых двух локальных минимумов, соседствующих с исключенными максимумами, и всех экстремумов перед наибольшим и после наименьшего максимумов, задают начальный уровень дискриминации локальных максимумов, последовательно уве личивая уровень дискриминации, формируютряд последовательностей отсчетов сигналов и ряд соответствующих модельных последовательностей, для каждой из которых формируют распределение числа интервалов между минимумами, максимумами от величины интервалов, а в полученных распределениях определяют величину, положение и число относительных максимумов, определяют перепады уровней дискриминации, соответствующие переходу к числу относительных 20 максимумов меньше на одно, причем порогфильтрации, определяют по формуле+ (С +2 СфСгде С+ - минимально допустимая величиналокального максимума;С - значение сигнала в локальном минимуме между максимумами каждой пары;Со - средний уровень шумов, нормиро ванный на интервале дискретизации сигнала;а - коэффициент, обеспечивающий требуемую вероятность соответствия максимума перепаду, обусловленному флуктуациями самого сиг нала,а начальный уровень дискриминации определяют по формуле- С+ 1/21597816 ВНИИПИ Государственного комитета по изобретениям и открытиям при ГКНТ СССР 113035, Москва, Ж - 35, Раушская наб., д. 4/5 Производственно-издательский комбинат Патент, г. Ужгород, ул. Гагарина, 1 ОРедактор С. ПатрушеваЗаказ 3053 Соста Техред А. Тираж 34 тель С. Сравцук ратов Корректор А. Обруца Подписное