EA043532B1 - Способ определения спектральных показателей приведенного рассеяния и поглощения светового излучения мелкодисперсными материалами и устройство для его осуществления - Google Patents

Способ определения спектральных показателей приведенного рассеяния и поглощения светового излучения мелкодисперсными материалами и устройство для его осуществления Download PDF

Info

Publication number
EA043532B1
EA043532B1 EA202292115 EA043532B1 EA 043532 B1 EA043532 B1 EA 043532B1 EA 202292115 EA202292115 EA 202292115 EA 043532 B1 EA043532 B1 EA 043532B1
Authority
EA
Eurasian Patent Office
Prior art keywords
spectral
radiation
dependences
optical
cor
Prior art date
Application number
EA202292115
Other languages
English (en)
Inventor
Владимир Александрович Фираго
Константин Игоревич Шулико
Original Assignee
Белорусский Государственный Университет (Бгу)
Filing date
Publication date
Application filed by Белорусский Государственный Университет (Бгу) filed Critical Белорусский Государственный Университет (Бгу)
Publication of EA043532B1 publication Critical patent/EA043532B1/ru

Links

Description

Изобретение относится к области спектроскопии диффузного отражения с пространственным разрешением и может использоваться в устройствах неразрушающего контроля оптических характеристик мелкодисперсных рассеивающих материалов. Эти характеристики используют для оценки отклонения свойств и состава контролируемых материалов от заданных или эталонных. Основными оптическими характеристиками мелкодисперсных рассеивающих сред являются спектральный приведенный или редуцированный показатель рассеяния μ's(λ) и спектральный показатель поглощения μa(λ). Из них μа(λ) имеет в основном резонансный характер, т.е. μa увеличивается в резонансных полосах поглощения молекул рассеивающей среды, которые для конденсированных сред достаточно широкие. Рассеяние же светового излучения определяется неоднородностью среды и зависит от ее микрофизических параметров, в том числе от показателя преломления дисперсионной среды nt и показателя преломления нерастворимых включений в нее или дисперсной фазы ndp. Поэтому размеры включений и отношение ndp/nt влияют на интенсивность и спектр рассеяния излучения. При этом спектральный приведенный показатель рассеяния или его спектральная зависимость μ's(λ) обычно растет при уменьшении длины волны λ. Поэтому спектральную зависимость p's в основном описывают степенной зависимостью μ's(λ)®а(λ/λ0)b, где а с размерностью см-1 определяет значение μ's0) на длине волны λ0, а b - наклон зависимости μs(λ) [1].
Известны способы [2, 3] вычисления зависимостей коэффициента отражения R(μ's, μa, р) или яркости L(μ's, μa, р) диффузно отраженного излучения полубесконечной рассеивающей средой от приведенного показателя рассеяния μ's0) и показателя поглощения μа0) на длине волны λ0, а также расстояния р между оптоволоконным зондом, через который излучение вводится в среду, и приемным зондом. Их называют диффузионными приближениями. Причем, предполагают, что оси этих оптоволоконных зондов расположены перпендикулярно поверхности исследуемого материала. Для определения искомых значений приведенного показателя рассеяния μ' s0) и показателя поглощения μ a0) на одной длине волны λ0, используют способ, основанный на составлении функции невязки ί[μ' s0), μ a0), р] и применении подгонки рассчитываемых зависимостей R[μ' s0), μ a0), р] или L[μ' s0), μ a0), р] к экспериментально полученным зависимостям Rm0, р) или Lm0, р) при варьировании значениями μ' s0) и μ a0) [4, 5]. Подгонку осуществляют путем минимизации функции невязки ί[μ' s0), μ a0), р] методом наименьших квадратов. В результате минимизации ί[μ' s0), μ a0), р] определяют искомые значения приведенного показателя рассеяния μ' s0) и показателя поглощения μ a0) на длине волны λ0.
Недостаток этого способа заключается в том, что сложная нелинейная связь между экспериментально получаемыми зависимостями L(λ0, рД или R(λ0, рк) и значениями μ's0), μa0), а также влияние состояния поверхности контролируемой рассеивающей среды, может приводить при минимизации функции невязки к неоднозначности получаемых решений μ' s0) и μ a0), так как эффективно разделить их влияние на экспериментально получаемые зависимости Rm0, р) или Lm0, р) не удается [6].
При наличии математической модели, адекватно описывающей оптические характеристики среды, отмеченные выше различия в поведении спектрального приведенного показателя рассеяния μ's(λ) и спектрального показателя поглощения μа(λ) позволяют при проведении подгонки параметров модели к экспериментально получаемым спектральным зависимостям Rm(λ, р) или Lm(λ, р) в широком диапазоне длин волн λ обеспечивать лучшую устойчивость решений, т.е. улучшить точность определения спектральных зависимостей μ' s(λ) и μ а(λ) [7-9]. Причем для устранения влияния на результаты определения μ' s(λ) и μ ^λ) спектральной чувствительности 8(λ) используемой аппаратуры в способе [8] при расчетах вместо экспериментально получаемых зависимостей R(λ0, ρk) используют их нормированный спектрально-пространственный профиль коэффициента диффузного отражения rnm(λ, ρk, c)=R(λ, ρk)/R(λ, ре), где рс - минимальное из используемых расстояний ρk. При выполнении операции нормировки R(λ, ρk)/R(λ, ре), спектральная чувствительность аппаратуры S(λ), которая неявно содержится в виде множителей в числителе и знаменателе отношения, сокращается. Применение этого профиля rnm(λ, ρk, c) позволяет исключить из действий, осуществляемых с экспериментально полученными зависимостями R(λ, ρk), операцию учета спектральной чувствительности применяемой аппаратуры S(λ), что укорачивает и упрощает последовательность проводимых действий. При проведении операции подгонки теоретические варьируемые зависимости rnv[μ' s(λ), μ a(λ), ρk, c] могут рассчитывать с помощью моделей диффузионного приближения [2-5] или заранее формировать путем моделирования процессов диффузного отражения излучения рассеивающими средами методом Монте-Карло [8].
Наиболее близким к заявляемому изобретению является способ [9] определения спектральных показателей приведенного рассеяния μ's(λ) и поглощения μа(λ) излучения мелкодисперсными материалами путем ввода широкополосного светового излучения в малый участок поверхности мелкодисперсного материала с известным отношением показателей преломления материала nt и среды над ним nm, т.е. nrel=nt/nm, измерения интенсивности диффузно отраженного излучения I(λ, ρk) при разных расстояниях ρk между точкой ввода излучения и регистрации отраженного излучения, определения зависимости его спектрального коэффициента диффузного отражения Rm(λ, ρk) и последующего численного определения искомых спектральных показателей μ' s(λ) и μ а(λ) при решении системы нелинейных уравнений в виде
- 1 043532 функции невязки, содержащей квадраты разностей измеренных Rm(λ, pk) и теоретических рассчитываемых с помощью диффузионного приближения [2, 3, 5] зависимостей Rv[μ' s(λ), μ a(λ), nrel, pk, aav], которую минимизируют методом наименьших квадратов при варьировании усредненным коэффициентом пропорциональности αav=Rm(λ, μk)/Rv(λ, μk) и параметрами а и b аппроксимации μ' s(λ) гиперболической функцией вида μ'*s(λ)=аλ’b.
В отличие от способа [4], использующего при определении показателей μ's и μ а только узкий участок спектра (например, при использовании в качестве источника излучения лазера с длиной волны генерации λ0, тогда μ' s= μ' s0), a μ aа0)), устойчивость получаемых решений улучшается.
Недостаток этого способа - отсутствие учета влияния состояния поверхности исследуемого образца мелкодисперсного материала на регистрируемые зависимости Rm(λ, pk, c), а также снижение вклада в общую сумму квадратов составляющих функции невязки с большими значениями pk, поскольку зависимости Rm(λ, pk, c) и Rv[μ' s(λ), μ а(λ), pk] быстро уменьшаются с ростом pk.
Известно устройство для измерения экспериментальных зависимостей коэффициента отражения R(λ, p) или яркости L(λ, p) мелкодисперсных материалов или рассеивающих сред, содержащее набор оптоволокон, которые подключают к многоканальному спектрометру [4, 9, 10, 11]. Многоканальные спектрометры используют вогнутые фокусирующие дифракционные решетки с плоским полем формируемого изображения и коррекцией аберраций типа IV [12, 13], принцип действия которых поясняет фиг. 1. В плоскости изображения, формируемого этой решеткой 2, размещают матрицу фотоприемников 3, что позволяет одновременно получать спектры от нескольких оптоволокон [9], размещаемых вдоль входной щели спектрометра 7, как показано на фиг. 1. По такой схеме выполнен многоканальный спектрометр S200-MF фирмы СоларЛС (Республика Беларусь) [14].
Недостатком таких устройств является ограниченное число используемых оптоволокон, а также влияние изгибов оптоволокон, которые неизбежно возникают в процессе подготовки к измерениям и их проведении, на значения спектральных потоков, поступающих на вход используемого спектрометра, что искажает получаемые спектральные зависимости.
Известно устройство для формирования гиперспектрального изображения [15, 16], содержащее объектив, входную щель, на которую объективом фокусируется изображение наблюдаемой поверхности, вогнутую фокусирующую дифракционную решетку с плоским полем формируемого изображения и коррекцией аберраций, матрицу фотоприемников с системой цифровой обработки получаемых изображений, на светочувствительную область которой дифракционная решетка фокусирует дифрагированные лучи, причем ориентация щели относительно матрицы фотоприемников выбрана так, что спектр формируется вдоль строк размещения светочувствительных элементов матрицы. Гиперспектральное изображение формируют путем смещения изображения, создаваемого объективом на входной щели, в направлении перпендикулярном ей. Это смещение обычно реализуют при линейном перемещении гиперспектральной видеокамеры в нужном направлении со скоростью, согласованной с необходимым спектральным разрешением устройства ц и требуемым временем экспозиции тэ, определяемом освещенностью элементов используемой матрицы фотоприемников дифрагированным решеткой излучением.
Недостаток этого устройства при его использовании в спектроскопии диффузного отражения с пространственным разрешением - ограниченные функциональные возможности, вследствие отсутствия источника излучения, вводимого в среду, и недостаточного динамического диапазона измеряемых спектральных зависимостей L(λ, p).
Наиболее близким техническим решением к заявляемому устройству является устройство [9], содержащее зонд из 10 оптоволокон, галогенную лампу, излучение которой фокусируют на торец первого оптоволокна зонда, служащего для освещения мелкодисперсного материала при непосредственном контакте всех торцов оптоволокон одной стороны зонда с поверхностью исследуемого мелкодисперсного материала, многоканальный спектрометр, в плоскости входной щели которого последовательно друг за другом вдоль щели располагают торцы 9 приемных оптоволокон второй стороны зонда, и компьютер, необходимый для управления работой многоканального спектрометра и проведения расчетов на основе использования диффузионного приближения и визуализации получаемых спектральных показателей μ'*s(λ) и μ*α(λ).
Недостатки этого устройства - необходимость перед измерениями с помощью эталонного рассеивающего образца контролировать коэффициенты пропускания излучения каждого приемного оптоволокна с последующим учетом их влияния на интенсивность регистрируемых сигналов, а также фиксированные значения расстояний pk между излучающим и приемными зондами.
Задачей изобретения является повышение точности способа определения спектральных показателей приведенного рассеяния μ's(λ) и поглощения μа(λ) мелкодисперсных рассеивающих сред, основанного на регистрации спектральной яркости диффузно отраженного излучения L(λ, pk) на разных расстояниях pk между точкой ввода зондирующего излучения в среду и точками регистрации, и создание компактного устройства для его осуществления, которое не имеет подвижных механических узлов и не требует использования набора оптоволокон, подключаемых к многоканальному спектрометру.
Сущность изобретения заключается в том, что при определении спектральных показателей приве- 2 043532 денного рассеяния μ's(λ) и поглощения μа(λ) светового излучения мелкодисперсными материалами используют способ, в котором вводят излучение широкополосного светового источника в малый участок поверхности мелкодисперсного материала с известным отношением показателей преломления материала nt и среды над ним nm, т.е. nrel=nt/nm, регистрируют зависимости интенсивности I(λi, pk) диффузно отраженного излучения от используемых длин волн λi и расстояний pk между точкой ввода и регистрации, рассчитывают нормированный спектрально-пространственный профиль интенсивности Inmi, pk, c)=I(λi, pk)/I(λi, pc) отраженного излучения, где pc - минимальное из используемых расстояний pk, а искомые спектральные зависимости μ' si) и μ ai) определяют при решении системы нелинейных уравнений в виде функции невязки
которую минимизируют методом наименьших квадратов, где kh - число используемых расстояний pk, kb - индекс минимального расстояния из набора расстояний pk, которые используют при расчете разностей искомых спектральных зависимостей μ аi, pk, c) и их усредненной по значениям pk зависимости (μ*а(λi)),
Inv[μ' s0), m λ, μ ai, pk, c), nrel] - варьируемый нормированный спектрально-пространственный профиль, который рассчитывают с использованием варьируемой степенной гиперболической зависимости , где μ'*s0) - варьируемое значение показателя приведенного рассеяния при длине волны λ0, m λ - варьируемый показатель степени, и диффузионного приближения, с помощью которого рассчитывают зависимость интенсивности диффузно отражаемого излучения от четырех пара метров I(μ'*s, μ*a, pk, c, nrel), а варьируемые спектральные зависимости μ*аi, pk, c) рассчитывают путем интерполяции на неравномерной сетке отсчетов семейства заранее рассчитываемых поверхностей Inr(μ' s, μ a, pk, c, nrel) = !(μ* s, μ a, pk, nrel)/I(μ’ s, μ a, pc, nrel) с учетом текущих значений варьируемой спектральной зависимости μ' si), причем варьируемую в процессе расчетов нормированную зависимость Inv[μ' s0), m λ, μ ai, pk, c), nrel] и заранее рассчитываемые плоскости Inr(μ' s, μ a, pk, c, nrel) делят на варьируемую корректирующую функцию вида
Лог (<ог > Р^е ) = 1 + (ДРА,С/Ру Г + (2 - )(ДРА,с/Ру ), (2) где Δpk, c = pk-pc, m*cor - варьируемый коэффициент, с помощью подбора которого учитывают влияние состояния поверхности мелкодисперсного материала на интенсивность диффузно отражаемого излу чения.
Предлагаемый способ может быть осуществлен с помощью устройства для определения спектральных показателей приведенного рассеяния μ's(λ) и поглощения μа(λ) светового излучения мелкодисперсными материалами, содержащее объектив, широкополосный источник светового излучения с микропроцессорным блоком его управления, которое вводится в материал непосредственно или с помощью оптоволоконного зонда, причем точка ввода излучения расположена в плоскости фокусировки объектива, щель, на которую объективом фокусируется исходное изображение поверхности материала в участке прилегающим к точке ввода излучения в материал, вогнутую фокусирующую дифракционную решетку с плоским полем формируемого изображения и коррекцией аберраций, матрицу фотоприемников с блоком цифровой регистрации получаемых изображений, на светочувствительную область которой дифракционная решетка фокусирует дифрагированные лучи, причем ориентация щели и дифракционной решетки относительно матрицы фотоприемников выбрана так, что спектр формируется вдоль строк матрицы, управляющий компьютер, который совместно с блоком цифровой регистрации получаемых изображений и микропроцессорным блоком управления широкополосным источником светового излучения осуществляет управление процессом регистрации и передачи семейства получаемых изображений с оптимальным выбором времени экспозиции тэ для каждого из кадров с целью обеспечения регистрации пространственной зависимости спектральной интенсивности диффузно отраженного излучения I(λi, p) в задаваемом диапазоне значений расстояний p от точки ввода излучения в материал до точки регистрации, а также обработку получаемых семейств спектральных зависимостей с оперативным расчетом искомых зависимостей μ' si) и μ αi).
Поставленная цель достигается путем ввода широкополосного излучения в малый участок поверхности исследуемого материала или среды, регистрации спектрально-пространственных зависимостей интенсивности I(λi, pk) диффузно отраженного материалом излучения, расчета нормированного спектрально-пространственного профиля интенсивности отраженного излучения Inmi, pk, c)=I(λi, pk)/I(λi, pc) и решения системы нелинейных уравнений в виде функции невязки (1), включающей две суммы. Первая
- 3 043532 из них содержит отношение варьируемого теоретически рассчитываемого с помощью диффузионного приближения спектрально-пространственного профиля Inv[p' s(λ), μ ai, pk, c), nrel]=Inv^' ^λ0), m λ, μ α(λ„ pk, c), nrel], где = M·' Л)(Ш) ' , к экспериментально определенному профилю Inmi, pk, c). Для учета влияния состояния поверхности мелкодисперсного материала на получаемые зависимости Inm(Xi, pk, c) моделируемую зависимость Inv[μ' s(λ), μ a(Xi, pk, c)], делят на корректирующую варьируемую функцию fcor(m cor, pk, c). Вторая сумма содержит разности искомых спектральных зависимостей μ аi, pk, c) и их усредненной по значениям μk зависимости (μ аi)). Причем варьируемые спектральные зависимости μ аi, pk, c) рассчитывают путем интерполяции на неравномерной сетке отсчетов семейства заранее рассчитанных поверхностей I(μ' s, μ a, pk, c, nrel) с учетом текущих значений варьируемых спектральных зависимостей μ' si), которые делят на значения варьируемой корректирующей функции fcor(m cor, pk, c). Техническая реализация заявляемого способа осуществляется при помощи устройства, содержащего источник светового излучения с микропроцессорным блоком его управления, которое вводится в исследуемый мелкодисперсный материал, объектив, щель, фокусирующую дифракционную решетку с плоским полем формируемого изображения и коррекцией аберраций, матрицу фотоприемников с блоком цифровой регистрации получаемых изображений, управляющий компьютер, который совместно с блоком цифровой регистрации получаемых изображений и микропроцессорным блоком управления широкополосным источником светового излучения обеспечивает регистрацию интенсивности диффузно отраженного излучения I(λi, p) в задаваемом диапазоне значений расстояний p от точки ввода излучения в материал до точки регистрации без необходимости применения подвижных оптико-механических узлов или набора оптоволокон с внешним многоканальным спектрометром.
Полезные свойства, появляющиеся у заявляемого технического решения:
улучшение устойчивости и точности получаемых решений;
отсутствие подвижных оптико-механических узлов или набора оптоволокон с многоканальным спектрометром, что уменьшает погрешности определения нормированных спектральнопространственных профилей интенсивности отраженного излучения Inm(Xi, pk, c).
Отличия заявляемого способа от известных решений:
учет влияния состояния поверхности исследуемого мелкодисперсного материала на варьируемые в процессе подгонки зависимости Inv[μ' s0), m λ, μ ai, pk, c), nrel] путем деления на значения варьируемой корректирующей функции (2);
использование в функции невязки (1) квадратичной суммы отношений {Inv[μ' s(λ), μ ai, pk, c), * nrel]/fcor(m cor, pk, c)}/Inmi, pk, c)-1, что позволяет с одинаковым весом учитывать при минимизации функции невязки зависимости, полученные для разных значений pk;
дополнение функции невязки (1) суммой квадратов разности искомых спектральных зависимостей μ аi, pk, c) и их усредненной по значениям pk зависимости (μ аi);
вычисление варьируемых значений μ аi, pk, c) с помощью двумерной интерполяции на неравномерной сетке заранее рассчитанных с использованием диффузионного приближения плоскостей Inr(μ' s, μ a, pk, c, nrei), деленных на текущие значения варьируемой корректирующей функции (2).
Заявляемое устройство содержит следующие признаки, отличные от известных решений:
обеспечение непосредственного ввода излучения в поверхность исследуемого образца мелкодисперсного материала с помощью широкополосного источника светового излучения, выполненного в виде сверхяркого светодиода с люминесцентным покрытием или галогенной лампы с оптоволокном, причем светодиод или выходной торец оптоволокна размещается в плоскости фокусировки объектива;
применение набора времен экспозиций Tm, обеспечивающего регистрацию гиперспектрального изображения щели в широком диапазоне изменений интенсивности I(Xi, pk, c), что позволяет проводить необходимую сшивку регистрируемых кадров и получать плавно изменяющиеся нормированные спектрально-пространственные профили Inm(Xi, pk, c)=I(Xi, pk)/I(Xi, pc) в необъходимом диапазоне значений p от pmin до pmax.
Сущность изобретения поясняется фиг. 1-14, где на фиг. 1 приведена схема, иллюстрирующая принцип действия многоканального спектрометра на основе матрицы фотоприемников и вогнутой фокусирующей дифракционной решетки с плоским полем формируемого изображения и коррекцией аберраций;
на фиг. 2 - схема устройства, осуществляющего предложенный способ определения спектральных зависимей приведенного рассеяния μ's(λ) и поглощения μа(λ);
на фиг. 3 - пример заранее рассчитываемых с использованием моделей диффузионного приближения [2, 3] нормированных плоскостей Im(^*s, μ*α, pk, c, nrei)=I(g’*s, μ*α, pk, nrei)/I(^*s, μ*α, pc, nrei) при pc=p1=2,3 мм;
на фиг. 4 - пример четырех экспериментально измеренных нормированных зависимостей Inm[μ' s(λ), μ а(λ), pk, c] для образца полиуретана при разных значениях pk и значении pc=2,3 мм, где связь между μ's и μ а представлена в параметрическом виде с целью иллюстрации, что каждая из полученных зависимо- 4 043532 стей расположена в пределах своей плоскости значений Inr(μ'*s, μ*α, pk, c, nrel);
на фиг. 5 - возможные изменения варьируемой спектральной зависимости ^(λ) = μν0)(λ0/λ) при разных значениях ш\ в диапазоне от -0,5 до 4, λ0= 500 нм и μ'*s0)=50 см-1;
на фиг. 6 - траектория минимума среднеквадратичного отклонения σ4(μ'°) функции невязки (1) при λχ=977,8 нм, fcor(m cor, pk, c)=1 и использовании модели диффузионного приближения [2] для образца полиуретана;
на фиг. 7 - траектория минимума среднеквадратичного отклонения σ% функции невязки (1) для образца белого полистирола при λχ=750 нм, fcor(m cor, pk, c)=1 и использовании модели диффузионного приближения [2], на фиг. 8 - зависимости Iχ, pk, c) полученные экспериментально для образца полиуретана при разных положении объектива слева (жирные сплошные линии) и справа (тонкие сплошные линии черного цвета) от точки ввода излучения в образец, а также подогнанные зависимости Ιην[μ' s0), μ λ, μ αχ, pk, c), nrel]/fcor(m cor, pk, c) при двух положениях объектива (штриховые линии) и использовании функции невязки (1);
на фиг. 9 - спектральные зависимости μ' s(λ) (а) и (μ а(λ)) (б) для образца полиуретана, полученные в результате подгонки, где штриховыми линиями представлены спектральные зависимости (μ aR(λ)) и (μ aL(λ)), рассчитанные для двух положений объектива;
на фиг. 10 - рассчитанные спектральные зависимости произведения μeff(λ)·pk для образца полиуретана, иллюстрирующие наличие участков спектра, в которых не выполняются условия применимости моделей диффузионных приближений [2, 3];
на фиг. 11 - зависимости значений корректирующей функции fcor от приращений расстояния pk, c=pk-pc, полученные в результате подгонки для образца полиуретана при двух положениях (слева fcorL и справа fcorR) объектива;
на фиг. 12 - зависимости Iχ, pk, c) полученные экспериментально для образца белого пищевого полистирола при разных положениях объектива слева (жирные сплошные линии) и справа (тонкие сплошные линии черного цвета) от точки ввода излучения в образец, а также подогнанные зависимости Ыр' s(λo), μ λ, μ αχ, pk, c), nrel]/fcor(m cor, pk, c) при двух положениях объектива (штриховые линии) и использовании функции невязки (1);
на фиг. 13 - спектральные зависимости μ' s(λ) (а) и (μ а(λ)) (б) для образца белого полистирола, полученные в результате подгонки, где штриховыми линиями представлены усредненные зависимости (μ aR(λ)} и (μ aL(λ)), рассчитанные для двух положений объектива;
на фиг. 14 - зависимости значений корректирующей функции fcor от приращений расстояния pk, c=pk-pc, полученные в результате подгонки для образца полистирола при двух положениях (слева fcorL, справа fcorR) объектива.
Устройство для определения спектральных показателей приведенного рассеяния μ's(λ) и поглощения μа(λ), функциональная схема которого показана на фиг. 2, имеет объектив 4, который сфокусирован на поверхность мелкодисперсного материала 5. Для ввода светового излучения в малый участок этого материала, расположенный в центре поля зрения объектива, служит широкополосный источник светового излучения. Он может быть выполнен в виде компактного сверхяркого светодиода с люминесцентным покрытием 6 или в виде оптоволоконного источника, выходной торец которого помещается на место светодиода 6, а на второй входной торец которого фокусируется излучение галогенной лампы или сверхяркого светодиода. Для устранения влияния фонового излучения и темновых токов используемой матрицы фотоприемников 7 применяют управление режимом работы (включение/выключение) этого широкополосного источника светового излучения с помощью микропроцессорного блока его управления 8 по командам управляющего компьютера 9. Это позволяет при выключении источника излучения или перекрытии светового потока от галогенной лампы с помощью управляемой шторки регистрировать фоновые сигналы вырабатываемые матрицей фотоприемников 7 и потом вычитать их из сигналов, регистрируемых при вводе светового излучения в рассеивающий материал.
При изотропном рассеянии вокруг точки ввода излучения формируется осесимметричный ореол диффузно отраженного излучения, показанный на фиг. 2 штриховыми изолиниями, интенсивность которого быстро спадает с ростом расстояния p. Объектив 4 формирует изображение этого ореола в плоскости размещения щели 10. Расположение объектива для схемы измерений, показанной на фиг. 2, выбрано таким, чтобы на щель попадало изображение ореола, находящегося справа от точки ввода излучения в материал.
Излучение, прошедшее через щель 10, попадает на вогнутую дифракционную решетку 11 с плоским полем формируемого изображения и коррекцией аберраций. Дифрагированные этой решеткой пучки излучения формируют на светочувствительной поверхности матрицы фотоприемников 7 гиперспектральное изображение щели. Причем, как показано на фиг. 2, расположение матрицы 7 выбрано таким, чтобы спектр формировался вдоль строк расположения элементов этой матрицы фотоприемников. Пара
- 5 043532 метры дифракционной решетки 11, расстояния от щели 10 до ее центра и от решетки до матрицы фотоприемников 7 выбирают такими, чтобы на светочувствительной поверхности матрицы фотоприемников формировался спектр в заданном диапазоне длин волн от коротковолновой Xsw до длинноволновой Xlw его границ. Для схемы измерений, приведенной на фиг. 2, направление спада интенсивности диффузно отраженного излучения и интенсивности потоков, регистрируемых светочувствительными элементами матрицы 7, наблюдается при возрастании ρ и ρ'. Для обеспечения регистрации гиперспектрального изображения щели в большом динамическом диапазоне изменения интенсивности дифрагированных потоков матрица фотоприемников 7 с блоком цифровой регистрации получаемых изображений 12 осуществляет последовательную регистрацию несколько кадров с нарастающим временем экспозиции Tm. Причем, для каждого значения Tm регистрируется и фоновые изображения при выключенном широкополосным источнике светового излучения на основе светодиода 6 или перекрытой шторке при использовании оптоволокна и источника на основе галогенной лампы. Значения Tm выбирают таким образом, чтобы обеспечить сшивку кадров и получение результирующего гиперспектрального изображения щели с плавно спадающей зависимостью Ι(λ^ pk) в большом динамическом диапазоне значений I. Получаемые кадры с гиперспектральными изображениями щели и фона пересылаются блоком 12 в управляющий компьютер 9. Там они накапливаются в оперативной памяти компьютера и могут сохраняться в архиве на жестком диске.
При обработке с помощью компьютера mh полученных кадров гиперспектральных изображений щели с разными временами экспозиции Tm вычитают фоновые изображения и формируют результирующую спектрально-пространственную плавно спадающую зависимость Ι(λ^ p'k). Потом с учетом масштабного коэффициента km, связывающего расстояния ρ и ρ', т.е. p=kmp', рассчитывают нормированный спектрально-пространственный профиль интенсивности Τ^(λι, pk, c)=I(X^ рД/ДХ^ pc) отраженного излучения, где рс - минимальное из используемых расстояний pk. Значение pc выбирают из критериев:
(йа + Ps) · Р 1 при использовании приближения [2, 4] или для приближения [3]. Эти критерии следуют из условий применимости используемого диффузного приближения, которое в аналитическом виде описывает зависимость коэффициента отражения R(μ's, μ a, pk, nrel) или интенсивности Ι(μ's, μ a, pk, nrel) диффузно отражаемого излучения образцом материала 5 от трех переменных и отношения nrel показателей преломления материала nt и среды над ним nm.
Аналитическое выражение для описания зависимости интенсивности I(p) диффузно отражаемого материалом излучения при использовании модели диффузионного приближения с одним дипольным источником, содержащим положительный и отрицательный источники излучения [2], имеет вид
I (ρ) = fV ’ Ιόπ2^ гг
L-iWi 1 zo + Ί·ζι ,(3) где D=l/[3(μ'sα)] - оптическая диффузионная длина;
ФS- мощность излучения источника 6;
η = 7р2 + Ц; г2 = /р2 +(го + *)2; ^^/(р^+Ря/з);
Zb=2A-D; А=(1+Ггь)/(1-Ггь);
rrh=0,668+0,0636nrel+0,710/nrel-1,440/n2 rel - коэффициент внутреннего отражения в заднюю полусферу излучения, попадающего на границу раздела из среды, полученный с учетом формул Френеля [5]; а эффективный показатель ослабления излучения средой μ^ определяется выражением Peff = y^/D = ^(Ра+цД.
Недостатки аналитического выражения (3) - занижение значений зависимости I(p) при Р < Weff - (μα+ Щ) и неполный учет влияния коэффициента преломления материала nt и состояния его поверхности [5, 17].
Модификация диффузионного приближения (3) путем использования второго вспомогательного дипольного источника [3] позволяет смоделировать более острый пик диффузно отраженного излучения в окрестности точки ввода излучения в ткань. Мощность излучения вспомогательного источника ФS* определяется мощностью главного источника OS, умноженной на два множителя * (4) где z*0=(1-η)(1/μ's)2/z0 - глубина размещения положительного вспомогательного источника, а η - коэффициент пропорциональности для вспомогательного источника зависящий от фактора анизотропии g рассеивающего материала и подгоночного параметра γ. Авторы этого приближения [3] протестировали его на ограниченных комбинациях оптических свойств рассеивающей полубесконечной среды и рекомендуют использовать значение γ=10 для расстояний между зондами p>0,1/μ's. Для изотропно рассеивающих сред фактор анизотропии g мал, поэтому значение η близко к нулю [3] и приближения [2] и [3] практически совпадают.
- 6 043532
Влияние границы раздела между материалом и воздухом учитывается так же, как и в приближении с одним дипольным источником [2] введением экстраполированной границы на высоте zb=2AD. Поэтому второй негативный вспомогательный источник располагается на высоте z 0+2zb. Поскольку интенсивность диффузно отражаемого излучения в модифицированной модели формируется при воздействии двух дипольных источников, она определяется суммой [3] где
Для формирования функции невязки (1) при используемом наборе значений pk заранее рассчитывают поверхности Ir(μ' s, μ а, μk, nrei) на заданной сетке значений μ's и μ a, которые необходимы для создания нормированных плоскостей Inr(μ' s, μ a, pk, c, nrel)=I(μ' s, μ a, pk, nrel)/I(μ' s, μ a, pc, nrel), используемых при вычислении теоретического нормированного спектрально-пространственного профиля Inv[μ' s0), μ λ, μ ax, pk, c), nrel], варьируемого в процессе минимизации (1). Пример трехмерного изображения четырех рассчитанных плоскостей Inr(μ' s, μ a, pk, c) приведен на фиг. 3. Монотонное изменение и отсутствие пересечения плоскостей Inr(μ' s, μ a, pk, c, nrel) значительно упрощает процесс подбора значений μ а и μ's при подгонке параметров используемой модели диффузного приближения. Действительно, каждой точке плоскости μ a, μ's будут соответствовать свои значения Inra, μ's, pk, c), как показано на фиг. 3. Это указывает, что при выполнении диффузионного приближения, каждая из получаемых при измерениях зависимостей Inmx, pk, c)=Inmax), μ'sx), pk, c, nrel], при параметрическом представлении μ' sа) будет в виде линии лежать в пределах соответствующей плоскости со значением pk, что иллюстрирует фиг. 4. Причем, влияние спектральной чувствительности S(λ) применяемой аппаратуры на результаты подгонки нормированных варьируемых зависимостей Invx, pk, c) к экспериментально измеряемым Inmx, pk, c), как уже указывалось выше, исключается.
Для вычисления варьируемых спектрально-пространственных зависимостей Inv^x, pk, c), входящих в (1), используют различие в поведении спектральных зависимостей μ'sx) и μаx). Значения μ'sx) обычно плавно возрастают при уменьшении длины волны, a μ^λχ) увеличивается в основном в полосах поглощения молекул, входящих в состав исследуемого мелкодисперсного материала, т.е. μax) имеет резонансный характер. Поэтому при формировании варьируемой спектральной зависимости μ' sx) используют гиперболическую степенную зависимость * μΓ(λ/) = μΓ(λο)·(λΟ/λ/Γ^ ПО) * где μ' s0) - варьируемое значение показателя приведенного рассеяния при длине волны λ0, а m λ варьируемый показатель степени гиперболы λ0χ. Как видно из фиг. 5, выражение (10), зависящее от трех параметров, достаточно гибко описывает возможные спектральные зависимости μ' sx).
Варьируемые спектральные зависимости μ ax, pk, c) при задаваемых в процессе подгонки параметрах зависимости μ' sx) рассчитывают, используя параметрическую связь μ a(μ's, pk, c)=F[μ' s, Inm(μ' s, pk, c)], путем двумерной интерполяции на неравномерной сетке отсчетов семейства заранее рассчитанных плоскостей Inr(μ' s, μ a, pk, c, nrel)=I(μ' s, μ a, pk, nrel)/I(μ' s, μ a, pc, nrel) с учетом текущих значений варьируемой спектральной зависимости μ' sx) и корректирующей функции fcor(m cor, pk, c). Правомерность такого подхода подтверждают зависимости на фиг. 3 и 4. Действительно, как видно из зависимостей на фиг. 3, по заданным значениям μ' sx, pk, c), pk, c и экспериментально измеренным значениям Inmx, pk, c), можно при выполнении диффузионного приближения определить соответствующие значения μ ax, pk, c). Это позволяет рассчитывать семейство варьируемых нормированных спектральных зависимостей Inv[μ' s0), μ λ, μ ax, pk, c), nrel], которые используют при минимизации функции невязки (1).
Поскольку существующие диффузионные приближения [2, 3] не в полной мере учитывают показа тель преломления исследуемого мелкодисперсного материала nt, а также влияние состояния его поверхности на результаты подгонки [5, 17], варьируемую в процессе расчетов нормированную зависимость Inv[μ' s0), μ λ, μ ax, pk, c), nrel] делят на эмпирическую варьируемую корректирующую функцию (2)
- 7 043532
где m*cor - варьируемый корректирующий коэффициент, с помощью которого обеспечивают максимально возможное совпадение в процессе подгонки семейства варьируемых расчетных зависимостей
Invi, pk, c) с экспериментально полученным семейством Inmi, pk, c). Выбор такой формы эмпирической корректирующей функции основан на анализе отклонений семейства подогнанных при минимизации (1) ** зависимостей Invi, pk, c) без учета fcor(m cor, pk, c), т.е. когда fcor(m cor, pk, c)=1, от экспериментально полученных Inmi, pk, c). Для ряда мелкодисперсных материалов. Получаемые значения корректирующей функции отличаются от единицы не более чем на 10%, что соответствует существующим оценкам относительных погрешностей определения μ's и μ a с помощью диффузионного приближения [4, 5, 17].
Поскольку значения спектрально-пространственные профилей Inmi, pk, c) и Inv[μ' s0), μ λ, μ ai, pk, c), nrel] быстро падают при увеличении pk, при составлении функции невязки в первой сумме (1) используют их отношение, из которого вычитается единица
Это позволяет с одинаковым весом учитывать в функции невязки зависимости, полученные для разных значений pk.
Как видно из поведения траектории минимума среднеквадратичного отклонения (СКО) для первой суммы (11) функции невязки (1), приведенной на фиг. 6 для полиуретана при fcor(m cor, pk, c)=1, наблюдается монотонное падение СКО невязки при уменьшении μ a вплоть до точки глобального минимума СКО, показанного маркером в виде кружка. Выявлено, что по мере роста среднего по спектру значения приведенного показателя рассеяния μ's мелкодисперсного материала спад СКО для первой суммы функции невязки (11) становится более пологим и начинает наблюдаться отсутствие глобального минимума, что иллюстрирует зависимость на фиг. 7 для образца белого пищевого полистирола. Это еще раз подтверждает, что модели [2, 3] диффузионного приближения не полностью учитывают механизмы процессов подповерхностного переноса излучения в мелкодисперсных средах.
Поэтому функцию невязки (1) предложено дополнять второй квадратичной суммой разностей искомых спектральных зависимостей μ ai, pk, c) и их усредненной по значениям pk зависимости (μ ai))
При вычислении варьируемых спектральных зависимостей μ*ai, pk, c), входящих в (1) и (12), ис* пользуют интерполяцию на неравномерной сетке плоскостей нормированных расчетных значений Inr( μ' s, μ a, pk, c, nrel), которые делят на варьируемые значения корректирующей функции fcor(m cor, pk, c), что позволяет при варьировании значения m* cor дополнительно учесть как влияние процессов подповерхностного переноса излучения, так . и отклонения от диффузионного приближения, вызываемые нарушением условий его применимости при использовании приближения [2] или °·1 для приближения [3]. Введение в функцию невязки (1) второй суммы (12) обеспечивает выявление глобального минимума СКО функции невязки, поскольку рассчитываемые в процессе подгонки спектральной зависимости μ ai, pk, c) будут сближаться друг с другом, что приведет к формированию своего минимум СКО квадратичной суммы (12) и, следовательно, к возникновению глобального минимума функции невязки (1).
Результаты минимизации (1), т.е. подгонки расчетных варьируемых зависимостей Invi, pk, c)/fcor(m*cor, pk, c)=Inv[μ'*s(λo), μ*λ, μ*a(λi, pk, c, nrel)]/fcor(m*cor, pk, c) к экспериментальным зависимостям Inm(Xi, pk, c) для образца полиуретана в широком участке спектра, приведены на фиг. 8, а полученные при этом спектральные зависимости μ' s(λ) и μ α(λ) - на фиг. 9. Вследствие хаотичной ориентации частиц дисперсной фазы используемого образца полиуретана логично полагать, что процессы рассеяния вводимого в него излучения, носят изотропный характер, т.е. фактор анизотропии g~0. При этом результаты расчетов при использовании моделей диффузионного приближения [2, 3] совпадают. Формируя произведение можно приблизительно выявить участки спектра, в которых условия при менимости диффузионного приближения нарушаются. Для используемого образца полиуретана на фиг. 10 приведены полученные в результате подгонки зависимости μeff(λ)·pk, показывающие, что в участке спектра примерно от 500 до 900 нм условия применимости диффузионного приближения нарушены.
Неплохое сближение подогнанных зависимостей Invi, pk, c)/fcor(m cor, pk, c) и полученных экспериментально Inmi, pk, c) при таких условиях обеспечивается корректирующей функцией, вид которой для образца полиуретана показан на фиг. 11. Внимательно анализируя зависимости на фиг. 8, можно убедиться, что некоторое отклонение подогнанных зависимостей Invi, pk, c)/fcor(m cor, pk, c) от экспериментально измеренных Inmi, pk, c) наблюдаются в участке спектра, где условия применимости диффузионного приближения нарушены. Корректирующая функция в таких случаях уменьшает вогнутость резуль* тирующих зависимостей Invi, pk, c)/fcor(m cor, pk, c) вдоль координаты р, поскольку коэффициент
- 8 043532 m*cor=1,07 оказался больше единицы.
Как уже отмечалось выше, при увеличении среднего по спектру значения приведенного показателя рассеяния μ's материала крутизна спада СКО для первой суммы функции невязки (1) падает, что будет снижать точность определения искомых спектральных зависимостей μ' si) и μ ai). Зависимости, приведенные на фиг. 12-13, показывают, что, несмотря на выполнение условий применимости диффузионного приближения и неплохое совпадение подогнанных зависимостей Invi, pk, c)/fcor(m cor, pk, c) с экспериментально полученными Inmi, pk, c), усредненные спектральные зависимости (μ aR(λ)) и (μ aL(λ)) для образца белого полистирола, полученные при регистрации ореола отраженного излучения справа и слева от точки ввода излучения в полистирол и представленные на фиг. 13б, имеют смещение относительно друг друга. Причина этого смещения - малая крутизна спада СКО первой суммы в (1).
Поскольку условие применимости диффузионного приближения для образца полистирола соблюдается, полученная корректирующая функция fcor(m cor, pk, c) с коэффициентом m cor=0,812<1, вид которой представлен на фиг. 14, увеличивает вогнутость результирующих зависимостей Inv[μ' s0), μ λ, μ ai, pk, c, nrel)]fcor(m cor, pk, c) вдоль координаты p, учитывая таким образом влияние показателя преломления полистирола nt и состояния его поверхности на процессы подповерхностного переноса излучения.
Приведенные примеры для двух образцов материалов с разной удельной плотностью дисперсной фазы и, следовательно, усредненных по спектру приведенных показателей рассеяния p's иллюстрируют, что предложенные решения позволяют достичь поставленных в изобретении целей.
Таким образом, предлагаемые способ определения спектральных показателей приведенного рассеяния и поглощения излучения мелкодисперсными материалами и устройство для его осуществления позволяют за счет использования широкополосного источника светового излучения, вводимого в образец, регистрации нормированного спектрально-пространственного профиля интенсивности Inmi, pk, c)=I(λi, pk)/I(λi, pc) отраженного излучения, решения системы нелинейных уравнений путем минимизации функции невязки в виде (1), в состав которой входит варьируемый спектрально-пространственный профиль Inv[μ' 80), μ λ, μ ai, pk, c), nrel], деленный на корректирующую функцию (2), а также варьируемые спектральные зависимости μ ai, pk, c), рассчитываемые путем интерполяции на неравномерной сетке отсчетов семейства заранее рассчитанных с помощью диффузионного приближения поверхностей Inr(μ' s, μ a, pk, c)=I(μ' s, μ a, pk, nrel)/I(μ' s, μ a, pc, nrel) с учетом текущих значений варьируемой зависимости μ' si) и текущих значений корректирующей функции (2), повысить устойчивость определения искомых спектральных зависимостей μ' si) и μ ai) при варьировании в процессе минимизации (1) всего тремя параметрами μ'*s(λo), μ*λ и m*cor.
Определенные при минимизации (1) спектральные зависимости μ' si) и μ ai) далее можно использовать как для оперативной проверки отклонения контролируемых параметров и характеристик исследуемого мелкодисперсного материала, в частности μ'8(λ) и μa(λ), от аналогичных свойств эталонного образца, так и оперативного определения концентрации оптически активных компонент, входящих в состав этого материала. Определение концентрации компонент осуществляют путем решения обратной спектрометрической задачи по определению состава многокомпонентных сред по измеренной или определенной экспериментально в широком участке спектра суммарной спектральной зависимости показате^1, μ»Σ(λ) = Σε™μ™(λ)’ ля поглощения от длины волны где cm - неизвестные значения концентраций m-ых компонент исследуемого образца материала, μam(λ) - заранее известные спектральные показатели поглощения светового излучения оптически активными m-ыми компонентами материала.
Предлагаемый способ и устройство является основой для развития перспективного направления спектроскопии диффузного рассеяния с пространственным разрешением для оперативного контроля оптических свойств и компонентного состава различных мелкодисперсных материалов и веществ.
Широкое применение это направление найдет при оперативном контроле свойств ряда пищевых продуктов (молоко, сыры, масло, фрукты, овощи, корнеплоды и т.д.), а также при исследованиях изменений свойств этих продуктов в процессе вызревания и хранения. Описанный способ и устройство имеет перспективы широкого применения и в медицине при контроле состояния микроциркуляторного русла кровеносной системы человека (контроль процессов доставки и потребления кислорода тканями на основе определения средней оксигенации крови в капиллярах и венулах) и гидратации тканей, что важно при лечении их отеков и возрастных изменений, а также при контроле процессов восстановления тканей при лечении травм, ранений и ожогов.
Совокупность предлагаемых решений позволяет улучшить устойчивость определения искомых спектральных показателей по сравнению с известными способами и упростить конструкцию устройства для его осуществления за счет исключения подвижных оптико-механических узлов или набора оптоволокон, подключаемых к многоканальному спектрометру.
-
EA202292115 2022-06-30 Способ определения спектральных показателей приведенного рассеяния и поглощения светового излучения мелкодисперсными материалами и устройство для его осуществления EA043532B1 (ru)

