JPH09270004A - X線画像デジタル処理装置 - Google Patents
X線画像デジタル処理装置Info
- Publication number
- JPH09270004A JPH09270004A JP8101887A JP10188796A JPH09270004A JP H09270004 A JPH09270004 A JP H09270004A JP 8101887 A JP8101887 A JP 8101887A JP 10188796 A JP10188796 A JP 10188796A JP H09270004 A JPH09270004 A JP H09270004A
- Authority
- JP
- Japan
- Prior art keywords
- frequency component
- image
- average luminance
- low frequency
- influence
- 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
Landscapes
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Closed-Circuit Television Systems (AREA)
Abstract
ントラスト・シャープさを再現し、より見易い画像を得
る。 【解決手段】 元の画像データから低周波成分を抽出す
る回路21と、着目画素の小近傍および周辺部の平均輝
度をそれぞれ算出する回路22、23と、この小近傍お
よび周辺部の平均輝度の関係から散乱線影響度Kを求め
る回路24と、Kを上記の低周波成分に乗じる乗算器2
5と、その乗算結果を元の画像データから差し引く減算
器26とからなる。
Description
て用いられるX線画像信号をデジタル画像データに変換
して処理する装置に関する。
の身体など)のX線透視像を電気的な画像信号として撮
像し、これをモニター装置に表示して観察することが広
く行なわれている。さらに、この画像信号をデジタル画
像データに変換して種々のデジタル画像処理を行なうこ
とも普及している。
にその物質内部で散乱するため、その散乱線の影響によ
り、X線透視像のコントラストが悪くなったり、あるい
はシャープさ(尖鋭度)が鈍くなったりするという、画
質低下が生じる。そのため、現状では、X線の線質を変
えて散乱しにくくしたり、あるいはX線コリメータの形
状を工夫して、コリメータでの散乱線を抑制することが
試みられている。
ようにX線の線質の調整やコリメータの工夫によって
も、被写体内の散乱線の影響による画質低下は本質的に
除去することは不可能である。
たX線透視像を処理することにより、散乱線の影響を除
去するように改善した、X線画像デジタル処理装置を提
供することを目的とする。
め、この発明によるX線画像デジタル処理装置において
は、入力されたアナログX線画像信号をA/D変換する
A/D変換手段と、デジタル化された画像信号から低周
波成分を抽出する手段と、デジタル化された画像信号か
ら、着目する画素の小近傍およびその周辺部における平
均輝度をそれぞれ求める手段と、該小近傍の平均輝度と
周辺部の平均輝度との関係を求める手段と、この関係に
応じて上記の低周波成分を調整した上で元のデジタル化
された画像信号から差し引く手段とを有することが特徴
となっている。
部分から低い部分へとにじみだすようなものとなる。そ
こで、画像信号から低周波成分を抽出すれば、散乱線に
よる影響分を取り出すことができる。ところが、この低
周波成分には真の画像の低周波成分も含まれている。そ
こで、この取り出された低周波成分を、元の画像信号か
ら差し引くだけでは、散乱線の影響を除くことについて
は不十分である。上記のように、散乱線の影響は、画像
の輝度の高い部分から低い部分へとにじみだすようなも
のとなっているため、小領域の輝度が周辺に対して低い
ものであれば、その小領域に含まれる画素については散
乱線によって輝度が高くなり、逆に小領域の輝度が周辺
に対して高いものであれば、その小領域に含まれる画素
については散乱によって輝度が低くなるものと考えられ
る。そこで、この小領域の平均輝度と周辺部の平均輝度
との関係を求めて、これに応じて上記の低周波成分を調
整すれば、散乱線の影響をより正確に反映した信号を得
ることができ、これを元の画像信号から差し引くことに
より、散乱線の影響分を適切に除去することが可能とな
る。
ついて図面を参照しながら詳細に説明する。図1におい
て、X線管11から発射されたX線が被写体10を透過
してイメージインテンシファイア12に入射し、X線透
過像が光学像に変換される。イメージインテンシファイ
ア12には光学系13を介してTVカメラ14が結合さ
れており、イメージインテンシファイア12の出力光学
像の画像信号が得られる。この画像信号はA/D変換器
15によりデジタル画像データに変換された後、デジタ
ル画像処理装置16に送られて種々のデジタル画像処理
を受け、その後D/A変換器17でアナログの画像信号
に戻され、TVモニター装置18に送られる。
抽出回路21と、小近傍平均輝度算出回路22と、周辺
部平均輝度算出回路23と、散乱線影響度算出回路24
と、乗算器25と、減算器26とからなる散乱線除去部
20と、他の画像処理を行う画像処理部27とから構成
されている。
回路21に送られ、画像の低周波成分ISが抽出され
る。この低周波成分抽出回路21は、たとえば図2の
(a)で示すようなNF×NFのマトリクスのテンプレ
ートを用いて畳み込み演算を行なうローパス空間フィル
タにより構成される。また、デジタル画像データは小近
傍平均輝度算出回路22および周辺部平均輝度算出回路
23にも送られる。
周辺部平均輝度算出回路23は、いずれもたとえばテン
プレートによる畳み込み演算で構成された空間フィルタ
である。小近傍平均輝度算出回路22は着目する画素の
周囲の小さい領域の平均輝度MSを求めるものであり、
周辺部平均輝度算出回路23はその周辺の比較的広い領
域での平均輝度MLを求めるものである。そのため、前
者で用いるテンプレートは図2の(b)に示すようにN
S×NSの小さいマトリクスで構成され、後者で用いる
テンプレートはそれを囲むような比較的大きなNL×N
Lで構成される。これらNS×NSおよびNL×NLの
テンプレートはすべて「1」の重みを正規化して用いて
いる。NS×NSのマトリクスは着目画素のごく近傍の
平均輝度を求めるため小さく設定し、NL×NLのマト
リクスは比較的広い範囲の平均輝度を求めるために大き
く設定することが望ましい。
たMSとMLとから各画素ごとに散乱線影響度Kを算出
する。この例では、図3に示すようなグラフを用いてM
S−MLの値からKを求める。そして、乗算器25にお
いて低周波成分ISに対して影響度Kを乗じ(IS×
K)、これを補正値ICとして減算器26に送り、元の
デジタル画像データから補正値ICを減算する。これに
より散乱線の影響が除去された画像データは画像処理部
27で、ウインドウ変換、エッジ強調などの各種のデジ
タル画像処理を受けた後、D/A変換器17に送られ、
アナログの画像信号に戻されてTVモニター装置18に
送られて表示される。
ロフィールが図4の(a)のようになる場合を想定す
る。輝度が急峻に落ち込んでいる狭い部分は、たとえば
血管造影の場合の血管部分に相当する。ところが散乱線
の影響があるということは、輝度の高い部分の周囲に輝
度の低い部分があると、その輝度の高い部分から低い部
分へとX線がにじみ出すようなものであるから、その輝
度の高い部分における輝度は低くなり、逆に輝度の低い
部分の輝度は高くなる。そこで、この散乱線の影響によ
り、実際に得られる画像データの1ライン分のプロフィ
ールは、図4の(b)のようになる。この画像データか
ら低周波成分ISを取り出すと図4の(c)のようにな
るが、これには真のデータの低周波成分と散乱の影響分
とが含まれている。
とは、周辺部より小近傍の画素値が大きいということで
あり、そのような場合は散乱によって画素値が真の値よ
りも低くなっていると考えられる。逆にMS−MLの値
が負であるということは、周辺部より小近傍の画素値が
小さいということであり、そのような場合は散乱によっ
て画素値が真の値よりも高くなっていると考えられる。
そのため、MS−MLの値が正のときは、画素値に対し
て散乱の影響がマイナス側に働き、逆にMS−MLの値
が負のときは、画素値に対して散乱の影響がプラス側に
働く。そこで、図3に示すようなカーブ(たとえば太
線)によりMS−MLの値に応じてKを求め、これをI
Sに乗じた上で元の画像データから差し引けばよいこと
が分かる。この図3においてMS−MLの値が0付近で
K値が0としたのは、ノイズによる影響を避けるための
不感帯を設けるためである。この図3のカーブは、X線
条件(管電圧、管電流、イメージインテンシファイアの
種類・視野等)に応じて図3の細線や点線のように変え
ることもできる。最適なカーブは、実験やコンピュータ
シュミレーションなどで求める。また、ここでは、MS
とMLとの差を用いたが、MSとMLとの比を用いて同
様のカーブからKを求めるようにしてもよい。
Cを図4の(d)に示すようになり、このICを元の画
像データから差し引くことにより、図4の(e)に示す
ような画像データを得ることができる。補正値ICは、
上記の通り、散乱による影響をより反映したものとなっ
ているので、散乱の影響を除いたデータを得ることが可
能となる。図4の(b)で示した元の画像データと図4
の(e)で示した補正後のデータとの比較からも分かる
通り、血管部などの急峻な細い落ち込み部のコントラス
ト・シャープさが再現され、血管部などがより見易くな
る。
いるが、小領域ごとにKを算出し、補間によって各画素
ごとのKを求めるようにしてもよい。また、小近傍平均
輝度MSは、上記では低周波成分の抽出処理とは別個の
処理によって求めているが、小近傍平均輝度MSは一種
のローパスフィルタ出力であるから、低周波成分抽出回
路21で得た値ISを兼用してもよい。さらに、上記で
は、低周波成分抽出回路21、小近傍平均輝度算出回路
22、周辺部平均輝度算出回路23、散乱線影響度算出
回路24等をハードウェアで構成しているかのように説
明したが、ソフトウェアによる処理で求めることももち
ろん可能である。
像デジタル処理装置によれば、散乱線の影響を適切に除
去して、元来の被写体のコントラスト・シャープさを再
現し、より見易い、優れた画質の画像を得ることができ
る。また、量子ノイズなどの高周波ノイズ成分が強調さ
れることなく、S/N比の改善が可能である。
Claims (1)
- 【請求項1】 入力されたアナログX線画像信号をA/
D変換するA/D変換手段と、デジタル化された画像信
号から低周波成分を抽出する手段と、デジタル化された
画像信号から、着目する画素の小近傍およびその周辺部
における平均輝度をそれぞれ求める手段と、該小近傍の
平均輝度と周辺部の平均輝度との関係を求める手段と、
この関係に応じて上記の低周波成分を調整した上で元の
デジタル化された画像信号から差し引く手段とを備える
ことを特徴とするX線画像デジタル処理装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10188796A JP3721632B2 (ja) | 1996-03-31 | 1996-03-31 | X線画像デジタル処理装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP10188796A JP3721632B2 (ja) | 1996-03-31 | 1996-03-31 | X線画像デジタル処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH09270004A true JPH09270004A (ja) | 1997-10-14 |
| JP3721632B2 JP3721632B2 (ja) | 2005-11-30 |
Family
ID=14312452
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP10188796A Expired - Fee Related JP3721632B2 (ja) | 1996-03-31 | 1996-03-31 | X線画像デジタル処理装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3721632B2 (ja) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005005981A (ja) * | 2003-06-11 | 2005-01-06 | Shimadzu Corp | 画像処理装置 |
| CN100411445C (zh) * | 2003-07-18 | 2008-08-13 | 佳能株式会社 | 校正图像亮度分布的图像处理方法及装置 |
| CN105011955A (zh) * | 2014-04-30 | 2015-11-04 | 株式会社东芝 | X射线诊断装置及图像处理装置 |
| US9648212B2 (en) | 2015-02-24 | 2017-05-09 | Samsung Display Co., Ltd. | Image processing device, image processing method, and computer program |
| US9866761B2 (en) | 2015-02-24 | 2018-01-09 | Samsung Display Co., Ltd. | Image dynamic range compensation device and method |
| US10786219B2 (en) | 2018-01-24 | 2020-09-29 | Konica Minolta, Inc. | Radiographic image processing apparatus, scattered radiation correction method, and computer readable storage medium |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105427355B (zh) * | 2015-11-26 | 2018-04-03 | 上海联影医疗科技有限公司 | 罩壳、x射线图像的散射成分计算、重建的方法及装置 |
-
1996
- 1996-03-31 JP JP10188796A patent/JP3721632B2/ja not_active Expired - Fee Related
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005005981A (ja) * | 2003-06-11 | 2005-01-06 | Shimadzu Corp | 画像処理装置 |
| CN100411445C (zh) * | 2003-07-18 | 2008-08-13 | 佳能株式会社 | 校正图像亮度分布的图像处理方法及装置 |
| CN105011955A (zh) * | 2014-04-30 | 2015-11-04 | 株式会社东芝 | X射线诊断装置及图像处理装置 |
| CN105011955B (zh) * | 2014-04-30 | 2019-03-01 | 东芝医疗系统株式会社 | X射线诊断装置及图像处理装置 |
| US9648212B2 (en) | 2015-02-24 | 2017-05-09 | Samsung Display Co., Ltd. | Image processing device, image processing method, and computer program |
| US9866761B2 (en) | 2015-02-24 | 2018-01-09 | Samsung Display Co., Ltd. | Image dynamic range compensation device and method |
| US10786219B2 (en) | 2018-01-24 | 2020-09-29 | Konica Minolta, Inc. | Radiographic image processing apparatus, scattered radiation correction method, and computer readable storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3721632B2 (ja) | 2005-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3976337B2 (ja) | ノイズ低減用画像処理 | |
| US6766064B1 (en) | Method and apparatus for performing a contrast based dynamic range management algorithm | |
| JP3690882B2 (ja) | 画像の強調処理方法および装置 | |
| US4618990A (en) | Edge enhancement filtering for digital fluorography images | |
| EP0866607A2 (en) | Noise filter for digital x-ray imaging system | |
| US6760401B2 (en) | Apparatus and method for processing of digital images | |
| US9311695B2 (en) | Image processing method and radiographic apparatus using the same | |
| JP2003244542A (ja) | 2度露光二重エネルギーx線撮影のためのモーション・アーチファクト減少アルゴリズム | |
| US9922409B2 (en) | Edge emphasis in processing images based on radiation images | |
| US7149358B2 (en) | Method and system for improving contrast using multi-resolution contrast based dynamic range management | |
| CN103561654A (zh) | 医用图像处理装置、x 射线诊断装置、医用图像处理方法及x 射线诊断方法 | |
| JP3721632B2 (ja) | X線画像デジタル処理装置 | |
| US8374417B2 (en) | Image processing method and radiographic apparatus using the same | |
| US20130308841A1 (en) | Method and apparatus for image processing | |
| JPS6340533A (ja) | X線診断装置 | |
| US20110305382A1 (en) | Image processing method and radiographic apparatus using the same | |
| JPH10105701A (ja) | 放射線画像強調処理方法および装置 | |
| US6956977B2 (en) | Methods for improving contrast based dynamic range management | |
| JP2009054013A (ja) | 画像処理装置 | |
| JPH10229519A (ja) | X線画像デジタル処理装置 | |
| JP2001178686A (ja) | 画像表示方法および画像生成装置 | |
| JP3629875B2 (ja) | X線画像デジタル処理装置 | |
| JP4006083B2 (ja) | 画像処理装置、方法及びコンピュータ読み取り可能な記憶媒体 | |
| US20200305829A1 (en) | Radiographic image processing device, method of operating radiographic image processing device, and radiographic image processing program | |
| JPH0991423A (ja) | 画像処理方法および装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20041108 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20041124 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050124 |
|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20050331 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20050331 |
|
| 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: 20050823 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050905 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080922 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090922 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090922 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100922 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110922 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110922 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120922 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120922 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130922 Year of fee payment: 8 |
|
| LAPS | Cancellation because of no payment of annual fees |