JP6006454B2 - 放射線画像処理装置および方法並びにプログラム - Google Patents
放射線画像処理装置および方法並びにプログラム Download PDFInfo
- Publication number
- JP6006454B2 JP6006454B2 JP2016132257A JP2016132257A JP6006454B2 JP 6006454 B2 JP6006454 B2 JP 6006454B2 JP 2016132257 A JP2016132257 A JP 2016132257A JP 2016132257 A JP2016132257 A JP 2016132257A JP 6006454 B2 JP6006454 B2 JP 6006454B2
- Authority
- JP
- Japan
- Prior art keywords
- grid
- radiation
- scattered
- image
- scattered radiation
- 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.)
- Active
Links
- 230000005855 radiation Effects 0.000 title claims description 380
- 238000000034 method Methods 0.000 title claims description 65
- 238000003384 imaging method Methods 0.000 claims description 74
- 238000002834 transmittance Methods 0.000 claims description 66
- 238000003672 processing method Methods 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 2
- 230000035699 permeability Effects 0.000 claims 2
- 239000012466 permeate Substances 0.000 claims 1
- 238000006243 chemical reaction Methods 0.000 description 31
- 230000000694 effects Effects 0.000 description 8
- 210000000038 chest Anatomy 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 6
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 3
- 229910052782 aluminium Inorganic materials 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000001678 irradiating effect Effects 0.000 description 3
- 210000004072 lung Anatomy 0.000 description 3
- 210000001370 mediastinum Anatomy 0.000 description 3
- 238000003786 synthesis reaction Methods 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 2
- 210000001015 abdomen Anatomy 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 210000003128 head Anatomy 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 229920001342 Bakelite® Polymers 0.000 description 1
- 102100035353 Cyclin-dependent kinase 2-associated protein 1 Human genes 0.000 description 1
- 101000585359 Homo sapiens Suppressor of tumorigenicity 20 protein Proteins 0.000 description 1
- 102100029860 Suppressor of tumorigenicity 20 protein Human genes 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000004637 bakelite Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000035876 healing Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
放射線画像に対して使用したグリッドに起因する縞模様を除去する処理を行って、グリッド使用放射線画像を取得する手段と、
使用したグリッドとは異なる仮想的なグリッドの特性である仮想グリッド特性を取得する特性取得手段と、
仮想グリッド特性に基づいて、グリッド使用放射線画像の散乱線除去処理を行う散乱線除去手段とを備えたことを特徴とするものである。
放射線画像の撮影時に、散乱線を除去するために使用が想定される、仮想的なグリッドの特性である仮想グリッド特性を取得する特性取得手段と、
仮想グリッド特性に基づいて、放射線画像の散乱線除去処理を行う散乱線除去手段とを備えたことを特徴とするものである。
散乱線除去手段を、散乱成分情報にも基づいて散乱線除去処理を行う手段としてもよい。
放射線画像に対して使用したグリッドに起因する縞模様を除去する処理を行って、グリッド使用放射線画像を取得し、
使用したグリッドとは異なる仮想的なグリッドの特性である仮想グリッド特性を取得し、
仮想グリッド特性に基づいて、グリッド使用放射線画像の散乱線除去処理を行うことを特徴とするものである。
放射線画像の撮影時に、散乱線を除去するために使用が想定される、仮想的なグリッドの特性である仮想グリッド特性を取得し、
仮想グリッド特性に基づいて、放射線画像の散乱線除去処理を行うことを特徴とするものである。
Is(x,y) = Io(x,y)*Sσ(T(x,y)) …(2)
S(x,y) = Is(x,y)/(Is(x,y)+Ip(x,y)) …(3)
ここで、(x,y)は放射線画像の画素位置の座標、Ip(x,y)は画素位置(x,y)における一次線像、Is(x,y)は画素位置(x,y)における散乱線像、Io(x,y)は画素位置(x,y)における被写体表面への入射線量、μは被写体の線減弱係数、Sσ(T(x,y))は画素位置(x,y)における被写体厚に応じた散乱の特性を表す畳みこみカーネルである。式(1)は公知の指数減弱則に基づく式であり、式(2)は「J M Boon et al, An analytical model of the scattered radiation distribution in diagnostic radiolog, Med. Phys. 15(5), Sep/Oct 1988」(参考文献1)に記載された手法に基づく式である。なお、被写体表面への入射線量Io(x,y)は、どのような値を定義してもS(x,y)を算出する際に除算によってキャンセルされるため、例えば値を1とする等、任意の値とすればよい。
散乱線透過率Tsおよび一次線透過率Tp、並びに散乱線含有率分布S(x,y)は0〜1の間の値となるため、変換係数R(x,y)も0〜1の間の値となる。散乱線除去部36は、変換係数R(x,y)を複数の周波数帯域のそれぞれについて算出する。
=Σr{I(x,y,r)×(S(x,y)×Ts(r)+(1-S(x,y))×Tp(r))} …(5)
ここで、本実施形態においては、放射線画像を6つの周波数帯域に周波数分解するものとし、散乱線透過率Tsおよび一次線透過率Tpは6つの周波数帯域について取得されるものとする。この場合、散乱線透過率Tsおよび一次線透過率Tpは、例えば下記式(6)に示す値となる。なお、式(6)では右側ほど低周波数帯域の値を表すものとする。
Tp={0.7, 0.7, 0.7, 0.7, 0.7, 0.7} …(6)
式(6)に示すように、散乱線透過率Tsおよび一次線透過率Tpは、高周波数帯域(r=1〜4)では同一の値であるが、低周波数帯域(r=5〜6)においては、散乱線透過率Tsの方が低い値となる。これはグリッドは散乱線の周波数成分が支配的である低周波帯域ほどその除去率が高いが、一次線については除去率の周波数依存性が小さいからである。
Ip(x,y,r)=(1-S(x,y))×I(x,y,r) …(7)
さらに散乱線除去部36は、下記の式(8)により、散乱線成分Is(x,y,r)および一次線成分Ip(x,y,r)のそれぞれに対して、仮想グリッド特性である散乱線透過率Ts(r)および一次線透過率Tp(r)を適用して画像変換し、変換された散乱線成分Is′(x,y,r)および一次線成分Ip′(x,y,r)を算出する。
Ip′(x,y,r)=Ip(x,y,r)×Tp(r)=(1-S(x,y))×I(x,y,r)×Tp(r) …(8)
そして下記の式(9)により、Is′(x,y,r)および一次線成分Ip′(x,y,r)を周波数合成して、処理済みの放射線画像I(x,y)′を算出する。
=Σr{S(x,y)×I(x,y,r)×Ts(r)+(1-S(x,y))×I(x,y,r)×Tp(r)}
=Σr{I(x,y,r)×(S(x,y)×Ts(r)+(1-S(x,y))×Tp(r))} …(9)
次いで、第1の実施形態において行われる処理について説明する。図4は第1の実施形態において行われる処理を示すフローチャートである。撮影装置10において取得された放射線画像が画像処理装置30に入力されると(ステップST1)、特性取得部32がグリッド情報、被写体情報および撮影条件の少なくとも1つの入力部38からの入力を受け付け(グリッド情報等入力、ステップST2)、仮想グリッド特性、すなわち散乱線透過率Tsおよび一次線透過率Tpを取得する(ステップST3)。
そして、下記の式(11)により、周波数帯域毎の除去成分ΔI(x,y,r)を算出する。
=I(x,y,r)×{S(x,y)×(1-Ts(r))+(1-S(x,y))×(1-Tp(r))} …(11)
そして、除去成分ΔI(x,y,r)を周波数合成し、放射線画像I(x,y)から周波数合成した除去成分ΣrΔI(x,y,r)を減算して、処理済みの放射線画像I′(x,y)を取得する。
=I(x,y)-Σr{I(x,y,r)×{S(x,y)×(1-Ts(r))+(1-S(x,y))×(1-Tp(r))}}
…(12)
次いで、第2の実施形態において行われる処理について説明する。図5は第2の実施形態において行われる処理を示すフローチャートである。なお、ステップST11〜ステップST16までの処理は第1の実施形態におけるステップST1〜ステップST6の処理と同一であるため、ここでは詳細な説明は省略する。
ΔIp(x,y,r)=Ip(x,y)×(1-Tp(r))=(1-S(x,y))×I(x,y)×(1-Tp(r)) …(13)
そして下記の式(14)により、散乱線除去成分ΔIs(x,y,r)および一次線除去成分ΔIp(x,y,r)を周波数合成し、周波数合成した散乱線除去成分ΣrΔIs(x,y,r)および一次線除去成分ΣrΔIp(x,y,r)を放射線画像から減算して、処理済みの放射線画像I′(x,y)を算出する。
なお、上記第1および第2の実施形態において、放射線画像が放射線検出器への入射線量に比例した画素値を持つ、放射線量に対してリニアな空間において散乱線除去処理を行い、その後に対数変換を行って、人間の視覚と比例する対数リニアな空間に変換することが好ましい。
12 X線源
14 放射線検出器
20 制御装置
22 線源駆動制御部
24 検出器制御部
30 放射線画像処理装置
32 特性取得部
34 散乱情報取得部
36 散乱線除去部
38 入力部
40 表示部
42 記憶部
Claims (15)
- 被写体を透過した散乱線を除去するためのグリッドを使用して、前記被写体を撮影することにより取得された放射線画像に対して画像処理を行う放射線画像処理装置において、
前記放射線画像に対して前記使用したグリッドに起因する縞模様を除去する処理を行って、グリッド使用放射線画像を取得する手段と、
前記使用したグリッドとは異なる仮想的なグリッドの特性である仮想グリッド特性を取得する特性取得手段と、
前記仮想グリッド特性に基づいて、前記グリッド使用放射線画像の散乱線除去処理を行う散乱線除去手段とを備えたことを特徴とする放射線画像処理装置。 - 前記散乱線除去手段は、前記グリッド使用放射線画像中の散乱線成分および一次線成分が、前記仮想グリッド特性に対応するグリッドを用いて撮影を行うことにより取得した放射線画像に含まれる散乱線成分および一次線成分となるように、前記グリッド使用放射線画像を変換して前記散乱線除去処理を行う請求項1記載の放射線画像処理装置。
- 前記グリッド使用放射線画像に含まれる放射線の散乱成分を表す散乱成分情報を取得する散乱情報取得手段をさらに備え、
前記散乱線除去手段は、前記散乱成分情報にも基づいて前記散乱線除去処理を行う手段である請求項1または2記載の放射線画像処理装置。 - 前記散乱成分情報は、前記グリッド使用放射線画像内の散乱線含有率分布である請求項3記載の放射線画像処理装置。
- 前記散乱情報取得手段は、前記グリッド使用放射線画像から前記散乱線含有率分布を取得する手段である請求項4記載の放射線画像処理装置。
- 前記散乱情報取得手段は、前記グリッド使用放射線画像内の照射野の情報、前記被写体についての情報、および前記グリッド使用放射線画像の取得時の撮影条件の少なくとも1つに基づいて、前記散乱線含有率分布を取得する手段である請求項5記載の放射線画像処理装置。
- 前記特性取得手段は、前記散乱線に対する透過率および前記被写体を透過した放射線に含まれる一次線に対する透過率の少なくとも一方を、前記仮想グリッド特性として取得する手段である請求項1から6のいずれか1項記載の放射線画像処理装置。
- 前記特性取得手段は、前記仮想グリッド特性の指定を受け付けることにより、該仮想グリッド特性を取得する手段である請求項7記載の放射線画像処理装置。
- 前記特性取得手段は、グリッドの種類を表すグリッド情報、前記被写体についての情報、および前記グリッド使用放射線画像の取得時の撮影条件の少なくとも1つにより、前記仮想グリッド特性の指定を受け付ける手段である請求項8記載の放射線画像処理装置。
- 前記散乱線除去手段は、前記グリッド使用放射線画像に含まれる複数の周波数成分のうちの少なくとも1つの周波数成分に対して、前記散乱線除去処理を行う手段である請求項1から9のいずれか1項記載の放射線画像処理装置。
- 前記散乱線除去手段は、前記グリッド使用放射線画像を周波数分解して複数の周波数帯域毎の周波数成分を取得し、前記複数の周波数成分のうちの少なくとも1つの周波数成分に対して前記散乱線除去処理を行い、前記散乱線除去処理後の周波数成分およびこれ以外の周波数成分を合成して、処理済みの放射線画像を取得する手段である請求項10記載の放射線画像処理装置。
- 前記散乱線除去手段は、前記グリッド使用放射線画像から少なくとも1つの周波数帯域の周波数成分を抽出し、該抽出された周波数成分に対して前記散乱線除去処理を行い、処理後の周波数成分を前記グリッド使用放射線画像に加減算することにより、処理済みの放射線画像を取得する手段である請求項10記載の放射線画像処理装置。
- 前記特性取得手段は、前記抽出された周波数成分毎に前記仮想グリッド特性を取得する手段である請求項11または12記載の放射線画像処理装置。
- 被写体を透過した散乱線を除去するためのグリッドを使用して、前記被写体を撮影することにより取得された放射線画像に対して画像処理を行う放射線画像処理方法において、
前記放射線画像に対して前記使用したグリッドに起因する縞模様を除去する処理を行って、グリッド使用放射線画像を取得し、
前記使用したグリッドとは異なる仮想的なグリッドの特性である仮想グリッド特性を取得し、
前記仮想グリッド特性に基づいて、前記グリッド使用放射線画像の散乱線除去処理を行うことを特徴とする放射線画像処理方法。 - 被写体を透過した散乱線を除去するためのグリッドを使用して、前記被写体を撮影することにより取得された放射線画像に対して画像処理を行う放射線画像処理方法をコンピュータに実行させるための放射線画像処理プログラムにおいて、
前記放射線画像に対して前記使用したグリッドに起因する縞模様を除去する処理を行って、グリッド使用放射線画像を取得する手順と、
前記使用したグリッドとは異なる仮想的なグリッドの特性である仮想グリッド特性を取得する手順と、
前記仮想グリッド特性に基づいて、前記グリッド使用放射線画像の散乱線除去処理を行う手順とを有することを特徴とする放射線画像処理プログラム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013067794 | 2013-03-28 | ||
JP2013067794 | 2013-03-28 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013229940A Division JP6006193B2 (ja) | 2013-03-28 | 2013-11-06 | 放射線画像処理装置および方法並びにプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016172098A JP2016172098A (ja) | 2016-09-29 |
JP6006454B2 true JP6006454B2 (ja) | 2016-10-12 |
Family
ID=57008696
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016132257A Active JP6006454B2 (ja) | 2013-03-28 | 2016-07-04 | 放射線画像処理装置および方法並びにプログラム |
JP2016175200A Active JP6174217B2 (ja) | 2013-03-28 | 2016-09-08 | 放射線画像処理装置および方法並びにプログラム |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016175200A Active JP6174217B2 (ja) | 2013-03-28 | 2016-09-08 | 放射線画像処理装置および方法並びにプログラム |
Country Status (1)
Country | Link |
---|---|
JP (2) | JP6006454B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11317885B2 (en) | 2019-03-29 | 2022-05-03 | Fujifilm Corporation | Radiographic image processing device, method of operating radiographic image processing device, and radiographic image processing program |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3042285B1 (fr) * | 2015-10-09 | 2017-12-01 | Commissariat Energie Atomique | Procede de correction d’un spectre |
JP6776171B2 (ja) * | 2017-04-13 | 2020-10-28 | 株式会社日立製作所 | 放射線撮像装置、画像処理装置および画像処理用プログラム |
JP6743969B2 (ja) * | 2017-04-20 | 2020-08-19 | 株式会社島津製作所 | 放射線画像処理装置および放射線画像処理方法 |
JP7190344B2 (ja) * | 2018-12-11 | 2022-12-15 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
JP7069078B2 (ja) | 2019-03-29 | 2022-05-17 | 富士フイルム株式会社 | 放射線画像処理装置及び放射線画像処理方法 |
WO2021157403A1 (ja) * | 2020-02-07 | 2021-08-12 | キヤノン株式会社 | 画像処理装置および方法、プログラム |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61133847A (ja) * | 1984-12-04 | 1986-06-21 | Shoichi Yasui | X線画像における散乱線除去画像製造法 |
JP3423828B2 (ja) * | 1995-11-30 | 2003-07-07 | 株式会社日立メディコ | X線画像作成方法およびその装置 |
CN100510725C (zh) * | 2006-11-14 | 2009-07-08 | 北京国药恒瑞美联信息技术有限公司 | 用于消除散射辐射影响的虚拟滤线栅成像方法及其系统 |
CN102770075B (zh) * | 2010-02-24 | 2015-01-28 | 株式会社日立医疗器械 | X射线图像诊断装置和医用图像处理方法 |
-
2016
- 2016-07-04 JP JP2016132257A patent/JP6006454B2/ja active Active
- 2016-09-08 JP JP2016175200A patent/JP6174217B2/ja active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11317885B2 (en) | 2019-03-29 | 2022-05-03 | Fujifilm Corporation | Radiographic image processing device, method of operating radiographic image processing device, and radiographic image processing program |
Also Published As
Publication number | Publication date |
---|---|
JP2016198665A (ja) | 2016-12-01 |
JP6174217B2 (ja) | 2017-08-02 |
JP2016172098A (ja) | 2016-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6006193B2 (ja) | 放射線画像処理装置および方法並びにプログラム | |
JP6006454B2 (ja) | 放射線画像処理装置および方法並びにプログラム | |
JP6071144B2 (ja) | 放射線画像解析装置および方法並びにプログラム | |
JP6223781B2 (ja) | 放射線画像処理装置および方法並びにプログラム | |
JP6128463B2 (ja) | 放射線画像処理装置および方法並びにプログラム | |
JP6653629B2 (ja) | 放射線画像処理装置、方法およびプログラム | |
JP6042855B2 (ja) | 放射線画像撮影装置、放射線画像撮影方法、及び放射線画像撮影プログラム | |
JP6071853B2 (ja) | 放射線画像処理装置、方法およびプログラム | |
JP6301439B2 (ja) | 放射線画像解析装置および方法並びにプログラム | |
US10349913B2 (en) | Radiography apparatus, radiography method, and radiography program | |
JP6325141B2 (ja) | 放射線画像撮影装置、放射線画像撮影方法、及び放射線画像撮影プログラム | |
JP6313402B2 (ja) | 放射線画像撮影装置、放射線画像撮影方法、及び放射線画像撮影プログラム | |
JP6321267B2 (ja) | 画像解析装置および方法並びにプログラム | |
JP6285006B2 (ja) | 放射線画像処理装置、方法およびプログラム | |
Konstantinidis et al. | Digital Breast Tomosynthesis: Systems, Characterization and Simulation | |
JP2023068496A (ja) | 放射線画像処理装置、方法およびプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160704 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20160715 |
|
TRDD | Decision of grant or rejection written | ||
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20160816 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20160823 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160908 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6006454 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
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 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |