RU2701880C1 - Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью - Google Patents

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

Info

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
Application number
RU2018142117A
Other languages
English (en)
Inventor
Владимир Всеволодович Разевиг
Андрей Викторович Журавлев
Сергей Иванович Ивашов
Александр Иванович Ивашов
Original Assignee
Владимир Всеволодович Разевиг
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Владимир Всеволодович Разевиг filed Critical Владимир Всеволодович Разевиг
Priority to RU2018142117A priority Critical patent/RU2701880C1/ru
Application granted granted Critical
Publication of RU2701880C1 publication Critical patent/RU2701880C1/ru

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric 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

Использование: для подповерхностной радиолокации. Сущность изобретения заключается в том, что восстановление радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью, включает в себя ступенчатое изменение сигнала в заданном диапазоне частот с равномерным шагом в диапазоне
Figure 00000011
Figure 00000012
где 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) со ступенчатым изменением сигнала в заданном диапазоне частот, при этом частоты зондирующего сигнала выбираются с равномерным шагом в диапазоне
от
Figure 00000001
;
до
Figure 00000002
,
где:
kmin=0,72;
kmax=0,81;
D - диаметр антенны;
с - скорость света, при этом количество отдельных частот в диапазоне от fmin до fmax должно быть не менее пяти.
Недостатком аналогов и прототипа является то, что при работе на средах с неровной поверхностью данный способ не позволяет получить детализированные радиоизображения подповерхностных объектов, так как сильные отражения от неровной поверхности среды маскируют более слабые отражения от подповерхностных объектов, расположенных дальше от плоскости регистрации радиоголограммы.
Перечень фигур, чертежей и иных материалов.
На фиг. 1 представлена гипсовая панель с заложенным под нее объектом.
На фиг. 2 представлен стенд для проведения экспериментов.
На фиг. 3 показана действительная часть радиоголограммы.
На фиг. 4 показан результат восстановления без учета рельефа поверхности.
На фиг. 5 представлена оцифровка рельефа гипсовой панели.
На фиг. 6. представлена карта глубины.
На фиг. 7 представлен рельеф гипсовой панели, полученный по данным видеосенсора Kinect v2.
Фиг. 8. Действительная часть рассчитанной радиоголограммы.
Фиг. 9. Действительная часть разностной радиоголограммы.
Фиг. 10. Результат восстановления разностной радиоголограммы.
Задача.
Техническая задача состоит в устранении указанного недостатка за счет создания цифровой карты рельефа неровной поверхности среды с использованием датчика глубины, позволяющего измерять расстояние, соответствующее каждому пикселю получаемого датчиком глубины оптического изображения, рассчитывается радиоголограмма поверхности с использованием полученного рельефа неровной поверхности среды, рассчитывается разностная радиоголограмма (разность между экспериментальной и расчетной радиоголограммами), по разностной радиоголограмме вычисляется радиоизображение подповерхностного объекта методом обратного распространения.
Отличительные признаки.
1. В отличие от известного способа, включающего в себя ступенчатое изменение сигнала в заданном диапазоне частот с равномерным шагом в диапазоне
от
Figure 00000003
;
до
Figure 00000004
,
где:
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), вычислялась по формуле
Figure 00000005
,
где Е0 - амплитуда возбуждения передатчика; σ - эффективная площадь рассеяния, введенная в формулу для соблюдения размерности; NPS - количество описывающих поверхность точечных рассеивателей; r=[x, y, z0]T - вектор координат приемо-передающей антенны; z0=35 мм; rj=[xj, yj, zj]T - вектор координат j-го точечного рассеивателя; g(f, θ) - диаграмма направленности антенны в направлении на точечный рассеиватель, имеющая для использующейся антенны в виде открытого конца круглого волновода вид
Figure 00000006
; θ(r, rj) - угол между вектором нормали к плоскости регистрации и направлением из точки приемо-передатчика на j-й точечный рассеиватель; а - радиус волновода; k=2πf/c - волновое число; f - частота сигнала; c - скорость света.
При регистрации данных радиолокатором возникает набег фазы, складывающийся из некоторой начальной фазы и набега фазы в антенне и в разъемах. При вычитании экспериментальной и смоделированной радиоголограмм необходимо обеспечить, чтобы этот набег был одинаков в обеих радиогологрммах. Кроме того, необходимо учесть, что в плоскости регистрации между радиоголограммами может присутствовать остаточный сдвиг, вызванный ошибками при разметке области сканирования и при оцифровке поверхности. Обе эти задачи решались одновременно. Для каждого значения набега фазы ϕ в диапазоне от 0 до 2π с шагом π/100 вычислялась двумерная взаимная корреляционная функция между нормализованными экспериментальной радиоголограммой и смоделированной радиоголограммой, сдвинутой по фазе на ϕ:
Figure 00000007
.
Максимальное значение C(ϕ, u, ν) соответствует как искомому набегу фазы ϕ0, который нужно добавить к смоделированной радиоголограмме, так и сдвигу между радиоголограммами в плоскости регистрации (u0, ν0). После этого можно вычислить разностную радиоголограмму, являющуюся радиоголограммой подповерхностного объекта:
Figure 00000008
.
Смоделированная радиоголограмма гипсовой панели приведена на фиг. 8. Разностная радиоголограмма, полученная в результате вычитания экспериментальной и смоделированной радиоголограмм, приведена на фиг. 9. Результат ее восстановления приведен на фиг. 10.
Предложенный способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью, позволяет скомпенсировать отражения от поверхности и получить детализированные радиоизображения подповерхностных объектов.

Claims (9)

  1. Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью, включающий в себя ступенчатое изменение сигнала в заданном диапазоне частот с равномерным шагом в диапазоне
  2. Figure 00000009
  3. Figure 00000010
  4. где:
  5. kmin=0,72;
  6. kmax=0,81;
  7. D - диаметр антенны;
  8. c - скорость света,
  9. количество отдельных частот в диапазоне от fmin до fmax не менее пяти, автоматическое выравнивание амплитудно-частотной характеристики, отличающийся тем, что создается цифровая карта рельефа неровной поверхности среды с использованием датчика глубины, позволяющего измерять расстояние, соответствующее каждому пикселю получаемого датчиком глубины оптического изображения, рассчитывается радиоголограмма поверхности с использованием полученного рельефа неровной поверхности среды, рассчитывается разностная радиоголограмма (разность между экспериментальной и расчетной радиоголограммами), по разностной радиоголограмме вычисляется радиоизображение подповерхностного объекта методом обратного распространения.
RU2018142117A 2018-11-29 2018-11-29 Способ восстановления радиоголограмм подповерхностных объектов, находящихся в средах с неровной поверхностью RU2701880C1 (ru)

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)

* Cited by examiner, † Cited by third party
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 Владимир Всеволодович Разевиг Способ получения радиоголограмм подповерхностных объектов

Patent Citations (4)

* Cited by examiner, † Cited by third party
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