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
Links
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000001678 irradiating effect Effects 0.000 claims description 2
- 230000007423 decrease Effects 0.000 claims 1
- 238000000034 method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000003745 diagnosis Methods 0.000 description 3
- 210000000038 chest Anatomy 0.000 description 2
- 238000011976 chest X-ray Methods 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 206010002329 Aneurysm Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 210000001370 mediastinum Anatomy 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
置、X線量子計数器あるいは胸部デジタルX線画像診断
装置等のX線画像のデータ処理に利用される画像処理装
置に関する。
ミックレンジを持ち、その表示が困難である。そこで、
このような胸部X線画像を表示するため、CR(Compute
r Ra-diography)などでは、空間周波数フィルタ処理
(ぼけマスク処理など)が行われている。これらの処理
は、特定の空間周波数帯域(高周波数帯域)を強調した
デジタル処理で、線状の陰影を強調する画像処理(エッ
ジ強調処理)である。
ものとして、感度特性の異なる2種類の材料を組み合わ
せたX線フィルムも開発されている。
ジタル画像処理では、線状の陰影が強調されるので読影
はし易くなるが、しゅ瘤などの大きな陰影の読影は難し
く、また、従来のレントゲン写真に比べて画像品位が大
きく相違する画像(不自然な画像)となる等の欠点があ
る。そのため、この種のデジタル画像処理を採用したX
画像撮像装置(FCR;富士コンピュータラジオグラフ
ィ)などでは、医師らに診断時には必ず2枚の画像を読
影することを提示するするようにしている。しかし、こ
れは必ずしも良い処置とは言えない。
もので、被検体にX線を一様に照射して検出された透過
X線画像を表示するにあたり、自然で読影し易い画像を
得ることのできる画像処理装置の提供を目的とする。
め、本発明は、実施例に対応する図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画素分の画像データとし
て用いるように構成したことによって特徴づけられる。
ある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 を得る。
の積を積算器3で求め、このような計算をX線画像Pの
全ての画素について行うことにより画像全体の処理デー
タを得る。
Pのある画素の元の透過線量値が、目標検出線量;Cに
近い場合には、その値は処理後でも殆ど変化せずC付近
の値となるが、元の透過線量値がCに対して十分に低い
場合には、その処理後の値は積算器3での倍数gの掛け
算により元の値よりも大きくなる(図3参照)。従っ
て、積算器3は、空間的な低周波成分に対しては、あた
かもX線量を増やしたかのように機能することになり、
これによりダイナミックレンジの幅を所定範囲内に抑え
ることが可能になる。
の重み付け計算により、一つの画素の周辺画素に対する
減弱度が抽出されているので、小さな構造も読影するこ
とができるX線画像を得ることも可能となる。
明する。図1は本発明実施例の構成図である。この図で
は、処理の説明を簡単とするために、ある一つの画素p
i,j の処理結果Qi,j を得るときの処理を行う構成のみ
を示している。また、同様に処理説明を簡単とするた
め、この例では透過X線画像Sに対して縦方向の処理の
みを行うものとする。
計算器1、減弱変換器2及び積算器3によって構成され
ている。荷重平均計算器1は、X線画像Sの画素pi,j
を中央とする画素区間pi,j-L・・pi,j+L の各画素値
(透過線量)を用いて平均透過線量を重み付け計算によ
り求めるための計算器である。その計算式は
ための荷重係数で、この例では、次の(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として、被写体がないと
きの検出値、また肺野と縦隔との境界値あるいは肺野の
画素データの最高値等を用いる。
gi,j と、画素pi,j の元の透過線量値Pi,j とが掛け
算され、その演算結果;Qi,j =Pi,j *gi,j が画像
処理結果として出力される。
素について行うことにより、X線画像の全体の処理デー
タを得ることができ、このような処理において、荷重平
均計算器1は画像データ群に対して空間的なロウパスフ
ィルタとして作用することになり、これにより空間的な
低周波成分をある程度の範囲内で一定にしてダイナミッ
クレンジの幅を所定の範囲内に抑えることができる。
ば図3に示すように、X線画像の元の透過線量データP
が広いダイナミックレンジを持つ場合、このデータを、
本発明実施例の荷重平均計算器1、減弱変換器2及び積
算器3を用いて処理すると、その処理結果Qは、定数C
に近い透過線量ではあまり変化しないが、定数Cに比べ
少ない透過線量では変化が大きく、定数Cに近づく。そ
の結果、ダイナミックレンジの幅が、元のX線画像より
も狭くなる。しかも、この例での画像処理ではX線画像
の全体的な変化も、元の画像の状態が反映されるので、
自然で読影し易い画像を得ることができる。
Sとしてもよいが、この場合、低周波成分がほぼ一定と
なって小さな構造は明瞭になるものの、図3の模式的線
図のQ′に示すように、空間的な低周波成分がほば同一
のレベルとなってしまい、得られる画像がフィルムに比
べて不自然になるといった問題が発生する。そこで、自
然で読影し易い画像を得ることを目的とする場合には、
減弱変換器2の関数として、例えば、先に述べたg=√
(C/S)を用いることが好ましい。
な画像データを参照して説明する。まず、図に示すデー
タでは、元のX線画像より荷重平均計算器1の出力がS
=P(透過線量の計数値)、減弱変換器2により倍率g
=√(C/P)、積算器3の出力はQ=g・P=√(C
/P)・P=√(P)・√(C)となる。
慮すると計数Pの値はP±√(P)程ゆらいでいるが、
上記の積算を行った場合のゆらぎを求めると、 となり、目標計数値Cとしたときのゆらぎは、C±√
(C)であるから、g=√(C/P)で減弱変換を行っ
ても、ゆらぎは等しくなる。これにより、先の図3の模
式的線図に示したように画像の不自然さがなくなる。
する処理の方向を縦方向のみとしているが、その方向は
横方向であってもよいし、また、縦方向の処理を行った
後、横方向の処理を行うといった手法を採用してもよ
い。さらに、空間的に拡張した画像処理つまり一つの画
素についてその縦と横の双方の周辺画素の画素値を用い
たり、領域内の平均透過線量の計算を行うといった処理
手法を採用してもよい。
つの画素pi,j についての処理のみを示しているが、他
の画素についても画素pi,j と同等な構成により処理を
行うことは言うまでもない。
平均計算器の計算に適用する重み付けの係数としては、
先の(2) 式の関数のほか、ハミング窓関数あるいはハニ
ング窓関数が挙げられる。ただし、(2) 式の関数を採用
した場合、荷重平均計算器での演算速度が高速になり、
処理速度が速くなるといった利点がある。
線形の関数(図2参照)を採用すると、重み付け計算を
行う際の掛け算(Wi,j *Pi,J+JS)は、X線画像のあ
る画素区間についてのみ行い、この代表区間の計算結果
を基準として加減算を行うことにより画像の全ての画素
についての平均透過線量を求めるといった手法が可能と
なり、その掛け算を削減できる分だけ演算時間を大幅に
短縮できる。
換器に適用する関数としては、先に述べたg=C/Sや
g=1/S、さらに√(C/S)を一般化して(C/
S)Kあるいはシグモイド関数を用いた1+EXP(−k
1(S−k2))などが挙げられる。
処理装置によれば、X線画像の空間的な低周波成分をあ
る程度の範囲内で一定にして、ダイナミックレンジの幅
を所定の範囲内に抑えることが可能となるので、自然で
読影し易い画像を得ることができる。これにより、1枚
の画像で診断に必要な読影を十分に行うことが可能とな
り、その読影の能率が向上する。
いる関数を適宜に変更するだけで、ダイナミックレンジ
の特性を簡単かつ任意に変えることができる。これによ
り、例えば読影を行いたい部位の陰影を強調した画像を
得ることができる等、診断の目的に合った画像を作成す
ることも可能になる。
荷重係数を示す図
いての説明を行うための図
Claims (1)
- 【請求項1】 X線源から被写体に一様に照射され、そ
の被写体を通過した透過X線像を、複数の画素が配列さ
れた検出器で検出することにより得られるX線画像のデ
ータ処理を行う装置であって、X線画像の一つの画素と
その周辺の画素の各透過線量値を用い、その中央画素の
重み係数を最大として荷重平均計算を行って透過線量の
平均値を求める荷重平均計算器と、その平均値を変数と
する関数で、当該平均値の小→大の変化に応じて減少す
る関数を用いて減弱変換倍率を求める減弱変換器と、こ
の変換倍率と上記中央画素の透過線量値との掛け算を行
う積算器を備え、その演算結果を当該X線画像の1画素
分の画像データとして用いるように構成されていること
を特徴とするX線画像処理装置。
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) |
-
1994
- 1994-02-28 JP JP02956294A patent/JP3303504B2/ja not_active Expired - Lifetime
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 |