JP2021530059A - 高ダイナミックソナー又はレーダーデータを表示するための方法及び装置 - Google Patents
高ダイナミックソナー又はレーダーデータを表示するための方法及び装置 Download PDFInfo
- Publication number
- JP2021530059A JP2021530059A JP2021500236A JP2021500236A JP2021530059A JP 2021530059 A JP2021530059 A JP 2021530059A JP 2021500236 A JP2021500236 A JP 2021500236A JP 2021500236 A JP2021500236 A JP 2021500236A JP 2021530059 A JP2021530059 A JP 2021530059A
- Authority
- JP
- Japan
- Prior art keywords
- image
- intensity
- color
- matrix
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 86
- 238000009877 rendering Methods 0.000 claims abstract description 50
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000001914 filtration Methods 0.000 claims description 26
- 230000004044 response Effects 0.000 claims description 26
- 230000008859 change Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 5
- 230000002194 synthesizing effect Effects 0.000 claims description 5
- 238000007639 printing Methods 0.000 claims description 3
- 239000003086 colorant Substances 0.000 abstract description 4
- 238000012935 Averaging Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 29
- 238000012800 visualization Methods 0.000 description 19
- 125000001475 halogen functional group Chemical group 0.000 description 18
- 238000002592 echocardiography Methods 0.000 description 16
- 230000000694 effects Effects 0.000 description 16
- 235000019557 luminance Nutrition 0.000 description 11
- 238000013507 mapping Methods 0.000 description 9
- 238000010606 normalization Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004313 glare Effects 0.000 description 3
- 230000008447 perception Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000002829 reductive effect Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 230000000873 masking effect Effects 0.000 description 2
- 238000012634 optical imaging Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 210000001525 retina Anatomy 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 210000004087 cornea Anatomy 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/001—Texturing; Colouring; Generation of texture or colour
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20208—High dynamic range [HDR] image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Processing Or Creating Images (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
Description
− 色による単純な符号化のいわゆるナイーブな解決策において、色による情報の符号化は、ある強度I(x,y)の色を、低い強度が所与の階調(例えば、赤階調)に、中間強度が他の階調(例えば、緑階調)に、且つ最後に高強度が更に他の階調(例えば、青階調)にあるように適宜選択されたカラーマップに関連付けることを含む。出力データO(x,y)は、次式で表すことができる。
O(x,y)=c(I(x,y))
ここで、「c」は、カラーマップを表す。
(1)基準値に対するデータN(x,y)の正規化
を含み、後者は、ソナー又はレーダー領域において、すなわち、
1.a.時間可変利得(TVG)、換言すれば対象のジオメトリに依存する所定の法則f(ρ(x,y)):
N(x,y)=(I(x,y)−A)×f(ρ(x,y))
1.b.低速で可変な平均成分:
N(x,y)=(I(x,y)−A)*h(x,y)
を除去するため、畳み込み又は再帰的フィルタリングにより、高域通過フィルタhを信号に適用することになる局所自動利得制御(AGC)
により行うことができる。
(2)データE(x,y)をクランピングして区間[0,T](m<T)内に強制的に留めること:
E(x,y)=min(N(x,y),T)、Tの値が対象のパラメータである。
(3)データをカラーマップにトランスコーディングすること:
O(x,y)=c(E(x,y))(データを最適に知覚できるように、
が適宜選択される)。関数cは、データのレンダリングに用いる「カラーマップ」である。
L(c)=0.2126r+0.7152g+0.0722b
ここで、c=[r;g;b]Tであり、r、g、bは、(0〜1の範囲で符号化された値の)赤、緑、青色の三つ組である。以下では、L(x,y)を用いて、座標(x,y)のピクセルに対して符号化された色に関連付けられた輝度を表記する。また、S(c)を用いて、0〜1の範囲における色cの彩度を表記し、H(c)を用いて、cの色相を表記する。
L(x,y)=E(x,y)γ
ここで、パラメータγは、コントラストを調整する。γの値が小さいと、コントラストが低く、従って「マット」画像が増えるが、信号I(x)が低い撮像領域の可視性が増す。
L(x,y)=K1.logL(x,y)+K2
ここで、K1及びK2は、法則のパラメータである。
ここで、γは、E(x,y)=1の周りの関数Lの傾きを調整するパラメータである(この傾きは、値γ/4を有し、知覚的にコントラストに対応するため、一般に操作者が設定するパラメータとして残される。γの値が小さいほど、弱い信号ゾーンI(x,y)を読み取り易い)。本法則の形式は、視力を含む多くの刺激に対する人間の知覚が刺激の強度に対して対数的に増大すると述べるウェーバーフェヒナーの法則により物理的に正当化され、例えば星の等級を対数スケールで示す理由を説明する。本法則は、従って、1に近い値では対数法則に近づき、値が大きいほど対数法則から外れて行く。
− 極めて短い、従ってIk(x,y)が適用される区間内で時間的に連続して取得が実行され、センサの設定が輝度区間[Ak−1,Ak](A0=A且つAn=B)に設定されている、センサの異なる設定(露光、開口、利得)、又は
− 単一の画像I(x,y)を取得し、次いで区間[A,B]を複数の輝度区間[A0=A,A1],[A2,A3],...,[Anー2,Anー1],[Anー1,An=B]に再分割し、次いで、
を構築すること
のいずれかにより、同一シーンの複数の画像I1(x,y),...In(x,y)を取得することを含む。
ここで、
であり、wk(x,y)は、データIkを(x,y)で優先させるために適宜選択された重みであり、kは、Ik(x,y)が区間[Ak−1,Akの中央になるべく近いように選択され、関数fは、画像の「トーンマッピング」を実行して、その画像を区間[0,1]内の成分にマッピングする。具体的には、この定式化により、座標(x,y)のピクセルにおいて、画像Ik(x,y)が過小又は過剰に露出されない強度を与える値kを優先させることが可能になる。しかし、基本的な短所として、物体の輝度が必ずしも生データの強度に比例しなくなる。
− Jaroslaw Duda,“Zastosowanie techniki HDR do obrazow pozyskiwanych z sektorowego sonaru skanujacego”,Biuletyn WAT vol. LX,No.3,2011;
− Fredo Durand and Julie Dorsey,“Fast Bilateral Filtering for the Display of High−Dynamic−Range Images”.ACM Transactions on Graphics, 2002,21,3,257−266;
− P.E.Debevec,J.Malik,Recovering High Dynamic Range Radiance Maps from Photographs,International Conference on Computer Graphics and Interactive Techniques,Addison−Wesley Publishing Co.,New York,USA,1997;
− Martin Lambers,Andreas Kolb,“Adaptive Dynamic Range Reduction for SAR /mages”,7th European conference on synthetic aperture radar(EUSAR),2008;
− Satoshi Hisanaga,Koji Wakimoto and Koji Okamura,“Compression Method for High Dynamic Range lntensity to lmprove SAR Image Visibility”,Proc.International MultiConference of Engineers and Computer Scientists 2011 vol.1,lMECS 2011,March 16−18,2011,Hong Kong
− 西安電子科技大学から出願された中国特許第1004408752B号明細書は、混合トーンマッピングアルゴリズムに基づいて高ダイナミックレンジを用いる画像圧縮方法を記述しており、及び
− 国際公開第2007113824A2号パンフレットは、高ダイナミックレンジを用いてレーダーデータからの合成画像のボリューム可視化を容易にする方法を記述されている。
− ナイーブ法の場合:問題に適しており、且つ本質的に輝度により符号化された振幅データ(信号の振幅が大きいほど画面上の表現の輝度が高くなる)を視認するのに慣れた操作者にとって読み易いカラーマップを見つけることは、困難である。
− 従来のLDR方法(正規化、クランピング、トランスコーディング)の場合:後者は、有用な低ダイナミックレンジを有する画像のみを表示できる。このダイナミックレンジは、トーンマッピング技術を用いて高めることができるが、満足なレベルまでダイナミックレンジを高めることができない。
− HDRレンダリングの場合:満足できるレベルまでダイナミックレンジを高めることができるが、本方式は、(幅が1又は2ピクセルの)極小空間大きさのピーク強度を示すデータには限界がある。そのデータは、適切な方法により異なる色でレンダリングするにもかかわらず、これらのピークが極めて小さいため、多くの場合、全解像度を見ることさえ困難であり、且つ画像がスケールダウン(「ズームアウト」)された場合に更に困難である。
− 低及び高強度のデータの同時視覚化、
− 高強度の信号(多くの場合、目標に物理的に関連付けられた)に対応する点の空間配置の検出及び解釈(操作者によるこの検出の確率は、非線形低域通過フィルタリングステップによる空間大きさの増大により顕著に増大し、次いで以前には小さ過ぎるために目視できなかった高強度のゾーンを見えるようにする)、
− 色及び輝度並びに非線形フィルタリングステップにより得られたエコーの周りのハローの大きさによる点の強度のレベルの読み取り
をより容易に実行できるように適宜構成されている。
− 行列データから、色が入力信号の低強度を色で表す、「低強度画像」と称される第1のカラー画像を生成するステップと、
− 行列データから、色が入力信号の高強度を表す、「高強度画像」と称される第2のカラー画像を生成するステップであって、画像の各点(x,y)において、座標(0,0)の点に関して中心対称である単峰性インパルス応答の非線形低域通過フィルタを適用することを含む少なくとも1つのフィルタリングステップを含むステップと、
− 低強度画像と高強度画像とを、定数係数によって重み付けされている、色成分の各々の平均となるように一点ずつ合成して、併合カラー画像を生成するステップと、
− 各点において、得られた画像の輝度が低強度画像の輝度に等しいように、合成された画像を正規化するステップと、
− 得られた画像を行列の少なくとも1つの軸上で再寸法設定するステップと
を含むステップを含む。
− 低強度画像の生成は、第1の一時行列であって、原画像Iからその最小値Aが減算されたものに等しい第1の一時行列を構築し、その後、所定の閾値T1よりも高い値を有する入力を除いて、前記入力がT1に設定されている状態で、第1のものに等しい第2の一時行列を構築し、その後、第1のカラーパレットを適用することにより、第2の一時行列から低強度カラー画像を構築することを含むステップを含む。
− 高強度画像の生成は、
− 所定の閾値T1よりも低い値を有する入力を除いて、前記入力がゼロに設定されており、且つ閾値T2よりも高い値を有する入力を除いて、前記入力がT2に設定されている状態で、全ての点で第1の一時行列に等しい第3の一時行列を計算することと、
− 第3の一時行列の最大値を計算することと、
− 座標(0,0)の点に関して中心対称である単峰性インパルス応答の非線形低域通過フィルタリングを第3の一時行列に適用することと、
− 非線形低域通過フィルタリングの結果を一点ずつ結果の全体的な最大値で除算し、且つ第3の一時行列の最大値で乗算することと、
− 除算の結果に第2のカラーパレットを適用することにより、高強度カラー画像を構築することと
を含むステップを含む。
− 得られた画像の再寸法設定は、事前に設定されているか又は操作者の制御に委ねられているかのいずれかである要因による、2つの軸の少なくとも1つ上におけるアップサンプリング又はダウンサンプリングを含む。
− 画像の行列の2つの水平及び垂直軸上における非線形フィルタのインパルス応答の典型的な幅は、2つの軸に沿った入力データの空間分解能に依存する。
− 画像の行列の2つの水平及び垂直軸上における非線形フィルタのインパルス応答の典型的な幅は、2つの軸上における画像のアップサンプリングレート又はダウンサンプリングレートに依存する。
− 低強度画像を生成するために使用されるカラーパレットは、本質的に一定の色相及び彩度のものであり、且つデータの符号化は、信号の強度と共に増大する輝度の変化によって実行される。
− 入力信号は、局所水平軸上のその平均が目標値mに制約されるように事前調整され、低強度画像のために使用されるカラーパレットの輝度は、対数シグモイド則に従う。より具体的には、信号の強度をxと表記すれば、輝度は、法則(x/m)γ/((x/m)γ+1)に従い、ここで、γは、ハードコードされているか又は設定パラメータとして操作者によって提供される定数パラメータである。
− 高強度画像を生成するために使用されるカラーパレットは、一方ではI(x,y)≦A+T1について本質的にゼロ輝度のものであり、且つI(x,y)と共に増大し、及び他方では可変の色相及び彩度のものであり、色相及び彩度の変化は、信号強度情報を符号化する。
− 非線形フィルタのインパルス応答は、2Dガウス関数であって、その標準偏差は、「典型的な幅」パラメータである、2Dガウス関数である。
この関数は、2D画像である。値I(x,y)は、例えば、音響又は高周波波の振幅に対応する、センサにより記録されたスカラー値(例えば、ボルト単位の電圧に比例)である。Iに用いる値の単位を[U]と表記する。特定の実装形態において、2D画像は、行列の各入力がスカラーを含む2D行列である。
この関数は、表示対象の画像に対応するカラー2D画像である。特定の実装形態において、2D画像は、各入力において、表示対象の色成分を含む3次元ベクトルを含む2Dテーブルである。
−
から取得されて単位[U]で表される2つの閾値T1>0及びT2>T1、
− ピクセルで表わされ、画像Iの解像度及び前記画像の表示解像度の関数として可変な2つの半径r1及びr2、
− 所定の又はパラメータ化可能な定数重み付け係数α∈[0,1]、
− 第1のカラーマップ
これにより、Clo(0)=[0;0;0]T(黒)且つClo(v)=[1;1;1]T(v=T1の場合)となる。十分に大きい局所水平軸にわたり画像の平均が実質的に一定であるように事前調整されたソナー又はレーダーデータの場合、色Clo(v)は、必ずしもではないが、好適には、図6aに示す以下の3つの特徴を有する。
− 輝度L(Clo(v))は、vと共に増大する(但し必ずしも線形にではない)。
− 輝度L(Clo(v))は、更に好適には、対数シグモイド的に、すなわちlogvに対してシグモイド的に増大する。
− 色の色相及び彩度は、殆ど変化せず、理想的には一定である。
− 第2のカラーマップ
これにより、n個の値vi,i∈0..n(v0=T1,vi<vi+1且つvn=T2)及び値viに関連付けられたn個の色Ciが与えられ、それにより、
− vi≦v<vi+1の場合、
−
(色成分の加重平均ベクトル)となり;
− C0=[0;0;0]T(黒)である。
− 画像をディスプレイに合わせて再寸法設定するために、x及びy軸にそれぞれ適用する2つのスケール係数
及び
。fx(又はfy)が1よりも大きい場合、画像を拡大(「ズームイン」効果)し、さもなければ縮小(「ズームアウト」)する。これらのスケール係数は、事前に設定されているか、又は例えばキーボードキー、ボタン、サムホイール若しくはポテンショメーター等のグラフィック若しくは物理インターフェースを介して操作者/ユーザーにより制御される。
による画像データの取得及びフォーマッティングの第1のステップ(300)後、本方法は、Iのデータをオフセッティング及びクランピングすることにより強制的に[0,T1]内に収めることができる(302)。
Elo(x,y)=min(I(x,y)−A,T1)
mhi=max(x,y)Ehi
により、Ehiのフィルタリングされたバージョンを表記する。
ここで、フィルタの局所2Dインパルス応答hは、有限又は無限であるかに依らず、以下の特徴を有する。
− インパルス応答は、分割可能であり、すなわち、2Dフィルタリングは、x軸上の非線形フィルタhx、続いて他軸上で別のフィルタリングhyを用いて1Dフィルタリングと同様に実行することができる。
− インパルス応答hx(又はhy)は、u(又はv)が値0であれば値1を有し、u(又はv)の関数として0の周りに対称的に減少し、単峰性且つ中心対称である。
− これらの各々のフィルタの典型的な空間有効範囲、すなわちフィルタのインパルス応答がゼロよりも顕著に大きい領域の幅は、半径rx(x軸上における)及びry(y軸上における)に比例する。
− これらの半径rx(x軸上における)及びry(y軸上における)は、フィルタを非線形にレンダリングする座標(x,y)の関数として可変であり、以下の関係に従う。
rx(Ehi(x,y),r1,r2,fx,ρx,bf,bρ)=kx(ρx,fx,fy,bf,bρ)r(Ehi(x,y),r1,r2)
ry(Ehi(x,y),r1,r2,fy,ρy,bf,bρ)=ky(ρy,fy,bf,bρ)r(Ehi(x,y),r1,r2)
ここで、r(0,r1,r2,fx,ρx,bf,bρ)=r1,r(T2−T1,r1,r2,fx,ρx,,bf,bρ)=r2であり、関数rは、(場合によりレベルによって)Ehi(x,y)の値と共に増大する。
− kx=ky=kにおいてピクセル単位で表される定数kは、例えば、パラメータ化ファイルに事前に設定されているか又はパラメータ設定として操作者に委ねられている。この場合、フィルタは、入力画像の解像度から独立しており、Bのインパルス応答は、両軸上で同じ大きさである。この選択肢が用いられるのは、操作者がボタン等のグラフィック又は物理インターフェースを介してブール変数を「偽」にトグル切替を行った場合である。
− 又はkx=ρx×k且つky=ρy×kにおいて、メートル単位で表される定数kは、例えば、パラメータ化ファイルに事前に設定されているか又はパラメータ設定として操作者に委ねられている。この場合、フィルタは、入力画像の解像度に依存しており、Bのインパルス応答は、もはや両軸上の大きさが同じではない。この選択肢が用いられるのは、操作者がボタン等のグラフィック又は物理インターフェースを介してブール変数を「真」にトグル切替を行った場合である。
− kx、kyをスケール係数fx、fyに依存させることも可能である。例えば、スケール係数が増大する(「ズームイン」)場合、Bにより誘発されたハローの半径を減らすために、kx、kyを減らすことが有利であり、逆に「ズームアウト」の場合又はハローの半径を増大させたい場合にはこの逆である。
の最大値は、Ehiの最大値と異なる。本方法は、画像の最大値が再びmhiになるように(線形に)再正規化可能(312)にし、
次いでデータを第2のカラーマップにトランスコーディング可能にする(314)。
但し、上述のように、
である。値Ihi(x,y)は、従って、3成分(赤、緑、青)のベクトルとして記述された色である。
− 従来のHDRレンダリングと異なり、画像の「高強度」部分(データEhiに含まれる)に非線形低域通過フィルタリングを適用するという事実により、高強度のピクセルの空間有効範囲を広げることができる。従って、例えばレーダー又はソナー分野において、定時的な強いエコーが存在する場合、その空間有効範囲が広がるため強調され、従ってステップ(320)において画像が再サンプリング(「ズームアウト」)された場合でも知覚され易くなる。
− 典型的に、シーンの大部分を構成する弱いエコーに対する従来のレンダリングを近似するために、カラーパレットCloを一定の色相及び彩度の色の勾配として選択すること。対照的に、カラーパレットChiの階調においてレンダリングされた少数の強いエコーの場合、濃く且つ変化した色相を選択するという事実により、強いエコーの強度が直接符号化することができるため、操作者は、強度及び色の関係に応じてエコーの力を認識する動作を実行することができる。エコーの力の印象も、ハローの空間有効範囲がこれらのエコーを囲んでいるという事実により補強される。エコーの強調は、シーンの解析に際して操作者の認識負荷を軽減することに役立ち得る。
− ある選択肢を起動する簡単な動作により、得られたハローの大きさ(定数kx、ky)を視覚化スケールパラメータfx、fyに紐付けることができるという事実により、有利には、操作者は、視覚化スケールfx、fyを用いて、2種類の作業を、典型的に周辺機器(キーボードのキー、マウスのサムホイール等)を介して入力された単一コマンドとして実行できる。これらの2つの作業は、以下の通りである。
− 「ズームアウト」が顕著であり、操作者がシーンを全体として又は少なくとも広い水平軸にわたり知覚した場合、潜在的に目標である強いエコーの視覚的検出。この機能が可能なのは、「ハロー」効果が顕著であり、「遠いエコーを見る」ことができる場合であるのに対し、フィルタリングステップがなければ、これらのエコーは、スケーリング中に消失するであろう。
− 係数により符号化された「ズーム」レベルが顕著であり、目標に沿った顕著な信号ピクセルの分布(「明点」として当業者に知られる)に操作者がより注目しており、且つ高水準の空間精度が重要であるエコーの構造の視覚的解析。この場合、非線形低域通過フィルタリングステップで得られる「曖昧な」効果を減少させなければならないが、非線形フィルタリングステップなしに従来のHDR方法を用いて得られるものと比較して、シーンが分かり易くなる程度に十分強くなければならず、(全ての強いピクセルを単一の値にまとめるクランピングステップにより、目標が典型的に均一な色のピクセルのブロブ(「ピクセルブロブ」)である)LDR法ではなおさらそうである。
− 色相及び彩度が本質的に一定であるカラーパレットClo並びに輝度が増大する濃いカラーパレットChiを用いる基本要因は、加重平均ステップ(316)及びこれに続くステップ(318)で実行される再正規化を介しても注目すべき効果が生じ得る。実際、加重平均により、再正規化が継続するため、「低強度」画像と同じ輝度でレンダリングされ続ける近傍の弱い信号ピクセルの存在を一切マスキングすることなく、強いエコーの存在を認識することができる。この効果が重要である理由は、レーダー及びソナーでは、弱い信号が、シーンに関する情報を一般に影を通して提供するからである。このようなパレットの適切な選択がなければ、影は、消える場合がある。
ここで、δは、ディラック分布であり、aは、拡散されなかった光の一部であり、Fは、拡散を決定する角度θの関数である(人間観察者の場合、∫k/F(θ)は、典型的に10%の値を有する)。Fにはいくつかのモデルが存在し、特にGreg Spencer;Peter Shirley;Kurt Zimmerman;Donald P.Greenberg(1995),“Physically−based glare effects for digital images”.Siggraph:325.doi:10.1145/218380.218466.ISBN 0897917014に記述されているVosモデルが含まれる。ハローの知覚される半径は、観察者の位置における物体がなす角度θ又は同様に観察者から(センサから)光源までの距離に依存することを強調しておく。ソナー及びレーダーでは、この現象は生じない。従って、現在文献から知られるように、ブルーミングは、ソナー又はレーダー画像の表現に加わる効果ではない。
Claims (13)
- 高ダイナミックレンジを有する入力信号からの行列データIのための画像レンダリング方法であって、コンピュータによって実施され、
− 前記初期行列データから、色が前記入力信号の低強度を表す、「低強度画像」と称される第1のカラー画像を生成するステップと、
− 前記初期行列データから、色が前記入力信号の高強度を表す、「高強度画像」と称される第2のカラー画像を生成するステップであって、前記画像の各点(x,y)において、座標(0,0)の点に関して中心対称である単峰性インパルス応答の非線形低域通過フィルタを適用することを含む少なくとも1つのフィルタリングステップを含むステップと、
− 前記低強度画像と前記高強度画像とを、定数係数によって重み付けされている、色成分の各々の平均となるように一点ずつ合成して、併合カラー画像を生成するステップと、
− 各点において、前記得られた画像の輝度が前記低強度画像の輝度に等しいように、前記合成された画像を正規化するステップと、
− 前記得られた画像を行列の少なくとも1つの軸上で再寸法設定するステップと
を含むステップを含む方法。 - 低強度画像を生成する前記ステップは、第1の一時行列であって、前記原画像Iからその最小値Aが減算されたものに等しい第1の一時行列を構築し、その後、所定の閾値T1よりも高い値を有する入力を除いて、前記入力がT1に設定されている状態で、前記第1のものに等しい第2の一時行列を構築し、その後、第1のカラーパレットを適用することにより、前記第2の一時行列から低強度カラー画像を構築することを含むステップを含む、請求項1に記載の方法。
- 高強度画像を生成する前記ステップは、
− 前記所定の閾値T1よりも低い値を有する入力を除いて、前記入力がゼロに設定されており、且つ閾値T2よりも高い値を有する入力を除いて、前記入力がT2に設定されている状態で、全ての点で前記第1の一時行列に等しい第3の一時行列を計算することと、
− 前記第3の一時行列の最大値を計算することと、
− 前記座標(0,0)の点に関して中心対称である単峰性インパルス応答の非線形低域通過フィルタリングを前記第3の一時行列に適用することと、
− 前記非線形低域通過フィルタリングの結果を一点ずつ前記結果の全体的な最大値で除算し、且つ前記第3の一時行列の前記最大値で乗算することと、
− 前記除算の結果に第2のカラーパレットを適用することにより、高強度カラー画像を構築することと
を含むステップを含む、請求項2に記載の方法。 - 前記得られた画像の前記再寸法設定は、事前に設定されているか又は操作者の制御に委ねられているかのいずれかである要因による、2つの軸の少なくとも1つ上におけるアップサンプリング又はダウンサンプリングを含む、請求項1〜3のいずれか一項に記載の方法。
- 前記画像の前記行列の前記2つの水平及び垂直軸上における前記非線形フィルタの前記インパルス応答の典型的な幅は、前記2つの軸に沿った前記入力データの空間分解能に依存する、請求項1〜4のいずれか一項に記載の方法。
- 前記画像の前記行列の前記2つの水平及び垂直軸上における前記非線形フィルタの前記インパルス応答の前記典型的な幅は、前記2つの軸上における前記画像のアップサンプリング又はダウンサンプリングレートに依存する、請求項1〜5のいずれか一項に記載の方法。
- 前記低強度画像を生成するために使用される前記カラーパレットは、本質的に一定の色相及び彩度のものであり、且つ前記データの符号化は、前記信号の前記強度と共に増大する輝度の変化によって実行される、請求項1〜6のいずれか一項に記載の方法。
- 前記入力信号は、局所水平軸上のその平均が目標値mに制約されるように事前調整され、前記低強度画像のために使用される前記カラーパレットの輝度は、対数シグモイド則に従う、請求項1〜7のいずれか一項に記載の方法。
- 前記高強度画像を生成するために使用される前記カラーパレットは、一方ではI(x,y)≦A+T1について本質的にゼロ輝度のものであり、且つI(x,y)と共に増大し、及び他方では可変の色相及び彩度のものであり、色相及び彩度の変化は、信号強度情報を符号化する、請求項1〜8のいずれか一項に記載の方法。
- 前記非線形フィルタの前記インパルス応答は、2Dガウス関数であって、その標準偏差は、「典型的な幅」パラメータである、2Dガウス関数である、請求項1〜9のいずれか一項に記載の方法。
- 前記再寸法設定された画像を表示又は印刷するステップを更に含む、請求項1〜10のいずれか一項に記載の方法。
- 請求項1〜11のいずれか一項に記載の方法の前記ステップを実施するための手段を含む、高ダイナミックレンジを有する入力信号からの行列データのための画像レンダリング装置。
- コンピュータプログラム製品であって、コンピュータプログラムは、前記プログラムがコンピュータ上で実行されるとき、請求項1〜11のいずれか一項に記載の方法の前記ステップの実行を可能にする符号命令を含む、コンピュータプログラム製品。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1800714 | 2018-07-05 | ||
FR1800714A FR3083633B1 (fr) | 2018-07-05 | 2018-07-05 | Procede et dispositif d'affichage de donnees sonar ou radar a haute dynamique |
PCT/EP2019/067995 WO2020007981A1 (fr) | 2018-07-05 | 2019-07-04 | Procede et dispositif d'affichage de donnees sonar ou radar a haute dynamique |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021530059A true JP2021530059A (ja) | 2021-11-04 |
JP7220770B2 JP7220770B2 (ja) | 2023-02-10 |
Family
ID=65031144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021500236A Active JP7220770B2 (ja) | 2018-07-05 | 2019-07-04 | 高ダイナミックソナー又はレーダーデータを表示するための方法及び装置 |
Country Status (8)
Country | Link |
---|---|
US (1) | US11328393B2 (ja) |
EP (1) | EP3818495B1 (ja) |
JP (1) | JP7220770B2 (ja) |
AU (1) | AU2019298637B2 (ja) |
CA (1) | CA3104996A1 (ja) |
FR (1) | FR3083633B1 (ja) |
SG (1) | SG11202100061XA (ja) |
WO (1) | WO2020007981A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3872749A1 (en) * | 2020-02-25 | 2021-09-01 | STMicroelectronics (Research & Development) Limited | Local tone mapping for hdr video |
CN115685170B (zh) * | 2023-01-04 | 2023-05-09 | 青岛国数信息科技有限公司 | 基于强化学习的主动声纳目标回波检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014510339A (ja) * | 2011-03-02 | 2014-04-24 | ドルビー ラボラトリーズ ライセンシング コーポレイション | 局所マルチスケールトーンマッピングオペレータ |
JP2016517637A (ja) * | 2013-02-21 | 2016-06-16 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 改善されたhdrイメージ符号化及び復号化方法、装置 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59176767A (ja) | 1983-03-25 | 1984-10-06 | 株式会社デンソー | 車両用表示装置 |
FR2821163B1 (fr) | 2001-02-16 | 2003-05-09 | Thomson Marconi Sonar Sas | Procede de detection d'objets mobiles au moyen de bouees sous-marines passives |
US6617996B2 (en) * | 2001-06-14 | 2003-09-09 | Mala Geoscience Ab | Ground penetrating radar with audible output |
IL174737A (en) | 2006-04-03 | 2014-12-31 | Camero Tech Ltd | Device and method for visualizing volume in broadband radar |
GB2453927A (en) * | 2007-10-12 | 2009-04-29 | Curtiss Wright Controls Embedded Computing | Method for improving the representation of targets within radar images |
WO2011012646A2 (en) * | 2009-07-28 | 2011-02-03 | F. Hoffmann-La Roche Ag | Non-invasive in vivo optical imaging method |
US10787234B2 (en) * | 2014-05-23 | 2020-09-29 | Hegel Industrial Solutions Pty Ltd | Data capture device and system |
US9460499B2 (en) * | 2014-05-30 | 2016-10-04 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Systems and methods for selective enhancement of a region of interest in an image |
CN104408752B (zh) | 2014-12-11 | 2017-07-04 | 西安电子科技大学 | 一种基于混合色调映射算法的高动态范围图像压缩方法 |
EP3043289B1 (en) * | 2015-01-07 | 2023-04-19 | Honda Research Institute Europe GmbH | Control system for an autonomous vehicle and a method for generating a control signal and autonomous vehicle equipped with such control system |
JP6910068B2 (ja) * | 2015-06-11 | 2021-07-28 | ユニバーシティ オブ ピッツバーグ−オブ ザ コモンウェルス システム オブ ハイヤー エデュケーションUniversity Of Pittsburgh Of The Commonwealth System Of Higher Education | ヘマトキシリン・エオシン(h&e)染色組織画像における関心領域を調べて、多重化/高多重化蛍光組織画像で腫瘍内細胞空間的不均一性を定量化するシステム及び方法 |
EP3340167A1 (en) * | 2016-12-22 | 2018-06-27 | Thomson Licensing | Displaying method for displaying an image on a first display device and at least a portion of the same image on a second display device |
US11497389B2 (en) * | 2018-10-03 | 2022-11-15 | Karl Storz Imaging, Inc. | Attachment system for conditioning light between endoscope and camera |
-
2018
- 2018-07-05 FR FR1800714A patent/FR3083633B1/fr not_active Expired - Fee Related
-
2019
- 2019-07-04 US US17/257,845 patent/US11328393B2/en active Active
- 2019-07-04 JP JP2021500236A patent/JP7220770B2/ja active Active
- 2019-07-04 WO PCT/EP2019/067995 patent/WO2020007981A1/fr active Application Filing
- 2019-07-04 EP EP19734430.2A patent/EP3818495B1/fr active Active
- 2019-07-04 SG SG11202100061XA patent/SG11202100061XA/en unknown
- 2019-07-04 AU AU2019298637A patent/AU2019298637B2/en active Active
- 2019-07-04 CA CA3104996A patent/CA3104996A1/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014510339A (ja) * | 2011-03-02 | 2014-04-24 | ドルビー ラボラトリーズ ライセンシング コーポレイション | 局所マルチスケールトーンマッピングオペレータ |
JP2016517637A (ja) * | 2013-02-21 | 2016-06-16 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 改善されたhdrイメージ符号化及び復号化方法、装置 |
Also Published As
Publication number | Publication date |
---|---|
FR3083633A1 (fr) | 2020-01-10 |
SG11202100061XA (en) | 2021-02-25 |
JP7220770B2 (ja) | 2023-02-10 |
EP3818495A1 (fr) | 2021-05-12 |
US11328393B2 (en) | 2022-05-10 |
US20210295479A1 (en) | 2021-09-23 |
AU2019298637A1 (en) | 2021-02-25 |
FR3083633B1 (fr) | 2020-05-29 |
AU2019298637B2 (en) | 2023-12-07 |
CA3104996A1 (en) | 2020-01-09 |
WO2020007981A1 (fr) | 2020-01-09 |
EP3818495B1 (fr) | 2022-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9934562B2 (en) | Method for dynamic range editing | |
US10410327B2 (en) | Shallow depth of field rendering | |
US8477135B2 (en) | Method and apparatus for volume rendering using depth weighted colorization | |
US20060034536A1 (en) | Systems and methods relating to magnitude enhancement analysis suitable for high bit level displays on low bit level systems, determining the material thickness, and 3D visualization of color space dimensions | |
CA2306593A1 (en) | Smart erasure brush | |
WO2006079115A2 (en) | Spatial standard observer | |
CN110248242B (zh) | 一种图像处理和直播方法、装置、设备和存储介质 | |
JP6550138B2 (ja) | 映像処理装置 | |
JP7220770B2 (ja) | 高ダイナミックソナー又はレーダーデータを表示するための方法及び装置 | |
JPH06510133A (ja) | エイリアス除去画像の表示装置 | |
EP1481632B1 (en) | Ophtalmologic image processing apparatus | |
Zhang et al. | Adaptive local contrast enhancement for the visualization of high dynamic range images | |
Zhou et al. | Photographic high-dynamic-range scalar visualization | |
JP4359662B2 (ja) | カラー画像の露出補正方法 | |
US11521306B2 (en) | Image processing apparatus, image processing method, and storage medium | |
US20220012521A1 (en) | System for luminance qualified chromaticity | |
JP2009219140A (ja) | カラー画像の露出評価方法 | |
WO2006002327A2 (en) | Systems and methods relating to magnitude enhancement analysis suitable for high bit level displays on low bit level systems, determining the material thickness, and 3d visualization of color space dimensions | |
Biryukov et al. | Scientific Visualization for Interrogation of Results in Lighting Simulation System | |
HU196007B (hu) | Szabadalmi leírás szolgálati találmány b |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210305 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220531 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220804 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20221206 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20221221 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20230110 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230131 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7220770 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |