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

Авторы патента:

Шиховцев Артем Юрьевич (RU)

Ковадло Павел Гаврилович (RU)

G01W1/00 - Метеорология

G01N21/17 - системы, в которых на падающий свет влияют свойства исследуемого материала (в которых исследуемый материал оптически возбуждается с тем, чтобы вызвать изменение длины волны падающего света G01N 21/63)

G01N2021/1765 - Исследование или анализ материалов путем определения их химических или физических свойств (разделение материалов вообще B01D,B01J,B03,B07; аппараты, полностью охватываемые каким-либо подклассом, см. в соответствующем подклассе, например B01L; измерение или испытание с помощью ферментов или микроорганизмов C12M,C12Q; исследование грунта основания на стройплощадке E02D 1/00;мониторинговые или диагностические устройства для оборудования для обработки выхлопных газов F01N 11/00; определение изменений влажности при компенсационных измерениях других переменных величин или для коррекции показаний приборов при изменении влажности, см. G01D или соответствующий подкласс, относящийся к измеряемой величине; испытание

G01J9/00 - Измерение оптической разности фаз (управление фазой световых лучей G02F 1/01); определение степени когерентности; измерение оптической длины волны (спектрометрия G02F 3/00)

G01J2009/002 - Измерение интенсивности, скорости или спектрального состава, поляризации, фазы или импульсных характеристик инфракрасных, видимых или ультрафиолетовых лучей; колориметрия; радиационная пирометрия (источники света F21,H01J,H01K,H05B; определение свойств веществ оптическими средствами G01N)

Владельцы патента RU 2712464:

Федеральное государственное бюджетное учреждение науки Ордена Трудового Красного Знамени Институт солнечно-земной физики Сибирского отделения Российской академии наук (ИСЗФ СО РАН) (RU)

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

где C(δi, δj, δt) - пространственно-временная кросскорреляционная функция локальных наклонов волновых фронтов, A₁(δi, δj, δt) - взаимнокорреляционная функция, δi и δj - расстояния между субапертурами, δi и δj - расстояния между субапертурами (участков на апертуре телескопа), индексами i, j обозначены позиции субапертур на апертуре вдоль произвольно выбранных ортогональных осей координат, δt - временной промежуток между отсчетами, F - знак фурье-преобразования, - знак обратного фурье-преобразования, K - матричная весовая функция, имеющая размерность м^-2/3, определение матричной весовой функции производится по данным интегральных значений амплитуды искажений волнового фронта, высоты атмосферных слоев определяются по формулам где α - зенитный угол Солнца. Технический результат заключается в повышении пространственного разрешения и расширении области применимости. 3 ил.

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

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

Наиболее близким к заявляемому способу измерения вертикальных профилей показателя преломления воздуха для коррекции солнечных изображений в реальном времени является способ [1], заключающийся том, что вертикальные профили показателя преломления воздуха измеряют по локальным наклонам волновых фронтов от двух разнесенных на известный угол светящихся объектов (фрагментов на солнечном изображении), находящихся в поле зрения одного телескопа, с помощью двух одинаковых датчиков волнового фронта. Рассчитывают усредненные во времени пространственные кросскорреляционные функции локальных наклонов волновых фронтов от двух источников света, которые регистрируют с помощью двух датчиков от скрещенных оптических пучков в разных частях поля зрения астрономического телескопа. Определяют автокорреляционную функцию локальных наклонов волнового фронта на одном из датчиков. По усредненным по времени кросскорреляционным функциям локальных наклонов двух волновых фронтов и автокорреляционной функции локальных наклонов на одном из датчиков волнового фронта оценивают характеристику показателя преломления воздуха.

Формула для расчета вертикального профиля характеристики показателя преломления воздуха по этому способу выглядит следующим образом:

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

взаимная корреляционная функция для выбранного датчика

δi и δj - расстояния между субапертурами (участков на апертуре телескопа), индексами i, j обозначены позиции субапертур на апертуре вдоль произвольно выбранных ортогональных осей координат, - локальные наклоны волнового фронта на субапертуре i,j камеры (1), - локальные наклоны волнового фронта на субапертуре i+δi, j+δj камеры (2), O(δ(δδj) - число пересекающихся освещенных субапертур на одном из датчиков волнового фронта, <> - знак усреднения по временной реализации, F - знак Фурье преобразования, - знак обратного Фурье преобразования, р - матричная весовая функция, имеющая размерность м^-2/3, N и М - общее количество субапертур вдоль каждой из произвольно выбранных перпендикулярных осей координат.

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

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

Измерения локальных наклонов волновых фронтов на каждой субапертуре выполняют с помощью одного датчика волнового фронта, работающего по одному источнику света произвольной формы и размера, позиция которого изменяется за счет суточного вращения Земли. Пересечения световых пучков формируют между положением пучка в начальный момент времени и последующим положением пучка вследствие суточного движения объекта по небосводу. На втором шаге расчетов положение опорного пучка смещено на угол (на один шаг по времени). Далее, по данным измерений локальных наклонов волнового фронта рассчитывают усредненные по отдельным временным интервалам пространственно-временные кросскореляционные функции локальных наклонов, формируемых в скрещенных оптических пучках, в начальный момент и последующие моменты времени. Из зарегистрированного набора оптических пучков выбирают все скрещенные пучки, «точки пересечения» которых располагаются на разных высотных уровнях по лучу зрения. Длину временного интервала определяют исходя из минимальной высоты измерений профиля характеристик показателя преломления воздуха. Количество пересечений пучков выбирают в соответствии с заданным пространственным разрешением положения расчетных уровней по вертикали. Нижняя высота определения единичного профиля - 20 м, что по времени измерений занимает до 20 с. Пространственное разрешение способа возрастает за счет высокой дискретизации углов во времени и ограничивается скоростью регистрации датчика волнового фронта и угловым смещением наблюдаемого объекта во времени. Ранжирование «точек пересечения» световых пучков, разнесенных по апертуре телескопа и времени регистрации, связанного с углом между последовательными положениями пучков выполняют согласно соотношению:

где D - диаметр входного зрачка телескопа, h - высота слоя, n - число выбранных субапертур датчика волнового фронта, θ(t) - переменный угол между источниками световых пучков в разные моменты времени, α - зенитный угол Солнца. Полученные амплитуды пространственно-временных кросскореляционных функции показателя преломления воздуха, привязанные к высотным уровням, позволяют рассчитать вертикальный профиль характеристики показателя преломления воздуха по формуле (1):

где пространственно-временную кросскорреляционную функцию локальных наклонов волновых фронтов C(δi,δj,δt) находят по формуле (2):

взаимнокорреляционную функцию A₁((δi,δj,δt) определяют по формуле (3)

δi и δj - расстояния между субапертурами, индексами i, j обозначены позиции субапертур на апертуре вдоль произвольно выбранных ортогональных осей координат, s_0,0 - локальные наклоны волнового фронта на опорной субапертуре с координатами (0,0) в момент времени t=0, s_δi,δj,δt - локальные наклоны волнового фронта на субапертуре сдвинутой на шаг относительно исходной, в том числе, на субапертурах в последующие моменты времени O(δ(δδj) - число пересекающихся субапертур на датчике волнового фронта, <> - знак усреднения по временной реализации, F - знак Фурье преобразования, - знак обратного Фурье преобразования, N и M - общее количество субапертур вдоль каждой из произвольно выбранных перпендикулярных осей координат, K - матричная весовая функция, имеющая размерность м^-2/3. Определение матричной весовой функции производится по данным интегральных значений амплитуды искажений волнового фронта. Высоты атмосферных слоев определяются по формулам

Способ был реализован на Большом солнечном вакуумном телескопе (БСВТ) Байкальской астрофизической обсерватории ИСЗФ СО РАН. Датчик волнового фронта был установлен в параллельном пучке в плоскости, сопряженной с плоскостью апертуры телескопа. Датчик имел 8×8 субапертур и цифровую камеру с размером матрицы 7×7 мм, 740×480 пиксел, частота 100 кадров в секунду. В качестве источника света использовался край солнечного изображения. Полные серии измерений занимали ~ 45 сек, при этом на восстановление единичного вертикального профиля затрачивалось ~ 15 сек. Использование одного датчика волнового фронта позволило сохранить плотность светового потока на каждой субапертуре и, как следствие, повысить точность измерений локальных наклонов волнового фронта. По сравнению с аналогом при прочих равных условиях точность измерений возросла более чем в 2 раза. Кроме этого, что важно, увеличилось число высотных уровней более чем в 10 раз.

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

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

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

На рисунке 3 (верхняя фигура) показаны изменения приведенных к зениту средних значений радиуса Фрида r₀ в дневное время, рассчитанные по измерениям датчика волнового фронта Шака-Гартмана. По оси абсцисс показано время наблюдений, по оси ординат - значения радиуса Фрида. Непрерывной кривой показано изменение средних значений радиуса Фрида, определенных по дрожанию субизображений по 52 парам разнесенных субапертур. Заливкой показаны средние квадратические отклонения значений радиуса Фрида. Оценка радиуса Фрида производилась по данным измерений дрожания края солнечного изображения, выполненных за 28.06.2018.

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

Анализ рис. 2 и 3 показывает, что значения радиуса Фрида, определенные в слое от 0 до 20 км, имеют тенденцию к уменьшению при возрастании величины k, характеризующей степень развитости турбулентности и наоборот. Характер изменений этих величин указывает на реалистичность результатов измерений применяемым способом.

Характер распределения характеристик показателя преломления воздуха по вертикали согласуется с полученными данными выполненных лидарных наблюдений вертикальной структуры воздушных течений. Экспериментальные исследования в месте расположения БСВТ, выполненные с помощью Stream Line lidar, позволили выявить особенности течений в нижнем слое атмосферы. Был подтвержден двухслойный характер изменений скорости и направления ветра с высотой. Полученные из измерений оценки высот повышенной интенсивности турбулентности заявляемым способом оказались в согласии с характерными высотами, определяемыми из лидарных измерений.

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

1. Wilson R.W. "SLODAR: measuring optical turbulence altitude with a Shack-Hartmann wavefront sensor" / MNRAS Volume 337, Issue 1, 2002. pp. 103-108.

2. Townson M.J. "Correlation wavefront sensing and turbulence for solar adaptive optics / A thesis presented for the degree of Doctor of Philosophy" / University of Durham. United Kingdom. 2016. P. 139.

3. Hickson P. "Fundamentals of atmospheric and adaptive optics" / The University of British Columbia, Department of Physics and Astronomy. 2008. P. 68.

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

где C(δi, δj, δt) - пространственно-временная кросскорреляционная функция локальных наклонов волновых фронтов, A₁(δi, δj, δt) - взаимнокорреляционная функция, δi и δj - расстояния между субапертурами, δi и δj - расстояния между субапертурами (участков на апертуре телескопа), индексами i, j обозначены позиции субапертур на апертуре вдоль произвольно выбранных ортогональных осей координат, δt - временной промежуток между отсчетами, F - знак Фурье-преобразования, - знак обратного Фурье-преобразования, K - матричная весовая функция, имеющая размерность м^-2/3, определение матричной весовой функции производится по данным интегральных значений амплитуды искажений волнового фронта, высоты атмосферных слоев определяются по формулам где α - зенитный угол Солнца.

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

Способ определения структурно-стратиграфических особенностей строения снежного покрова // 2704432

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

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

Изобретение относится к способам экологического мониторинга. Сущность: обрабатывают материалы аэрокосмической съемки.

Система для мониторинга искусственных сооружений высокоскоростной магистрали // 2698419

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

Автоматизированная система экологического мониторинга и прогнозирования загрязнения атмосферного воздуха промышленного региона // 2697571

Изобретение относится к области экологии и может быть использовано для экологического мониторинга и прогнозирования загрязнения атмосферы промышленного региона.

Способ автоматического определения параметров ионосферных слоев по ионограммам // 2697433

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

Способ определения влагозапаса в снежном покрове // 2695949

Изобретение относится к области метеорологии и может быть использовано для определения влагозапаса снежного покрова. Сущность: проводят три измерения мощности дозы гамма-излучения в приземной атмосфере.

Способ разведки ледовой обстановки на северном морском пути // 2694085

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

Способ определения влагозапаса в снежном покрове // 2694080

Изобретение относится к области метеорологии и может быть использовано для определения влагозапаса снежного покрова. Сущность: измеряют плотность потока бета-излучения над снежным покровом в период перед началом таяния снега в дневное время суток не менее чем через 3,5 часа после выпадения снега.

Способ проведения метеомониторинга для прогнозирования вероятности повреждения электросетевого оборудования // 2675655

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

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

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

Устройство для измерения длины распространения инфракрасной поверхностной электромагнитной волны // 2703772

Изобретение относится к области исследования поверхности металлов и полупроводников оптическими методами и касается устройства для измерения длины распространения инфракрасной поверхностной электромагнитной волны (ПЭВ).

Устройство для определения длины распространения поверхностной электромагнитной волны инфракрасного диапазона за время одного импульса излучения // 2699304

Изобретение относится к области исследования поверхности материалов оптическими методами и касается устройства для определения длины распространения поверхностной электромагнитной волны (ПЭВ) инфракрасного диапазона за время одного импульса излучения.

Способ спектральной диагностики оптических осей и типов колебательных центров в кристаллах с водородными связями // 2697425

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

Способ и анализатор для осуществления способа определения содержания оксидов азота в ракетных окислителях // 2695816

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

Способ спектрального лазерного сканирования композитных материалов в соответствии с оптической плотностью его матрикса и составных компонентов // 2692825

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

Способ определения времени релаксации плёночного просветляющегося поглотителя с помощью фемтосекундного волоконного лазера в режиме генерации // 2687991

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

Устройство для определения коэффициента затухания поверхностной электромагнитной волны инфракрасного диапазона за время одного импульса излучения // 2681658

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

Устройство для измерения длины распространения инфракрасной поверхностной электромагнитной волны // 2681427

Изобретение относится к области бесконтактного исследования поверхности металлов и полупроводников и касается устройства для измерения длины распространения инфракрасной поверхностной электромагнитной волны (ПЭВ).

Неинвазивный анализ материала // 2681260

Изобретение относится к измерительной технике и может быть использовано для неинвазивного анализа материала. Раскрыты способ и система для анализа материала (100).

Способы обнаружения маркера для активного туберкулеза // 2712270

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