JP2892676B2 - 画像データ圧縮方法 - Google Patents
画像データ圧縮方法Info
- Publication number
- JP2892676B2 JP2892676B2 JP1090727A JP9072789A JP2892676B2 JP 2892676 B2 JP2892676 B2 JP 2892676B2 JP 1090727 A JP1090727 A JP 1090727A JP 9072789 A JP9072789 A JP 9072789A JP 2892676 B2 JP2892676 B2 JP 2892676B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- image data
- data
- compression
- gradations
- 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 - Lifetime
Links
- 238000000034 method Methods 0.000 title claims description 39
- 238000013144 data compression Methods 0.000 title claims description 12
- 238000007906 compression Methods 0.000 description 36
- 230000006835 compression Effects 0.000 description 36
- 238000012545 processing Methods 0.000 description 19
- 230000015654 memory Effects 0.000 description 15
- 230000008569 process Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 9
- 238000010606 normalization Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000007726 management method Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 238000013139 quantization Methods 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000013480 data collection Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- XQTLDIFVVHJORV-UHFFFAOYSA-N tecnazene Chemical compound [O-][N+](=O)C1=C(Cl)C(Cl)=CC(Cl)=C1Cl XQTLDIFVVHJORV-UHFFFAOYSA-N 0.000 description 1
- 230000003936 working memory Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【発明の詳細な説明】 〔産業上の利用分野〕 本発明は、医用画像管理システムの画像ワークステー
シヨンのような、複数の種類の多値画像データを取り扱
う装置における画像データのデータ圧縮を統一的に行な
うのに好適な画像データ圧縮方法に関する。
シヨンのような、複数の種類の多値画像データを取り扱
う装置における画像データのデータ圧縮を統一的に行な
うのに好適な画像データ圧縮方法に関する。
近年、多種類のデイジタル医用画像データを集中的に
保管・管理して、高速の画像検索と高精細の画像表示機
能により、医師による画像診断を支援する等の目的で、
医用画像管理システムが開発された。第9図は医用画像
管理システムの一例を示すブロツク図である。第9図の
一点鎖線の内側は、画像ワークステーシヨン1である。
この画像ワークステーシヨン1は、インタフエース2,
インタフエース3,インタフエース4により、3種類
のデイジタルラジオグラフイ装置(以下DR装置と略記す
る。)と接続している。DR装置は、X線による人体の透
過画像をデイジタル画像データとして収集する装置の総
称である。第9図の5はフイルム方式DR装置であり、通
常のX線フイルム画像から光学的・電気的手段を用いて
画像データを読み出す装置である。6はイメージインテ
ンシフアイア方式DR装置であり、X線イメージインテイ
シフアイアとテレビカメラを用いたDR装置である。7は
輝尽性蛍光体方式DR装置であり、X線により形成された
輝尽性蛍光体上の潜像を光学的・電気的手段で読み出
し、デイジタル化する装置である。こうしたデイジタル
画像収集装置で収集されるデイジタル画像データは、一
般に装置ごとに異なつた階調数を持つている。デイジタ
ルデータの階調数は、2進表現が一般に用いられている
ために、2nで表わされることが多い。ここでnは自然数
である。第9図のシステムでは、5の装置は212=4098
の、6の装置は210=1024の、7の装置は28=256の階調
数を持つている。画像ワークステーシヨン1の主たる機
能は、(1)上記DR装置より収集されたデイジタル画像
データをシステムバス8を介して外部記憶装置9に書き
込む。(2)外部記憶装置9より操作者の所望する画像
データを読み出し、システムバス8を介して画像メモリ
10に送り、該メモリに結合した高精細CRTモニタを用い
た画像表示装置11により画像として表示する。上記機能
を遂行するために、システムバス8には、マイクロコン
ピユータを搭載したCPU12が接続され、更に、操作者の
マイクロコンピユータに対する命令を受付るための操作
卓13が操作卓インタフエース14を介して接続されてい
る。この画像ワークステーシヨン1により画像データを
外部記憶装置9に蓄積する際、画像データの容量が大き
な問題となる。外部記憶装置9としては高速検索性を考
慮して、高速アクセスが可能で大容量である数十メガバ
イトの磁気デイスクや、数ギガバイトの光デイスクが用
いられるが、DR装置による画像データは画像1枚当り数
メガバイトになるので、光デイスクを用いる場合でも、
光デイスク1枚当り高々1000枚程度の画像しか記録でき
ない。そこで第9図の画像ワークステーシヨンでは、専
用の圧縮/再生モジユール15を用いて画像データのデー
タ圧縮が可能なように構成されている。圧縮/再生モジ
ユール15には高速演算装置16が搭載されており、これを
用いて特願昭62-297556号記載の方法による符号化・復
号化処理が行なわれる。該方法のデータ圧縮を用いる
と、X線画像に対して10:1程度の圧縮比が得られるの
で、1個の外部記憶装置に蓄積できる画像の枚数を、デ
ータ圧縮を用いない場合の10倍程度にすることができ
る。上記方法のデータ圧縮を用いる場合、上記(1),
(2)の手続きの代りに、次の(3)〜(8)の手続き
が行なわれる。即ち、(3)DR装置2〜4より収集され
たデイジタル画像データをシステムバス8を介して、画
像メモリに書き込む。(4)圧縮/再生モジユール15の
圧縮処理を起動すると、圧縮/再生モジユール15は、デ
ータ転送線17を用いて、画像メモリ10より画像データ
を読み出し、高速演算装置16を用いて上記方法による符
号化処理を行い、形成された圧縮データをデータ転送線
18を用いて、バツフアメモリ19に書き込む。(5)バ
ツフアメモリ19上の圧縮データをシステムバス8を介し
て外部記憶装置9に書き込む。(6)外部記憶装置9よ
り、操作者の所望する画像に対応する圧縮データを読み
出し、システムバス8を介してバツフアメモリ19に書き
込む。(7)圧縮/再生モジユール15の再生処理を起動
すると、圧縮/再生モジユール15は、データ転送線18
を用いて、バツフアメモリ19より圧縮データを読み出
し、高速演算装置16を用いて上記方法による復号化処理
を行ない、形成された再生画像データを、データ転送線
17を用いて、画像メモリ10に書き込む。(8)画像メ
モリ10上の再生画像データを、該メモリに結合した画像
表示装置11により画像として表示する。
保管・管理して、高速の画像検索と高精細の画像表示機
能により、医師による画像診断を支援する等の目的で、
医用画像管理システムが開発された。第9図は医用画像
管理システムの一例を示すブロツク図である。第9図の
一点鎖線の内側は、画像ワークステーシヨン1である。
この画像ワークステーシヨン1は、インタフエース2,
インタフエース3,インタフエース4により、3種類
のデイジタルラジオグラフイ装置(以下DR装置と略記す
る。)と接続している。DR装置は、X線による人体の透
過画像をデイジタル画像データとして収集する装置の総
称である。第9図の5はフイルム方式DR装置であり、通
常のX線フイルム画像から光学的・電気的手段を用いて
画像データを読み出す装置である。6はイメージインテ
ンシフアイア方式DR装置であり、X線イメージインテイ
シフアイアとテレビカメラを用いたDR装置である。7は
輝尽性蛍光体方式DR装置であり、X線により形成された
輝尽性蛍光体上の潜像を光学的・電気的手段で読み出
し、デイジタル化する装置である。こうしたデイジタル
画像収集装置で収集されるデイジタル画像データは、一
般に装置ごとに異なつた階調数を持つている。デイジタ
ルデータの階調数は、2進表現が一般に用いられている
ために、2nで表わされることが多い。ここでnは自然数
である。第9図のシステムでは、5の装置は212=4098
の、6の装置は210=1024の、7の装置は28=256の階調
数を持つている。画像ワークステーシヨン1の主たる機
能は、(1)上記DR装置より収集されたデイジタル画像
データをシステムバス8を介して外部記憶装置9に書き
込む。(2)外部記憶装置9より操作者の所望する画像
データを読み出し、システムバス8を介して画像メモリ
10に送り、該メモリに結合した高精細CRTモニタを用い
た画像表示装置11により画像として表示する。上記機能
を遂行するために、システムバス8には、マイクロコン
ピユータを搭載したCPU12が接続され、更に、操作者の
マイクロコンピユータに対する命令を受付るための操作
卓13が操作卓インタフエース14を介して接続されてい
る。この画像ワークステーシヨン1により画像データを
外部記憶装置9に蓄積する際、画像データの容量が大き
な問題となる。外部記憶装置9としては高速検索性を考
慮して、高速アクセスが可能で大容量である数十メガバ
イトの磁気デイスクや、数ギガバイトの光デイスクが用
いられるが、DR装置による画像データは画像1枚当り数
メガバイトになるので、光デイスクを用いる場合でも、
光デイスク1枚当り高々1000枚程度の画像しか記録でき
ない。そこで第9図の画像ワークステーシヨンでは、専
用の圧縮/再生モジユール15を用いて画像データのデー
タ圧縮が可能なように構成されている。圧縮/再生モジ
ユール15には高速演算装置16が搭載されており、これを
用いて特願昭62-297556号記載の方法による符号化・復
号化処理が行なわれる。該方法のデータ圧縮を用いる
と、X線画像に対して10:1程度の圧縮比が得られるの
で、1個の外部記憶装置に蓄積できる画像の枚数を、デ
ータ圧縮を用いない場合の10倍程度にすることができ
る。上記方法のデータ圧縮を用いる場合、上記(1),
(2)の手続きの代りに、次の(3)〜(8)の手続き
が行なわれる。即ち、(3)DR装置2〜4より収集され
たデイジタル画像データをシステムバス8を介して、画
像メモリに書き込む。(4)圧縮/再生モジユール15の
圧縮処理を起動すると、圧縮/再生モジユール15は、デ
ータ転送線17を用いて、画像メモリ10より画像データ
を読み出し、高速演算装置16を用いて上記方法による符
号化処理を行い、形成された圧縮データをデータ転送線
18を用いて、バツフアメモリ19に書き込む。(5)バ
ツフアメモリ19上の圧縮データをシステムバス8を介し
て外部記憶装置9に書き込む。(6)外部記憶装置9よ
り、操作者の所望する画像に対応する圧縮データを読み
出し、システムバス8を介してバツフアメモリ19に書き
込む。(7)圧縮/再生モジユール15の再生処理を起動
すると、圧縮/再生モジユール15は、データ転送線18
を用いて、バツフアメモリ19より圧縮データを読み出
し、高速演算装置16を用いて上記方法による復号化処理
を行ない、形成された再生画像データを、データ転送線
17を用いて、画像メモリ10に書き込む。(8)画像メ
モリ10上の再生画像データを、該メモリに結合した画像
表示装置11により画像として表示する。
上記従来の手続きによつて画像データ圧縮を行なう
際、操作者は画像ごとに特願昭62-297556号記載の符号
化方法に含まれる再生画像の原画像に対する忠実度に影
響する2種類のパラメータを操作卓13より入力する必要
があつた。上記符号化方法において、画像は第10図に示
すように一辺16画素の正方形のサブ画像20に分割して取
り扱われる。そして、各サブ画像20の画像データは第11
図の流れ図に従う符号化処理を受けて圧縮される。第11
図は上記文献の第1図に対応する。第11図において、前
処理21では、上記文献の(1),(2)式による直流成
分の抑圧を行なう。双線形関数近似22では、上記文献の
(4)〜(24)式を用いて、入力されたサブ画像のブロ
ツク分割と、ブロツク内の画像データ値の変化に対する
双線形近似を行なう。この過程で、記号dで表わす許容
歪と呼ぶパラメータが用いられる。これは、上記文献の
(24)式に示されており、双線形近似による誤差の許容
程度に関するパラメータである。即ち、dが大きい程、
近似の程度が進み圧縮比が大きくなる反面、再生画像の
原画像に対する忠実度が低下する。次に、量子化23で
は、まず双線形近似の4種類の係数A〜Dを、上記文献
の(25)〜(28)式を用いて正規化する。ここで、第2
のパラメータである4個の正規化係数SDi(i=1〜
4)が用いられる。SDiは量子化の精度に関係してお
り、SDiを小さくする程精度が上がり、再生画像の忠実
度が向上するが、小さくし過ぎると、上記文献の第10図
に示す拡張モードの発生頻度が多くなり、圧縮比が低下
する。正規化された4種類の係数は、それらの統計的性
質に適合するように設計された量子化器により、数ビツ
トに量子化される。次に、係数パツキング24では、量子
化された係数データをブロツク分割形態を記述するQuad
-Treeデータ等と共に上記文献の第7図,第8図に示す
ように圧縮データとしてメモリに格納する。上記第11図
の処理で用いたdとSDi(i=1〜4)の圧縮パラメー
タ25は、第12図に示すように、圧縮データのデータ量を
示す情報等と共に圧縮データ付属情報26となり、圧縮デ
ータ27と圧縮データ付属情報26とを併せて全圧縮データ
28が外部記憶装置に書き込まれる。再生時の復号化処理
においても、第10図のサブ画像20ごとに処理が行なわれ
る。第13図はサブ画像20の復号化処理の流れ図であり、
上記文献の第11図に対応する。第13図においてRM,RF解
析29では、圧縮データを解析して、そこから、Quad-Tre
eデータや、量子化された係数データを取り出す。取り
出された係数データは、逆量子化30で逆量子化された
後、逆正規化31で上記文献の(29)〜(32)式によつて
逆正規化される。ここで第12図の圧縮パラメータ25に含
まれる正規化係数SDi(i=1〜4)を使用する。最後
に双線形補間32において、このようにして再生された係
数データを用い、Quad-Treeデータの示す分割ブロツク
内の画像データ値を、上記文献の(32)〜(35)式に従
つて計算し、各サブ画像の再生が終了する。上述のよう
に、上記文献の方法で画像データを圧縮/再生する場
合、圧縮時の符号化処理において、dとSDi(i=1〜
4)の2種類のパラメータを与える必要がある。第9図
の画像ワークステーシヨン1において、これらのパラメ
ータの値は、画像ごとに変えなければならない。例え
ば、第9図の3種類のDR装置5〜7で同一の被写体のX
線画像の3種類の画像データを得た場合を考えてみる。
いずれの画像データも、透過X線量の大小を表わすもの
に相当するので、画像データの2次元的変動の様子は相
互に相似したものになることは明らかである。しかし、
それらの画像データの2次元的変動の振幅(以下、これ
をコントラストと呼ぶ。)は、それらの画像データの階
調数の比程度に異なつている。それら3種類の画像デー
タを圧縮し再生した時、同様の画質の再生画像を得よう
とする際、元来大きなコントラストを有する画像データ
に対しては、比較的大きな誤差が許されるし、逆に小さ
いなコントラストしか持たない画像データの誤差の大き
さは小さく抑えられなければならない。従つて、圧縮パ
ラメータdとSDi(i=1〜4)は、212の階調数を有す
るフイルム方式DR装置5より得られた画像データに対し
ては、大きく選ぶべきであり、逆に、28の階調数しかな
い輝尽性蛍光体方式DR装置7より得られた画像データに
対しては、小さく選ぶべきである。更に、同一のDR装置
で同一の被写体の画像データを得たとしても、画像を得
る際のX線条件の違いによつてコントラストが異なつて
くる。一般に、X線管電圧が高い場合、画像のコントラ
ストは小さくなる。従つて、そのような場合には、圧縮
パラメータ、特に再生画像の忠実度に直接に関係するd
を小さくするべきである。上述のように、圧縮パラメー
タの値は画像ごとに変えるべきであるが、従来の手続で
は、この要請に応えるために、操作者による手入力を必
要としていた。即ち、第9図において圧縮/再生モジユ
ール15の圧縮処理を起動する前に、操作者13は、圧縮対
象の画像に適すると考えられる圧縮パラメータの値を、
操作卓13から入力していた。
際、操作者は画像ごとに特願昭62-297556号記載の符号
化方法に含まれる再生画像の原画像に対する忠実度に影
響する2種類のパラメータを操作卓13より入力する必要
があつた。上記符号化方法において、画像は第10図に示
すように一辺16画素の正方形のサブ画像20に分割して取
り扱われる。そして、各サブ画像20の画像データは第11
図の流れ図に従う符号化処理を受けて圧縮される。第11
図は上記文献の第1図に対応する。第11図において、前
処理21では、上記文献の(1),(2)式による直流成
分の抑圧を行なう。双線形関数近似22では、上記文献の
(4)〜(24)式を用いて、入力されたサブ画像のブロ
ツク分割と、ブロツク内の画像データ値の変化に対する
双線形近似を行なう。この過程で、記号dで表わす許容
歪と呼ぶパラメータが用いられる。これは、上記文献の
(24)式に示されており、双線形近似による誤差の許容
程度に関するパラメータである。即ち、dが大きい程、
近似の程度が進み圧縮比が大きくなる反面、再生画像の
原画像に対する忠実度が低下する。次に、量子化23で
は、まず双線形近似の4種類の係数A〜Dを、上記文献
の(25)〜(28)式を用いて正規化する。ここで、第2
のパラメータである4個の正規化係数SDi(i=1〜
4)が用いられる。SDiは量子化の精度に関係してお
り、SDiを小さくする程精度が上がり、再生画像の忠実
度が向上するが、小さくし過ぎると、上記文献の第10図
に示す拡張モードの発生頻度が多くなり、圧縮比が低下
する。正規化された4種類の係数は、それらの統計的性
質に適合するように設計された量子化器により、数ビツ
トに量子化される。次に、係数パツキング24では、量子
化された係数データをブロツク分割形態を記述するQuad
-Treeデータ等と共に上記文献の第7図,第8図に示す
ように圧縮データとしてメモリに格納する。上記第11図
の処理で用いたdとSDi(i=1〜4)の圧縮パラメー
タ25は、第12図に示すように、圧縮データのデータ量を
示す情報等と共に圧縮データ付属情報26となり、圧縮デ
ータ27と圧縮データ付属情報26とを併せて全圧縮データ
28が外部記憶装置に書き込まれる。再生時の復号化処理
においても、第10図のサブ画像20ごとに処理が行なわれ
る。第13図はサブ画像20の復号化処理の流れ図であり、
上記文献の第11図に対応する。第13図においてRM,RF解
析29では、圧縮データを解析して、そこから、Quad-Tre
eデータや、量子化された係数データを取り出す。取り
出された係数データは、逆量子化30で逆量子化された
後、逆正規化31で上記文献の(29)〜(32)式によつて
逆正規化される。ここで第12図の圧縮パラメータ25に含
まれる正規化係数SDi(i=1〜4)を使用する。最後
に双線形補間32において、このようにして再生された係
数データを用い、Quad-Treeデータの示す分割ブロツク
内の画像データ値を、上記文献の(32)〜(35)式に従
つて計算し、各サブ画像の再生が終了する。上述のよう
に、上記文献の方法で画像データを圧縮/再生する場
合、圧縮時の符号化処理において、dとSDi(i=1〜
4)の2種類のパラメータを与える必要がある。第9図
の画像ワークステーシヨン1において、これらのパラメ
ータの値は、画像ごとに変えなければならない。例え
ば、第9図の3種類のDR装置5〜7で同一の被写体のX
線画像の3種類の画像データを得た場合を考えてみる。
いずれの画像データも、透過X線量の大小を表わすもの
に相当するので、画像データの2次元的変動の様子は相
互に相似したものになることは明らかである。しかし、
それらの画像データの2次元的変動の振幅(以下、これ
をコントラストと呼ぶ。)は、それらの画像データの階
調数の比程度に異なつている。それら3種類の画像デー
タを圧縮し再生した時、同様の画質の再生画像を得よう
とする際、元来大きなコントラストを有する画像データ
に対しては、比較的大きな誤差が許されるし、逆に小さ
いなコントラストしか持たない画像データの誤差の大き
さは小さく抑えられなければならない。従つて、圧縮パ
ラメータdとSDi(i=1〜4)は、212の階調数を有す
るフイルム方式DR装置5より得られた画像データに対し
ては、大きく選ぶべきであり、逆に、28の階調数しかな
い輝尽性蛍光体方式DR装置7より得られた画像データに
対しては、小さく選ぶべきである。更に、同一のDR装置
で同一の被写体の画像データを得たとしても、画像を得
る際のX線条件の違いによつてコントラストが異なつて
くる。一般に、X線管電圧が高い場合、画像のコントラ
ストは小さくなる。従つて、そのような場合には、圧縮
パラメータ、特に再生画像の忠実度に直接に関係するd
を小さくするべきである。上述のように、圧縮パラメー
タの値は画像ごとに変えるべきであるが、従来の手続で
は、この要請に応えるために、操作者による手入力を必
要としていた。即ち、第9図において圧縮/再生モジユ
ール15の圧縮処理を起動する前に、操作者13は、圧縮対
象の画像に適すると考えられる圧縮パラメータの値を、
操作卓13から入力していた。
本発明の目的は、従来技術の欠点である、上記の圧縮
パラメータの手入力を除去し、第9図の画像ワークステ
ーシヨン1のような、複数の多値画像データを取り扱う
装置において、圧縮パラメータを統一的かつ自動的に決
定する画像データ圧縮方法を提供する事にある。
パラメータの手入力を除去し、第9図の画像ワークステ
ーシヨン1のような、複数の多値画像データを取り扱う
装置において、圧縮パラメータを統一的かつ自動的に決
定する画像データ圧縮方法を提供する事にある。
上記目的を達成するために、本発明では、画像データ
を収集した装置によって定まる1枚の画像の画像データ
全体の階調数と、該画像データ全体ごとに定められる標
準表示階調数とに着目し、圧縮パラメータをそれらの関
数として自動的に決定する。ここで、標準表示階調数と
は次のようなものである。第9図の画像ワークステーシ
ヨン1において、画像メモリ10上の画像データを画像表
示装置11で表示する際に、画像表示装置11に表示階調数
Wと表示階調中央値Lを与える必要がある。画像表示装
置11は第14図に示すように、階調数2nのデイジタルの画
像データ33を最大振幅Vのアナログのビデオ信号34に変
換するデイジタルアナログ変換器35と高精細CRTモニタ3
6よりなる。WとLは第14図中にグラフで示したよう
に、画像データからビデオ信号への変換時に用いられ
る。即ち、2nの幅のデータの内Lを中心としたWの幅の
階調をVの幅のビデオ信号に変換する。このWの標準値
が標準表示階調数である。標準表示階調数は通常画像デ
ータ階調数と同一に設定されるが、異なる場合もある。
第9図の医用画像管理システムにおいて、3種類のDR装
置5〜7は、それぞれ収集した画像データを検査するた
めの画像表示装置を備えている。そして、収集した画像
データのコントラストが小さい場合、表示階調数を小さ
くして、表示画像に充分大きなコントラストを持たせら
れるようになつている。このようにして選ばれた表示階
調数は標準表示階調数として、画像データと共にDR装置
5〜7から画像ワークステーシヨン1に送られ、画像ワ
ークステーシヨン1の画像表示装置11で画像データを表
示するときに使われる。
を収集した装置によって定まる1枚の画像の画像データ
全体の階調数と、該画像データ全体ごとに定められる標
準表示階調数とに着目し、圧縮パラメータをそれらの関
数として自動的に決定する。ここで、標準表示階調数と
は次のようなものである。第9図の画像ワークステーシ
ヨン1において、画像メモリ10上の画像データを画像表
示装置11で表示する際に、画像表示装置11に表示階調数
Wと表示階調中央値Lを与える必要がある。画像表示装
置11は第14図に示すように、階調数2nのデイジタルの画
像データ33を最大振幅Vのアナログのビデオ信号34に変
換するデイジタルアナログ変換器35と高精細CRTモニタ3
6よりなる。WとLは第14図中にグラフで示したよう
に、画像データからビデオ信号への変換時に用いられ
る。即ち、2nの幅のデータの内Lを中心としたWの幅の
階調をVの幅のビデオ信号に変換する。このWの標準値
が標準表示階調数である。標準表示階調数は通常画像デ
ータ階調数と同一に設定されるが、異なる場合もある。
第9図の医用画像管理システムにおいて、3種類のDR装
置5〜7は、それぞれ収集した画像データを検査するた
めの画像表示装置を備えている。そして、収集した画像
データのコントラストが小さい場合、表示階調数を小さ
くして、表示画像に充分大きなコントラストを持たせら
れるようになつている。このようにして選ばれた表示階
調数は標準表示階調数として、画像データと共にDR装置
5〜7から画像ワークステーシヨン1に送られ、画像ワ
ークステーシヨン1の画像表示装置11で画像データを表
示するときに使われる。
第9図のDR装置5〜7のような画像データ収集装置
は、対象画像の最大コントラストを表現するのに必要か
つ充分な大きさに、その画像データ階調数を設計するの
が通例である。従つて、画像データ階調数の比は、画像
データ収集装置に依存する画像データのコントラストの
比をほぼ反映している。先に述べたように、元来大きな
コントラストを持つ画像データに対しては、再生画像デ
ータの原画像データに対する誤差は比較的大きくてよい
ので、再生画像の忠実度に関する圧縮パラメータをデー
タ階調数の関数として決めるようにすれば、多種の画像
データ収集装置を取り扱う場合も、統一的,自動的に圧
縮パラメータが決定できる。又、標準表示階調数は画像
ごとのコントラストの違いを反映するものなので、圧縮
パラメータを標準表示階調数の関数とすれば、よりきめ
細かい取り扱いとなる。
は、対象画像の最大コントラストを表現するのに必要か
つ充分な大きさに、その画像データ階調数を設計するの
が通例である。従つて、画像データ階調数の比は、画像
データ収集装置に依存する画像データのコントラストの
比をほぼ反映している。先に述べたように、元来大きな
コントラストを持つ画像データに対しては、再生画像デ
ータの原画像データに対する誤差は比較的大きくてよい
ので、再生画像の忠実度に関する圧縮パラメータをデー
タ階調数の関数として決めるようにすれば、多種の画像
データ収集装置を取り扱う場合も、統一的,自動的に圧
縮パラメータが決定できる。又、標準表示階調数は画像
ごとのコントラストの違いを反映するものなので、圧縮
パラメータを標準表示階調数の関数とすれば、よりきめ
細かい取り扱いとなる。
以下、本発明の一実施例を添付図面を参照して説明す
る。
る。
本実施例の画像データ圧縮は、第9図の画像ワークス
テーシヨン1で、CPU12によつて、第1図の流れ図に従
つて行なわれる。第1図の処理に先立ち、画像ワークス
テーシヨン1には、対象画像に関して第2図に示す画像
付属情報37と画像データ38よりなる全画像データ39が存
在している。画像付属情報37には、患者名,撮影日付等
と共に、画像データ38が収集されたDR装置によつて定ま
る画像データ階調数40と、そのDR装置より画像データ38
と共に送られて来た標準表示階調数41が含まれている。
全画像データ39の内、画像付属情報37はCPU12の作業用
メモリに、画像データ38は画像メモリ10に書き込まれて
いる。第1図の処理において、ステツプ42では、画像
付属情報37中の画像データ階調数40と標準表示階調数41
より第11図の圧縮パラメータ25を求め、圧縮/再生モジ
ユール15にそのパラメータを入力する。本実施例におい
ては、第11図の圧縮パラメータ25の内、許容歪dは第3
図に示すように、標準表示階調数WSに比例するものとし
て計算される。ここで比例係数はそれぞれの画像によら
ず定められるパラメータである。もう1個の圧縮パラメ
ータである正規化係数SDi(i=1〜4)は第4図に示
すように、画像データ階調数2nに比例するものとして計
算される。ここでも、比例係数は画像によらず定められ
るパラメータである。次にステツプ43では、圧縮/再
生モジユールの圧縮処理を起動する。画像メモリ10上の
画像データは、第10図に示したサブ画像ごとに、ステツ
プで与えられた圧縮パラメータ25を使つた、第11図の
符号化処理を受け、圧縮データとなつてバツフアメモリ
19に書き込まれる。最後にステツプ44では、バツフア
メモリ19上の圧縮データ27を第12図に示すように圧縮デ
ータ付属情報26と共に全圧縮データ28として外部記憶装
置9に記録する。本実施例によれば、コントラストの小
さな画像データに対しては、小さな標準表示階調数WSが
付属するので、第3図より許容歪dも小さくなり、再生
画像の原画像に対する誤差が小さく抑えられる。又、画
像データ階調数2nが小さなDR装置より得られたコントラ
ストの小さな画像データに対しては、第4図より小さな
正規化係数SDi(i=1〜4)が与えられる。このよう
に、画像ごとに適合した圧縮パラメータを自動的に得る
ことができる。本実施例で、正規化係数SDiを標準表示
階調数WSに比例するように変更しても同様な効果が得ら
れる。しかし、正規化係数SDi(i=1〜4)は、先に
述べたように、圧縮データの精度に関係するものなの
で、画像データの濃淡表現の尺度に関係する画像データ
階調数2nの関数として取り扱つた方がふさわしいと考え
る。
テーシヨン1で、CPU12によつて、第1図の流れ図に従
つて行なわれる。第1図の処理に先立ち、画像ワークス
テーシヨン1には、対象画像に関して第2図に示す画像
付属情報37と画像データ38よりなる全画像データ39が存
在している。画像付属情報37には、患者名,撮影日付等
と共に、画像データ38が収集されたDR装置によつて定ま
る画像データ階調数40と、そのDR装置より画像データ38
と共に送られて来た標準表示階調数41が含まれている。
全画像データ39の内、画像付属情報37はCPU12の作業用
メモリに、画像データ38は画像メモリ10に書き込まれて
いる。第1図の処理において、ステツプ42では、画像
付属情報37中の画像データ階調数40と標準表示階調数41
より第11図の圧縮パラメータ25を求め、圧縮/再生モジ
ユール15にそのパラメータを入力する。本実施例におい
ては、第11図の圧縮パラメータ25の内、許容歪dは第3
図に示すように、標準表示階調数WSに比例するものとし
て計算される。ここで比例係数はそれぞれの画像によら
ず定められるパラメータである。もう1個の圧縮パラメ
ータである正規化係数SDi(i=1〜4)は第4図に示
すように、画像データ階調数2nに比例するものとして計
算される。ここでも、比例係数は画像によらず定められ
るパラメータである。次にステツプ43では、圧縮/再
生モジユールの圧縮処理を起動する。画像メモリ10上の
画像データは、第10図に示したサブ画像ごとに、ステツ
プで与えられた圧縮パラメータ25を使つた、第11図の
符号化処理を受け、圧縮データとなつてバツフアメモリ
19に書き込まれる。最後にステツプ44では、バツフア
メモリ19上の圧縮データ27を第12図に示すように圧縮デ
ータ付属情報26と共に全圧縮データ28として外部記憶装
置9に記録する。本実施例によれば、コントラストの小
さな画像データに対しては、小さな標準表示階調数WSが
付属するので、第3図より許容歪dも小さくなり、再生
画像の原画像に対する誤差が小さく抑えられる。又、画
像データ階調数2nが小さなDR装置より得られたコントラ
ストの小さな画像データに対しては、第4図より小さな
正規化係数SDi(i=1〜4)が与えられる。このよう
に、画像ごとに適合した圧縮パラメータを自動的に得る
ことができる。本実施例で、正規化係数SDiを標準表示
階調数WSに比例するように変更しても同様な効果が得ら
れる。しかし、正規化係数SDi(i=1〜4)は、先に
述べたように、圧縮データの精度に関係するものなの
で、画像データの濃淡表現の尺度に関係する画像データ
階調数2nの関数として取り扱つた方がふさわしいと考え
る。
次に、本発明の第2の実施例について述べる。本実施
例も上記実施例と同様に、第9図のワークステーシヨン
1で、CPU12によつて第1図の流れ図に従つて行なわれ
る。但し、圧縮/再生モジユール15におけるサブ画像ご
との符号化処理の内容が第11図のものから若干異なる第
5図のものに変更されており、これに伴い、圧縮パラメ
ータ25にスケーリングフアクタSFが追加され、これを含
めた圧縮パラメータの決定法が変更されている。即ち、
第5図の流れ図では、第11図と較べて、前処理21の前に
乗算45が挿入されている。乗算45では、画像データに新
しいパラメータであるスケーリングフアクタSFを乗ず
る。前処理21以降では、SFを乗じた画像データに対して
第11図と同じ処理が行なわれる。スケーリングフアクタ
SFは第1図のステツプ42において第6図に示すように
画像データ階調数2nに反比例するように決定される。こ
こで比例定数は画像によらず定められたパラメータであ
る。このようにすると、第5図の乗算45でSFを乗じるこ
とにより、DR装置5〜7による画像データ階調数2nの違
いがキヤンセルされるので、前処理21以降の処理では、
DR装置の違いによるコントラストの違いを考えなくても
良いことになる。これに対応して、正規化係数SDi(i
=1〜4)は、画像によらず一定の値に固定することが
できる。又、許容歪dは第7図に示すように、標準表示
階調数WSの画像データ階調数2nに対する比WS/2nに比例
して決定される。ここでも、比例定数は画像によらず定
められたパラメータである。本実施例の方法で作成した
圧縮データから画像を再生する際、第13図に示したサブ
画像ごとの復号化処理の代りに、第8図の流れ図による
復号化処理が用いられる。即ち、第8図では、第13図と
較べて、双線形補間32の後に除算46が挿入されている。
除算46では、双線形補間32までで計算された再生画像デ
ータをスケーリングフアクタSFで除している。本実施例
の方法による圧縮データから双線形補間32までの処理で
再生される画像データは原画像データにSFを乗じたもの
に対応している。従つて、除算46でSFで除することによ
り、原画像データと同じ階調数の画像データとなる。本
実施例は、前記実施例と比較して、サブ画像ごとの符号
化・復号化処理に乗算45,除算46の処理が付加されてい
るため、その分処理時間が増す欠点がある。しかし、本
実施例によれば、符号化処理の本体である前処理21以降
の処理及び復号化処理の本体である双線形補間32以前の
処理において、画像データの階調数が一定になる。この
性質は、特に処理を行なう高速演算装置における演算情
報の一定化や、データのオーバフロー、アンダーフロー
の防止法が容易になる利点がある。
例も上記実施例と同様に、第9図のワークステーシヨン
1で、CPU12によつて第1図の流れ図に従つて行なわれ
る。但し、圧縮/再生モジユール15におけるサブ画像ご
との符号化処理の内容が第11図のものから若干異なる第
5図のものに変更されており、これに伴い、圧縮パラメ
ータ25にスケーリングフアクタSFが追加され、これを含
めた圧縮パラメータの決定法が変更されている。即ち、
第5図の流れ図では、第11図と較べて、前処理21の前に
乗算45が挿入されている。乗算45では、画像データに新
しいパラメータであるスケーリングフアクタSFを乗ず
る。前処理21以降では、SFを乗じた画像データに対して
第11図と同じ処理が行なわれる。スケーリングフアクタ
SFは第1図のステツプ42において第6図に示すように
画像データ階調数2nに反比例するように決定される。こ
こで比例定数は画像によらず定められたパラメータであ
る。このようにすると、第5図の乗算45でSFを乗じるこ
とにより、DR装置5〜7による画像データ階調数2nの違
いがキヤンセルされるので、前処理21以降の処理では、
DR装置の違いによるコントラストの違いを考えなくても
良いことになる。これに対応して、正規化係数SDi(i
=1〜4)は、画像によらず一定の値に固定することが
できる。又、許容歪dは第7図に示すように、標準表示
階調数WSの画像データ階調数2nに対する比WS/2nに比例
して決定される。ここでも、比例定数は画像によらず定
められたパラメータである。本実施例の方法で作成した
圧縮データから画像を再生する際、第13図に示したサブ
画像ごとの復号化処理の代りに、第8図の流れ図による
復号化処理が用いられる。即ち、第8図では、第13図と
較べて、双線形補間32の後に除算46が挿入されている。
除算46では、双線形補間32までで計算された再生画像デ
ータをスケーリングフアクタSFで除している。本実施例
の方法による圧縮データから双線形補間32までの処理で
再生される画像データは原画像データにSFを乗じたもの
に対応している。従つて、除算46でSFで除することによ
り、原画像データと同じ階調数の画像データとなる。本
実施例は、前記実施例と比較して、サブ画像ごとの符号
化・復号化処理に乗算45,除算46の処理が付加されてい
るため、その分処理時間が増す欠点がある。しかし、本
実施例によれば、符号化処理の本体である前処理21以降
の処理及び復号化処理の本体である双線形補間32以前の
処理において、画像データの階調数が一定になる。この
性質は、特に処理を行なう高速演算装置における演算情
報の一定化や、データのオーバフロー、アンダーフロー
の防止法が容易になる利点がある。
以上、本発明の実施例について述べてきたが、本発明
は上記実施例に限定されることなく種々の変形が可能な
ことは言うまでもない。例えば、第3図,第7図で設定
された比例関係は横軸にオフセツト値を持たせて一次関
数形に置き換えることも可能であり、場合によつてはも
つと一般的な広義単調増加関数形の方がよい結果となる
可能性がある。又、本発明の考え方は、実施例で用いた
特願昭62-297556号記載の符号化方法以外の符号化方法
にも適用できる。例えば、NHK技術月報、昭和60年5月
(1985年)第164頁から第171頁で論じられている2階層
符号化では、該文献の第168頁の第6図中の「モード選
択1」及び「モード選択2」で用いる閾値Th0及びTh1が
忠実度に影響するパラメータであるが、これらは上記実
施例におけるdやSDi(i=1〜4)と同様に、画像デ
ータ階調数2nや標準表示階調数WSに比例して決定するこ
とにより、適切な値が自動的に得られる。又、映像情報
(M)、18、(1986年)第719頁から第724頁で論じられ
ている離散コサイン変換法でも、該文献の第720頁の第
1図中の「係数の切捨て」で用いる閾値Th及び「係数の
量子化」で用いる量子化精度qは、上記実施例のdやSD
i(i=1〜4)と同様な性質を備えているので、上記
実施例と同様の取り扱いを行なうことにより、適切な値
が自動的に求められる。
は上記実施例に限定されることなく種々の変形が可能な
ことは言うまでもない。例えば、第3図,第7図で設定
された比例関係は横軸にオフセツト値を持たせて一次関
数形に置き換えることも可能であり、場合によつてはも
つと一般的な広義単調増加関数形の方がよい結果となる
可能性がある。又、本発明の考え方は、実施例で用いた
特願昭62-297556号記載の符号化方法以外の符号化方法
にも適用できる。例えば、NHK技術月報、昭和60年5月
(1985年)第164頁から第171頁で論じられている2階層
符号化では、該文献の第168頁の第6図中の「モード選
択1」及び「モード選択2」で用いる閾値Th0及びTh1が
忠実度に影響するパラメータであるが、これらは上記実
施例におけるdやSDi(i=1〜4)と同様に、画像デ
ータ階調数2nや標準表示階調数WSに比例して決定するこ
とにより、適切な値が自動的に得られる。又、映像情報
(M)、18、(1986年)第719頁から第724頁で論じられ
ている離散コサイン変換法でも、該文献の第720頁の第
1図中の「係数の切捨て」で用いる閾値Th及び「係数の
量子化」で用いる量子化精度qは、上記実施例のdやSD
i(i=1〜4)と同様な性質を備えているので、上記
実施例と同様の取り扱いを行なうことにより、適切な値
が自動的に求められる。
本発明によれば、再生画質の原画像に対する忠実度に
影響するパラメータの適切な値を自動的に決定すること
ができるので、従来の方法で行なわれていたパラメータ
の手入力が不要となり、操作性が向上すると共に、入力
誤りによる事故も防止できる。
影響するパラメータの適切な値を自動的に決定すること
ができるので、従来の方法で行なわれていたパラメータ
の手入力が不要となり、操作性が向上すると共に、入力
誤りによる事故も防止できる。
第1図は本発明の実施例における画像データ圧縮処理の
全体の手順を示す流れ図、第2図は実施例における対象
画像に関する情報の格納形態を示す説明図、第3図は第
1の実施例における許容歪と標準表示階調数の関係を示
すグラフ、第4図は第1の実施例における正規化係数と
画像データ階調数の関係を示すグラフ、第5図は第2の
実施例におけるサブ画像の符号化処理手順を示す流れ
図、第6図は第2の実施例におけるスケーリングフアク
タと画像データ階調数の関係を示すグラフ、第7図は第
2の実施例における許容歪と標準表示階調数/画像デー
タ階調数の関係を示すグラフ、第8図は第2の実施例に
おけるサブ画像の復号化処理手順を示す流れ図、第9図
は医用画像管理システムの一例を示すブロツク図、第10
図はサブ画像の概念を示す説明図、第11図は特願昭62-2
97556号の画像データ符号化方式によるサブ画像の符号
化手順を示す流れ図、第12図は全圧縮データの構成を示
す模式図、第13図は特願昭62-297556号の画像データ符
号化方式によるサブ画像の復号化手順を示す流れ図、第
14図は画像表示処理の構成と機能を示すブロツク図であ
る。 1……画像ワークステーシヨン、11……画像表示装置、
15……圧縮/再生モジユール、20……サブ画像、25……
圧縮パラメータ、40……画像データ階調数、41……標準
表示階調数、45……乗算、46……除算。
全体の手順を示す流れ図、第2図は実施例における対象
画像に関する情報の格納形態を示す説明図、第3図は第
1の実施例における許容歪と標準表示階調数の関係を示
すグラフ、第4図は第1の実施例における正規化係数と
画像データ階調数の関係を示すグラフ、第5図は第2の
実施例におけるサブ画像の符号化処理手順を示す流れ
図、第6図は第2の実施例におけるスケーリングフアク
タと画像データ階調数の関係を示すグラフ、第7図は第
2の実施例における許容歪と標準表示階調数/画像デー
タ階調数の関係を示すグラフ、第8図は第2の実施例に
おけるサブ画像の復号化処理手順を示す流れ図、第9図
は医用画像管理システムの一例を示すブロツク図、第10
図はサブ画像の概念を示す説明図、第11図は特願昭62-2
97556号の画像データ符号化方式によるサブ画像の符号
化手順を示す流れ図、第12図は全圧縮データの構成を示
す模式図、第13図は特願昭62-297556号の画像データ符
号化方式によるサブ画像の復号化手順を示す流れ図、第
14図は画像表示処理の構成と機能を示すブロツク図であ
る。 1……画像ワークステーシヨン、11……画像表示装置、
15……圧縮/再生モジユール、20……サブ画像、25……
圧縮パラメータ、40……画像データ階調数、41……標準
表示階調数、45……乗算、46……除算。
フロントページの続き (58)調査した分野(Int.Cl.6,DB名) G06T 1/00 G06T 5/00 - 5/50 H04N 1/40 - 1/419 A61B 5/00 A61B 6/00
Claims (1)
- 【請求項1】画像データの階調数と該画像データの標準
表示階調数の双方もしくは一方が既知である画像データ
を圧縮する際に、再生画像の原画像に対する忠実度に影
響する少なくとも1個のパラメータを含む画像データ圧
縮方法において、1枚の画像の画像データ全体の階調数
と該画像データ全体の標準表示階調数の双方もしくは一
方を前記パラメータの関数とすることを特徴とする画像
データ圧縮方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP1090727A JP2892676B2 (ja) | 1989-04-12 | 1989-04-12 | 画像データ圧縮方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP1090727A JP2892676B2 (ja) | 1989-04-12 | 1989-04-12 | 画像データ圧縮方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH02270078A JPH02270078A (ja) | 1990-11-05 |
JP2892676B2 true JP2892676B2 (ja) | 1999-05-17 |
Family
ID=14006588
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP1090727A Expired - Lifetime JP2892676B2 (ja) | 1989-04-12 | 1989-04-12 | 画像データ圧縮方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2892676B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100943273B1 (ko) * | 2003-05-07 | 2010-02-23 | 삼성전자주식회사 | 4-컬러 변환 방법 및 그 장치와 이를 이용한 유기전계발광표시장치 |
JP5159063B2 (ja) * | 2006-08-10 | 2013-03-06 | 株式会社東芝 | 医用画像保管装置 |
-
1989
- 1989-04-12 JP JP1090727A patent/JP2892676B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JPH02270078A (ja) | 1990-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6115486A (en) | Teleradiology system for the storage and transmission of angiographic and related image sequences | |
US20020090140A1 (en) | Method and apparatus for providing clinically adaptive compression of imaging data | |
JP2847525B2 (ja) | 画像データ圧縮方法 | |
US6144772A (en) | Variable compression encoding of digitized images | |
JP2001525622A (ja) | 画像圧縮方法 | |
US6909436B1 (en) | Radiologist workstation | |
US5884005A (en) | System for recording and displaying images on a record carrier in at least a first and a second data stream | |
US6215900B1 (en) | Image data processing for digital x-ray detectors | |
US20030123719A1 (en) | Medical image processing apparatus, medical image processing method, program and a recording medium thereof | |
JPH11262482A (ja) | 信号依存性のノイズを含む画像の処理方法 | |
US6468217B1 (en) | Method and apparatus for performing real-time storage of ultrasound video image information | |
JP3343070B2 (ja) | 画像出力システム | |
JP2892676B2 (ja) | 画像データ圧縮方法 | |
CN113538253B (zh) | 一种医疗影像的处理方法、介质及ct装置 | |
US4412248A (en) | Ultrasonic image storage device and method | |
JP4124380B2 (ja) | 画像処理再生装置 | |
JP2001520818A (ja) | Jpeg圧縮を使用して画像シーケンスの所望の平均圧縮比を発生する量子化係数値を自動的に決定するための方法および装置 | |
EP0083985B1 (en) | Improved ultrasonic image storage system | |
KR0155803B1 (ko) | 원격 의료진단시스템 | |
Cohen et al. | Image compression for radiology and telemedicine | |
JP2791558B2 (ja) | データ圧縮方法 | |
JPS62150475A (ja) | 画像処理装置 | |
JPH05249948A (ja) | 画像データ表示装置 | |
Skov et al. | Prototype system for archiving and transmitting digitized radiographs | |
JP2946683B2 (ja) | フィルムプリンタ装置 |