Способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора



Способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора
Способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора
Способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора
Способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора
Способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора
Способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора
Способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора

 


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

Открытое акционерное общество "Государственный научно-исследовательский навигационно-гидрографический институт" (ОАО "ГНИНГИ") (RU)

Изобретение относится к области гидроакустики. Сущность: способ определения глубин в реальном масштабе времени при обследовании рельефа дна гидролокатором бокового обзора с последующим его восстановлением, включающий измерения времени задержки синфазных сигналов донной реверберации, принимаемых двумя антеннами, разнесенными по вертикали на несколько длин волн упругих колебаний, и разрешение неоднозначности измерений, вычисление глубин, в котором для достижения технического результата при каждом совпадении фаз интерферирующих сигналов регистрируют мгновенное значение частоты сигнала в нижнем канале, измеряют время запаздывания появления сигнала в верхнем канале с тем же значением мгновенной частоты, измеренное значение времени запаздывания умножают на значение рабочей частоты интерферометра, определяют порядковую нумерацию ряда измерений задержки прихода синфазных сигналов в период каждого зондирования в реальном масштабе времени, глубины вычисляют, соответствующие каждой интерференционной полосе, а при последующем восстановлении рельефа дна по измеренным глубинам выполняют оценку репрезентативности (значимости) критических точек рельефа путем представления гладкой непрерывной поверхности рельефа дна деревом Кронрода-Риба. Гидролокатор бокового обзора включает генератор 1 зондирующих импульсов, приемо-передающую верхнюю антенну 2, схему 3 измерения времени tn, приемо-передающую нижнюю антенну 4, интерферометр 5, схему 6 измерения времени Δtn, селектор 7, частотный детектор 8, опорный генератор 9, запоминающее устройство 10, схему 11 сравнения амплитуд, частотный детектор 12, вычислитель 13, антенный коммутатор 14, устройство 15 отображения и документирования. Технический результат: повышение достоверности определения глубин посредством гидролокатора бокового обзора и последующего восстановления рельефа дна по измеренным глубинам посредством гидролокатора бокового обзора. 7 ил., 2 табл.

 

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

Известен способ определения глубин в стороне от судна посредством гидролокатора бокового обзора (интерферометра) путем измерения дальностей и направлений на элементы отражения донного рельефа на основе использования явления интерференции сигналов донной реверберации, принимаемых на две антенны, разнесенные по вертикали на величину нескольких длин волн упругих колебаний, возбуждаемых одной из антенн интерферометра. В результате выделения и фиксации времени задержки прихода на обе приемные антенны ряда синфазных сигналов, например, на самописце, формируется семейство интерференционных полос, содержащих информацию о направлениях прихода ряда синфазных сигналов и соответствующих им наклонных дальностей до дна (О некоторых особенностях обработки интерференционных эхограмм//Записки по гидрографии. Л.: ГУНиО МО СССР, №203, 1979, с.11-16, [1]. Тарасюк Ю.Ф. Дистанционный промер глубин//Судостроение за рубежом, 1975, №1, с.104-105, [2]. Наумов Е.А. Использование явления интерференции звуковых волн для определения угла наклона дна//Судовождение. Сборник научных трудов ЛВИМУ им. С.О.Макарова, 1974, вып.14, с.252-255, [3]).

Глубины Zn и соответствующие им горизонтальные отстояния Yn от антенн интерферометра вычисляются по формулам [1]:

Z n = D n sin Θ n = C t n n / 2 N = D n n / N ,                                      ( 1 )

Y n = D n cos Θ n = D n 1 ( n / N ) 2 ,                                             ( 2 )

где Θn - углы от горизонта до направления прихода синфазных сигналов;

Dn - наклонные дальности до дна в направлении прихода ряда синфазных сигналов;

С - скорость распространения звука в воде;

n - порядковый номер интерференционной полосы;

N=b/λ - максимальное число ряда регистрируемых интерференционных полос;

λ - длина волны акустических колебаний;

b - величина базового разноса антенн по вертикали;

tn - время задержки прихода ряда синфазных сигналов.

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

Известен также способ определения номеров интерференционных полос, основанный на визуальном дешифровании предвычисленных совпадений полос с заданными номерами из двух семейств интерференционных полос, полученных, например, при использовании двух разных, но близких по величине антенных баз (О некоторых особенностях обработки интерференционных эхограмм//Записки по гидрографии. Л.: ГУНиО МО СССР, №203, 1979, с.11-16, [1]).

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

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

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

Новые отличительные признаки, заключающиеся в том, что при каждом совпадении фаз интерферирующих сигналов регистрируют мгновенное значение частоты сигнала в нижнем канале, измеряют время запаздывания появления сигнала в верхнем канале с тем же значением мгновенной частоты, измеренное значение времени запаздывания умножают на значение рабочей частоты интерферометра, определяют порядковую нумерацию ряда измерений задержки прихода синфазных сигналов в период каждого зондирования в реальном масштабе времени, глубины вычисляют, соответствующие каждой интерференционной полосе, а при последующем восстановлении рельефа дна по измеренным глубинам выполняют оценку репрезентативности (значимости) критических точек рельефа путем представления гладкой непрерывной поверхности рельефа дна деревом Кронрода-Риба, позволяют достичь поставленного технического результата путем выполнения дополнительных измерений времени относительного запаздывания прихода сигналов с одинаковыми мгновенными частотами на разнесенные антенны (Δtn), при этом относительное запаздывание сигналов всегда будет наблюдаться в канале верхней антенны.

Для повышения достоверности измерений величины Δtn в условиях помех они производятся в пределах строба длительностью τ=b/с, а результаты измерений по каждой интерференционной полосе накапливаются за несколько циклов зондирования и осредняются. Повышение точности измерений величины Δtn обеспечивается тем, что в момент выделения синфазных сигналов фиксируется мгновенная частота сигнала на нижней антенне в момент формирования n-й полосы и относительное запаздывание определяется по времени запаздывания прихода сигнала той же частоты на верхнюю антенну.

Измерение времени запаздывания осуществляется по каждой полосе. Основанием к этому служит известное теоретическое положение о том, что сигнал донной реверберации на входе антенны представляет собой узкополосный гауссовский шум с нормальным законом распределения мгновенной частоты в пределах спектра сигнала. Ширина спектра сигнала донной реверберации примерно равна ширине спектра излученного импульса. Плотность вероятности мгновенных значений частоты такого сигнала на выходе предварительного усилителя приемника с высокой степенью точности соответствует также нормальному закону распределения (Ольшевский В.В. Статистические методы в гидролокации. Л.: Судостроение, 1983, с.100-130).

При этом ширина спектра у сигналов донной реверберации, принимаемых интерферометром, зависит от длительности излучаемых импульсов и лежит в пределах 1500÷500 Гц. Мгновенная частота сигнала принимает значение fT±(0÷750)·10-3 кГц, распределенное по нормальному закону. Дисперсия процесса определяется соотношением сигнал/шум.

При условии, что расстояние между антеннами значительно меньше интервала корреляции сигналов (обычно b<50λ), можно утверждать, что разность хода синфазных сигналов пропорциональна времени запаздывания появления в верхнем канале сигнала с мгновенной частотой, фиксированной в нижнем канале на момент стробирования. Очевидно, что Δ t n = b / c sin Θ n .                                                       ( 3 )

Подставим в (3) значение sin Θn из формулы (1), получим Δ t n = b n / c N .                                                                                           ( 4 )

Отсюда n = F 0 Δ t n ,                                                                      ( 5 )

где F0 - рабочая частота интерферометра.

В случае получения дробных значений n из-за погрешности измерений полученные числа округляются до ближайшего целого.

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

Сущность предлагаемого технического решения поясняется чертежами.

Фиг.1. Блок-схема гидролокатора бокового обзора. Гидролокатор бокового обзора включает генератор 1 зондирующих импульсов, приемо-передающую верхнюю антенну 2, схему 3 измерения времени tn, приемо-передающую нижнюю антенну 4, интерферометр 5, схему 6 измерения времени Δtn, селектор 7, частотный детектор 8, опорный генератор 9, запоминающее устройство 10, схему 11 сравнения амплитуд, частотный детектор 12, вычислитель 13, антенный коммутатор 14, устройство 15 отображения и документирования.

Фиг.2. Пример получения интервалов устойчивости (persistence intervals) для одномерной кривой. 16 - уровни сечения, 17 - отрезки сечения, 18 - интервалы устойчивости.

Фиг.3. График поверхности: а) - соответствующие изолинии высот и дерева Кронрода-Риба; б) - 18 критических точек.

