RU2701880C1 - Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью - Google Patents
Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью Download PDFInfo
- Publication number
- RU2701880C1 RU2701880C1 RU2018142117A RU2018142117A RU2701880C1 RU 2701880 C1 RU2701880 C1 RU 2701880C1 RU 2018142117 A RU2018142117 A RU 2018142117A RU 2018142117 A RU2018142117 A RU 2018142117A RU 2701880 C1 RU2701880 C1 RU 2701880C1
- Authority
- RU
- Russia
- Prior art keywords
- radio
- calculated
- hologram
- uneven surface
- subsurface
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Использование: для подповерхностной радиолокации. Сущность изобретения заключается в том, что восстановление радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью, включает в себя ступенчатое изменение сигнала в заданном диапазоне частот с равномерным шагом в диапазоне где kmin=0,72; kmax=0,81; D - диаметр антенны; c - скорость света, количество отдельных частот в диапазоне от fmin до fmax не менее пяти, автоматическое выравнивание амплитудно-частотной характеристики, при этом создается цифровая карта рельефа неровной поверхности среды с использованием датчика глубины, позволяющего измерять расстояние, соответствующее каждому пикселю получаемого датчиком глубины оптического изображения, рассчитывается радиоголограмма поверхности с использованием полученного рельефа неровной поверхности среды, рассчитывается разностная радиоголограмма (разность между экспериментальной и расчетной радиоголограммами), по разностной радиоголограмме вычисляется радиоизображение подповерхностного объекта методом обратного распространения. Технический результат: обеспечение возможности скомпенсировать отражения от поверхности и получить детализированные радиоизображения подповерхностных объектов. 10 ил.
Description
Область техники
Изобретение относится к области подповерхностной радиолокации, а именно, к способам определения расположения и формы неоднородностей и включений в конденсированных средах.
Уровень техники.
Известен способ подповерхностного зондирования, (Финкельштейн М.И., Кутев В.А., Золотарев В.П. Применение радиолокационного подповерхностного зондирования в инженерной геологии, М., Недра, 1986, с. 46). Он основано на использовании непрерывного сигнала с изменением частоты по симметричному или несимметричному пилообразному закону. Частота биений между опорным / прямым / и отраженным сигналами, является функцией расстояния до объекта.
Известен также способ зондирования конденсированных сред (Journal of Applied Physics, v.56, №9, 1984, p.2575) со ступенчатым изменением частоты в заданном диапазоне.
Наиболее близким аналогом (прототипом) является способ зондирования конденсированных сред (RU 2000103678 A, 10.01.2002, G01V 3/12) со ступенчатым изменением сигнала в заданном диапазоне частот, при этом частоты зондирующего сигнала выбираются с равномерным шагом в диапазоне
где:
kmin=0,72;
kmax=0,81;
D - диаметр антенны;
с - скорость света, при этом количество отдельных частот в диапазоне от fmin до fmax должно быть не менее пяти.
Недостатком аналогов и прототипа является то, что при работе на средах с неровной поверхностью данный способ не позволяет получить детализированные радиоизображения подповерхностных объектов, так как сильные отражения от неровной поверхности среды маскируют более слабые отражения от подповерхностных объектов, расположенных дальше от плоскости регистрации радиоголограммы.
Перечень фигур, чертежей и иных материалов.
На фиг. 1 представлена гипсовая панель с заложенным под нее объектом.
На фиг. 2 представлен стенд для проведения экспериментов.
На фиг. 3 показана действительная часть радиоголограммы.
На фиг. 4 показан результат восстановления без учета рельефа поверхности.
На фиг. 5 представлена оцифровка рельефа гипсовой панели.
На фиг. 6. представлена карта глубины.
На фиг. 7 представлен рельеф гипсовой панели, полученный по данным видеосенсора Kinect v2.
Фиг. 8. Действительная часть рассчитанной радиоголограммы.
Фиг. 9. Действительная часть разностной радиоголограммы.
Фиг. 10. Результат восстановления разностной радиоголограммы.
Задача.
Техническая задача состоит в устранении указанного недостатка за счет создания цифровой карты рельефа неровной поверхности среды с использованием датчика глубины, позволяющего измерять расстояние, соответствующее каждому пикселю получаемого датчиком глубины оптического изображения, рассчитывается радиоголограмма поверхности с использованием полученного рельефа неровной поверхности среды, рассчитывается разностная радиоголограмма (разность между экспериментальной и расчетной радиоголограммами), по разностной радиоголограмме вычисляется радиоизображение подповерхностного объекта методом обратного распространения.
Отличительные признаки.
1. В отличие от известного способа, включающего в себя ступенчатое изменение сигнала в заданном диапазоне частот с равномерным шагом в диапазоне
где:
kmin=0,72;
kmax=0,81;
D - диаметр антенны;
с - скорость света, количество отдельных частот в диапазоне от fmin до fmax не меньше пяти и осуществляется выравнивание амплитудно-частотной характеристики, дополнительно создается цифровая карта рельефа неровной поверхности среды с использованием датчика глубины, позволяющего измерять расстояние, соответствующее каждому пикселю получаемого датчиком глубины оптического изображения, рассчитывается радиоголограмма поверхности с использованием полученного рельефа неровной поверхности среды, рассчитывается разностная радиоголограмма (разность между экспериментальной и расчетной радиоголограммами), по разностной радиоголограмме вычисляется радиоизображение подповерхностного объекта методом обратного распространения.
Сведения, подтверждающие возможность осуществления изобретения.
Одним из методов получения изображений объектов, находящихся в оптически непрозрачных средах, является радиовидение, когда радиоизображение подповерхностного объекта восстанавливается по рассеянному им полю (радиоголограмме), регистрируемому с помощью приемо-передатчика, перемещаемого по плоскости регистрации радиоголограммы, расположенной на небольшом расстоянии от поверхности среды.
Обычно при восстановлении радиоголограмм используется модель однородного полупространства с плоской границей раздела. При работе на средах с неровной поверхностью данная модель не позволяет получать детализированные радиоизображения подповерхностных объектов, так как сильные отражения от поверхности маскируют более слабые отражения от подповерхностных объектов, расположенных дальше от плоскости регистрации радиоголограммы.
Для получения радиоизображения подповерхностного объекта, расположенного в среде с неровной верхней поверхностью, предложен способ восстановления радиоизображений объектов, находящихся в однородных средах с неровной верхней поверхностью, который заключается в создании цифровой карты рельефа поверхности с использованием датчика глубины (RGB-D видеосенсора, позволяющего регистрировать карту глубины, то есть измерять расстояние, соответствующее каждому пикселю получаемого оптического изображения), моделировании радиоголограммы поверхности с использованием полученного рельефа среды, вычислении разностной радиоголограммы с предварительной подгонкой амплитуды и фазы экспериментальной и расчетной радиоголограмм, получении радиоизображения подповерхностного объекта путем применения метода обратного распространенения к разностной радиоголограмме.
Эксперименты проводились с использованием голографического подповерхностного радиолокатора РАСКАН с непрерывным монохроматическим сигналом частотой 6.8 ГГц. Объектом выступала вырезанная из фольги буква «R», заложенная под стеновую гипсовую 3D панель размером 600×600×38 мм (Фиг. 1). Панель располагалась на стопке листов из сухой штукатурки для имитации однородного полупространства (на фиг. 1 панель приподнята, чтобы объект был виден). Ручное сканирование осуществлялось по листу оргстекла толщиной 10 мм, расстояние от верхней поверхности листа до вершин пирамид гипсовой панели составляло 35 мм (фиг. 2). Размер области сканирования W×Н составлял 400×400 мм с шагом между выборками радиолокационного сигнала 5 мм по обеим осям. Оргстекло было размечено так, чтобы область регистрации радиоголограммы располагалась относительно панели симметрично. На фиг. 3 приведена зарегистрированная радиоголограмма, а на фиг. 4 - результат ее восстановления без учета рельефа поверхности, в предположении, что поверхность является плоской. В данном случае регулярная структура рельефа позволяет предположить наличие в правом нижнем углу постороннего объекта (но его форму определить нельзя), но в случае хаотичного рельефа даже наличие объекта определить будет невозможно.
Для оцифровки рельефа поверхности использовался RGB-D видеосенсор Kinect v2, позволяющий измерять расстояние с точностью до нескольких миллиметров. На фиг. 5 показан Kinect v2 во время получения рельефа гипсовой панели, а на фиг. 6 - данные канала дальности в виде полутонового изображения (карта глубины), где на цветовой шкале расстояние, измеряемое от видеосенсора, указано в метрах, белый цвет соответствует минимальному расстоянию, а черный - максимальному.
При работе с видеосенсором было невозможно обеспечить абсолютную параллельность оптической оси видеосенсора и нормали к поверхности панели (имеется в виду нормаль к панели как к целому объекту). На фиг. 5 видно, что левый верхний угол панели оказался несколько ближе к видеосенсору, чем правый нижний угол. Поэтому после оцифровки рельефа точки, соответствующие вершинам пирамид на гипсовой панели, были аппроксимированы плоскостью, был вычислен угол β между этой плоскостью и плоскостью z=0, и все точки оцифрованной поверхности были повернуты на угол β, а также смещены таким образом, чтобы вершины пирамид лежали в плоскости z=0, а начало координат по осям x и y находилось в центре панели. В результате этого был получен рельеф, представленный на фиг. 7, где для улучшения восприятия изображения шаг между точками увеличен в три раза по сравнению с исходным шагом. Моделирование радиоголограммы поверхности осуществлялось в соответствии с методикой, заключающейся в представлении поверхности набором плотно расположенных всенаправленных точечных рассеивателей. В данном случае рассеиватели располагались в узлах сетки с шагом 2 мм.
Комплексная амплитуда отраженного от поверхности сигнала, зарегистрированного приемником в точке плоскости регистрации с координатами (x, y, z0), вычислялась по формуле
где Е0 - амплитуда возбуждения передатчика; σ - эффективная площадь рассеяния, введенная в формулу для соблюдения размерности; NPS - количество описывающих поверхность точечных рассеивателей; r=[x, y, z0]T - вектор координат приемо-передающей антенны; z0=35 мм; rj=[xj, yj, zj]T - вектор координат j-го точечного рассеивателя; g(f, θ) - диаграмма направленности антенны в направлении на точечный рассеиватель, имеющая для использующейся антенны в виде открытого конца круглого волновода вид ; θ(r, rj) - угол между вектором нормали к плоскости регистрации и направлением из точки приемо-передатчика на j-й точечный рассеиватель; а - радиус волновода; k=2πf/c - волновое число; f - частота сигнала; c - скорость света.
При регистрации данных радиолокатором возникает набег фазы, складывающийся из некоторой начальной фазы и набега фазы в антенне и в разъемах. При вычитании экспериментальной и смоделированной радиоголограмм необходимо обеспечить, чтобы этот набег был одинаков в обеих радиогологрммах. Кроме того, необходимо учесть, что в плоскости регистрации между радиоголограммами может присутствовать остаточный сдвиг, вызванный ошибками при разметке области сканирования и при оцифровке поверхности. Обе эти задачи решались одновременно. Для каждого значения набега фазы ϕ в диапазоне от 0 до 2π с шагом π/100 вычислялась двумерная взаимная корреляционная функция между нормализованными экспериментальной радиоголограммой и смоделированной радиоголограммой, сдвинутой по фазе на ϕ:
Максимальное значение C(ϕ, u, ν) соответствует как искомому набегу фазы ϕ0, который нужно добавить к смоделированной радиоголограмме, так и сдвигу между радиоголограммами в плоскости регистрации (u0, ν0). После этого можно вычислить разностную радиоголограмму, являющуюся радиоголограммой подповерхностного объекта:
Смоделированная радиоголограмма гипсовой панели приведена на фиг. 8. Разностная радиоголограмма, полученная в результате вычитания экспериментальной и смоделированной радиоголограмм, приведена на фиг. 9. Результат ее восстановления приведен на фиг. 10.
Предложенный способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью, позволяет скомпенсировать отражения от поверхности и получить детализированные радиоизображения подповерхностных объектов.
Claims (9)
- Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью, включающий в себя ступенчатое изменение сигнала в заданном диапазоне частот с равномерным шагом в диапазоне
- где:
- kmin=0,72;
- kmax=0,81;
- D - диаметр антенны;
- c - скорость света,
- количество отдельных частот в диапазоне от fmin до fmax не менее пяти, автоматическое выравнивание амплитудно-частотной характеристики, отличающийся тем, что создается цифровая карта рельефа неровной поверхности среды с использованием датчика глубины, позволяющего измерять расстояние, соответствующее каждому пикселю получаемого датчиком глубины оптического изображения, рассчитывается радиоголограмма поверхности с использованием полученного рельефа неровной поверхности среды, рассчитывается разностная радиоголограмма (разность между экспериментальной и расчетной радиоголограммами), по разностной радиоголограмме вычисляется радиоизображение подповерхностного объекта методом обратного распространения.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2018142117A RU2701880C1 (ru) | 2018-11-29 | 2018-11-29 | Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2018142117A RU2701880C1 (ru) | 2018-11-29 | 2018-11-29 | Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2701880C1 true RU2701880C1 (ru) | 2019-10-02 |
Family
ID=68171041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2018142117A RU2701880C1 (ru) | 2018-11-29 | 2018-11-29 | Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2701880C1 (ru) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2000103678A (ru) * | 2000-02-17 | 2002-01-10 | ОАО "Центральный научно-исследовательский институт радиоэлектронных систем" | Способ зондирования конденсированных сред |
US6657577B1 (en) * | 1997-07-02 | 2003-12-02 | Malaa Geoscience Forvaltning Ab | Radar plant and measurement technique for determination of the orientation and the depth of buried objects |
JP2012137461A (ja) * | 2010-12-28 | 2012-07-19 | Oyo Corp | 車載型地中レーダ計測装置のシステム同期方法 |
RU2482518C1 (ru) * | 2011-11-01 | 2013-05-20 | Владимир Всеволодович Разевиг | Способ получения радиоголограмм подповерхностных объектов |
-
2018
- 2018-11-29 RU RU2018142117A patent/RU2701880C1/ru not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6657577B1 (en) * | 1997-07-02 | 2003-12-02 | Malaa Geoscience Forvaltning Ab | Radar plant and measurement technique for determination of the orientation and the depth of buried objects |
RU2000103678A (ru) * | 2000-02-17 | 2002-01-10 | ОАО "Центральный научно-исследовательский институт радиоэлектронных систем" | Способ зондирования конденсированных сред |
JP2012137461A (ja) * | 2010-12-28 | 2012-07-19 | Oyo Corp | 車載型地中レーダ計測装置のシステム同期方法 |
RU2482518C1 (ru) * | 2011-11-01 | 2013-05-20 | Владимир Всеволодович Разевиг | Способ получения радиоголограмм подповерхностных объектов |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Johnson et al. | The geological interpretation of side‐scan sonar | |
Godin et al. | An assessment of laser range measurement of marble surfaces | |
US7466628B2 (en) | Method of constructing mathematical representations of objects from reflected sonar signals | |
CN102298156A (zh) | 用于反虚反射地震数据的方法和装置 | |
CN101655357B (zh) | 一种用于二维相位展开的相位梯度相关质量图获取方法 | |
CN110050205B (zh) | 使用偏移道集的潜波照明 | |
Huthwaite et al. | A new regularization technique for limited-view sound-speed imaging | |
KR20200037729A (ko) | 3d gpr을 이용한 터널 콘크리트 구조물 안전성 분석 장치 및 이를 이용한 안전성 분석 방법 | |
Sac et al. | 2D high-frequency forward-looking sonar simulator based on continuous surfaces approach. | |
MX2013015016A (es) | Indicador de repetibilidad 4d basado en iluminacion por disparo par adquision sismica. | |
Choi et al. | Array type miniaturized ultrasonic sensors to detect urban sinkholes | |
Leckebusch et al. | Investigating the true resolution and three‐dimensional capabilities of ground‐penetrating radar data in archaeological surveys: measurements in a sand box | |
RU2701880C1 (ru) | Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью | |
Bell | Application of optical ray tracking techniques to the simulation of sonar images | |
Drira et al. | Image source detection for geoacoustic inversion by the Teager-Kaiser energy operator | |
Richter et al. | Numerical simulation and experimental validation of wave pattern induced coordinate errors in airborne LiDAR bathymetry | |
Goff et al. | Comparison of a stochastic seafloor model with SeaMARC II bathymetry and Sea Beam data near the East Pacific Rise 13°–15° N | |
RU2635332C1 (ru) | Способ определения состояния ледяного покрова | |
Müller et al. | Mapping of bridge pier scour with projection moiré | |
CN114494603A (zh) | 基于Unity3D的仿真声呐图像数据生成方法 | |
Isidoro et al. | Design and characterisation of customised-roughness beds for open-channel flow experiments | |
Moya et al. | Side scan sonar images attributes characterization for seabed mapping | |
Ogushwitz | Realistic computer simulation of echo-sounding | |
Smith Neto et al. | The Shadow Moiré method using the phase shifting technique and digital image processing: computational implementation and application to the 3D-reconstruction of a buckled plate | |
Walker et al. | Noniterative techniques for GPR imaging through a nonplanar air-ground interface |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20201130 |