RU2059280C1 - Способ определения геометрических характеристик объекта многоапертурной оптической системой - Google Patents
Способ определения геометрических характеристик объекта многоапертурной оптической системой Download PDFInfo
- Publication number
- RU2059280C1 RU2059280C1 SU3186760A RU2059280C1 RU 2059280 C1 RU2059280 C1 RU 2059280C1 SU 3186760 A SU3186760 A SU 3186760A RU 2059280 C1 RU2059280 C1 RU 2059280C1
- Authority
- RU
- Russia
- Prior art keywords
- image
- averaging
- average energy
- optical system
- geometric characteristics
- Prior art date
Links
Images
Landscapes
- Image Analysis (AREA)
Abstract
Использование: техническая физика. Сущность изобретения: в способе измеряют среднюю энергию изображения E, выделяют из изображения со средней энергией, превышающей уровень E в пределах дифракционного разрешения, а усреднение проводят с весами, пропорциональными измеренным средним энергиям в выделенных спеклах. 3 ил.
Description
Изобретение относится к технической физике, и может быть использовано в оптической астрономии для определения геометрических характеристик удаленных объектов.
В современной астрономии чрезвычайно важной является задача получения информации об удаленных объектах с высоким угловым разрешением. Решение многих проблем в астрофизике связано с определением геометрических характеристик (например, диаметров) небесных тел с угловым разрешением порядка 0,001". Хорошо известно, что угловое разрешение телескопа принципиально ограничено конечным размером D его апертуры. Существующие крупные телескопы имеют апертуры с диаметрами 4-6 м, однако, дальнейшее увеличение их размеров связано с резким увеличением стоимости, которая возрастает приблизительно пропорционально D3. В то же время известен альтернативный подход к проблеме достижения высокого углового разрешения, восходящий еще к пионерским работам Майкельсона и Пиза по определению угловых размеров ряда крупных звезд (Борн М. и Вольф Э. Основы оптики. М. Мир, 1973).
Известный способ звездного интерферометра Майкельсона основан на формировании интерференционного изображения наблюдаемого объекта путем когерентного сложения световых пучков от разнесенных на большое расстояние апертур, измерении видности интерференционных полос и определении по значениям видности при различных базовых расстояниях между апертурами геометрических характеристик объекта.
Недостатком способа является необходимость ограничения размеров апертур характерным размером корреляции rс атмосферных искажений фазы принимаемого излучения (≈10 см) с целью устранения влияния турбулентной атмосферы Земли, что резко ограничивает количества света в интерференционном изображении и, в конечном итоге, обуславливает низкую точность метода. Данное ограничение привело к тому, что, войдя во все учебники по оптике, метод звездного интерферометра Майкельсона фактически не применялся на практике в течение нескольких десятилетий.
Возрождение интереса к многоапертурным оптическим системам связано с появлением способа, предложенного Лабейри (A.Labeyrie, Ann. Rev. Astron. Astrophys. 1978, v.16, р.77-102), который и взят за прототип. Данный способ основан на короткоэкспозиционной регистрации в многоапертурной системе серии М искаженных атмосферой изображений, формировании их Фурье-полей, усреднении их интенсивностей путем накопления и оценке по результатам усреднения геометрических характеристик объекта. Данный способ применим к многоапертурным системам телескопов любого диаметра.
Однако для достижения удовлетворительной точности оценки необходимо зарегистрировать достаточно большое число М изображений (М 1000). На практике реальное число изображений в серии составляет около 100. При этом метод дает недостаточную точность. Кроме того, чрезвычайно велики требования к детектору изображения: в системах с большими (D >> ro) телескопами, разнесенными на расстояние L, превосходящее диаметр апертур D (L >> D) общий размер искаженного изображения определяется выражением F, где F фокусное расстояние системы; λ длина волны света, в то время, как ширина интерференционной полосы составляет F. Число элементов разрешения детектора изображения вдоль одной координатной оси, таким образом, должно быть порядка 2 , а в двумерной системе 2 . Для многоапертурной системы с L ≈ 100 м (напомним, что в экспериментах Лабейри L превышало 50 м) и ro ≈ 10 см необходим приемник с разрешением 2000х2000. Создание таких высокочувствительных детекторов (например, ПЗС-камер) представляет собой сложную техническую задачу.
Цель изобретения заключается в повышении точности определения геометрических характеристик объекта при одновременном понижении требований к устройству реализации способа.
Для достижения цели перед регистрацией изображения измеряют среднюю энергию изображения Е, выделяют области размера, не превосходящего дифракционный элемент разрешения отдельной апертуры, средняя энергия в пределах которых превышает Е, регистрируют части формируемого изображения в этих областях, формируют Фурье-поля зарегистрированных частей изображения, измеряют их интенсивности и усредняют с весами, пропорциональными измерениям средним энергиям в выделенных областях, по усредненным интенсивностям определяют геометрические характеристики объекта.
Сравнительный анализ известных технических решений и предложенного технического решения позволяет выделить в качестве характерных следующие отличительные признаки:
выделение областей изображения размера, не превосходящего дифракционный элемент апертуры, средняя энергия в пределах которых превышает среднюю энергию во всем изображении;
измерение интенсивностей Фурье-полей по отдельным областям;
усреднение измеренных интенсивностей.
выделение областей изображения размера, не превосходящего дифракционный элемент апертуры, средняя энергия в пределах которых превышает среднюю энергию во всем изображении;
измерение интенсивностей Фурье-полей по отдельным областям;
усреднение измеренных интенсивностей.
Кратко рассмотрим математическую сторону способа. Остановимся вначале на модели искаженного атмосферой изображения объекта, формируемого многоапертурной системой. В качестве простейшей многоапертурной системы будем рассматривать звездный интерферометр Майкельсона с большими (D >> ro) апертурами.
Изображение точечного источника, формируемое каждой апертурой, представляет собой картину случайно расположенных пятен спеклов, причем размер каждого близок к дифракционному элементу разрешения апертуры, а число их оценивается как (D/ro)2. Интерференционное изображение аналогичным образом представимо в виде набора пятен (фиг.2), однако каждое пятно размера F промодулировано тонкой интерференционной структурой с периодом F, причем ее фаза случайно меняется от пятна к пятну. При наблюдении объекта, не разрешаемого отдельной апертурой, пятна размера F практически не расплываются, однако меняется видность тонкой интерференционной структуры в пределах каждого пятна (по видности или, что математически эквивалентно, по интенсивности Фурье-поля на соответствующей пространственой частоте / можем определять геометрические характеристики объекта). На основании данной модели легко определить точность оценки интенсивности Фурье-поля по Лабейри: обозначим распределения интенсивности в каждом i-м пятне изображения Ii(x), их Фурье-поля (f). Тогда все изображение имеет вид
I(x)Ii(x); (f)(f)
Формируемая в методе Лабейри величина
= (fL)(fL) (fL)+ (1) распадается на 2 суммы, первая из которых состоит из N членов и дает искомую величину, а вторая сумма N2-N членов со случайной фазой, величина которой ≈ N, представляет собой ошибку метода. Таким образом, отношение сигнал-шум метода Лабейри при усреднении (1) по М изображениям составит при ≫ 1 величину (она совпадает с оценкой, полученной более строгими методами).
I(x)Ii(x); (f)(f)
Формируемая в методе Лабейри величина
= (fL)(fL) (fL)+ (1) распадается на 2 суммы, первая из которых состоит из N членов и дает искомую величину, а вторая сумма N2-N членов со случайной фазой, величина которой ≈ N, представляет собой ошибку метода. Таким образом, отношение сигнал-шум метода Лабейри при усреднении (1) по М изображениям составит при ≫ 1 величину (она совпадает с оценкой, полученной более строгими методами).
Если же в искаженном изображении, обработав каждый спекл, сформировать (fL) и сложить их, то точность оценки <(fL)> будет определяться отношением сигнал-шум для (fL) (оно порядка единицы) и оценивается, как , а при усреднении по М искаженным изображениям улучшается до . Обрабатывая 100 независимых искаженных атмосферой изображений, сформированных системой с метровыми телескопами, может получить оценку величин интенсивностей Фуpье-спектров распределения интенсивности наблюдаемого объекта с хорошим отношением сигнал-шум (порядка 100).
Спеклы в изображении могут быть выделены как области дифракционного размера F, энергия в которых превышает среднюю по изображению Е (см.фиг.3).
Регистрируемые части изображения неизбежно содержат аддитивный шум, который можно считать однородным по изображению. Чтобы простейшим образом учесть точность регистрации, достаточно складывать интенсивности Фуpье-полей зарегистрированных частей изображения с весами, пропорциональными средним энергиям в этих частях, что повышает относительный вклад более ярких и, следовательно, зарегистрированных с большей точностью пятен.
Реализация предложенного способа устройством, схема которого представлена на фиг.1, осуществляется следующим образом.
Искаженное атмосферой световое излучение от принимаемого объекта 1 принимают и фокусируют многоапертурной системой 2, разделяют сфокусированный пучок светоделительным зеркалом 3,формируют на матричном детекторе 5 изображение объекта и регистрируют его с разрешением , определяют с помощью компаратора 6 среднюю энергию в изображении и определяют положение областей размера спеклов по отношению энергии в них над средней энергией, с помощью устройства 7 управления светоделительным зеркалом путем поворота зеркала 3 направляют и фокусируют с изменением масштаба микрообъективом 4 на матричном детекторе 8 выделенные области изображения, последовательно регистрируют их с разрешением определяют интенсивность их Фурье-спектров с помощью спектрального анализатора 9, усредняют их в сумматоре 10 с весами, определенными в компараторе 6 и по величинам, полученным в результате усреднения с помощью цифрового процессора 11 определяют геометрические характеристики объекта. Следует отметить, что устройства 6, 9 и 10 легко могут быть реализованы на базе простейших аналого-цифровых устройств.
Положительный эффект от использования предлагаемого способа заключается в повышении точности оценки геометрических характеристик наблюдаемого объекта при одновременном понижении требований к устройству реализации метода.
При анализе предложенного способа было проведено его экспериментальное моделирование на ЭВМ СМ-1420 с использованием М=30 искаженных короткоэкспозиционных изображений объекта двойной звезды, смоделированных для 8-апертурной системы с 10. Оказалось, что при использовании для определения величин квадрата модуля Фурье-спектра объекта (по которым затем определялась его автокорреляция, несущая информацию о геометрических характеристиках расстоянии между звездами) на 8-ми формируемых пространственных частотах способа-прототипа относительная среднеквадратичная ошибка составила 39% а при определении указанных величин предложенным способом 8% что соответствует повышению точности в 4,9 раз.
В настоящее время на предприятии разрабатывается конструкторская документация на устройство, реализующее предложенный способ.
Claims (1)
- СПОСОБ ОПРЕДЕЛЕНИЯ ГЕОМЕТРИЧЕСКИХ ХАРАКТЕРИСТИК ОБЪЕКТА МНОГОАПЕРТУРНОЙ ОПТИЧЕСКОЙ СИСТЕМОЙ, заключающийся в короткоэкспозиционной регистрации серии изображений, формировании их Фурье-полей, усреднении их интенсивностей путем накопления и оценке по результатам усреднения геометрических характеристик объекта, отличающийся тем, что, с целью повышения точности определения, до регистрации измеряют среднюю энергию изображения Е, выделяют из изображения со средней энергией, превышающей уровень Е в пределах дифракционного элемента разрешения, а усреднение проводят с весами, пропорциональными измеренным средним энергиям в выделенных спеках.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SU3186760 RU2059280C1 (ru) | 1987-11-30 | 1987-11-30 | Способ определения геометрических характеристик объекта многоапертурной оптической системой |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SU3186760 RU2059280C1 (ru) | 1987-11-30 | 1987-11-30 | Способ определения геометрических характеристик объекта многоапертурной оптической системой |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2059280C1 true RU2059280C1 (ru) | 1996-04-27 |
Family
ID=20928834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
SU3186760 RU2059280C1 (ru) | 1987-11-30 | 1987-11-30 | Способ определения геометрических характеристик объекта многоапертурной оптической системой |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2059280C1 (ru) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2527252C1 (ru) * | 2013-06-07 | 2014-08-27 | Федеральное государственное автономное образовательное учреждение высшего профессионального образования "Волгоградский государственный университет" | Способ предотвращения угрозы для планеты путем оценки размеров пассивных космических объектов |
-
1987
- 1987-11-30 RU SU3186760 patent/RU2059280C1/ru active
Non-Patent Citations (1)
Title |
---|
Rev. Astron. Astrophys. 1978, v.16, p.77-102. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2527252C1 (ru) * | 2013-06-07 | 2014-08-27 | Федеральное государственное автономное образовательное учреждение высшего профессионального образования "Волгоградский государственный университет" | Способ предотвращения угрозы для планеты путем оценки размеров пассивных космических объектов |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
West et al. | A survey and examination of subpixel measurement techniques | |
Andersson | Long-range three-dimensional imaging using range-gated laser radar images | |
CN103592108A (zh) | Ccd芯片调制传递函数测试装置及方法 | |
CN102713503B (zh) | 干涉测距仪中的散斑减轻方法和对应的测距仪 | |
US10113945B2 (en) | Method and apparatus for combining measurements of particle characteristics using light scattering and optical imaging | |
CN109186759B (zh) | 一种光栅光谱仪像质测量方法和装置 | |
Sensiper et al. | Modulation transfer function testing of detector arrays using narrow-band laser speckle | |
US11953381B2 (en) | Wavefront sensors and wavefront shape determination using related imaging sensors | |
CN106570902B (zh) | 基于探测器像素响应频谱获取的psf相对质心计算方法 | |
CN108510480B (zh) | 基于辐射对比度的卫星探测性能评估方法、装置及存储器 | |
RU2059280C1 (ru) | Способ определения геометрических характеристик объекта многоапертурной оптической системой | |
CN112165570A (zh) | 一种基于计算鬼成像的多深度目标对焦方法 | |
Neilsen | Signal processing on digitized LADAR waveforms for enhanced resolution on surface edges | |
Guo et al. | Differential optical-path approach to improving the performance of Multiple-Slit Streak Tube Imaging LiDAR | |
Schiller et al. | Using vicarious calibration to evaluate small target radiometry | |
US4606067A (en) | Device for recognizing far-away electro-magnetically radiating objects | |
US7391519B1 (en) | Generation of spatially distributed spectral data using a multi-aperture system | |
Daniels et al. | Random targets for modulation transfer function testing | |
Berdja et al. | Simulation of pupil-plane observation of angle-of-arrival fluctuations in daytime turbulence | |
RU2669262C1 (ru) | Способ оценки и максимизации предельного инструментального разрешения космического аппарата дистанционного зондирования Земли на местности | |
Hawkins | Fourier deconvolution of electronographic images | |
Barnes et al. | Measuring aerosol optical depth (AOD) and aerosol profiles simultaneously with a camera lidar | |
Introne | Enhanced spectral modeling of sparse aperture imaging systems | |
Abdelsalam et al. | An automatic processing technique for accurate surface form measurement | |
Walczykowski et al. | Using XEVA video sensors in acquiring spectral reflectance coefficients |