JPH06133171A - 画像データ圧縮方法および装置 - Google Patents

Abstract

(57)【要約】 【目的】 放射線画像を担持した画像データを圧縮、伸
長し、この伸長された画像データに基づいて放射線画像
を再生する際に、再生画像におけるモアレ発生を防止す
る一方、圧縮処理の高速化および圧縮率の向上を実現
し、画質（鮮鋭度）の低下も防止する。 【構成】 圧縮器33により圧縮処理を施す前の画像デー
タＳd1あるいはＳd2’の生成条件に応じて、圧縮方法を
変える。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、放射線画像を担持した
画像データを圧縮する方法、および装置に関するもので
ある。

【０００２】

【従来の技術】従来より、Ｘ線フイルムに撮影記録され
ている放射線画像をフイルムディジタイザ等により読み
取って電気的画像信号を得、この画像信号に基づいて放
射線画像をＣＲＴ等の表示装置に出力させることが行な
われている。

【０００３】一方、ある種の蛍光体に放射線（Ｘ線、α
線、β線、γ線、電子線、紫外線等）を照射すると、こ
の放射線エネルギーの一部が蛍光体中に蓄積され、この
蛍光体に可視光等の励起光を照射すると、蓄積されたエ
ネルギーに応じて蛍光体が輝尽発光を示すことが知られ
ており、このような性質を示す蛍光体は蓄積性蛍光体
（輝尽性蛍光体）と呼ばれる。

【０００４】この蓄積性蛍光体を利用して、人体等の放
射線画像情報を一旦蓄積性蛍光体の層を有するシートに
記録し、この蓄積性蛍光体シートをレーザ光等の励起光
で走査して輝尽発光光を生ぜしめ、得られた輝尽発光光
を光電的に読み取って画像信号を得、この画像信号に基
づき写真感光材料等の記録材料、ＣＲＴ等の表示装置に
放射線画像を可視像として出力させる放射線画像情報記
録再生システムが本出願人によりすでに提案されている
（特開昭55-12492号、同56-11395号等）。

【０００５】ところで、上記Ｘ線フイルムや蓄積性蛍光
体シート等の記録シートに被写体の放射線画像を撮影記
録する際には、被写体で散乱された放射線が記録シート
に照射されないように、被写体と記録シートとの間にグ
リッドを配置することが多い。このグリッドは、4.0 本
／mm程度の細かなピッチで例えば鉛とアルミニウム等が
交互に配置されてなるものである。このようなグリッド
を用いて撮影を行なうと、被写体により散乱された放射
線が記録シートに照射され難くなるので、放射線画像の
コントラストを向上させることができるが、記録シート
には、被写体像とともに細かな縞模様状のグリッド像が
記録される。

【０００６】一方、放射線画像が記録された記録シート
から画像信号を得る放射線画像読取装置は、通常、放射
線画像情報を担持した光を光電的に読み取り、画像情報
として必要な最高の空間周波数（以下、この必要な最高
の空間周波数をｆssで表わす。）に対応したサンプリン
グ間隔Δｘ＝１／（２・ｆss）でサンプリングしてディ
ジタルの画像信号を得るように構成される。このように
して得られた画像信号には、被写体の放射線画像を表わ
す有用な情報だけでなく、このグリッド像に起因するノ
イズが含まれることがある。

【０００７】以下、この点について説明する。図６(a)
は記録シートに記録された放射線画像（グリッド像を含
む）の、該グリッド像の縞模様の並び方向の空間周波数
特性を表わすグラフである。ここでは4.0 本／mmのグリ
ッドを用いて撮影を行ない、グリッド像の空間周波数が
４cycle ／mmにあるものとし、被写体の放射線画像の再
現に必要な最高の空間周波数ｆssは2.5cycle／mmである
として説明する。

