Способ нейроподобного снижения размерности оптических спектров Российский патент 2017 года по МПК G01J3/28 G06N3/02 

Область техники, к которой относится изобретение

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

Предшествующий уровень техники

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

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

Другим важным направлением защиты промышленных товаров, в частности товаров, основанных на различных жидких средах, является спектральная компьютерная квалиметрия жидких сред, позволяющая проводить идентификацию, распознавание и оценивание параметров продукции пищевой и нефтехимической промышленностей с помощью методов оптической спектрометрии в широком диапазоне длин электромагнитных волн (от ультрафиолетовых до инфракрасных) [3].

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

Известен способ снижения размерности сигналов, основанный на модели нервной клетки по МакКаллоку-Питсу [МакКалок Дж., Питтс У. Логические исчисления идей, относящихся к нервной деятельности. // Автоматы. М.: ил, 1956].

Известный способ снижения размерности сигналов, состоящих из N компонент, заключается в мультиплексировании каждой компоненты х_i (i=1, 2, …, N) сигналов на М (М<N) каналов, увеличении/ уменьшении в каждом канале значений x_i, на весовой коэффициент w_ji (j=1, 2, …, М), последующего пространственного накопления в каждом канале всех результирующих сигналов w_jix_i,. Итоговая размерность М исходных данных определяется числом возможных редуцированных сигналов u_j=Σw_jix_i (j=1, 2, …, М).

Недостаток данного способа заключается в том, что само по себе описанное выше снижение размерности сигналов не имеет сколь-нибудь значимой интерпретации. Значимая интерпретация появляется лишь при дополнительном введении нелинейной функции активации y_j=ƒ(Σw_jix_i), где j=1, 2, …, М). Однако нелинейная функция активации ƒ и последующая настройка весов увеличения/уменьшения на ее основе целесообразна лишь для задач классификации, распознавания и прогнозирования данных. В этих случаях выходные величины y_j принимают значения 0 и 1, что не соответствует задачам снижения размерности сигналов с сохранение их топологических свойств, например, взаимных соотношений их сходства/ различия.

Известен способ снижения размерности сигналов, основанный на модели нервной клетки по Видроу Б. и Хоффу М. [Widrow В., Hoff М. Е. Jr., Adaptive switching circuits, Western Conf. Rec, IRE, 1960, cz. 4, s. 94-104].

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

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

Наиболее близким к заявляемому из известных по своему назначению, технической сущности и достигаемому результату является выбранный в качестве прототипа способ снижения размерности оптических спектров для адаптивного управления плазменным процессом на основе использования нейронных сетей, формирующих управляющие параметры по оптическому спектру, генерируемому лазерной плазмой [Chi Yung Fu. Closed loop adaptive control of spectrum-producing step using neural networks. US 5841651 A. 1998].

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

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

Раскрытие изобретения

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

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

Целью изобретения является расширение функциональных возможностей способа путем ортогонального отображения оптических спектров в пространство меньшей размерности с сохранением топологических свойств, связанных с соотношениями сходства/различия спектров как для их наглядной визуализации, так и дальнейшей автоматической обработки (идентификации и распознавания, оценивания параметров источников оптических спектров [3]).

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

Краткое описание чертежей

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

Фиг. 1 иллюстрирует реализацию предлагаемого способа при мультиплексировании компонент оптических спектров на 3 канала;

Фиг. 2 иллюстрирует люминесцентные спектры перед снижением их размерности;

Фиг. 3 иллюстрирует редуцированные трехмерные образы люминесцентных спектров;

Фиг. 4 иллюстрирует инфракрасные спектры перед снижением их размерности;

Фиг. 5 иллюстрирует редуцированные трехмерные образы инфракрасных спектров;

Фиг. 6 иллюстрирует множество редуцированных двухмерных образов инфракрасных спектров.

Лучшие варианты осуществления изобретения

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

Согласно предлагаемому способу снижение размерности оптических спектров производится для их совокупности - множества {S}_k (k=1, 2, …, К), состоящего из K элементов S_k=(S^k₁, S^k₂, …, S^k_N)^T, каждый из которых включает дискретные компоненты с амплитудами S^k_n(n=1, 2, …, N).

Каждая компонента S^k_n любого к-го оптического спектра мультиплексируется на несколько каналов, число которых значительно меньше числа компонент спектров (на Фиг. 1 показано три канала X, У, Z).

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

В качестве примера, вариации регистрируют в виде величин Δ_n=|S^max_n-S^min_n|, определяемых максимальными S^max_n и минимальными S^min_n значениями разброса значений амплитуд всех компонент (n=1, 2, …, N).

Далее, все вариации Δ_n (n=1, 2, …, N) сравнивают с первоначально назначенным порогом Δ_1thr, и при превышении этого порога формируют путем первой режекторной фильтрации отсчеты опорного дискретного сигнала W_X=(w_X1, w_X2, …, w_XN)^T первого канала в виде единичных орт:

По этим отсчетам настраивают величины коэффициентов w_Xn усиления/ослабления значений амплитуд компонент оптических спектров в первом канале X.

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

Затем, во втором канале Y, например с помощью Гильберт-фильтрации, формируют первый вспомогательный дискретный сигнал W₁=(w₁₁, W₁₂, …, w_1N)^T:

