Nothing Special   »   [go: up one dir, main page]

JP3039879B2 - 放射線画像の周波数処理方法 - Google Patents

放射線画像の周波数処理方法

Info

Publication number
JP3039879B2
JP3039879B2 JP3085583A JP8558391A JP3039879B2 JP 3039879 B2 JP3039879 B2 JP 3039879B2 JP 3085583 A JP3085583 A JP 3085583A JP 8558391 A JP8558391 A JP 8558391A JP 3039879 B2 JP3039879 B2 JP 3039879B2
Authority
JP
Japan
Prior art keywords
image
frequency processing
processing method
frequency
radiation image
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
JP3085583A
Other languages
English (en)
Other versions
JPH04318540A (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.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
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 Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP3085583A priority Critical patent/JP3039879B2/ja
Priority to US07/866,478 priority patent/US5272626A/en
Publication of JPH04318540A publication Critical patent/JPH04318540A/ja
Application granted granted Critical
Publication of JP3039879B2 publication Critical patent/JP3039879B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • G06T5/75Unsharp masking

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Conversion Of X-Rays Into Visible Images (AREA)
  • Radiography Using Non-Light Waves (AREA)
  • Closed-Circuit Television Systems (AREA)
  • Image Analysis (AREA)

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、医療用診断に用いる放
射線画像の周波数処理方法、さらに詳しくは放射線画像
の診断性能を向上させるための周波数処理において画像
周辺部におけるアーティファクトを低減させる方法に関
するものである。
【0002】
【従来の技術】ここで周波数処理とは、放射線画像の診
断性能を向上させるために行なわれるもので、例えば特
開昭55−163472号、特開昭55−87953号
等に開示されたボケマスク処理等の周波数強調処理を意
味するものである。すなわち、読取出力信号Dorg
に、該読取出力信号Dorgからボケマスク信号Dus
を引いたものに強調度βを掛けたものを加える処理を施
すもので、これにより、画像において所定の空間周波数
成分を強調することができる。これを式で表すと下記の
式のようになる。◆ D=Dorg+β(Dorg−Dus) (D:周波数処理された信号、Dorg:読取出力信
号、◆ Dus:ボケマスク信号、β:強調度)
【0003】
【発明が解決しようとする課題】上記のような周波数処
理は、放射線画像の診断性能を向上させるために行なわ
れるものであるが、この周波数処理には次のような問題
がある。すなわち読取光の2次元走査によって得た放射
線画像情報には、図1に示す原画像信号に対して図2お
よび図3に示すように、読み取られる画像シート1等の
画像2の端部に近い周辺部においてアーティファクト3
が目立ち、これをそのままにすると診断の妨げとなると
いう難点がある。
【0004】周波数処理における強調度βは、撮影部
位、撮影方法(造影、断層等)等に応じて自動的に設定
されたり、外部入力によりマニュアルで設定されるもの
であるが、単にそれだけによって設定したのでは上記の
アーティファクトに対処できないため、この点の改良が
望まれるところである。
【0005】本発明は、上記問題に鑑み、画像の周辺部
におけるアーティファクトが目立たないように周波数処
理をするようにした放射線画像の周波数処理方法を提供
することを目的とするものである。
【0006】
【課題を解決するための手段】前記課題を解決するた
め、本発明による放射線画像の周波数処理方法は、放射
線画像を光電的に読み取って得た画像信号を所定の強調
度βに従って周波数処理する放射線画像情報の周波数処
理方法において、前記強調度βを前記放射線画像の端部
において所定幅に亘って画像中心側から外方に向って単
調減少させることを特徴とするものである。
【0007】なお、本発明による放射線画像の周波数処
理方法は、前記放射線画像の形成を、被写体を透過した
放射線を蓄積性螢光体シートに蓄積記録せしめることに
より行ない、前記光電的読取りを前記蓄積記録された放
射線画像に励起光を照射することにより発光される輝尽
発光光を光電的に読み取ることによって行なう、蓄積性
螢光体シ―トを使用した放射線画像の記録読取方法に適
用されるものである。
【0008】
【作用および効果】本発明による周波数処理方法は、前
記周波数処理において、強調度βを放射線画像の端部に
おいて所定幅すなわちアーティファクトの生じやすい幅
に亘って画像中心側から外方に向って単調減少させるこ
とによって、アーティファクトの発生を抑制することが
できる。
【0009】
【実施例】以下、本発明による周波数処理方法の実施例
を図面に基づいて説明する。
【0010】図4は、本発明による周波数処理方法の一
実施例を説明する、画像の大きさを画素単位の数で示す
図である。
【0011】周波数処理とは、放射線画像の診断性能を
向上させるために行なわれるもので、例えば特開昭55
−163472号に開示されたボケマスク処理等の周波
数強調処理においては、下記の式により所定の空間周波
数成分を強調する。◆ D=Dorg+β(Dorg)(Dorg−Dus)◆ ここでDusは、画像中、縦N、横Nの画素に対応する
大きさの正方形の中の画素の信号Dorgの平均、すな
わちΣDorg/N2 である。また、強調度βは画像信
号の値(Dorg)によって決まるものである。
【0012】本発明による周波数処理方法の実施例で
は、上記式は◆ D=Dorg+α×β(Dorg)(Dorg−Du
s)◆ となり、βには別のファクターαが掛けられ、強調度は
実質的にαに応じて変化するものとなる。その変化は、
前記アーティファクト3を生じやすい画像2の周辺部に
おいて、所定幅に亘って画像中心側から外方に向ってβ
を単調減少させることである。
【0013】すなわち具体的には、図4に示すように画
像のサイズ(p×q)と周辺部の幅(m)を設定したと
き、横方向の画素の位置(i)と縦方向の画素の位置
(j)に応じてαを画像中心側から外方に向って単調減
少させるように変化させるものである。その例として
は、以下に示す(1)から(4)のうち少なくともいず
れか1つと(5)との組み合せが考えられる。
【0014】 α=i/m (i<m) (1) α=(p−i)/m (i>p−m) (2) α=j/m (j<m) (3) α=(q−j)/m (j>q−m) (4) α=1 (otherwise :その他) (5) すなわち、αを、前記アーティファクトを生じやすい幅
mの周辺部において、iまたはjの値、すなわちその画
素の画像端部からの距離に応じて、1/mからm/m=
1まで単調に増加するように、換言すれば画像中心側か
ら外方に向って単調減少するように変化させるものであ
る。これにより、前記アーティファクトを抑制すること
ができる。
【0015】なお、上記実施例におけるαの式の(1)
から(4)を二乗して下記のようにしてもよい。
【0016】 α=(i/m)2 (i<m) (1′) α={(p−i)/m}2 (i>p−m) (2′) α=(j/m)2 (j<m) (3′) α={(q−j)/m}2 (j>q−m) (4′) また、上記実施例では、アーティファクトの発生する周
辺部の幅を縦も横もmとしたが、これは走査方向によっ
て変えてもよい。
【0017】さらに、上記実施例におけるαをαpとα
q の積とし、それぞれ下記のようにしてもよい。
【0018】 αp=i/m (i<m) (1″) αp=(p−i)/m (i>p−m) (2″) αp=1 (otherwise :その他) (5′) αq =j/m (j<m) (3″) αq =(q−j)/m (j>q−m) (4″) αq =1 (otherwise :その他) (5″) なお、上記実施例では、βは原画像信号Dorgに応じ
て変化するようにしているため、特公昭62-62373号に開
示されているようなβを可変にした場合の効果をも併わ
せ奏するものであるが、これは変化させず一定としても
よい。βを原画像信号Dorgに応じて変化させる方法
としては、βを原画像信号の大きさに応じて単調に増加
させる方法や、原画像信号の大きさが小さいときと大き
いときは単調に増加させ、原画像信号の大きさが中程度
のときは一定とする方法がある。また放射線画像におい
て差信号の小さい部分では非鮮鋭マスク処理による周波
数強調を普通に行ない、差信号の大きい部分(たとえば
骨と筋肉の境界、軟部とガス部の境界、胃のBa 充填部
とその周囲との境界血管造影における血管陰影等)では
周波数強調の程度の増加を押えて偽画像の発生を防止す
る方法を採用してもよい。
【0019】なお、照射線量が高くノイズの小さいとき
や、複数画像を重ねてノイズを小さくしたときは、周波
数処理の強調度を大きくして診断性能をさらに向上させ
ることができる。
【図面の簡単な説明】
【図1】本発明の対象とする画像信号の原画像信号の状
態を示すグラフ
【図2】本発明の対象とするアーティファクトを有する
画像信号を示すグラフ
【図3】本発明の対象とするアーティファクトを有する
画像の一部を示す平面図
【図4】本発明による周波数処理方法の一実施例を説明
する、画像の大きさを画素単位の数で示す図
【符号の説明】
1 画像シート 2 画像 3 アーティファクト
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭56−11038(JP,A) 特開 昭61−115540(JP,A) 特開 平1−238638(JP,A)

Claims (1)

    (57)【特許請求の範囲】
  1. 【請求項1】 被写体を透過した放射線を蓄積性螢光体
    シートに蓄積記録せしめることにより放射線画像を形成
    し、 この蓄積記録された放射線画像に励起光を照射すること
    により発光される輝尽発光光を光電的に読み取って画像
    信号を得る放射線画像の記録読取方法において、 前記画像信号を所定の強調度βに従って周波数処理する
    放射線画像情報の周波数処理方法であって、 前記強調度βを前記放射線画像の端部において所定幅に
    亘って画像中心側から外方に向って単調減少させること
    を特徴とする放射線画像の周波数処理方法。
JP3085583A 1991-04-17 1991-04-17 放射線画像の周波数処理方法 Expired - Lifetime JP3039879B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP3085583A JP3039879B2 (ja) 1991-04-17 1991-04-17 放射線画像の周波数処理方法
US07/866,478 US5272626A (en) 1991-04-17 1992-04-10 Frequency processing method for radiation images

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP3085583A JP3039879B2 (ja) 1991-04-17 1991-04-17 放射線画像の周波数処理方法

Publications (2)

Publication Number Publication Date
JPH04318540A JPH04318540A (ja) 1992-11-10
JP3039879B2 true JP3039879B2 (ja) 2000-05-08

Family

ID=13862836

Family Applications (1)

Application Number Title Priority Date Filing Date
JP3085583A Expired - Lifetime JP3039879B2 (ja) 1991-04-17 1991-04-17 放射線画像の周波数処理方法

Country Status (2)

Country Link
US (1) US5272626A (ja)
JP (1) JP3039879B2 (ja)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2445536B1 (fr) * 1978-12-26 1989-12-22 Fuji Photo Film Co Ltd Procede et appareil pour traiter une image radiographique
JPS5611038A (en) * 1979-07-11 1981-02-04 Fuji Photo Film Co Ltd Method and device for treating radiation picture
US4315318A (en) * 1978-12-26 1982-02-09 Fuji Photo Film Co., Ltd. Method and apparatus for processing a radiation image
JPS61115540A (ja) * 1984-11-13 1986-06-03 株式会社日立製作所 フイルム処理方式
US4747052A (en) * 1984-11-14 1988-05-24 Philips Medical Systems, Inc. Radiation image processing
CA1316591C (en) * 1987-10-20 1993-04-20 Kazuhiro Hishinuma Method and apparatus for radiation image processing and x-ray image processing
JPH0789199B2 (ja) * 1988-03-19 1995-09-27 富士写真フイルム株式会社 放射線照射野輪郭候補点検出方法
JPH07101287B2 (ja) * 1988-04-20 1995-11-01 富士写真フイルム株式会社 放射線照射野認識方法

Also Published As

Publication number Publication date
US5272626A (en) 1993-12-21
JPH04318540A (ja) 1992-11-10

Similar Documents

Publication Publication Date Title
JP3002937B2 (ja) エネルギーサブトラクション画像処理方法
EP1909227B1 (en) Method of and apparatus for minimizing ringing artifacts in an input image
US7912316B2 (en) Method, apparatus and program for noise suppression
JP3196864B2 (ja) 画像のダイナミックレンジ圧縮処理方法
JPH0614168B2 (ja) 放射線画像の周波数処理方法および装置
JPH04156689A (ja) 放射線画像処理方法および装置
JP3039879B2 (ja) 放射線画像の周波数処理方法
JPH04318775A (ja) エネルギーサブトラクション画像生成方法
JP3783116B2 (ja) 放射線画像強調処理方法および装置
JPH08294006A (ja) 画像のダイナミックレンジ圧縮処理方法および装置
JP2002125153A (ja) ノイズ抑制処理装置並びに記録媒体
JP4083251B2 (ja) 照射野外画像処理方法および装置
JPH01237887A (ja) 放射線画像の粒状雑音低減方法
JP4323708B2 (ja) 画像処理方法、画像処理装置および記録媒体
JP3937616B2 (ja) 画像処理装置、方法及びコンピュータ読み取り可能な記憶媒体
JP2000232611A (ja) エネルギーサブトラクション画像生成方法および生成装置
JP3399684B2 (ja) 画像処理方法および装置
JPH0572624B2 (ja)
JP2739372B2 (ja) 画像平滑化処理方法および装置
JP2000235079A (ja) 放射線画像撮像装置
JP2867096B2 (ja) 画像平滑化処理方法
JPH0481982A (ja) デジタル放射線画像信号の処理装置
JPH0743765B2 (ja) 放射線画像処理方法および装置
JP2542183B2 (ja) X線画像処理方法および装置
JPS61119250A (ja) 放射線画像処理方法および装置

Legal Events

Date Code Title Description
A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 19981020

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

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

Free format text: PAYMENT UNTIL: 20080303

Year of fee payment: 8

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

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

Free format text: PAYMENT UNTIL: 20080303

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20090303

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20090303

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20100303

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20100303

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20110303

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20110303

Year of fee payment: 11

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

Free format text: PAYMENT UNTIL: 20120303

Year of fee payment: 12

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120303

Year of fee payment: 12