JPS642902B2 - - Google Patents

Info

Publication number
JPS642902B2
JPS642902B2 JP5595882A JP5595882A JPS642902B2 JP S642902 B2 JPS642902 B2 JP S642902B2 JP 5595882 A JP5595882 A JP 5595882A JP 5595882 A JP5595882 A JP 5595882A JP S642902 B2 JPS642902 B2 JP S642902B2
Authority
JP
Japan
Prior art keywords
data
value
frequency
signal strength
output
Prior art date
Legal status (The legal status 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 status listed.)
Expired
Application number
JP5595882A
Other languages
English (en)
Other versions
JPS58173473A (ja
Inventor
Shunichi Kohama
Hisazumi Fukuoka
Toshio Chiba
Takashi Matsumoto
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Oki Electric Industry Co Ltd
Original Assignee
Oki Electric Industry Co Ltd
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 Oki Electric Industry Co Ltd filed Critical Oki Electric Industry Co Ltd
Priority to JP5595882A priority Critical patent/JPS58173473A/ja
Publication of JPS58173473A publication Critical patent/JPS58173473A/ja
Publication of JPS642902B2 publication Critical patent/JPS642902B2/ja
Granted legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/56—Display arrangements
    • G01S7/60—Display arrangements for providing a permanent recording

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Description

【発明の詳細な説明】
本発明はソーナー装置から出力される周波数ご
との信号強度データを効果的に圧縮して記録紙等
に濃淡表示するためのソーナー信号表示処理方式
に関するものである。 一般にソーナー装置におけるローフアーグラム
とは、第1図に示すような周波数毎の受信信号強
度の分析結果を、第2図のように記録紙1上で周
波数を横軸、時間を縦軸にとつて継続的に濃淡表
示するものをいう。この表示面からスペクトル2
の存在とその特徴を見いだすのが目的であるが、
濃淡表示のダイナミツクレンジが分析結果のそれ
より小さいため、適切なダイナミツクレンジ圧縮
処理が必要とされる。従来はこの処理を周波数に
関する算術移動平均を用いた次のアルゴリズムを
用いて実行していた。 ここでGはゲイン定数、xo周波数nに対する信
号強度データ、1/2K+1o+K 〓i=n-K xiは算術平均値、yo はxoに対応する表示出力である。また、Kは整定
数であつて移動平均をとる範囲を決めるパラメー
タである。第3図は、上記式1の演算処理方法の
説明図で、周波数分析器3から出力される信号強
度データ列のうち、周波数nを中心とするxo-Kか
らxo+Kまでのデータがそれぞれ遅延回路41、42、
……42K+1によつて所定時間Tだけ遅延して加算
され、更に1/2K+1が乗算器5で乗算され周波数 nを中心とする算術平均値1/2K+1o+K 〓i=n-K xiが求め られ、別に遅延回路6で時間KTでけ遅延された
データxoを割算器7にて上記算術平均値で割り、
ゲイン定数Gを乗算器8で乗じて出力yoを得る如
くなつている。この方法によれば雑音レベルが帯
域によつて大きく異なる場合でも、その帯域の平
均値で各データを規格化しているため、出力yoは
全帯域にわたつて一様な分散をもつ背景雑音とし
て表示されるのである。これを示したのが第4図
であつて、各データxoは同図aに示す如くその周
辺のデータの平均値によつて規格化されるので、
同図bに示す如く全帯域にわたつて一様な濃度分
散をもつ表示出力が得られている。 ところが上記従来の方法は以下に述べるような
欠点をもつていた。 (1) 強いスペクトルが1本あると、その両側では
白抜け現象が起こり、そのため強いスペクトル
に近接して弱いスペクトルが(分析結果とし
て)得られていても表示されないことがある。 第5図がこの例を示したもので、同図aに示
す如く強いスペクトル11の近くに弱いスペク
トル12が検出されているが、これを規格化す
る移動平均値がスペクトル11のために大きく
持上り、そのため弱いスペクトル12の出力y
は非常に小さくなつて同図bに示す如く表示さ
れないのである。これは移動平均値が、本来そ
の帯域の雑音レベルを表わさなければいけない
のに、大きな信号スペクトルを含むことによつ
て雑音レベルよりはるかに大きな値をとつてし
まうことに起因する。従つて強いスペクトル1
1そのものについても対応する移動平均値が雑
音レベルより大きいため本来の強さが強調され
にくくなつている。 (2) 規格化された値 に対して出力yoは yo=Gy′o ……(式3) で示されるように線形である。従つて弱いスペ
クトルを明確に表示することができなかつた。
ここで何らかの非線形操作、たとえば yo=logy′o ……(式4) のような変換をすると、第6図に示すように
y′oに対して出力yoが上に凸である特性となり、
データの小さい所での変動部分は強調されるの
である。 (3) 出力yoを式(1)にもとずいてデイジタル処理に
て求める場合、yoをひとつ求めるのに割算を2
度行なわなければならず、この割算は有限ビツ
ト長で行なわれるためデータが小さいときは非
常に精度が悪かつた。 本発明は上記従来の欠点を除去するため、デー
タの対数値の算術平均が幾何平均の対数値に等し
いことを利用し、移動平均値として幾何平均を用
いて平均値を雑音レベルにより近似し、割算を対
数値の引き算に変換して処理するようにしたもの
で、その目的とするところは記録紙上の白抜け現
象を極力防止するとともに小さいデータに対する
演算精度を向上することにある。以下、図面につ
いて詳細に説明する。 第7図は本発明の第1の実施例の構成図であつ
て、この構成で示される演算を高速プロセツサ
(計算機)によつて実施する場合のデータ処理の
流れが第8図に示されている。さて、まず周波数
分析器(FFTが用いられるデイジタルプロセツ
サー)13が入力信号を分析して信号強度データ
列{xo}を出力する。このデータ列は高速プロセ
ツサがもつメモリ領域Mに入力される(第8図
1)。メモリ領域Mに格納された各データxoに対
し、対数増幅器14で10logxoを計算し同じメモ
リ領域Mに格納する。(第8図2)。この対数値の
データ列のうち、周波数nを中心とする
10logxo-Kから10logxo+Kまでのデータとそれぞれ
遅延回路151,152……152K+1、によつて所
定時間Tだけ遅延して加算し、更に乗算器16で
1/2K+1を乗算してデータの対数値の算術平均 10/2K+1o+K 〓i=n-K logxiすなわち幾何平均の対数値Zn
【式】を計算しメモ リの他の領域Nに格納する(第8図3)。この後
遅延回路17で時間KTだけ遅延させたデータ
10logxoから平均値Zoを引きかつ定数B(=
10logb)を加えて結果をメモリーMに入れる(第
8図4)。このメモリーMの値の中で負であるも
のはすべて0でおきかえ(第8図5)、乗算器1
8でゲイン数Gをかけた後、出力し表示器に送ら
れ記録紙上に表示される。 第9図は上記処理におけるメモリ内へ対数の格
納方法を示すもので、整数の入力xに対しx≧1
の場合、16ビツトのメモリの上位8ビツトに
10logxの整数部を、下位8ビツトに10logxの小
数部を格納する。 以上説明したように第1の実施例では移動平均
として幾何平均が用いられているため、1本の強
いスペクトルがあつても移動平均値は雑音レベル
から従来ほど大きくずれることがなくなつた。こ
の効果は雑音スペクトルのようにランダムに変動
しているデータ列に対して顕著にあらわれ、従来
のような白抜現象がほとんど見られなくなつた。
第10図aは周波数毎の信号強度データの分析結
果、第10図bは第10図aの分析結果に対する
従来の方式による表示結果、第10図cは本発明
方式による表示結果である。第10図aにあらわ
れているスペクトル21,22,23,24,2
5の第10図b,cにおけるあらわれ方及び白抜
現象の差異に注目すると、第10図bではスペク
トル23及び25があらわれておらず明らかに本
発明方式による結果が優れていることがわかる。
また第9図に示したような対数のとり方によつ
て、小さいデータについては有効ビツト数がふえ
た形になるので平均値(幾何平均の対数値)も精
度良く求められている。さらに対数関数は上に凸
な曲線であり、データの小さい範囲における変動
が強調されているため、弱いスペクトルに対する
S/N改善がなされている。 第1の実施例は演算をデイジタルプロセツサで
行なつたが第7図の構成にもとずいてこれをアナ
ログ系で実現することもできる。第11図はその
実施例を示すもので図中、31は周波数分析器、
32は対数増幅器、33は平均化回路、34は遅
延回路、35は固定バイアス回路、36は検波回
路、37は増幅器であり、第1の実施例と同様に
割算回路を必要としていない。図中、固定バイア
ス回路中35で固定バイアスBが加算されている
が、これは第1の実施例におけるB=10logbに相
当する電圧を印加する如くなつている。また、平
均化回路33は普通の積分回路すなわちローパス
フイルタで良いが、遅延回路34とゲインの相異
があつてはならない。なお、上記回路の動作は第
1の実施例のときと同様である。 更に第7図の構成にもとずき、移動平均をとる
回路を第7図の構成図そのままにトランスバーサ
ルフイルタを用いても実施できる。 以上説明したように本発明によれば、データの
対数値の算術平均をとることにより平均値として
幾何平均を用いることができるため、信号中に強
いスペクトルが存在していても平均値が雑音レベ
ルから大きくはずれることがなく、従つて、強い
スペクトルに近接した弱いスペクトルを欠落させ
ることなく信号強度データのダイナミツクレンジ
を圧縮でき、これを記録紙等に濃淡表示させた時
の白抜現象を抑えることができ、また、信号強度
データに対して対数圧縮を行なうことができるた
め、弱いスペクトルのレベルが上がり、その分、
明確に表示させることができる。また表示媒体と
しては記録紙上の濃淡表示に限らずブラウン管に
よるカラー表示にも適用できる。
【図面の簡単な説明】
図面は本発明の説明に供するもので、第1図は
ソーナー装置より出力される周波数に対する信号
強度を表わすグラフ、第2図は第1図のグラフに
対応するローフアグラム、第3図は従来のダイナ
ミツクレンジの圧縮処理方法の説明図、第4図
a,bは雑音レベルの濃淡表示の説明図、第5図
a,bは白抜現象の説明図、第6図は入出力特性
の比較の説明図、第7図は本発明の第1の実施例
の構成図、第8図はデータ処理の流れを表わす説
明図、第9図はメモリ内への対数の格納方法の説
明図、第10図aは周波数に対する信号強度デー
タの一例を示すグラフ、第10図bは第10図a
のグラフに対する従来の方式によるローフアグラ
ム、第10図cは第10図aのグラフに対する本
発明方式によるローフアグラム、第11図は本発
明の第2の実施例を示す構成図である。 1……記録紙、2,11,12,21,22,
23,24,25……スペクトル、13……周波
数分析器、14……対数増幅器、151,152,
…152K+1,17……遅延回路、16,18……
乗算器。

Claims (1)

    【特許請求の範囲】
  1. 1 ソーナー装置から出力される周波数ごとの信
    号強度データを対数圧縮し、該対数圧縮したデー
    タを周波数軸上の所定の範囲で算術平均し、その
    範囲内の上記対数圧縮した信号強度データから該
    平均値を差し引き、更に該差し引きした値に定数
    を加算し、その結果が正の値の場合はそのままと
    し負の値の場合は0に置き換えたことを特徴とす
    るソーナー信号表示処理方式。
JP5595882A 1982-04-06 1982-04-06 ソ−ナ−信号表示処理方式 Granted JPS58173473A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP5595882A JPS58173473A (ja) 1982-04-06 1982-04-06 ソ−ナ−信号表示処理方式

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP5595882A JPS58173473A (ja) 1982-04-06 1982-04-06 ソ−ナ−信号表示処理方式

Publications (2)

Publication Number Publication Date
JPS58173473A JPS58173473A (ja) 1983-10-12
JPS642902B2 true JPS642902B2 (ja) 1989-01-19

Family

ID=13013581

Family Applications (1)

Application Number Title Priority Date Filing Date
JP5595882A Granted JPS58173473A (ja) 1982-04-06 1982-04-06 ソ−ナ−信号表示処理方式

Country Status (1)

Country Link
JP (1) JPS58173473A (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61235778A (ja) * 1985-04-11 1986-10-21 Suzuki Giyotan Kk 魚群探知機における魚群表示方法
CN114458306B (zh) * 2020-11-06 2024-10-29 中国石油天然气集团有限公司 基于噪声测井的流体流量的确定方法、装置、设备及介质

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS557554A (en) * 1978-06-30 1980-01-19 Murata Mfg Co Ltd Dielectric thin film

Also Published As

Publication number Publication date
JPS58173473A (ja) 1983-10-12

Similar Documents

Publication Publication Date Title
US7386183B2 (en) Method and apparatus for image processing
JP3498842B2 (ja) 実時間人間視覚システム行動モデル化方法
KR0176601B1 (ko) 저역 필터링과 히스토그램 등화를 이용한 화질개선 방법 및 그 회로
KR960006477A (ko) 에지강조를 이루는 오차확산방식 이치화 방법 및 장치
KR20020026119A (ko) 영상신호의 콘트라스트 향상 장치
US5930402A (en) Method and device for local contrast enhancement of video signal
US6289367B1 (en) Digital signal processing circuits, systems, and method implementing approximations for logarithm and inverse logarithm
US8565549B2 (en) Image contrast enhancement
Chen et al. Two-pass bilateral smooth filtering for remote sensing imagery
JPH05176263A (ja) 階調補正装置
US20060039624A1 (en) System and method for fuzzy filtering images
Trivedi et al. A new contrast measurement index based on logarithmic image processing model
US6700626B2 (en) Video-apparatus with peaking filter
US12033649B2 (en) Noise floor estimation and noise reduction
US7668898B2 (en) Calculating circuit and method for computing an N-th root and a reciprocal of a number
Jernigan et al. Nonlinear lateral inhibition and image processing
CN112671412A (zh) 一种雷达回波数据压缩方法及装置
JPS58173473A (ja) ソ−ナ−信号表示処理方式
US6882735B2 (en) Dynamic range compression of an audio signal
JP2008500757A (ja) 映像信号の鮮鋭度を強調するための方法およびシステム
CN100473175C (zh) 确定调制信号的包络线的方法
US3717812A (en) Real time analysis of waves
US20220051384A1 (en) Scaled psnr for image quality assessment
JPH07131758A (ja) ビデオ信号の雑音低減方法
US20230305590A1 (en) A data processing method