Publications (1)

Publication Number Publication Date
EA043532B1 true EA043532B1 (ru) 2023-05-30

Family

ID=

Similar Documents

Publication Publication Date Title
CN102058393B (zh) 基于反射光谱测量的皮肤生理参数与光学特性参数的测量方法
US8634077B2 (en) Methods and systems for optically characterizing a turbid material using a structured incident beam
US20160146730A1 (en) Systems and methods for diagnosis of epithelial lesions
Pfefer et al. Reflectance-based determination of optical properties in highly attenuating tissue
Martelli et al. Phantoms for diffuse optical imaging based on totally absorbing objects, part 1: Basic concepts
Eisel et al. Investigation of optical properties of dissected and homogenized biological tissue
Rogelj et al. Effect of curvature correction on parameters extracted from hyperspectral images
Hohmann et al. Measurement of optical properties of pig esophagus by using a modified spectrometer set‐up
Vincely et al. Accuracy of retrieving optical properties from liquid tissue phantoms using a single integrating sphere
Tian et al. Quantitative characterization of turbidity by radiative transfer based reflectance imaging
EA043532B1 (ru) Способ определения спектральных показателей приведенного рассеяния и поглощения светового излучения мелкодисперсными материалами и устройство для его осуществления
JP2012063333A (ja) 光学特性値計測装置、光学特性値計測方法及び光学特性値計測プログラム
Lim et al. Hyperspectral photoacoustic spectroscopy of highly-absorbing samples for diagnostic ocular imaging applications
Diffey et al. Ultraviolet erythema: dose response and mediator diffusion
Jütte et al. Mueller matrix-based approach for the ex vivo detection of riboflavin-treated transparent biotissue
Kang et al. Polarization color imaging system for on‐line quantitative evaluation of facial skin lesions
Lee et al. Performance estimation of optical skin probe in short wavelength infrared spectroscopy based on Monte-Carlo simulation
CN113208562A (zh) 一种基于光检测技术的皮肤含水量检测系统及方法
Pham et al. Quantifying the optical properties and chromophore concentrations of turbid media by chemometric analysis of hyperspectral diffuse reflectance data collected using a Fourier interferometric imaging system
Wares et al. Noninvasive evaluation of hemodynamics and light scattering property during two-stage mouse cutaneous carcinogenesis based on multispectral diffuse reflectance images at isosbestic wavelengths of hemoglobin
Niedorf et al. Noninvasive reflection spectra provide quantitative information about the spatial distribution of skin chromophores
Lisenko et al. Method for determining skin pigment concentrations from multispectral images of the skin
He et al. Phantom preparation and optical property determination
JP2012063332A (ja) 表示装置、表示方法、プログラム、および、光学特性表示システム
Lisenko et al. Systems for real-time optical diagnostics of biological objects