JP3303504B2 - X線画像処理装置 - Google Patents

X線画像処理装置

Info

Publication number
JP3303504B2
JP3303504B2 JP02956294A JP2956294A JP3303504B2 JP 3303504 B2 JP3303504 B2 JP 3303504B2 JP 02956294 A JP02956294 A JP 02956294A JP 2956294 A JP2956294 A JP 2956294A JP 3303504 B2 JP3303504 B2 JP 3303504B2
Authority
JP
Japan
Prior art keywords
ray image
image
pixel
ray
value
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
Application number
JP02956294A
Other languages
English (en)
Other versions
JPH07236626A (ja
Inventor
昌昭 浮田
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.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
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 Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP02956294A priority Critical patent/JP3303504B2/ja
Publication of JPH07236626A publication Critical patent/JPH07236626A/ja
Application granted granted Critical
Publication of JP3303504B2 publication Critical patent/JP3303504B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、例えばX線画像表示装
置、X線量子計数器あるいは胸部デジタルX線画像診断
装置等のX線画像のデータ処理に利用される画像処理装
置に関する。
【0002】
【従来の技術】胸部X線画像は、3桁程度の広いダイナ
ミックレンジを持ち、その表示が困難である。そこで、
このような胸部X線画像を表示するため、CR(Compute
r Ra-diography)などでは、空間周波数フィルタ処理
(ぼけマスク処理など)が行われている。これらの処理
は、特定の空間周波数帯域(高周波数帯域)を強調した
デジタル処理で、線状の陰影を強調する画像処理(エッ
ジ強調処理)である。
【0003】なお、このほか、胸部X線画像を表示する
ものとして、感度特性の異なる2種類の材料を組み合わ
せたX線フィルムも開発されている。
【0004】
【発明が解決しようとする課題】ところで、上記したデ
ジタル画像処理では、線状の陰影が強調されるので読影
はし易くなるが、しゅ瘤などの大きな陰影の読影は難し
く、また、従来のレントゲン写真に比べて画像品位が大
きく相違する画像(不自然な画像)となる等の欠点があ
る。そのため、この種のデジタル画像処理を採用したX
画像撮像装置(FCR;富士コンピュータラジオグラフ
ィ)などでは、医師らに診断時には必ず2枚の画像を読
影することを提示するするようにしている。しかし、こ
れは必ずしも良い処置とは言えない。
【0005】本発明はそのような事情に鑑みてなされた
もので、被検体にX線を一様に照射して検出された透過
X線画像を表示するにあたり、自然で読影し易い画像を
得ることのできる画像処理装置の提供を目的とする。
【0006】
【課題を解決するための手段】上記の目的を達成するた
め、本発明は、実施例に対応する図1に示すように、X
線画像Pの一つの画素pi,j とその周辺の画素pi,j-L
・・pi,j+L の各透過線量値を用い、その中央画素pi,j
の重み係数を最大として荷重平均計算を行って透過線量
の平均値Si,j を求める荷重平均計算器1と、その平均
値Si,j を変数とする関数で、当該平均値Si,j の小→
大の変化に応じて減少する関数(例えば√(C/Si,j
))を用いて減弱変換倍率gi,j を求める減弱変換器
2と、この変換倍率gi,j と上記の中央画素pi,j の透
過線量値Pi,j との掛け算を行う積算器3を備え、その
演算結果を当該X線画像Pの1画素分の画像データとし
て用いるように構成したことによって特徴づけられる。
【0007】
【作用】荷重平均計算器1では、例えば、X線画像Pの
ある1画素pi,j について区間pi,j-L ・・pi,j ・・pi,
j+L の範囲の平均透過線量Si,j を、例えば図2に示す
荷重係数を用いた重み付け計算、すなわち処理対象の画
素pi,j を最大とする重み付け計算により求め、減弱変
換器2では、その平均値Si,j を、関数(g=√(C/
S))により変換して、元の透過線量値Pi,j を目標検
出線量;Cに近づけるのに必要な倍率gi,j を得る。
【0008】そして、倍率gと元の透過線量値Pi,j と
の積を積算器3で求め、このような計算をX線画像Pの
全ての画素について行うことにより画像全体の処理デー
タを得る。
【0009】ここで、本発明の画像処理では、X線画像
Pのある画素の元の透過線量値が、目標検出線量;Cに
近い場合には、その値は処理後でも殆ど変化せずC付近
の値となるが、元の透過線量値がCに対して十分に低い
場合には、その処理後の値は積算器3での倍数gの掛け
算により元の値よりも大きくなる(図3参照)。従っ
て、積算器3は、空間的な低周波成分に対しては、あた
かもX線量を増やしたかのように機能することになり、
これによりダイナミックレンジの幅を所定範囲内に抑え
ることが可能になる。
【0010】しかも、本発明では、荷重平均算出器1で
の重み付け計算により、一つの画素の周辺画素に対する
減弱度が抽出されているので、小さな構造も読影するこ
とができるX線画像を得ることも可能となる。
【0011】
【実施例】本発明の実施例を、以下、図面に基づいて説
明する。図1は本発明実施例の構成図である。この図で
は、処理の説明を簡単とするために、ある一つの画素p
i,j の処理結果Qi,j を得るときの処理を行う構成のみ
を示している。また、同様に処理説明を簡単とするた
め、この例では透過X線画像Sに対して縦方向の処理の
みを行うものとする。
【0012】まず、この例の画像処理装置は、荷重平均
計算器1、減弱変換器2及び積算器3によって構成され
ている。荷重平均計算器1は、X線画像Sの画素pi,j
を中央とする画素区間pi,j-L・・pi,j+L の各画素値
(透過線量)を用いて平均透過線量を重み付け計算によ
り求めるための計算器である。その計算式は
【0013】
【数1】
【0014】が用いられる。ここで、Wjsは重み付けの
ための荷重係数で、この例では、次の(2) 式及び図2の
グラフに示す関数を用いる。 Wjs=( js+L )/(L2),〔−L ≦js<0〕 (−js+L )/(L2),〔 0≦js≦L〕 ・・・・(2) 減弱変換器2は、荷重平均計算器1で計算された平均透
過線量Wjsを変換する倍率gi,j を 関数;gi,j =√(1/(Si,j /C)),C=定数 を用いて求める。なお、定数Cとして、被写体がないと
きの検出値、また肺野と縦隔との境界値あるいは肺野の
画素データの最高値等を用いる。
【0015】積算器3では、減弱変換器2で求めた倍率
gi,j と、画素pi,j の元の透過線量値Pi,j とが掛け
算され、その演算結果;Qi,j =Pi,j *gi,j が画像
処理結果として出力される。
【0016】そして、以上の計算をX線画像の全ての画
素について行うことにより、X線画像の全体の処理デー
タを得ることができ、このような処理において、荷重平
均計算器1は画像データ群に対して空間的なロウパスフ
ィルタとして作用することになり、これにより空間的な
低周波成分をある程度の範囲内で一定にしてダイナミッ
クレンジの幅を所定の範囲内に抑えることができる。
【0017】その作用効果を模式的に説明すると、例え
ば図3に示すように、X線画像の元の透過線量データP
が広いダイナミックレンジを持つ場合、このデータを、
本発明実施例の荷重平均計算器1、減弱変換器2及び積
算器3を用いて処理すると、その処理結果Qは、定数C
に近い透過線量ではあまり変化しないが、定数Cに比べ
少ない透過線量では変化が大きく、定数Cに近づく。そ
の結果、ダイナミックレンジの幅が、元のX線画像より
も狭くなる。しかも、この例での画像処理ではX線画像
の全体的な変化も、元の画像の状態が反映されるので、
自然で読影し易い画像を得ることができる。
【0018】ここで、減弱変換器2の関数は、g=C/
Sとしてもよいが、この場合、低周波成分がほぼ一定と
なって小さな構造は明瞭になるものの、図3の模式的線
図のQ′に示すように、空間的な低周波成分がほば同一
のレベルとなってしまい、得られる画像がフィルムに比
べて不自然になるといった問題が発生する。そこで、自
然で読影し易い画像を得ることを目的とする場合には、
減弱変換器2の関数として、例えば、先に述べたg=√
(C/S)を用いることが好ましい。
【0019】その理由を、以下、図4に示すような均一
な画像データを参照して説明する。まず、図に示すデー
タでは、元のX線画像より荷重平均計算器1の出力がS
=P(透過線量の計数値)、減弱変換器2により倍率g
=√(C/P)、積算器3の出力はQ=g・P=√(C
/P)・P=√(P)・√(C)となる。
【0020】さて、画素のゆらぎ(ポアソン分布)を考
慮すると計数Pの値はP±√(P)程ゆらいでいるが、
上記の積算を行った場合のゆらぎを求めると、 となり、目標計数値Cとしたときのゆらぎは、C±√
(C)であるから、g=√(C/P)で減弱変換を行っ
ても、ゆらぎは等しくなる。これにより、先の図3の模
式的線図に示したように画像の不自然さがなくなる。
【0021】なお、以上の本発明実施例では、画像に対
する処理の方向を縦方向のみとしているが、その方向は
横方向であってもよいし、また、縦方向の処理を行った
後、横方向の処理を行うといった手法を採用してもよ
い。さらに、空間的に拡張した画像処理つまり一つの画
素についてその縦と横の双方の周辺画素の画素値を用い
たり、領域内の平均透過線量の計算を行うといった処理
手法を採用してもよい。
【0022】また、以上の実施例では、X線画像Pの一
つの画素pi,j についての処理のみを示しているが、他
の画素についても画素pi,j と同等な構成により処理を
行うことは言うまでもない。
【0023】ここで、本発明の処理装置において、荷重
平均計算器の計算に適用する重み付けの係数としては、
先の(2) 式の関数のほか、ハミング窓関数あるいはハニ
ング窓関数が挙げられる。ただし、(2) 式の関数を採用
した場合、荷重平均計算器での演算速度が高速になり、
処理速度が速くなるといった利点がある。
【0024】すなわち、荷重係数として(2) 式のような
線形の関数(図2参照)を採用すると、重み付け計算を
行う際の掛け算(Wi,j *Pi,J+JS)は、X線画像のあ
る画素区間についてのみ行い、この代表区間の計算結果
を基準として加減算を行うことにより画像の全ての画素
についての平均透過線量を求めるといった手法が可能と
なり、その掛け算を削減できる分だけ演算時間を大幅に
短縮できる。
【0025】また、本発明の処理装置において、減弱変
換器に適用する関数としては、先に述べたg=C/Sや
g=1/S、さらに√(C/S)を一般化して(C/
S)Kあるいはシグモイド関数を用いた1+EXP(−k
1(S−k2))などが挙げられる。
【0026】
【発明の効果】以上説明したように、本発明のX線画像
処理装置によれば、X線画像の空間的な低周波成分をあ
る程度の範囲内で一定にして、ダイナミックレンジの幅
を所定の範囲内に抑えることが可能となるので、自然で
読影し易い画像を得ることができる。これにより、1枚
の画像で診断に必要な読影を十分に行うことが可能とな
り、その読影の能率が向上する。
【0027】さらに、本発明によると、減弱変換器で用
いる関数を適宜に変更するだけで、ダイナミックレンジ
の特性を簡単かつ任意に変えることができる。これによ
り、例えば読影を行いたい部位の陰影を強調した画像を
得ることができる等、診断の目的に合った画像を作成す
ることも可能になる。
【図面の簡単な説明】
【図1】本発明実施例の構成図
【図2】その実施例の荷重平均計算器1の計算に用いる
荷重係数を示す図
【図3】本発明実施例の作用説明図
【図4】本発明実施例の減弱変換器2で用いる関数につ
いての説明を行うための図
【符号の説明】
1 荷重平均計算器 2 減弱変換器 3 積算器 pi,j 画素 Pi,j 画素値(透過線量)

Claims (1)

    (57)【特許請求の範囲】
  1. 【請求項1】 X線源から被写体に一様に照射され、そ
    の被写体を通過した透過X線像を、複数の画素が配列さ
    れた検出器で検出することにより得られるX線画像のデ
    ータ処理を行う装置であって、X線画像の一つの画素と
    その周辺の画素の各透過線量値を用い、その中央画素の
    重み係数を最大として荷重平均計算を行って透過線量の
    平均値を求める荷重平均計算器と、その平均値を変数と
    する関数で、当該平均値の小→大の変化に応じて減少す
    る関数を用いて減弱変換倍率を求める減弱変換器と、こ
    の変換倍率と上記中央画素の透過線量値との掛け算を行
    う積算器を備え、その演算結果を当該X線画像の1画素
    分の画像データとして用いるように構成されていること
    を特徴とするX線画像処理装置。
JP02956294A 1994-02-28 1994-02-28 X線画像処理装置 Expired - Lifetime JP3303504B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP02956294A JP3303504B2 (ja) 1994-02-28 1994-02-28 X線画像処理装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP02956294A JP3303504B2 (ja) 1994-02-28 1994-02-28 X線画像処理装置

Publications (2)

Publication Number Publication Date
JPH07236626A JPH07236626A (ja) 1995-09-12
JP3303504B2 true JP3303504B2 (ja) 2002-07-22

Family

ID=12279580

Family Applications (1)

Application Number Title Priority Date Filing Date
JP02956294A Expired - Lifetime JP3303504B2 (ja) 1994-02-28 1994-02-28 X線画像処理装置

Country Status (1)

Country Link
JP (1) JP3303504B2 (ja)

Also Published As

Publication number Publication date
JPH07236626A (ja) 1995-09-12

Similar Documents

Publication Publication Date Title
EP0527525B1 (en) Method and apparatus for contrast enhancement
EP0032237B1 (en) Method of and apparatus for processing a radiographic image
US5805721A (en) Method and apparatus for contrast enhancement
JP2575612B2 (ja) 画像のコントラスト劣化を修正する方法及びその装置
JP3188491B2 (ja) X線記録のダイナミック圧縮方法及びその装置
US9996910B2 (en) Radiographic image processing device, method, and recording medium
US5550888A (en) Method of displaying notably a digital X-ray image as a visible image, and device for carrying out the method
JPS6262383B2 (ja)
JPH10171983A (ja) 画像処理方法および装置
CN116018611A (zh) 使用深度卷积网络的噪声抑制
JP3540914B2 (ja) X線撮影装置
EP4190243A1 (en) Image processing device, image processing method, learning device, learning method, and program
JPS58201169A (ja) 画像デ−タ処理装置
CN111127371B (zh) 图像增强参数自动优化方法、存储介质和x射线扫描装置
WO2020241110A1 (ja) 画像処理装置、画像処理方法及びプログラム
US8111887B2 (en) Image processing method, image processing apparatus and image processing program
JP3303504B2 (ja) X線画像処理装置
JP2013123148A (ja) 画像処理装置及び画像処理方法
Angelhed et al. Optimization and evaluation of an image intensifier TV system for digital chest imaging
Kawahara et al. Spatial frequency components of normal radiographic anatomical features on intraoral computed radiography
JP3443948B2 (ja) 放射線撮像装置
Flynn Application of the ILO International Classification of Radiographs of Pneumoconioses to Digital Chest Radiographic Images
JPS6262380B2 (ja)
Fujita et al. A simple method for evaluating image quality of screen-film system using a high-performance digital camera
Saunders Jr et al. Impact of resolution and noise characteristics of digital radiographic detectors on the detectability of lung nodules

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080510

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090510

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100510

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100510

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110510

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110510

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120510

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130510

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130510

Year of fee payment: 11

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140510

Year of fee payment: 12

EXPY Cancellation because of completion of term