где h_k - импульсная характеристика Гильберта (ИХГ) [4]:

Далее регистрируют абсолютные значения |w_1n| отсчетов первого вспомогательного дискретного сигнала (n=1, 2, …, N) и формируют путем второй режекторной фильтрации отсчеты второго вспомогательного дискретного сигнала W₂=(w₂₁, w₂₂, …, w_2N)^T в виде единичных орт:

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

Таким образом, если первая режекторная фильтрация выделяет компоненты, ответственные за различие, то вторая - за сходство оптических спектров.

Отсчеты опорного дискретного сигнала W_Y=(w_Y1, w_Y2, …, w_YN)^T второго канала формируют на основе процесса ортогонализации Грамма-Шмидта:

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

В последующих каналах формируют соответствующие опорные дискретные сигналы, ортогональные всем опорным сигналам предыдущих каналов. Так, например, для третьего канала Z третий опорный дискретный сигнал W_Z=(w_Z1, w_Z2, …, w_ZN)^T формируют как:

где W - любой начальный дискретный сигнал, например, W=(1, 2, 3, …, N).

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

Как показано на Фиг. 1, редуцированные образы всех K оптических спектров формируют в каждом канале с помощью пространственного накопления соответствующих отфильтрованных компонент:

При описанной процедуре снижения размерности K оптических спектров каждому спектру - многомерному элементу S_k будет соответствовать свой редуцированный трехмерный образ X_k, Y_k, Z_k (k=1, 2, …, K).

Для целей наглядной визуализации спектров путем увеличения различия редуцированных образов итерационно изменяют пороги Δ_1thr и Δ_2thr в первых двух каналах, увеличивая какой либо критерий различия (например, минимум минимумов расстояний между точкам X_k, Y_k, Z_k или множественный дискриминантный критерий Фишера [5]).

Приведем примеры применения рассмотренного способа для снижения размерности люминесцентных инфракрасных спектров.

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

На Фиг. 4 приведены примеры четырех различных инфракрасных спектров автомобильных бензинов (с октановыми числами 82, 92, 95, 98), а на Фиг. 5 - их редуцированные трехмерные образы.

На Фиг. 6 приведен пример множества редуцированных двухмерных образов инфракрасных спектров для разных классов бензинов (80_нормаль, 92_регуляр, 95_премиум, 98_супер), с 5-ю представителями в каждом классе.

Для наглядной визуализации спектров итерационное изменением порогов Δ_1thr и Δ_2thr в двух каналах было проведено на основании оптимизации множественного дискриминантного критерия Фишера («сумма квадратов расстояний между классами»/«сумма дисперсий классов» [5]).

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

Промышленная применимость

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

Источники информации

1. Maximovsky S., et all. Shattuckite synthesis and the pattern formation by the scanning laser beam. Tecnol. Metal. Mater. Miner. So’о Poulo. V. 13. N. 3. 2016. P. 248-251.

2. Краснов A.E., Максимовский C.H., Смирнов B.M., Ставцев А.Ю. Метрологическое обеспечение защиты ценных бумаг, документов и товаров на основе нанотехнологий. Технологии XXI века в легкой промышленности (электронное научное издание). №7. Часть I. Раздел №4. 2013.Статья №6.

3. Краснов А.Е. и др. Основы спектральной компьютерной квалиметрии жидких сред / Под ред. д.ф.-м.н., проф. Краснова А.Е. - М.: ИД «Юриспруденция», 2006. - 264 с.

4. Klingspor М. Hilbert transform: Mathematical theory and applications to signal processing. LiTH - MAT - EX, 2015. - 76 p.

5. Дуда P., Харт П. Распознавание образов и анализ сцен / Перевод с англ. - М.: Мир, 1976. - 511 с.

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

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

2. Способ нейроподобного снижения размерности оптических спектров по п. 1, отличающийся тем, что вариации компонент оптических спектров регистрируют в виде величин Δ_n=|S^max_n-S^min_n|, определяемых максимальными S^max_n и минимальными S^min_n значениями разброса значений амплитуд всех компонент (n=1, 2, …, N), а отсчеты первого дискретного опорного сигнала W_X=(w_X1, w_X2, …, w_XN)^T формируют путем первой режекторной фильтрации с помощью первого порога Δ_1thr в виде единичных орт:

3. Способ нейроподобного снижения размерности оптических спектров по п. 2, отличающийся тем, что с помощью Гильберт-фильтрации формируют первый вспомогательный дискретный сигнал W₁=(w₁₁, w₁₂ …, w_1N)^T:

где - импульсная характеристика фильтра Гильберта, регистрируют его абсолютные значения |w_1n|, формируют путем второй режекторной фильтрации с помощью второго порога Δ_2thr отсчеты второго вспомогательного дискретного сигнала W₂=(w₂₁, w₂₂, …, w_2N)^T в виде единичных орт:

а отсчеты опорных дискретных сигналов W_Y=(w_Y1, w_Y2, …, w_YN)^T и W_Z=(w_Z1, w_Z2, …, w_ZN)^T второго и третьего каналов формируют на основе процесса ортогонализации Грамма-Шмидта:

где W- любой начальный дискретный сигнал, например, W=(1, 2, 3, …, N).