【０００８】図６(b) は、被写体の放射線画像の再現に
必要なｆss＝2.5(cycle ／mm）以下の空間周波数帯の情
報を得るために、ｆss＝2.5(cycle ／mm）に対応したサ
ンプリング間隔Δｘ＝１／（２・ｆss）＝0.2 （mm）
で、すなわち１mmにつき５回サンプリングを行なった場
合のノイズの混入状態を説明する図である。この場合、
図に実線で示すグラフ（同図(a) と同一）がｆss＝2.5
(cycle ／mm) で折り返した位置に、いわゆるエリアジ
ングと呼ばれるノイズが混入する。したがってこの場
合、グリッド像の空間周波数４cycle ／mmのエリアジン
グノイズは、図中破線で示すように１cycle ／mmに生じ
る。

【０００９】図６(c) は、上記サンプリング間隔Δｘ＝
１／（２・ｆss）＝0.2 （mm）でサンプリングして得ら
れた画像信号が担持する放射線画像の空間周波数特性を
示している。図示されるように、１cycle ／mmの位置に
上記グリッド像に対応したノイズが混入するので、可視
画像を再生記録した際、該可視画像に１cycle ／mmの縞
模様として現われる結果となる。

【００１０】そこで例えば特開平3-114039号公報に示さ
れるように、Ｘ線フイルムや蓄積性蛍光体シート等の記
録シートを読取処理にかけて得られたアナログ画像信号
を、まず、所望とする空間周波数帯の最高空間周波数と
エリアジングの空間周波数とが一致するサンプリング間
隔よりも小さな間隔でサンプリングしてディジタルの初
期画像信号を得、次いでこの初期画像信号に対して、グ
リッドの空間周波数に対応する周波数成分を低減あるい
は除去するフィルタリング処理を施し、その後このフィ
ルタリング処理を受けた初期画像信号を、前記最高空間
周波数をナイキスト周波数とするサンプリング間隔で再
サンプリングして放射線画像再生用画像信号を生成する
ことも考えられている。

【００１１】他方、上記のような放射線画像を担持する
画像信号（画像データ）は、例えば光ディスク等の記録
媒体に格納（ファイリング）する際に一つの記録媒体に
より多くの画像を格納可能とするため、あるいは伝送効
率を高めるために、圧縮処理が施されることが多い。こ
の圧縮処理を受けた画像データは、記録媒体から読み出
された後、あるいは伝送された後に伸長処理が施された
上で、放射線画像再生のために供せられる。

【００１２】画像データの圧縮方法は従来より種々知ら
れているが、圧縮率の向上、および画質（鮮鋭度）低下
の防止を目的として、新しい圧縮方法が次々と提案され
ている。例えば特開昭62-181576 号公報には、放射線画
像を担持する原画像データを冗長度抑圧符号化処理する
放射線画像情報圧縮方法において、上記原画像データに
対して予め適当な方法で画像データ数を減少させる前処
理を施し、その前処理後の画像データを上記冗長度抑圧
符号化処理にかけるようにした方法が示されている。

【００１３】

【発明が解決しようとする課題】ところが従来より、圧
縮、伸長処理を受けた画像データから再生された放射線
画像においては、モアレが生じることがあるという問題
が認められている。このモアレは、前述したグリッドを
使用して撮影された放射線画像に多く発生する。

【００１４】そこで特開平2-280462号公報には、このモ
アレの発生を防止可能とした画像データ圧縮方法が提案
されている。この画像データ圧縮方法は、画像上に互い
に平行に設定された複数のブロックライン上の各々の上
に、１つの画素を含むブロックを所定間隔を置いて又は
複数の画素を含むブロックを連続的にもしくは所定間隔
を置いて設定すると共に該設定にあたっては各ブロック
ライン毎に位相をずらして設定し、各ブロック内の画素
の画像データに基づいて決定された各ブロックの代表画
像データのみを画像データとしてサンプリングするデー
タ数減少処理を原画像データに対して施し、上記データ
数減少処理が施された後の画像データを、適当な間隔で
サンプリングした主データと該主データ以外の補間デー
タとに区分し、該補間データに対し上記主データに基づ
く内挿予測符号化処理を施すことを特徴とするものであ
る。

【００１５】この画像データ圧縮方法は、モアレの発生
を防止する上では非常に効果的であるが、その半面、圧
縮処理に要する時間が長くなりがちで、圧縮率および画
質（鮮鋭度）の点でも改善の余地が残されている。

【００１６】本発明は上記のような事情に鑑みてなされ
たものであり、グリッド像によるモアレの発生を防止可
能で、その一方、圧縮処理の高速化および圧縮率の向上
を実現し、画質（鮮鋭度）の低下も防止できる画像デー
タ圧縮方法を提供することを目的とするものである。ま
た本発明は、その画像データ圧縮方法を実施する装置を
提供することも目的とする。

【００１７】

【課題を解決するための手段】本発明による画像データ
圧縮方法は、先に述べたように放射線画像を担持した画
像データを圧縮する方法において、圧縮処理を施す前の
画像データの生成条件に応じて、この圧縮処理の方法を
変えることを特徴とするものである。

【００１８】なお上記画像データの「生成条件」とは具
体的には、画像データがグリッド像を除去するための処
理、例えば、先に述べた特開平3-114039号公報に示され
るフィルタリング処理および再サンプリング処理等を受
けたか否か、ということを指すものである。

【００１９】一方上記方法を実施する本発明の画像デー
タ圧縮装置は、放射線画像を担持した画像データを圧縮
する装置において、相異なる複数の圧縮方法を実行可能
で、圧縮処理を施す前の画像データの生成条件に応じて
圧縮方法を変えるように形成されたことを特徴とするも
のである。

【００２０】

【作用および発明の効果】グリッド像を除去するための
処理を受けないで生成された画像データは、そのまま放
射線画像の再生のために供すれば、前述のモアレを発生
させる恐れがある。そこでこのような画像データに対し
ては、モアレ発生を抑止する効果のある圧縮方法、例え
ば前述の特開平2-280462号公報に示された圧縮方法を適
用すれば、伸長後の画像データに基づいて再生される放
射線画像にモアレが生じることを確実に防止できる。

【００２１】他方、グリッド像を除去するための処理を
受けて生成された画像データは、元よりモアレを発生さ
せる恐れがないものである。そこでこのような画像デー
タに対しては、モアレ発生を抑止する効果はなくても、
圧縮処理速度、圧縮率および画質（鮮鋭度）の点で優れ
ている圧縮方法、例えば前述の特開昭62-181576 号公報
に示された圧縮方法を適用すれば、圧縮処理の高速化、
圧縮率の向上および画質改善の効果が得られる。

【００２２】以上により本発明によれば、多数の放射線
画像を再生した場合の平均的な圧縮処理速度、圧縮率お
よび画質が向上し、その一方、再生放射線画像において
モアレが発生することも確実に防止可能となる。

【００２３】

【実施例】以下、図面に示す実施例に基づいて本発明を
詳細に説明する。図１は、本発明による画像データ圧縮
方法を実施する放射線画像情報読取再生システムの一例
を示している。このシステムにおいては、記録シートと
して前述の蓄積性蛍光体シートが用いられる。まず、図
２を参照して、蓄積性蛍光体シートへの放射線画像の撮
影記録について説明する。

【００２４】放射線源１から発せられた放射線２は、被
写体３を透過した後さらにグリッド４を経由して、蓄積
性蛍光体シート11に照射される。グリッド４は、４本／
mmのピッチで鉛４ａとアルミニウム４ｂとが交互に配置
されてなるものである。放射線２は鉛４ａには遮られ、
アルミニウム４ｂは透過して蓄積性蛍光体シート11に到
達する。このため蓄積性蛍光体シート11には、被写体像
とともに４本／mmの縞模様状のグリッド像が蓄積記録さ
れる。被写体３内で散乱された放射線２ａはグリッド４
に斜めに入射するので、そこで遮られ、あるいは反射し
て、蓄積性蛍光体シート11には多く照射されない。した
がって蓄積性蛍光体シート11には、散乱線の影響が少な
い鮮明な放射線画像が蓄積記録される。

【００２５】図３は、蓄積性蛍光体シート11に蓄積記録
された被写体像５（図の斜線部）と、縞模様状のグリッ
ド像６（図の縦縞）を概略的に示した図である。このよ
うに蓄積性蛍光体シート11には、被写体像５と縞模様６
とが重畳された放射線画像が記録される。

【００２６】なお、グリッド４を使用しないで放射線画
像が撮影記録されることもあり、その場合は当然、蓄積
性蛍光体シート11に上記縞模様６が記録されることはな
い。

【００２７】以上のようにして放射線画像が撮影記録さ
れた蓄積性蛍光体シート11は、図１のシステムの第１の
放射線画像情報読取装置31あるいは第２の放射線画像情
報読取装置32において、放射線画像情報の読取処理にか
けられる。ここでこれらの放射線画像情報読取装置31お
よび32には、グリッド４を使用して放射線画像が撮影記
録された蓄積性蛍光体シート11も、またグリッド４を使
用しないで放射線画像が撮影記録された蓄積性蛍光体シ
ート11も供給される。図４は第１の放射線画像情報読取
装置31の構成を詳しく示すものであり、以下この図４を
参照して、蓄積性蛍光体シート11からの放射線画像情報
読取りについて説明する。

【００２８】読取部10の所定位置にセットされた蓄積性
蛍光体シート11は、図示しない駆動手段により駆動され
るエンドレスベルト等のシート搬送手段15により、矢印
Ｙ方向に搬送（副走査）される。一方、レーザー光源16
から発せられた光ビーム17はモータ24により高速回転す
る回転多面鏡18によって反射偏向され、ｆθレンズ等の
集束レンズ19を通過した後、ミラー20により光路を変え
て蓄積性蛍光体シート11に入射し、このシート11上を副
走査方向（矢印Ｙ方向）と略垂直な矢印Ｘ方向に主走査
する。

【００２９】光ビーム17が照射された蓄積性蛍光体シー
ト11の部分からは、蓄積記録されている放射線画像情報
に応じた光量の輝尽発光光21が発散され、この輝尽発光
光２１は光ガイド２２によって導かれて、フォトマルチ
プライヤ（光電子増倍管）23によって光電的に検出され
る。上記光ガイド22はアクリル板等の導光性材料を成形
して作られたものであり、直線状をなす入射端面22ａが
蓄積性蛍光体シート11上の主走査線に沿って延びるよう
に配され、円環状に形成された射出端面22ｂにフォトマ
ルチプライヤ23の受光面が結合されている。入射端面22
ａから光ガイド22内に入射した輝尽発光光21は、この光
ガイド22の内部を全反射を繰り返して進み、射出端面22
ｂから射出してフォトマルチプライヤ23に受光され、輝
尽発光光21の光量がフォトマルチプライヤ23によって電
気信号に変換される。

【００３０】フォトマルチプライヤ23から出力されたア
ナログ出力信号Ｓには、良好な放射線画像を再生出力す
るために必要な、所望とする範囲の空間周波数帯のうち
の最高の空間周波数ｆss＝2.5cycle／mmより高い空間周
波数帯域の情報、特にこの信号Ｓがグリッド４を使用し
て放射線画像が撮影記録された蓄積性蛍光体シート11に
関するものである場合には、図３に示したグリッド像６
に関する情報も含まれている。このグリッド像６に関す
る情報は、可視像化された被写体３の放射線画像を見に
くくするものであり、取り除く必要のある情報である。
なお本例において、グリッド像の空間周波数は４cycle
／mmであるとする。

【００３１】上記アナログ出力信号Ｓはログアンプ26で
対数的に増幅された後、Ａ／Ｄ変換器28で所定のサンプ
リング間隔でサンプリングされてディジタル化され、デ
ィジタルの画像信号Ｓd1が得られる。この第１の放射線
画像情報読取装置31においては、Ａ／Ｄ変換器28のサン
プリング間隔Δｘ＝１／（２・ｆss）＝0.2 （mm）、す
なわち１mmにつき５回サンプリングとされており、そし
て特に上記グリッド像６を除去するような対策は施され
ないので、この画像信号Ｓd1にはグリッド像６に関する
情報が含まれることもある。この画像信号（画像デー
タ）Ｓd1は、第１の読取装置31で得られたことを示すコ
ードが付加された上で、圧縮器33において後述する圧縮
処理にかけられ、圧縮後の画像データＳd1’は光ディス
ク装置等の画像ファイリング装置34に格納される。

【００３２】図１に示される第２の放射線画像情報読取
装置32も、その基本構成は上記放射線画像情報読取装置
31と同様であり、そこからはディジタル画像信号（初期
画像信号）Ｓd2が出力される。

【００３３】ここで、第２の放射線画像情報読取装置32
におけるＡ／Ｄ変換器28のサンプリング間隔Δｘは0.1
（mm）、すなわち１mmにつき10回サンプリングとされて
おり、それにより得られたディジタル画像信号（初期画
像信号）Ｓd2は記憶部41に一時記憶される。この画像信
号Ｓd2は記憶部41から読み出されて画像信号生成部42に
入力され、そこで、グリッド像６の空間周波数４cycle
／mmに対応する信号成分を除去するフィルタリング処理
にかけられる。図５には、このフィルタリング処理の特
性を示してある。このフィルタリング処理を受けた画像
信号Ｓd2は、次いでこの画像信号生成部42において、再
サンプリング処理にかけられる。この再サンプリング処
理のサンプリング間隔Δｘは0.2 （mm）、すなわち１mm
につき10回サンプリングとされている。

【００３４】なお以上のフィルタリング処理および再サ
ンプリング処理は、前述した特開平3-114039号公報に詳
しく示されているものである。これらの処理を行なう
と、初期画像信号Ｓd2がグリッド４を使用して放射線画
像が撮影記録された蓄積性蛍光体シート11からのもので
あっても、処理後のディジタル画像信号Ｓd2’には、基
本的に、グリッド像６を示す成分が含まれなくなる。こ
の画像信号（画像データ）Ｓd2’は画像信号生成部42の
処理を受けたことを示すコードが付加された上で、記憶
部41に一時記憶され、次いで圧縮器33に入力される。そ
こで圧縮された画像データＳd2”は、画像ファイリング
装置34に格納される。

【００３５】画像ファイリング装置34に格納された画像
データＳd1’およびＳd2”は適宜読み出され、伸長器35
において所定の伸長処理を受けて各々伸長画像データＳ
e1およびＳe2とされる。これらの伸長画像データＳe1お
よびＳe2は、ＣＲＴ表示装置、光走査記録装置等の画像
再生装置36に入力される。この画像再生装置36は、伸長
画像データＳe1およびＳe2が担持する放射線画像を再生
する。

【００３６】次に、上記圧縮器33による圧縮処理につい
て詳しく説明する。この圧縮器33は、第１の放射線画像
情報読取装置31から出力されたことを示すコードが付加
された画像データＳd1に対しては、前述のモアレ発生を
抑止できる圧縮処理を行ない、その一方画像信号生成部
42の処理を受けたことを示すコードが付加されたディジ
タル画像データＳd2’に対しては、圧縮処理速度、圧縮
率および画質（鮮鋭度）の点で優れている圧縮処理を行
なう。

【００３７】まず、画像データＳd1に対する圧縮処理に
ついて、この処理の流れを示す図７を参照して説明す
る。なおこの圧縮処理は、前述の特開平2-280462号公報
に示されたものである。ここでは、画像データ（原画像
データ）Ｓd1が2000×2000画素分のデータからなり、各
画素のデータは８ビット０〜255 レベルであるとする。
なお、蓄積性蛍光体シート11に対する光ビーム17の走査
方向を図８に示してあり、各画素の副走査方向位置を
ｉ、主走査方向位置をｊで示して、各画素の原画像デー
タをＡ_ijと示す。

【００３８】最初に、上記原画像データＡ_ijに対して画
像データ数を減少させるべくデータ数減少処理を施すと
共に該減少処理を施すにあたって各ライン毎に位相をず
らせてサンプリングする位相シフトサンプリングを適用
する。

【００３９】図９は本実施例の位相シフトサンプリング
を説明する図であり、図中小さなます目は１つの画素を
示す。

【００４０】図示の位相シフトサンプリングは、２×２
画素ブロックから１つの代表画像データをサンプリング
して画像データ数を１／４に減少させると共に該代表画
像データとして各ブロック内の４つの画素の画像データ
の平均値を用いる２×２画素ブロック平均縮小サンプリ
ングであり、まず、図示の如く、画像上に、副走査方向
に２画素分の幅を有して主走査方向に延びる互いに平行
な複数のブロックラインＬ₁ ，Ｌ₂ ，Ｌ₃ ，Ｌ₄ ……を
設定し、各ブロックライン上に２×２画素ブロックを連
続的にかつ各ブロックライン毎に半ブロック長（１画素
分）づつ位相をずらして設定し、次に各ブロック毎に各
ブロック内の４つの画素の画像データの平均値を求めて
それを各ブロックの代表画像データ（平均縮小データ）
とし、その代表画像データのみを画像データとしてサン
プリングするものである。

【００４１】図10は上記サンプリングによって得られた
上記代表画像データを示す図であり、図中のます目は１
つの２×２画素ブロックを示し、各ます目の中のａ_ij，
ｂ_ij，ｃ_ij，ｄ_ijは各２×２画素ブロックの代表画像デ
ータである。

【００４２】次に、上記サンプリングによって得られた
代表画像データに対して補間符号化を行なう。まず、上
記代表画像データを、適当な間隔でサンプリングした主
データと補間データとに区分する。本実施例では、縦横
１ブロック置きに主データをサンプリングし、従って図
10中において○で囲った代表画像データａ_ijが主データ
であり、その他の代表画像データｂ_ij，ｃ_ij，ｄ_ijが補
間データである。

【００４３】次に、上記主データａ_ijに対し濃度分解能
を１ビット分低下させるべく量子化を行なう。具体的に
は、上記主データａ_ijを１ビット下位方向にシフトして
最下位ビットを削ることにより行う。

【００４４】続いて、上記量子化された主データａ_ij′
を各ブロックライン毎に前値予測を行なって予測誤差Δ
ａ_ijに変換する。ただし、各ブロックラインの先頭の主
データａ_ij′（j ＝１）はそのままにしておく。図11
は、かかる前値予測を行なった後の予測誤差（ただし先
頭のブロックについては主データａ_ij′（j ＝１）その
もの）を示す図である。なお、上記前値予測誤差Δａ_ij
は、 Δａ_ij＝ａ_ij′−ａ_ij-1′ であり、具体的には、 Δａ₁₂＝ａ₁₂′−ａ₁₁′，Δａ₁₃＝ａ₁₃′−ａ₁₂′，…
… Δａ₂₂＝ａ₂₂′−ａ₂₁′，Δａ₂₃＝ａ₂₃′−ａ₂₂′，…
… Δａ₃₂＝ａ₃₂′−ａ₃₁′，Δａ₃₃＝ａ₃₃′−ａ₃₂′，…
… である。

【００４５】続いて、この様にして求めた主データの各
予測誤差Δａ_ijに対しハフマン符号化を行なう。使用す
るハフマン符号表の一例を下の表１に示す。

【００４６】

【表１】

【００４７】すなわち、符号化対象データである上記予
測誤差Δａ_ijをこのハフマン符号表に従ってハフマン符
号に変換する。ただし、予測誤差が所定のしきい値（表
１のハフマン符号表では±７）を越える場合には、その
予測誤差Δａ_ijを表１中のエクステンション（ＥＸＴ）
コードの後に続けて、つまり0000001 ＋予測誤差Δａ_ij
の形で符号化する。

【００４８】そして、上記のように圧縮処理を施した主
データを前述の光ディスク等に格納する場合には、上記
各ブロックラインの先頭ブロック（ＴＯＬ：ＴＯＰ Ｏ
ＦＬＩＮＥ）については量子化後の主データａ_ij′（j
＝１）をそのままの形で、他のブロックについては上記
ハフマン符号表に基づいて符号化したハフマン符号およ
びＥＸＴコード＋予測誤差Δａ_ijの形で格納する。

【００４９】次に、上記補間データについて説明する。
まず上記各補間データｂ_ij，ｃ_ij，ｄ_ijについて主デー
タに基づく内挿予測を行ない、それらの予測誤差Δ
ｂ_ij，Δｃ_ij，Δｄ_ijを求める。かかる内挿予測は種々
の方法で行なうことができるが、本実施例では、主デー
タとして上記量子化後の主データａ_ij′を用い、内挿予
測式として下記の式 Δｂ₁₁＝ｂ₁₁−（ａ₁₁′×２＋ａ₁₂′×２）／２ Δｃ₁₁＝ｃ₁₁−（ａ₁₁′×２×３＋ａ₁₂′×２ ＋ａ₂₁′×２×３＋ａ₂₂′×２）／８ Δｄ₁₁＝ｄ₁₁−（ａ₁₁′×２＋ａ₁₂′×２×３ ＋ａ₂₁′×２＋ａ₂₂′×２×３）／８ を用いて、それぞれの内挿予測誤差Δｂ_ij，Δｃ_ij，Δ
ｄ_ijを求める。ただし、ａ₁₁′，ａ₁₂′，ａ₂₁′，
ａ₂₂′とΔｂ₁₁，Δｃ₁₁，Δｄ₁₁との位置関係は図12に
示す通りである。勿論、他の補間データの予測誤差も、
これらと同様の式で求める。

【００５０】次に、このようにして求めた各補間データ
の予測誤差Δｂ_ij，Δｃ_ij，Δｄ_ｉｊをさらに量子化
し、ビット分解能を１ビット分落とす。この場合、下の
表２に示す量子化特性Ａと表３に示す量子化特性Ｂとの
２つを用い、補間データｂ_ｉｊの予測誤差Δｂ_ijについ
ては量子化特性Ａを適用し、補間データｃ_ij，ｄ_ijの予
測誤差Δｃ_ij，Δｄ_ijについては量子化特性Ｂを適用し
て量子化する。

【００５１】

【表２】

【００５２】

【表３】

【００５３】続いて、上記量子化した後の各予測誤差Δ
ｂ_ij′，Δｃ_ij′，Δｄ_ij′に対しハフマン符号化を行
なう。このハフマン符号化も、上記主データのハフマン
符号化と同様に表１に示すハフマン符号表を用い、符号
化対象データである予測誤差をこのハフマン符号表に従
ってハフマン符号に変換する。また、予測誤差が±７を
超える場合には、その予測誤差Δｂ_ij′，Δｃ_ij′，Δ
ｄ_ij′を表１のＥＸＴコード0000001 の後に続けた形で
表わす。

【００５４】そして、上記の圧縮処理を施した補間デー
タを光ディスク等に格納する場合には、上記ハフマン符
号表に基づいて符号化したハフマン符号およびＥＸＴコ
ード＋内挿予測誤差Δｂ_ij′，Δｃ_ij′，Δｄ_ij′の形
で格納する。

【００５５】次に格納された圧縮画像データＳd1’を読
み出し、伸長器35で伸長する処理について説明する。圧
縮画像データＳd1’には以上述べた圧縮処理を示すコー
ドが圧縮器35において付加され、伸長器35はこのコード
が付加されている圧縮画像データＳd1’に対しては、以
下に説明する伸長処理を施す。図13はこの伸長処理を示
すフローチャートであり、各々の処理は図７に示した圧
縮処理の各ステップを逆に行なうものである。

【００５６】まず、主データの伸長について図14を参照
しながら説明する。上記の通り、主データの圧縮画像デ
ータは、ハフマン符号と、ＥＸＴコード＋量子化された
予測誤差Δａ_ijと、ＴＯＬブロックの量子化された主デ
ータａ_ij′そのままの３種類から成っている。従って、
まず読み出した圧縮画像データがそれらの３種類のどれ
であるかを検討し、ＴＯＬブロックのデータ（この場合
には格納する際ＴＯＬブロックのデータである旨の識別
符号を付して格納しておく）である場合にはそれからそ
のまま主データａ_ij′を求める。

【００５７】一方、ハフマン符号の場合には、該ハフマ
ン符号を表１のハフマン符号表に基づいて復号化するこ
とにより前述の前値予測誤差Δａ_ijを求め、該前値予測
誤差Δａ_ijと上記の既に求められているＴＯＬブロック
の主データａ_ij′（j ＝１）およびＥＸＴコード付デー
タの場合の主データａ_ij′とから前値予測復元を行なっ
て主データａ_ij′を求める。なお、この前値予測復元
は、 ａ_ij′＝ａ_ij-1′＋Δａ_ｉｊ なる式を用いて行なわれ、より具体的には、例えば各Ｔ
ＯＬブロックの主データａ_ｉｊ′（j ＝１）は既に求め
られているので、 ａ₁₂′＝ａ₁₁′＋Δａ₁₂，ａ₁₃′＝ａ₁₂′＋Δａ₁₃，…
… ａ₂₂′＝ａ₂₁′＋Δａ₂₂，ａ₂₃′＝ａ₂₂′＋Δａ₂₃，…
… ａ₃₂′＝ａ₃₁′＋Δａ₃₂，ａ₃₃′＝ａ₃₂′＋Δａ₃₃，…
… という式を用いて行なわれる。ＥＸＴコード付データの
場合には、そのＥＸＴコードの後に続いている予測誤差
Δａ_ijからハフマン符号の場合と同様に主データａ_ij′
を求める。

【００５８】上述のようにして全ての圧縮画像データＳ
d1’について主データａ_ij′を求めたら、次にそれを逆
量子化して、即ち１ビット上位にシフトすることによっ
て主データａ_ij（代表画像データ）を求める。

【００５９】次に、補間データの伸長について図15を参
照しながら説明する。上記の如く補間データの圧縮画像
データは、ハフマン符号と、ＥＸＴコード＋内挿予測誤
差Δｂ_ij′，Δｃ_ij′，Δｄ_ij′との２つから成ってい
る。従って、まず読み出した圧縮画像データがそのどち
らであるかを検討し、ハフマン符号の場合には、該ハフ
マン符号を表１のハフマン符号表に基づいて復号化する
ことにより前述の内挿予測誤差Δｂ_ij′，Δｃ_ij′，Δ
ｄ_ij′を求め、それらを逆量子化して、内挿予測誤差Δ
ｂ_ij，Δｃ_ij，Δｄ_ijを求め、これらと上記した既に求
められている前値予測復元主データａ_ijとから内挿予測
復元を行なって補間データ（代表画像データ）ｂ_ij，ｃ
_ij，ｄ_ijを求める。この内挿予測復元は下式 ｂ₁₁＝Δｂ₁₁＋（ａ₁₁＋ａ₁₂）／２ ｃ₁₁＝Δｃ₁₁＋（ａ₁₁×３＋ａ₁₂＋ａ₂₁×３＋ａ₂₂）／８ ｄ₁₁＝Δｄ₁₁＋（ａ₁₁＋ａ₁₂×３＋ａ₂₁＋ａ₂₂×３）／８ に基づいて行ない、他の８ビット補間データｂ_ij，
ｃ_ij，ｄ_ijも同様の式に基づいて求める。一方ＥＸＴコ
ード付のデータである場合にはそのＥＸＴコードの後に
続いているデータからそのままΔｂ_ij′，Δｃ_ij′，Δ
ｄ_ij′を求め、ハフマン符号の場合と同様に補間データ
ｂ_ij，ｃ_ij，ｄ_ijを求める。

【００６０】そして、上記の如くして主データと補間デ
ータとから成る代表画像データａ_ij，ｂ_ij，ｃ_ij，ｄ_ij
の全てを求めたら、次に図13に示す如くそれらに基づい
て補間拡大を行ない画像データ数を４倍にして元通りに
復元して原画像データに対応する伸長画像データＡ_ijを
求める。この伸長画像データＡ_ijの集合である画像デー
タＳe1は、先に述べたように画像再生装置36に入力さ
れ、放射線画像の再生に供せられる。

【００６１】上記画像データの補間拡大（復元補間）
は、例えば図16（図中の示す目は１つの画素を示す）に
示す様に、上記復元した各代表画像データａ_ij，ｂ_ij，
ｃ_ij，ｄ_ijをそれぞれのブロックの左上の画素の伸長画
像データＡ_ijとし、それらに基づいて他の伸長画像デー
タＡ_ijを主走査方向および副走査方向の直線補間により
復元すれば良く、例えば Ａ₁₂＝（ａ₁₁＋ｂ₁₁）／２ Ａ₂₂＝（Ａ₁₂＋ｃ₁₁）／２＝｛（ａ₁₁＋ｂ₁₁）／２＋ｃ₁₁｝／２ Ａ₃₃＝（ｃ₁₁＋ｄ₁₁）／２ Ａ₂₃＝（ｂ₁₁＋Ａ₃₃）／２＝｛ｂ₁₁＋（ｃ₁₁＋ｄ₁₁）／２｝／２ のようにして各伸長画像データＡ₁₂，Ａ₂₂，Ａ₂₃，Ａ₃₃
を復元し、他の伸長画像データＡ_ijも同様にして復元す
れば良い。

【００６２】なお、上記実施例では行なっていなかった
が、画像データ数減少処理の前に平滑化処理（スムージ
ング処理）を施しても良い。かかる平滑化処理の具体例
としては、例えば３×３画素サイズのマスクを用意し、
該マスクの中央に注目画素を位置させ、マスク内の９画
素の画像データの平均値をその注目画素の新たな画像デ
ータとし、かかる処理を各画素について行なう移動平均
フィルタ処理を挙げることができるが、勿論その他の平
滑化処理であっても良い。上記の如き平滑化処理を行な
うことにより、近傍画素の画像データ同志の相関を向上
させることができ、それによって上記前値予測や内挿予
測の予測誤差の零近傍への集中度を高め、圧縮率の向上
を図ることができる。

【００６３】上記の如く各ライン毎に位相をずらしてブ
ロックを設定する位相シフトサンプリングによりデータ
数減少処理を行なった後補間符号化処理を行なうように
すると、位相をずらさない同位相サンプリングによりデ
ータ数減少処理を行なった後補間符号化処理を行なう場
合よりも、圧縮率の向上を図ることができる。これは、
同位相サンプリングによる主データに基づくよりも位相
シフトサンプリングによる主データに基づいて補間デー
タを内挿予測した方が予測精度が向上することによるも
のと思われる。

【００６４】上述の如く位相シフトサンプリングによる
データ数減少処理を採用するとデータ圧縮率の向上を図
ることができるが、該位相シフトサンプリングは同位相
サンプリングに比して人間の視覚上粗さが目立たず、見
かけ上の鮮鋭度を向上させることができると共に、画像
データＳd1がグリッド４を使用して放射線画像が撮影記
録された蓄積性蛍光体シート11からのものである場合
に、再生放射線画像に前述のモアレが発生することを防
止できる。

【００６５】上記実施例の位相シフトサンプリングは２
×２画素ブロックを各ブロックライン毎に位相をずらせ
て設定したものであったが、該ブロックはｎ×ｎ画素ブ
ロック（ｎ≧３）であっても良く、また各ブロックの代
表画像データとしても各ブロック内の画像データの平均
値に限らず、例えば各ブロック内の特定位置の画素の画
像データそのものを用いることも可能である。

【００６６】また、上記実施例では主データに対して量
子化および前値予測＋ハフマン符号化を行なっている
が、この主データはそのままの形にしておいても良いし
他の適当な圧縮処理を施してもよい。また、量子化も種
々の態様の量子化を用い得るし、前値予測による予測誤
差を量子化するようにしてもよい。

【００６７】また、上記実施例では補間データに関し内
挿予測誤差の量子化を行なっているが、該量子化は、省
略してもよいし、他の態様で行なってもよいし、補間デ
ータそのものに対して行なってもよい。また、内挿予測
も他の方法（内挿予測式）で行なってもよい。勿論予測
誤差の符号化もハフマン符号表以外の符号表を用いて行
なってもよい。

【００６８】次に、画像データＳd2’に対する圧縮処理
について詳しく説明する。なおこの圧縮処理は、前述の
特開昭62-181576 号公報に示されたものである。図17
（ａ）はこの圧縮処理のプロセスフローチャートであ
る。図示の如く、まず原画像データＳd2’に対してデー
タ数を減少させる前処理を施し、この前処理によって減
少せしめられた画像データに対して冗長度抑圧符号化処
理を施す。この処理によって符号化された圧縮データ
は、前述したように画像ファイリング装置34において光
ディスク等の記録媒体に格納される。

【００６９】なお、この様にして圧縮格納されたデータ
に基づいて画像を再生する際は、図17（ｂ）の伸長プロ
セスフローチャートに示す様に、まず画像ファイリング
装置34の記録媒体から圧縮データＳd2”を読み出し、そ
の圧縮データＳd2”に対して上記冗長度抑圧符号化処理
において用いた符号表に従って復号処理を施し、復号化
された画像データに基づいてデータ数復元後処理を行な
う、即ちデータの数が前記データ数減少前処理によって
減らされ前のデータ数になるように上記復号処理された
画像データに基づいて幾つかの画像データを算出復元
し、その様にして復元された復元データ（勿論、この中
にはデータ数復元後処理によって復元されたものに限ら
ず、前記復号処理によって復元された画像データも含ま
れる）を画像再生装置36に送ればよい。

【００７０】上記データ数減少前処理は、要するに記録
すべき画像データの数を減少させる処理であって、その
処理の方法としては種々のものが考えられる。図18はか
かるデータ数減少前処理の一実施態様を説明する図であ
って、図中の各ます目はそれぞれ放射線画像の各画素を
示し、かつ各ます目中のＡ，Ｂ，Ｃ，……はそれぞれそ
の画素の画像データ（原画像データ）を示す。

【００７１】本例では図18（ａ）に示す如く、放射線画
像を２×２画素づつブロック分けし、各ブロック（太線
で囲まれたます目）Ｘ₁₁，Ｘ₁₂，……，Ｘ₂₁，Ｘ₂₂，…
…の所定位置（図示のものは○印を付した左上位置）の
画素の画像データＡ，Ｅ，……，Ｉ，Ｍ，……をそれぞ
れのブロックＸ₁₁，Ｘ₁₂，……，Ｘ₂₁，Ｘ₂₂，……の代
表画像データａ₁₁，ａ₁₂，……，ａ₂₁，ａ₂₂，……（ａ
₁₁＝Ａ，ａ₁₂＝Ｅ，……，ａ₂₁＝Ｉ，ａ₂₂＝Ｍ，……）
とし、この代表画像データのみを記録すべき画像データ
としてサンプリングするものである。この方法によれ
ば、結局画像データ数を１／４に減少させることができ
る。

【００７２】次に、このようにして減少せしめられた画
像データａ₁₁，ａ₁₂，……，ａ₂₁，ａ₂₂，……に対し
て、前述の如く冗長度抑圧符号化処理を施す。この冗長
度抑圧符号化処理は、要するに画像データが有している
情報の冗長性を抑圧して符号化する処理であって、代表
的なものとして予測符号化や直交変換による符号化等を
挙げることができる。

【００７３】ここでは前値予測とハフマン符号とを用い
て行なう予測符号化の場合について説明する。前記デー
タ数減少前処理によって求められた各ブロックの代表画
像データａ₁₁，ａ₁₂，……，ａ₂₁，ａ₂₂，……は図19に
示す如く並んでいる。図中の各ます目は、図18における
各ブロックＸ₁₁，Ｘ₁₂，……，Ｘ₂₁，Ｘ₂₂，……を示
す。

【００７４】まず、各画像データを前値予測し、その予
測値と実際のデータ値との差、即ち予測誤差を求める。
即ち、今代表画像データａ₁₁，ａ₁₂，ａ₁₃，……を持つ
各ブロックＸ₁₁，Ｘ₁₂，Ｘ₁₃，……の並びを第１ブロッ
クラインＬ₁ とし、代表画像データａ₂₁，ａ₂₂，ａ₂₃…
…を持つ各ブロックの並びを第２ブロックラインＬ₂、
代表画像データａ₃₁，ａ₃₂，ａ₃₃を持つ各ブロックの並
びを第３ブロックラインＬ₃ とし、かつ上記代表画像デ
ータａ₁₁，ａ₂₁，ａ₃₁はそれぞれ各ブロックライン
Ｌ₁ ，Ｌ₂ ，Ｌ₃ の先頭のブロックの代表画像データで
あるという前提の下に、第１ブロックラインＬ₁ におい
て、ａ₁₁はそのままとし、ａ₂₁についてはその前のデー
タ値ａ₁₁であると予測して実際のデータ値ａ₁₂との予測
誤差△ａ₁₂＝ａ₁₂−ａ₁₁を求め、ａ₁₃についてはその前
のデータ値ａ₁₂であると予測して実際のデータ値ａ₁₃と
の予測誤差△ａ₁₃＝ａ₁₃−ａ₁₂を求め、以下同様にして
△ａ₁₄，△ａ₁₅，……を求める。また、第２ブロックラ
インＬ₂ 、第３ブロックラインＬ₃ についても同様の方
法で先頭のａ₂₁，ａ₃₁はそのままとし、以下のａ₂₂，ａ
₂₃，ａ₂₄，……およびａ₃₂，ａ₃₃，ａ₃₄，……について
△ａ₂₂，△ａ₂₃，△ａ₂₄，……および△ａ₃₂，△ａ₃₃，
△ａ₃₄，……を求める。

【００７５】次に、この様にして求められたａ₁₁，△ａ
₁₂，△ａ₁₃，……，ａ₂₁，△ａ₂₂，△ａ₂₃，……，
ａ₃₁，△ａ₃₂，△ａ₃₃，……（図20参照）をハフマン符
号により符号化する。

【００７６】一般に、ある画像データ（注目画像デー
タ）はその近傍の画像データとあまり変わらないデータ
値を有する、すなわち両データ値にはあまり大きな差は
ないと考えることができ、従って注目画像データを近傍
の画像データ（上記の場合は前に位置する画像データ）
と同じ値であると予測し、そのときの予測誤差（前述の
△ａ₂₁，△ａ₂₂等）を求めれば、その予測誤差分布は零
の近傍に集中することになる。本予測誤差符号化は、こ
の様に予測誤差が零近傍に集中するという特性を利用
し、頻度の高いものには短い符号を、頻度の低いものに
は長い符号を割り当てるハフマン符号によりその予測誤
差および各ブロックラインの先頭ブロックのデータ
ａ₁₁，ａ₂₁，ａ₃₁，……を符号化することによって、結
局画像データの冗長度を抑圧し、トータルとしてのデー
タ量を圧縮することができるものである。下の表４にハ
フマン符号の一例を示す。

【００７７】

【表４】

【００７８】なお、上記例では符号化対象データが±13
を超えるものについては、エクステンションコード（Ｅ
ＸＴ）＋生データ、例えば符号化対象データが13の場
合、

【００７９】

【数１】

【００８０】という形で表現するように構成されてい
る。符号化対象データである前記予測誤差が例えば− 2
55〜＋ 255レベルまで作られた場合、本来ならばその全
ての符号化対象データに対してハフマン符号を割り当て
るのが圧縮率を向上させる意味では望ましいが、上述の
如く例えば±13を超えるような頻度の低いものについて
はＥＸＴ＋生データという形で表現するようにすると、
元々頻度が非常に低いものであるから圧縮率には殆んど
悪影響は生じないと共に符号表がコンパクトになるとい
う利点がある。

【００８１】上述の方法で符号化された圧縮データは、
画像ファイリング装置34において光ディスク等の記録媒
体に記録格納される。格納例を図21に示す。なお、図21
中のａ₁₁′，△ａ₁₂′……は、それぞれａ₁₁，△ａ₁₂，
……をハフマン符号化した圧縮データであり、ＴＯＩは
画像の始めを示す識別符号、ＴＯＬは各走査線（ブロッ
クライン）の始めを示す識別符号である。

【００８２】上記記録媒体に記録格納された圧縮データ
は、必要に応じて読み出され、可視像再生に供される。
可視像再生について簡単に説明すると、前述の如くまず
記録媒体から圧縮データを読み出し、その圧縮データを
上記圧縮データ作成の際に使用したハフマン符号表に基
づいて圧縮符号化される前の予測誤差データ等（ａ₁₁，
△ａ₁₂，△ａ₁₃，……，ａ₂₁，△ａ₂₂，△ａ₂₃，……，
ａ₃₁，△ａ₃₂，△ａ₃₃，……）に直すと共に該予測誤差
データ等から各ブロックの代表画像データを下記の如く
して算出復元する復号処理を行なう。

【００８３】 ａ₁₁，ａ₁₂＝ａ₁₁＋△ａ₁₂，ａ₁₃＝ａ₁₂＋△ａ₁₃，…… ａ₂₁，ａ₂₂＝ａ₂₁＋△ａ₂₂，ａ₂₃＝ａ₂₂＋△ａ₂₃，…… ａ₃₁，ａ₃₂＝ａ₃₁＋△ａ₃₂，ａ₃₃＝ａ₃₂＋△ａ₃₃，…… ・ ・ ・ ・ ・ ・ ・ ・ ・ 次に、上記の如くして復号、復元されたブロック代表画
像データに対してデータ数復元後処理を施し、画像デー
タ数をもとの原画像データＳd2’の数に戻す。

【００８４】このデータ数復元後処理によって復元され
た伸長画像データＳe2を図18（ｂ）に示す。この図にお
ける各ブロックＸ₁₁，Ｘ₁₂，Ｘ₂₁，Ｘ₂₂中の左上画素の
復元画像データＡ′，Ｅ′，Ｉ′，Ｍ′は前述の如く既
に代表画像データとして復元されており、Ａ′＝ａ₁₁，
Ｅ′＝ａ₁₂，Ｉ′＝ａ₂₁，Ｍ′＝ａ₂₂である。これら以
外の画素の画像データの復元にはいろいろな方法が考え
られるが、例えば水平方向および垂直方向の直線補間に
より復元することが考えられる。

【００８５】この方法により図18（ｂ）における画像デ
ータＢ′，Ｃ′，Ｄ′，Ｇ′，Ｊ′を復号すると次の通
りである。

【００８６】 Ｂ′＝（Ａ′＋Ｅ′）／２＝（ａ₁₁＋ａ₁₂）／２ Ｃ′＝（Ａ′＋Ｉ′）／２＝（ａ₁₁＋ａ₂₁）／２ Ｇ′＝（Ｅ′＋Ｍ′）／２＝（ａ₁₂＋ａ₂₂）／２ Ｊ′＝（Ｉ′＋Ｍ′）／２＝（ａ₂₁＋ａ₂₂）／２ Ｄ′＝（Ｃ＋Ｇ）／２＝（Ｂ＋Ｊ）／２＝（Ａ＋Ｅ＋Ｉ＋Ｍ）／４ ＝（ａ₁₁＋ａ₁₂＋ａ₂₁＋ａ₂₂）／４ そして、この様にして復元した復元データ（勿論、この
場合の直線補間による復元データＢ′，Ｃ′，Ｄ′，…
…は元の原画像データＢ，Ｃ，Ｄ，……と完全に一致す
るものではない）に基づいてＣＲＴ等に可視像を再生す
ればよい。

【００８７】上記データ数減少前処理は、図22に示すよ
うな方法で行なってもよい。すなわち、図18の場合と同
様に放射線画像を２×２画素づつブロック分けし、各ブ
ロックＸ₁₁，Ｘ₁₂，……，Ｘ₂₁，Ｘ₂₂，……内の各画素
の画像データの平均値を各ブロックの代表画像データａ
₁₁，ａ₁₂，……ａ₂₁，ａ₂₂，……とし、 ａ₁₁＝（Ａ＋Ｂ＋Ｃ＋Ｄ）／４， ａ₁₂＝（Ｅ＋Ｆ＋Ｇ
＋Ｈ）／４，…… ａ₂₁＝（Ｉ＋Ｊ＋Ｋ＋Ｌ）／４， ａ₂₂＝（Ｍ＋Ｎ＋Ｏ
＋Ｐ）／４，…… これらの代表画像データのみをサンプリングする方法で
行なってもよい。なお、この方法で前処理を行なった場
合は、前記伸長プロセスにおけるデータ数復元後処理に
おいて、例えば復元された各ブロックの代表画像データ
ａ₁₁，ａ₁₂，……，ａ₂₁，ａ₂₂，……をそれぞれのブロ
ックの特定位置、例えば左上位置の画素の画像データと
して取扱い、それらに基づいて前述と同様に水平、垂直
直線補間により他の画素の画像データを復元すればよ
い。

【００８８】この方法によると、データ数復元後処理に
よって復元された図22（ｂ）に示す画像データＡ′，
Ｅ′，Ｉ′，Ｍ′はそれぞれ Ａ′＝ａ₁₁， Ｅ′＝ａ₁₂， Ｉ′＝ａ₂₁， Ｍ′＝ａ
₂₂ となり、またＢ′，Ｃ′，Ｄ′，Ｇ′，Ｊ′は前記図18
の場合と同様の形となる。

【００８９】上記図18および図22に示すデータ数減少前
処理は、いずれもブロックを各ブロックラインＬ₁ ，Ｌ
₂ ，……毎に同位相で設定し、その各ブロックの代表画
像データをサンプリングする方法であったが、このブロ
ックの設定は各ブロックライン毎に位相をずらしたもの
であればより好ましい。なぜならば、位相をずらしてサ
ンプリングすると位相をずらさないでサンプリングした
場合に比べて、人間の視覚上粗さが目立たず、見かけ上
の鮮鋭度を向上させることができると共に、例えば散乱
線除去用の静止グリッドを用いた場合等に同位相でサン
プリングして画像データ数を減少させると復元された画
像上にモアレパターンが発生し易くなるが、位相をずら
せてサンプリングするとモアレパターンの発生を抑圧す
ることができるからである。

【００９０】以下に、位相をずらせてサンプリングする
例、即ち位相シフトサンプリングの代表的な例を図23〜
図25を参照しながら説明する。

【００９１】図23に示す位相シフトサンプリングは、放
射線画像上に互いに平行に設定された複数のブロックラ
インＬ₁ ，Ｌ₂ ，Ｌ₃ ……の各々の上に、１つのます目
で示される１つの画素を含んで成るブロックＸ₁₁，
Ｘ₁₂，Ｘ₁₃，……，Ｘ₂₁，Ｘ₂₂，Ｘ₂₃，……，Ｘ₃₁，Ｘ
₃₂，Ｘ₃₃，……を所定間隔を置いて（この場合は、１ブ
ロックに対応する間隔を置いて）設定すると共に、該ブ
ロックの設定に当っては各ブロックライン毎に位相をず
らして設定し、その様に設定された各ブロック内の画素
の画像データに基づいて決定された代表画像データ
ａ₁₁，ａ₁₂，ａ₁₃，……，ａ₂₁，ａ₂₂，ａ₂₃，……，ａ
₃₁，ａ₃₂，ａ₃₃，……のみをサンプリングするものであ
る。なお、この場合の代表画像データａ₁₁，ａ₁₂，……
は、例えばそのブロック内画素の画像データそのものを
使用し、 ａ₁₁＝Ａ，ａ₁₂＝Ｂ，ａ₁₃＝Ｃ，…… ａ₂₁＝Ｄ，ａ₂₂＝Ｅ，ａ₂₃＝Ｆ，…… ａ₃₁＝Ｇ，ａ₃₂＝Ｈ，ａ₃₃＝Ｉ，…… とすればよい。

【００９２】なお、この場合の伸長プロセスでのデータ
数復元後処理におけるデータが記録されていない画素の
画像データ復元は、例えば周囲４画素の平均値を取るこ
とによって行なえばよい。つまり、図中の例えばブロッ
クＸ₁₂，Ｘ₂₁，Ｘ₂₂，Ｘ₃₂によって囲まれている部分の
画素の画像データａは、 ａ＝（Ｂ＋Ｄ＋Ｅ＋Ｈ）／４ によって求めればよい。

【００９３】図24に示す位相シフトサンプリングは、同
図（ａ）に示す如く設定するブロックが２×２画素を含
んで成るものであり、各ブロックラインＬ₁ ，Ｌ₂ ，Ｌ
₃ ，……に連続的にかつ各ブロックラインＬ₁ ，Ｌ₂ ，
Ｌ₃ 毎に半ブロック長づつ位相をずらしてブロック
Ｘ₁₁，Ｘ₁₂，Ｘ₁₃，……，Ｘ₂₁，Ｘ₂₂，……，Ｘ₃₁，Ｘ
₃₂，Ｘ₃₃，……を設定し、各ブロック内の４画素の画像
データの平均値をそれぞれのブロックの代表画像データ
ａ₁₁，ａ₁₂，ａ₁₃，……，ａ₂₁，ａ₂₂，……，ａ₃₁，ａ
₃₂，ａ₃₃，……として図24（ｂ）に示す如くサンプリン
グするものである。ここで、ａ₁₁，ａ₁₂について具体例
を示すと以下の通りである。 ａ₁₁＝（Ａ＋Ｂ＋Ｃ＋Ｄ）／４， ａ₁₂＝（Ｅ＋Ｆ＋Ｇ＋Ｈ）／４ なお、この場合の伸長プロセスにおける画像データ数復
元後処理は、例えば各ブロックの代表画像データａ₁₁，
ａ₁₂，……，ａ₂₁，ａ₂₂，……，ａ₃₁，ａ₃₂，……を図
24（ｃ）に示す如くそれぞれのブロックの特定位置、例
えば左上位置の画素の画像データとし、他の画素の画像
データはそれらの代表画像データから以下の様に水平、
垂直直線補間により算出復元すればよい。

【００９４】 Ａ′＝ａ₁₁，Ｅ′＝ａ₁₂，Ｉ′＝ａ₂₁，Ｋ′＝ａ₂₂， Ｂ′＝（Ａ′＋Ｅ′）／２＝（ａ₁₁＋ａ₁₂）／２ Ｄ′＝（Ｂ′＋Ｉ′）／２＝｛（ａ₁₁＋ａ₁₂）／２＋ａ₂₁｝／２ Ｊ′＝（Ｉ′＋Ｋ′）／２＝（ａ₂₁＋ａ₂₂）／２ Ｇ′＝（Ｅ′＋Ｊ′）／２＝｛ａ₁₂＋（ａ₂₁＋ａ₂₂）／２｝／２ 図25に示す位相シフトサンプリングは、同図（ａ）に示
す如く設定するブロックが２×２画素を含んで成るもの
であり、各ブロックラインＬ₁ ，Ｌ₂ ，Ｌ₃ ，……に１
ブロック分の所定間隔を置いてかつ各ブロックラインＬ
₁ ，Ｌ₂ ，Ｌ₃毎に１ブロック長づつ位相をずらしてブ
ロックＸ₁₁，Ｘ₁₂，……，Ｘ₂₁，Ｘ₂₂，……，Ｘ₃₁，Ｘ
₃₂，……を設定し、各ブロック内の４画素の画像データ
の平均値をそれぞれのブロックの代表画像データａ₁₁，
ａ₁₂，……，ａ₂₁，ａ₂₂，……，ａ₃₁，ａ₃₂，……とし
て図25（ｂ）に示す如くサンプリングするものである。 ａ₁₁＝（Ａ＋Ｂ＋Ｃ＋Ｄ）／４ なお、この場合の伸長プロセスでのデータ数復元後処理
は、例えば各ブロックの代表画像データａ₁₁，ａ₁₂，…
…，ａ₂₁，ａ₂₂，……，ａ₃₁，ａ₃₂，……を図25（ｃ）
に示す如くそれぞれのブロックの特定位置、例えば左上
位置の画素の画像データとし、次に４つのブロックで囲
まれている２×２画素部分の左上画素の画像データとし
て該囲っている４つのブロックの代表画像データの平均
値を割り当てる、例えばブロックＸ₁₂，Ｘ₂₁，Ｘ₂₂，Ｘ
₃₂で囲まれている２×２画素部分の左上画像の画像デー
タａとして ａ＝（ａ₁₂＋ａ₂₁＋ａ₂₂＋ａ₃₂）／４ を割り当て、その後その他の画素の画像データを上記の
各代表画像データａ₁₁，ａ₁₂，……とそれらに基づいて
求めた画像データａとから水平および垂直直線補間によ
って求めればよい。

【００９５】以上、種々のデータ数減少前処理の実施態
様について説明したが、この前処理はその他にも種々の
態様を取り得、例えば図18、図22、図24、図25の例では
２×２画素のブロックを設定しているが、これをｎ×ｎ
（ｎは３以上の整数）画素のブロックとしてもよいし、
あるいは図24、図25においてブロックの代表画像データ
をブロック内画素の画像データの平均値としているが、
例えば図18の例の場合の如くブロック内特定位置の画素
の画像データをブロック代表画像データとすることも可
能であり、また図23、図25の例では要するにブロックを
市松状に設定してその他の部分の画像データは捨てるよ
うな処理を行なうが、この画像上に部分的にブロックを
設定し、他の部分の画像データを捨てる場合のブロック
の設定態様も上記市松状以外の種々のものを採用し得
る。

【００９６】また、上記冗長度抑圧符号化処理は、前述
の如く予測符号化に限られるものではなく、他の処理例
えば直交変換による符号化処理であってもよい。また、
予測符号化においても、予測方法は前述の前値予測に限
らず近傍画素の画像データを用いた２次元予測等種々の
予測方法を採用することができ、予測誤差の符号化もハ
フマン符号に限らず、他の符号を使用することができ
る。

【００９７】上述した画像データ圧縮方法は、前述の如
く、冗長度抑圧符号化処理を行なう前に予め画像データ
数そのものを減らし、即ち空間分解能を減らし、その減
少せしめられた画像データのみを冗長度抑圧符号化処理
をするので、画像データ数減少前処理を行なうことなく
冗長度抑圧符号化処理する場合に比べて、終局的に大き
なデータ圧縮率を達成することができる。

【００９８】また、人体等の放射線画像は、他の一般的
な画像に比べ、診断対象等の重要部分は低空間周波数域
に存在し、高空間周波数成分はあまり持っていないとい
う本発明者の研究結果により、結局上記の如く空間分解
能を減少させて高空間周波数成分を消失させても、例え
ば、濃度分解能を低下させる場合に比べて診断能の低下
等の画質の劣化が非常に少ない。

【００９９】またこの画像データ圧縮方法は、先に述べ
たモアレ発生を抑止する作用も有する画像データ圧縮方
法に比べれば、より高速で画像データ圧縮処理を行なう
ことができる。

【０１００】そして画像データＳd2’にこのような画像
データ圧縮方法を適用しても、このディジタル画像デー
タＳd2’は前述のフィルタリング処理および再サンプリ
ング処理を受けて、元よりグリッド像６の成分を含まな
いものであるから、圧縮、伸長処理後の画像データＳe2
に基づいて再生される放射線画像にモアレが生じること
はない。

【０１０１】なお、上述したグリッド像除去のためのフ
ィルタリング処理および再サンプリング処理を行なう信
号生成系（上記実施例では第２の放射線画像情報読取装
置32を含む系）において、これらの処理を行なうモード
と行なわないモードとを適宜選択できるようにしてもよ
い。その場合は、それらのモードを選択するスイッチの
状態が画像信号の生成条件と対応するから、このスイッ
チと連動させて圧縮処理の方法を変えるようにすればよ
い。

【図面の簡単な説明】

【図１】本発明方法を実施する放射線画像情報読取再生
システムを示すブロック図

【図２】蓄積性蛍光体シートに放射線画像を撮影記録す
る様子を示す概略図

【図３】グリッドが記録された放射線画像を示す概略図

【図４】上記システムの放射線画像情報読取装置を示す
斜視図

【図５】グリッド像を除去するフィルタリング処理の特
性を示すグラフ

【図６】グリッドが記録された放射線画像の空間周波数
特性を示すグラフ

【図７】本発明に適用される一つの画像データ圧縮処理
のフローチャート

【図８】蓄積性蛍光体シートからの画像データの読み取
り方式を示す概略図

【図９】図７に示す画像データ圧縮処理のデータ数減少
処理（位相シフトサンプリング）の態様を示す説明図

【図１０】主データのサンプリング態様を示す説明図

【図１１】主データの前値予測を説明する説明図

【図１２】補間データの内挿予測を説明する説明図

【図１３】図７の方法によって圧縮した画像データの伸
長方法を示すフローチャート

【図１４】図１３のフローチャートのうち主データに関
する部分の詳細フローチャート

【図１５】図１３のフローチャートのうち補間データに
関する部分の詳細フローチャート

【図１６】図１３に示す伸長処理における補間拡大によ
るデータ数復元方法を説明するための図

【図１７】本発明に適用されるもう一つの画像データ圧
縮処理と、それに対応する伸長処理を示すプロセスフロ
ーチャート

【図１８】図１７に示した処理におけるデータ数減少前
処理と、この前処理に対応するデータ数復元後処理の一
例の説明図

【図１９】上記前処理によって減じられた後の画像デー
タを示す説明図

【図２０】図１９の画像データに関する前値予測誤差等
を示す図

【図２１】図２０に示す前値予測誤差等をハフマン符号
化したものの格納状態を示す図

【図２２】データ数減少前処理の別の例と、その前処理
に対応するデータ数復元後処理の一例の説明図

【図２３】位相シフトサンプリングタイプのデータ数減
少前処理の説明図

【図２４】位相シフトサンプリングタイプのデータ数減
少前処理の別の例と、それに対応する各ブロックの代表
画像データを示す説明図

【図２５】位相シフトサンプリングタイプのデータ数減
少前処理のさらに別の例と、それに対応する各ブロック
の代表画像データを示す説明図

【符号の説明】

31、32 放射線画像情報読取装置 33 圧縮器 34 画像ファイリング装置 35 伸長器 36 画像再生装置 41 記憶部 42 画像信号生成部

Claims (2)

【特許請求の範囲】
1. 【請求項１】 放射線画像を担持した画像データを圧縮
   する方法において、圧縮処理を施す前の画像データの生
   成条件に応じて、この圧縮処理の方法を変えることを特
   徴とする画像データ圧縮方法。
2. 【請求項２】 放射線画像を担持した画像データを圧縮
   する装置において、複数の圧縮方法を実行可能で、圧縮
   処理を施す前の画像データの生成条件に応じて圧縮方法
   を変えるように形成されたことを特徴とする画像データ
   圧縮装置。