Фиг.4. Вертикальное представление для поверхности дерева Кронрода-Риба фиг.2а. Первая цифра в прямоугольнике - номер критической точки на фиг.26, вторая - высота критической точки.

Фиг.5. Дерево Кронрода-Риба для репрезентативных точек поверхности, изображенной на фиг.2а.

Фиг.6. Дерево Кронрода-Риба для зашумленной поверхности, изображенной на фиг.2а. (88 критических точек).

Фиг.7. Дерево Кронрода-Риба для поверхности, изображенной на фиг.2а, и дерево Кронрода-Риба для зашумленной поверхности.

Предлагаемый способ реализуется посредством гидролокатора бокового обзора (фиг.1). С генератора 1 зондирующих импульсов в момент времени t0 импульс с несущей частотой F0 через антенный коммутатор 14 поступает на приемо-передающую антенну 2. Одновременно передним фронтом этого импульса запускается схема 3 измерения времени tn. За период приема сигналов донной реверберации антеннами 2 и 4, в моменты совпадения фаз сигналов на этих антеннах с выхода интерферометра 5 на схему 3 поступает импульс для снятия очередного отсчета времени tn, на запуск схемы 6 измерения времени Δtn и на селектор 7. В селекторе 7 в момент tn производится выборка значения напряжения сигнала, поступающего с частотного детектора 8, при этом значение выборки пропорционально уходу мгновенной частоты сигнала от частоты опорного генератора 9. Значение выборки напряжения сигнала мгновенной частоты фиксируется в запоминающем устройстве 10 в своей ячейке памяти для каждой интерференционной полосы и подается на схему 11 сравнения амплитуд. В момент появления на втором входе схемы 11 сравнения амплитуд такого же значения напряжения сигнала, вырабатываемого частотным детектором 12, с выхода схемы 11 сравнения амплитуд на схему 6 поступает импульс для снятия отсчета времени Δtn.

Измеренные значения tn и соответствующие им значения Δtn подаются на вычислитель 13, реализованный на базе микроЭВМ, где последовательно вычисляется заданный ряд величин Dn и n и соответствующие им значения Zn и Yn. Значения Zn и Yn, увязанные со временем, поступают на устройство отображения и документирования 15.

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

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

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

Так как величины абсолютной высоты (глубины) критических точек не позволяют судить, какие из них более репрезентативные, а какие нет, за исключением только точек с наибольшим значением минимума и максимума. Для определения понятия значимости критической точки поверхности в качестве подходящей основы воспользуемся широко используемым в зарубежной картографии понятием "топографической значимости" (topographic prominence) (Christopherson G.L. Using ARC/GRID to Calculate Topographic Prominence in an Archaeological Landscape. // Arc/INFO User Conference, 2003. - 15 pp. Podobnikar T. Method for Determination of the Mountain Peaks // 12th AGILE International Conference on Geographic Information Science, Leibniz Universitat Hannover, Germany, 2009, p.1-8).

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

Вместе с тем, в отечественной географии существует схожее по смыслу понятие - "относительная высота", однако оно не совпадает с понятием "топографическая значимость". Относительная высота - это топографическое превышение какой-либо точки земной поверхности относительно другой точки, отсчитываемое по вертикали, равное разности абсолютных высот этих точек (например, высота горной вершины над уровнем дна ближайшей долины); расстояние по вертикали от указанного исходного уровня до уровня, точки или объекта, принятого за точку. Следовательно, понятие относительная высота имеет разный смысл в зависимости от контекста, что не позволяет им воспользоваться. Эта проблема может быть решена, если обобщить понятие топографической значимости путем включения в него локальных минимумов рельефа путем построения однозначного метода выявления для каждой седловой точки соответствующей ей точки минимума или максимума. Главная проблема на этом пути состоит в определении алгоритма выявления соответствия между седловыми точками и точками экстремумов. Эффективный метод оценки значимости можно получить, если использовать инструменты описания гладких функций, к которым относится картографическое представление поверхности рельефа (Жуков Ю.Н Математические инструменты описания картографического отображения рельефа Земли // Навигация и океанография. 2011, №32, стр.60-69).

Картографическое представление рельефа есть аналог математического объекта - невырожденной функции Морса. Такая функция имеет только простые критические точки:

седла, минимумы и максимумы, ее топологические свойства описаны, например, в работе (Милнор Дж. Теория Морса. - М.: Издательство ЛКИ, 2011. - 184 с.). Для функций этого типа разработаны методы выявления и упорядочивания критических точек поверхности с использованием топологических характеристик. В вычислительной топологии этот метод носит название топологической значимости (topological persistence) (Bauer U. Persistence in discrete Morse theory. Dissertation zur Eriangung des mathematisch-naturwissenschaftlichen Doktorgrades Doctor rerum naturalium der Georg-August-Universitat GOttingen, 2011. - 109 p.).

Естественно воспользоваться аналогом этого метода для решения нашей задачи нахождения соответствия между седловыми точками и точками экстремумов поверхности.

Не загружая текст математическими подробностями и деталями, приведем неформальное описание метода, адаптированное к нашей задаче, на примере получения интервалов устойчивости (persistence intervals) для одномерной кривой (фиг.2). Пусть дана гладкая одномерная кривая. Возьмем горизонтальную прямую, секущую данную кривую. На некотором уровне сечение кривой представляет собой набор несвязных горизонтальных интервалов (отрезков). Грубо говоря, целое число интервалов определяет топологический тип сечения на данном уровне. Рассмотрим, как изменяется топологический тип сечений кривой на каждом уровне при перемещении секущей прямой от минимального уровня кривой до наибольшего. При увеличении уровня сечения число интервалов меняется только в момент прохождении секущей через точку, соответствующую либо максимуму, либо минимуму. В диапазоне уровней между двумя последовательными экстремумами число интервалов сечения не меняется. Поэтому динамика числа интервалов сечения в процессе увеличения уровня сечения изменяется в соответствии со следующим простым правилом: если встречается минимум, то число интервалов увеличивается на единицу, если встречается максимум, то число интервалов уменьшается на единицу. В этом процессе всегда есть минимум, который создает интервал, и существует некоторый максимум, который этот интервал уничтожает. Более конкретно: текущий максимум является парой последнему предыдущему встреченному минимуму. Последовательные максимум и минимум являются сопряженной парой экстремумов. Сопряжение пары критических точек обусловлено собственно топологическими свойствами рассматриваемой кривой.

Указанная процедура позволяет однозначно получить множество пар минимумов-максимумов, определяемых топологией кривой. Величина абсолютной разности высот критических точек Δ=|hmin-hmax|, составляющих сопряженную пару, определяет количественную оценку интервала устойчивости и соответствующих сопряженных точек. Полученные указанным методом интервалы устойчивости являются искомыми объектами для одномерного случая.

Приведенный пример для одномерного случая дает основание ввести понятие "значимость" критической точки поверхности рельефа как меры топологической устойчивости той сопряженной пары критических точек, в которую она входит. Количественный критерий значимости ε определим следующим образом. Множеству пар критических точек П сопоставим множество разностей высот критических точек П={Δi} (i=1,…,N, N - число пар). Значимость конкретной пары критических точек в данном множестве пар П определим как отношение εiimax - максимальное значение среди {Δi}. Таким образом, каждая точка в паре имеет одну и ту же величину значимости. Величина ε всегда нормирована и лежит в диапазоне 0÷1.

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

В двумерном случае пара критических точек, образующих интервал устойчивости, всегда состоит либо из седла и минимума, либо из седла и максимума. Каждой такой паре критических точек сопоставим величину абсолютной разности высот критических точек Δ=|hs-he|, составляющих пару. Здесь hs - высота седловой точки, he - высота точки экстремума, минимума или максимума.

Для вычисления множества пар критических точек для двумерных гладких поверхностей обычно применяют вспомогательный математический аппарат - дерево Кронрода-Риба (дКР). Картографическое описание рельефа можно представить математическим объектом - функцией Морса, для которой дКР является дискретным аналогом, однозначно описывающим положение, высоту точек седел и локальных экстремумов, и главное, дКР описывает связь между ними, которая определяется топологией рассматриваемой поверхности. Каждой поверхности однозначно соответствует некоторое дКР. К настоящему времени разработаны алгоритмы вычисления дКР для всех типов представления поверхности в ЭВМ (Doraiswamy H., Natarajan V. Efficient Algorithms for Computing Reeb Graphs // Computational Geometry: Theory and Applications, Volume 42, Issue 6-7, August, 2009, p.606-616. Kunii Т.L. Constructing a reeb graph automatically from cross sections. // IEEE Comput. Graph. Appl. 11,6 (1991), 44-51. Pascucci V. Loops in Reeb graphs of 2-manifolds. // Discrete and Computational Geometry 32, 2 (2004), 231-244), и поэтому на процедуре вычисления дКР останавливаться не будем. Отметим важное практическое обстоятельство: вершины дКР естественным образом оснащаются координатами местоположения и высотой соответствующих критических точек.

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

Алгоритм вычисления пар критических точек по дКР очевиден, но достаточно громоздок, чтобы его приводить полностью. В алгоритме используется представление, что точки дКР, соответствующие седлам, являются точками тройного ветвления, а точки, соответствующие минимумам и максимумам, являются точками с одним примыкающим ребром. Приведем только основную схему алгоритма. Ключ алгоритма состоит в использовании динамически меняющегося дКР в соответствии с последовательностью выявляемых пар критических точек. Последовательно просматривается снизу вверх список упорядоченных по высоте критических точек, начиная со второй снизу. Если точка является минимумом или максимумом, то ищется ближайшее к ней в дКР седло. Они образуют пару седло-экстремум. Вершины экстремумов, седел и связывающее их ребро удаляются из дКР, а нарушенные удалением связи восстанавливаются в соответствии с порядком вершин до удаления. Если вершина является седлом и имеет двух нижележащих потомков, то выбирается ближайший по высоте минимум, и они образуют пару критических точек. Соответствующие вершины и ребра удаляются из дКР с последующим восстановлением связей. В результате такого процесса окончательное дКР будет иметь только две вершины, соответствующие наибольшим минимуму и максимуму, и одно соединяющее их ребро. Эта пара критических точек образует наиболее значимую пару с Δmax.

Заметим, что приведенный метод оценки значимости форм рельефа включает понятие топографической значимости, но в отличие от последнего дает однозначный алгоритм получения пар точек седло-экстремум.

Для поверхности, представленной на фиг.3, в соответствии с ее дКР вычислены пары критических точек и их значимость. Результаты представлены в таблице 1.

Таблица 1

Таблица значимости критических точек для поверхности, изображенной на фиг.3.

Номер точки соответствует номеру критических точек на фиг.3б.

Порядковый номер появления в списке пар Тип точки Номер точки Номер точки Тип точки Величина интервала устойчивости Величина значимости пары Порядковый номер значимости пары
1 Седло 18 12 Макс. 1,28 0,08796 4
2 Седло 8 4 Макс. 2,78 0,19086 3
3 Мин. 1 3 Седло 0,0007 0,00003 7
4 Мин. 17 16 Седло 0,03 0,00218 6
5 Мин. 6 5 Седло 0,0001 0,000008 8
6 Седло 14 13 Макс. 0,0000 0,0000003 9
7 Мин. 15 11 Седло 0,34 0,02355 5
8 Мин. 7 2 Седло 3,03 0,20785 2
9 Мин. 10 9 Макс. 14,57 1 1

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

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

Введем случайный равномерный на интервале [-0.5 0.5] шум в высотные координаты функции, изображенной на фиг.3а. Результаты построения изолиний такой поверхности и вычисленное для нее дКР представлены на фиг.6.

Из последнего примера следует, что предложенный метод чувствителен к шумовой составляющей. Поэтому он может быть эффективным инструментом в алгоритме предварительного сглаживания поверхности рельефа дна, полученной по измеренным глубинам. Кроме того, приведенные примеры показывают, что подобие двух поверхностей одного участка рельефа можно в первом приближении установить, используя значимые критические точки. Если значимость и пространственные координаты имеют достаточно близкие (с точки зрения решаемой задачи) значения, то соответствующие поверхности подобны. Другие, более тонкие, методы сравнения следует применять после проверки этого соответствия. Для иллюстрации этого утверждения на фиг.6 отображены два дКР для репрезентативных критических точек поверхностей, изображенных на фиг.3 и фиг.6, а в таблице 2 представлены величины значимости для репрезентативных точек этих поверхностей. Последние незначительно меняются под действием шума (таблица 2).

Таблица 2
Величины интервала устойчивости и значимости для репрезентативных критических точек
Номер точки Величина интервала устойчивости критической точки Величина интервала устойчивости критической точки Величина значимости крити
ческой точки
Величина значимости крити
ческой точки
Без шума С шумом Без шума С шумом
1 14,57 14,64 1 1
2 14,57 14,64 1 1
3 3,03 3,03 0,20786 0,20691
4 3,03 3,03 0,20786 0,20691
5 2,78 2,67 0,19086 0,18200
6 2,78 2,67 0,19086 0,18200
7 1,28 1,30 0,08796 0,08897
8 1,280 1,30 0,08796 0,08897

Определение репрезентативных критических точек представляет практическую ценность с точки зрения выявления шумовых возмущений в поверхности и сравнения двух поверхностей между собой. Дерево Кронрода-Риба является основным инструментом и структурой представления топологических и геометрических особенностей репрезентативных критических точек поверхности рельефа. Проблема выявления репрезентативности (значимости) критических точек рельефа с помощью ЭВМ решается путем использования представления гладкой непрерывной поверхности рельефа дна деревом Кронрода-Риба. Такая параметризация поверхности позволяет выявить значимые критические точки поверхности с учетом ее топологической структуры. Данный метод позволяет выявить шумовые возмущения в данных о поверхности и осуществить сравнение топологических свойств двух или нескольких поверхностей, полученных по результатам гидроакустических измерений. Совокупность общих с прототипом и отличительных признаков обеспечивает появление новых свойств у предлагаемого способа, а именно:

- повышение достоверности определения порядковой нумерации ряда наклонных расстояний, соответствующих интерференционным максимумам;

- упрощен процесс разрешения неоднозначности измерений;

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

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

1. О некоторых особенностях обработки интерференционных эхограмм//Записки по гидрографии. Л.: ГУНиО МО СССР, №203,1979, с.11-16.

2. Тарасюк Ю.Ф. Дистанционный промер глубин//Судостроение за рубежом, 1975, №1, с.104-105.

3. Наумов Е.А. Использование явления интерференции звуковых волн для определения угла наклона дна//Судовождение. Сборник научных трудов ЛВИМУ им. С.О.Макарова, 1974, вып.14, с.252-255.

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



 

Похожие патенты:

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

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

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

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

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

Изобретение относится к области гидрографии и может быть использовано для стереосъемки рельефа дна акватории гидроакустическим средством (ГАС), а также поиска подводных объектов, расположенных на поверхности дна акватории.

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

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

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

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

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

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

Изобретение относится к области гидрографии и может быть использовано для изучения поля рельефа дна в Мировом океане в целях навигационного гидрографического обеспечения флота и народного хозяйства. Способ съемки рельефа дна акватории, согласно изобретению, включает вождение судна, на котором установлен эхолот, по запланированным галсам, излучение гидроакустических сигналов в сторону дна, прием отраженных от поверхности дна этих сигналов эхолотом, измерение расстояния гидроакустическим путем от приемоизлучающей антенны эхолота до места отражающей поверхности дна акватории, определение геодезических координат приемоизлучающей антенны эхолота, измерение параметров бортовой и вертикальной качек истинного курса, скорости движения судна, в начале каждого галса и в конце проходят судном через место расположения опорного гидрографического пункта (ОГП) с известными глубиной и его геодезическими координатами, измеренные в моменты прохода через ОГП, и ее геодезические координаты фиксируются, а истинные значения глубин и их геодезические координаты в пунктах, расположенных вдоль галса, определяют вычислительным путем по установленным новым формулам, приведенным в описании изобретения. Эхолот для осуществления данного способа съемки рельефа дна акватории содержит приемоизлучающую антенну, передающий блок, приемоизмерительный блок, блок управления, а также введенный блок автоматического определения истинных глубин и их геодезических координат, свободных от случайных и систематических погрешностей, вычислительным путем по установленным формулам, и регистратор, при этом выход определяющего и вход приемоизмерительного блоков соединены соответственно через блок управления с входом приемоизлучающей антенны, а выходы приемоизмерительного блока и навигационного комплекса через блок управления соединены с входом блока автоматического определения истинных глубин и их геодезических координат, свободных от случайных и систематических погрешностей, выход которого соединен с входом регистратора. Технический результат: обеспечение автоматического исключения случайных систематических погрешностей при съёмке дна акватории эхолотом. 1 з.п. ф-лы, 2 ил.

Изобретение относится к диагностическим ультразвуковым системам для трехмерной визуализации. Ультразвуковая диагностическая система визуализации содержит ультразвуковой датчик, выполненный с возможностью сбора набора данных 3-мерного изображения объемной области, блок мультипланарного переформатирования, реагирующий на набор данных 3-мерного изображения, выполненный с возможностью формирования множества 2-мерных изображений, блок задания последовательности изображений, реагирующий на 2-мерные изображения, выполненный с возможностью формирования последовательности 2-мерных изображений, которые могут быть воспроизведены в виде последовательности 2-мерных изображений стандартного формата, порт данных, связанный с блоком задания последовательности изображений, выполненный с возможностью передачи последовательности 2-мерных изображений в другую систему визуализации, и дисплей просмотра последовательностей 2-мерных изображений. Система визуализации дополнительно содержит пользовательский интерфейс управления для выбора нормального направления через набор 3-мерных данных, который содержит выбор плоскости 2-мерного изображения, проходящей через набор 3-мерных данных, причем изображения последовательности 2-мерных изображений, сформированных блоком переформатирования данных изображения, параллельны плоскости выбранной плоскости 2-мерного изображения. Использование изобретения позволяет облегчить перенос и использование данных 3-мерного изображения на других платформах для медицинских изображений. 11 з.п. ф-лы, 4 ил.
Наверх