JP6917356B2 - 主成分材料組み合わせによる分解方法および装置 - Google Patents
主成分材料組み合わせによる分解方法および装置 Download PDFInfo
- Publication number
- JP6917356B2 JP6917356B2 JP2018240727A JP2018240727A JP6917356B2 JP 6917356 B2 JP6917356 B2 JP 6917356B2 JP 2018240727 A JP2018240727 A JP 2018240727A JP 2018240727 A JP2018240727 A JP 2018240727A JP 6917356 B2 JP6917356 B2 JP 6917356B2
- Authority
- JP
- Japan
- Prior art keywords
- region
- decomposition
- point
- combination
- main component
- 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
- 238000000354 decomposition reaction Methods 0.000 title claims description 131
- 238000000034 method Methods 0.000 title claims description 51
- 239000000463 material Substances 0.000 claims description 104
- 239000000126 substance Substances 0.000 claims description 16
- 238000001228 spectrum Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 13
- 230000006870 function Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 7
- 238000005192 partition Methods 0.000 claims description 4
- 238000002591 computed tomography Methods 0.000 description 22
- 230000009977 dual effect Effects 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 5
- 238000013170 computed tomography imaging Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 238000013016 damping Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 2
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 2
- 229940039231 contrast media Drugs 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000699670 Mus sp. Species 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4266—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a plurality of detector units
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/488—Diagnostic techniques involving pre-scan acquisition
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/06—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
- G01N23/083—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays
- G01N23/087—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption the radiation being X-rays using polyenergetic X-rays
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/408—Dual energy
Description
［優先権の主張］
本開示は、２０１７年１２月２６日に出願された中国特許出願第２０１７１１４３２８４０．５号に基づき、その優先権を主張するものである。この出願の開示は、ここにその全体を参照援用する。
［技術分野］
本開示は放射線イメージングの分野に関し、特に、主成分材料組み合わせに基づく分解方法および分解装置に関する。
本開示は、２０１７年１２月２６日に出願された中国特許出願第２０１７１１４３２８４０．５号に基づき、その優先権を主張するものである。この出願の開示は、ここにその全体を参照援用する。
［技術分野］
本開示は放射線イメージングの分野に関し、特に、主成分材料組み合わせに基づく分解方法および分解装置に関する。
関連するＸ線デュアルエネルギーＣＴ（コンピュータ断層撮影）イメージング法は、物質の減衰係数関数を、２つの既知の主成分材料の減衰係数関数の線形結合に分解し、その線形結合の係数は解くべき未知パラメータである。スキャン対象の分解には、一定セットのダブル主成分材料と、対応する一定セットの減衰係数関数とを用いる。たとえば、水と骨が医療用デュアルエネルギーＣＴにおけるダブル主成分材料として選択され得る。
本開示の一態様は、主成分材料の組み合わせに基づく分解方法を提供する。該分解方法は、スキャン対象の各点の分解係数と等価原子番号とを取得する取得ステップと、前記スキャン対象の各点の等価原子番号にしたがって、前記スキャン対象を複数の領域に分割する領域分割ステップであって、分割領域の各々で用いられる主成分材料の組み合わせは互いに異なるステップと、各分割領域について、前記領域で用いられる前記主成分材料の組み合わせに基づいて、前記領域中の各点の分解係数を決定する分解係数決定ステップと、すべての分割領域のすべての点の分解係数に基づく分解係数の変化が所定閾値より小さくない場合、前記分割領域中のすべての点の決定された分解係数と、前記分割領域中で用いられる前記主成分材料の組み合わせのすべての主成分材料の原子情報とに基づいて、前記分割領域の各々について、前記分割領域中のすべての点の等価原子番号を決定する等価原子番号決定ステップと、全分割領域中のすべての点の分解係数に基づいて分解係数の変化が所定閾値より小さくなるまで、前記領域分割ステップと、前記分解係数決定ステップと、等価原子番号決定ステップとをループするステップと、を含む。
本開示の幾つかの実施形態では、前記領域分割ステップは、所定の原子番号関係が同じ副区間にあるすべての点を、原子番号関係チェーン中のスキャン対象の各点の等価原子番号の位置により１つの領域に分割し、副区間の両端の原子番号に対応する材料が、その領域の主成分材料の組み合わせを構成する主成分材料と決定されることを含む。
本開示の幾つかの実施形態では、前記分解係数決定ステップは、各分割領域について、様々なエネルギーの構成の投影データと、前記領域のエネルギースペクトルと、前記領域の主成分材料の組み合わせに対応する減衰係数関数とにしたがって、各光線が通過する各経路に沿って、前記領域の各点の分解係数の積分を計算することと、各光線が通過する各経路に沿って、前記領域の各点の分解係数の積分にしたがって、前記領域の各点の分解係数を決定することとを含む。
本開示の幾つかの実施形態では、領域のエネルギースペクトルは、光源のエネルギースペクトルと、前記領域の外の他の領域のすべての点の分解係数と、前記他の領域の主成分材料の組み合わせに対応する減衰係数関数とにしたがって決定される。
本開示の幾つかの実施形態では、前記等価原子番号決定ステップは、前記領域の各点の分解係数と、前記領域の主成分材料の組み合わせのすべての主成分材料の原子番号と原子量とにしたがって、分割領域の各点の等価電子密度を決定することと、前記領域の各点の分解係数及び等価電子密度と、前記領域の主成分材料の組み合わせのすべての主成分材料の原子番号及び原子量とにしたがって、前記領域の各点の等価原子番号を決定することとを含む。
本開示の幾つかの実施形態では、取得ステップは、所定の主成分材料組み合わせに応じて、スキャン対象の全ての点の分解係数と等価原子番号とを取得する。
本開示の幾つかの実施形態では、前記方法はさらに、分解係数の変化が所定閾値より小さい場合、前記領域の各点の等価原子番号を用いて、各分割領域の物質識別を実行する。
本開示の幾つかの実施形態では、点はピクセル点またはボクセル点である。
本開示の他の一態様は、主成分材料の組み合わせに基づく分解装置であって、メモリと、前記メモリに結合したプロセッサであって、前記メモリに記憶された命令に基づいて、実施形態の方法を実行するように構成されたプロセッサとを有する、
本開示のもう１つの態様は、コンピュータプログラムを記憶するコンピュータ可読記憶媒体であって、前記コンピュータプログラムがプロセッサにより実行され、実施形態の方法を実施する、コンピュータ可読記憶媒体を含む。
本開示のもう１つの態様は、コンピュータプログラムを記憶するコンピュータ可読記憶媒体であって、前記コンピュータプログラムがプロセッサにより実行され、実施形態の方法を実施する、コンピュータ可読記憶媒体を含む。
本開示において、スキャン対象は複数の領域に分割され、各分割領域で用いられる主成分材料の組み合わせは互いに異なり、分解係数の変化が一定条件を満たすまで、スキャン対象は、各点の再決定される等価原子番号にしたがって、再分割される。これにより、動的な主成分材料の組み合わせに基づく分解方法が実現され、主成分材料の組み合わせの妥当でない選択により生じる分解エラーが減少し、マルチエネルギーＣＴの分解と物質識別の精度が向上する。（ヨウ素などの物質の造影剤を用いる）医療用ＣＴエンハンスイメージング、（セキュリティ、航空荷物・コンテナＣＴイメージングなどの）高エネルギーまたはマルチエネルギーＣＴイメージングなどの、スキャン対象の原子番号のばらつき範囲が大きい場合、利点は特に明らかである。
実施形態または関連技術の説明で用いる図面を簡単に説明する。図面を参照して以下の詳細な説明を読めば、本開示を明確に理解することができる。
言うまでもなく、以下に説明する図面は本開示の実施形態の一部に過ぎない。当業者は添付した図面から発明的努力をせずにその他の図面を作成できるだろう。
本開示の主成分材料の組み合わせに基づく分解方法の一実施形態を示すフロー図である。
数式（６−１３）を計算に用いる動的主成分材料の組み合わせに基づく反復実行分解方法を示すフロー図である。
一セットのシミュレーション実験結果を示す図である。
本開示の主成分材料の組み合わせに基づく分解装置の一実施形態を示す概略的構造図である。
本開示の主成分材料の組み合わせに基づく分解装置の他の一実施形態を示す概略的構造図である。
発明者は、識別すべき物質の原子番号のばらつきが小さい範囲にある場合、一定の主成分材料の組み合わせに基づく分解方法はよい効果を有し、原子番号の範囲が大きい場合、大きな分解エラーが生じることを見いだした。例えば、原子番号の範囲が炭素（原子番号Ｚ＝６であり、有機物を表す）から鉛（Ｚ＝８２であり、重金属を表す）の物質の減衰係数関数は非常に異なる。一定セットの主成分材料が使われる場合、例えば、炭素（Ｚ＝６）とすず（Ｚ＝５０）が選択された場合、スキャン対象中のＺ＞５０の物質とＺ＝２８近くの物質との分解エラーは非常に大きい。炭素（Ｚ＝６）と鉛（Ｚ＝８２）が選択された場合、スキャン対象中の３０＜Ｚ＜６０の物質の分解エラーは非常に大きい。すず（Ｚ＝５０）と鉛（Ｚ＝８２）が選択された場合、スキャン対象中のＺ＜５０の物質の分解エラーは非常に大きい。
本開示により解決される１つの技術的問題は、主成分材料の組み合わせの妥当でない選択により生じる分解エラーを低減することである。
本開示の実施形態における技術的解決策は、本開示の実施形態における添付の図面を参照して、明確かつ完全な方法で説明される。
図１は、本開示の主成分材料の組み合わせに基づく分解方法の一実施形態を示すフロー図である。図１に示すように、本実施形態の方法はステップ１００〜１４０を含む。
取得ステップと呼ぶことができるステップ１００では、スキャン対象の各点の分解係数や等価原子番号などの初期情報を取得することができる。
取得ステップは、一定の主成分材料組み合わせに応じて、スキャン対象の全ての点の等価原子番号と分解係数とを取得する。すなわち、スキャン対象の全ての点の等価原子番号と分解係数とを、一定の主成分材料組み合わせに基づく分解結果に応じて決定することができる。一定の主成分材料組み合わせに基づく分解方法は、従来技術における一定の二重主成分材料に基づく分解法を参照することで実現することができる。
領域分割ステップと言えるステップ１１０において、スキャン対象は、その各点の等価原子番号にしたがって、複数の領域に分割され得る。分割される各領域で用いられる主成分材料の組み合わせは互いに異なり、１つの領域は同じ主成分材料の組み合わせを採用する。
点はピクセル点またはボクセル点である。スキャン対象は、単一材料でないとき、通常は複数の領域に分割される。コンピュータの処理速度が十分速ければ、スキャン対象の各ピクセル点またはボクセル点は、多くても１つの領域に分割され得る。コンピュータの処理速度が限られている場合、等価原子番号の相違が所定範囲より小さいスキャン対象の複数の点は、１つの領域に分割できる。
分解係数決定ステップまたは分解ステップと言えるステップ１２０において、分割される各領域に関して、その中の各点の分解係数を、その領域で用いられる主成分材料の組み合わせに基づき決定できる。
等価原子番号決定ステップと言えるステップ１３０において、すべての分割領域のすべての点の分解係数に基づく分解係数の変化が所定閾値より小さい場合、各分割領域について、分割領域のすべての点の決定された分解係数と、分割領域で用いられる主成分材料の組み合わせのすべての主成分材料の原子情報とに基づいて、分割領域のすべての点の原子番号を再決定することができる。
原子情報は例えば原子番号と原子量とを含む。
ループステップと言えるステップ１４０において、分解係数の変化が所定閾値より小さくなるまで、前記領域分割ステップと、前記分解係数決定ステップと、等価原子番号決定ステップとをループするステップとを、全分割領域中のすべての点の分解係数にしたがって、反復的に実行することができる。
例えば、ある領域のすべての点の分解係数の変化が所定閾値より小さい場合、その領域の分解係数は収束したと見なす。
各分割領域の物質識別は、収束後の各領域の各点の等価原子番号を用いて、実行できる。分解係数の変化が所定閾値より小さい場合、前記領域の各点の等価原子番号を用いて、各分割領域の物質識別を実行することができる。
上記の実施形態において、スキャン対象は複数の領域に分割され、各分割領域で用いられる主成分材料の組み合わせは互いに異なり、分解係数の変化が一定条件を満たすまで、スキャン対象は、各点の再決定される等価原子番号にしたがって、再分割される。これにより、動的な主成分材料の組み合わせに基づく分解方法が実現され、主成分材料の組み合わせの妥当でない選択により生じる分解エラーが減少し、マルチエネルギーＣＴの分解と物質識別の精度が向上する。
（ヨウ素などの物質の造影剤を用いる）医療用ＣＴエンハンスイメージング、（セキュリティ、航空荷物・コンテナＣＴイメージングなどの）高エネルギーまたはマルチエネルギーＣＴイメージングなどの、スキャン対象の原子番号のばらつき範囲が大きい場合、利点は特に明らかである。この方法は、臨床医療、セキュリティ、産業用非破壊試験その他の分野で重要な適用価値を有する。
本開示の主成分材料の組み合わせに基づく分解方法は、マルチエネルギーフォトンカウンティング検出器に基づく３エネルギーＣＴデータやエネルギースペクトルＣＴデータなど、より多くのエネルギースペクトルの投影データが取得できる場合、デュアルエネルギー（エネルギースペクトル）ＣＴに適用可能であるだけでなく、マルチエネルギーＣＴにも適用可能である。デュアルエネルギーＣＴでは、各領域は２つの主成分材料の組み合わせを用いる。マルチエネルギーＣＴでは、各領域は、エネルギースペクトルと同じ数の主成分材料の組み合わせを必要とする。例えば、３エネルギーＣＴでは、各領域は３つの主成分材料の組み合わせを必要とする。後述の実施形態の具体的な数式は、デュアルエネルギーＣＴを例として記述されている。本開示のデュアルエネルギーＣＴを例の２つの主成分材料に基づく分解方法により、当業者は、創作的努力を要せずに、マルチエネルギーＣＴの分解方法を得ることができる。
ステップ１２０の例示的な実施は、反復的方法を用いて実施される。反復的方法は、ｔ＝１から開始して、領域ｔの主成分材料の組み合わせに基づき、領域ｔの各点の分解係数を計算し、次いで、領域ｔの各点の分解係数を次の計算の既知パラメータとして、領域ｔ＋１の主成分材料の組み合わせに基づいて領域ｔ＋１の各点の分解係数を計算し、ｔがＮまで増加したら終了し、各領域の各点の分解係数を取得する。Ｎは領域数を示す。
領域ｔの主成分材料の組み合わせに基づき、ある領域ｔの各点の分解係数の計算方法は、様々なエネルギーの構成の投影データ（例えば、高エネルギーの投影データと、低エネルギーの投影データ）と、前記領域のエネルギースペクトルと、前記領域の主成分材料の組み合わせに対応する減衰係数関数とにしたがって、各光線が通過する各経路に沿って、前記領域の各点の分解係数の積分を計算することと、各光線が通過する各経路に沿って、前記領域の各点の分解係数の積分にしたがって、所定の再構成アルゴリズムで、前記領域の各点の分解係数を決定することとを含む。再構成アルゴリズムは、例えば、フィルタードバックプロジェクション（ｆｉｌｔｅｒｅｄ ｂａｃｋ ｐｒｏｊｅｃｔｉｏｎ（ＦＢＰ））であってもよい。
ある領域ｔの主成分材料の組み合わせに基づき、その領域ｔの各点の分解係数の計算方法は、本実施形態では２つの主成分材料の組み合わせを例に取って説明し、数式（１−１０）を参照して具体的に実施することができる。
スキャン対象がＮ個に分割され、
各領域で用いられる２つの主成分材料の減衰係数関数は
２つの主成分材料に基づく各点（ピクセル点またはボクセル点）の分解係数（または分解係数画像）は
Ｘ線デュアルエネルギーＣＴの高（Ｈ）／低（Ｌ）エネルギーを有するＸ線の投影データは、
「ｔ」は現在解こうとしている領域（ｃｕｒｒｅｎｔｌｙ ｓｏｌｖｅｄ ｒｅｇｉｏｎ）を表す。領域ｔ以外の領域の分解係数画像は既知であると仮定する。領域ｔの分解係数画像は、上記の分解方法で取得でき、デュアルエネルギー投影の数式（４）から、
領域ｔ以外の領域の分解係数画像ｂ１ ｊ、ｂ２ ｊ（ｊ＝１，２，・・・Ｎ；ｊ≠ｔ）は既知であり、ＳＬ／Ｈ（Ｅ）も既知であるから、領域ｔのエネルギースペクトルＳＬ／Ｈ（Ｅ）を計算でき、光線ｒａｙｉが通る経路に沿った領域ｔの分解係数画像の積分を、デュアルエネルギー投影の数式（６）により、すなわち
各ピクセル／ボクセル点の分解係数は、各光線が通る各経路に沿った領域ｔの分解係数画像の積分にしたがって、一般的なＣＴ再構成アルゴリズム（ＦＢＰアルゴリズムなど）で再構成して、すなわち
領域ｔに限定されているので、最終的な分解係数画像は
ステップ１３０において、ある領域の各点の等価原子番号を決定する例示的方法は、前記領域の各点の分解係数と、前記領域の主成分材料の組み合わせのすべての主成分材料の原子番号と原子量とにしたがって、分割領域の各点の等価電子密度を決定することと、前記領域の各点の分解係数及び等価電子密度と、前記領域の主成分材料の組み合わせのすべての主成分材料の原子番号及び原子量とにしたがって、前記領域の各点の等価原子番号を決定することとを含む。
２つの主成分材料の組み合わせを例に取ると、各点の等価電子密度ρｅｆｆと等価原子番号Ｚｅｆｆの計算式は、
パラメータαは経験的に
ステップ１１０におけるスキャン対象の例示的な領域分割方法は、所定の原子番号関係が同じ副区間にあるすべての点を、原子番号関係チェーン中のスキャン対象の各点の等価原子番号の位置により１つの領域に分割し、副区間の両端の原子番号に対応する材料が、その領域の主成分材料の組み合わせを構成する主成分材料と決定されることを含む。
原子番号関係チェーン（ａｔｏｍｉｃ ｎｕｍｂｅｒ ｒｅｌａｔｉｏｎｓｈｉｐ ｃｈａｉｎ）の長さ（すなわち、数式１３のＮの大きさ）は、物質識別の正確性に関係する。原子番号関係チェーンの長さが長いほど、オブジェクトはより多くの領域に分割され、物質識別の正確性は向上する。原子番号関係チェーンにおける領域境界を定めるのに用いられる原子番号について、一般的な材料に対応する原子番号を選択することができ、または検出対象に含まれる材料にしたがって設定することができる。
上記の領域分割方法の式は、
上記の通り、本開示の動的な主成分材料の組み合わせに基づく分解方法は、すべての領域の分解係数が収束する前に、反復的に実行される。図２は、数式（６−１３）を計算に用いる動的主成分材料の組み合わせに基づく反復実行分解方法を示すフロー図である。式（Ｅｑ）は数式を表し、例えば、式（１０）は数式１０を表す。
図３は、一セットのシミュレーション実験結果を示す。画像サイズは２５６×２５６であり、ピクセルサイズは１．６ｍｍ×１．６ｍｍであり、ディテクタユニット数は７３６であり、回転角の数は５４０であり、回転角は３６０°である。デュアルエネルギーＣＴ光源エネルギースペクトルの最高エネルギーは、それぞれ９ＭｅＶと６ＭｅＶである。モデルの大きいシリンダーの材料は水であり、他の小さいシリンダーは内側グループと外側グループとに分割される。内側グループと外側グループとは同じ材料を有するが、材料の原子番号は１３、２０、２６、３０、３４、４０、５０、６０、７０、８０（反時計回り、色が薄いものから濃いものまで）である。図３の３つの画像は、左から右へ、実際の原子番号画像と、固定した２つの主成分材料に基づいて従来の分解方法で取得した原子番号画像であって、原子番号が６と５０の材料を固定の２つの主成分材料として用いたものと、本開示の動的な主成分材料の組み合わせに基づく分解方法による取得された原子番号画像とである。シミュレーション結果から、本開示は、Ｘ線デュアルエネルギーＣＴの分解精度を効果的に改善でき、デュアルエネルギーＣＴイメージングモデルと物理的プロセスとの間の材料の非一貫性により生じる分解エラーを低減でき、デュアルエネルギーＣＴのイメージング品質と物質識別を改善する重要性と適用価値を有する。
本開示は主成分材料の組み合わせに基づく分解装置も提案する。該分解装置は、上記の分解方法を実行する一以上のモジュールを有する。
図４は、本開示による主成分材料の組み合わせに基づく分解装置の一実施形態を示す概略的構造図である。
図４に示すように、分解装置は主成分材料の組み合わせ４０に基づき、モジュール４００−４４０を含む。
取得モジュール４００は、スキャン対象の各点の分解係数と等価原子番号とを取得するように構成される。
領域分割モジュール４１０は、前記スキャン対象の各点の等価原子番号にしたがって、前記スキャン対象を複数の領域に分割し、分割領域の各々で用いられる主成分材料の組み合わせは互いに異なるように構成される。
分解係数決定モジュール４２０は、各分割領域について、前記領域で用いられる前記主成分材料の組み合わせに基づいて、前記領域中の各点の分解係数を決定するように構成される。
等価原子番号決定モジュール４３０は、すべての分割領域のすべての点の分解係数に基づく分解係数の変化が所定閾値より小さくない場合、前記分割領域中のすべての点の決定された分解係数と、前記分割領域中で用いられる前記主成分材料の組み合わせのすべての主成分材料の原子情報とに基づいて、前記分割領域の各々について、前記分割領域中のすべての点の等価原子番号を決定するように構成される。
ループモジュール４４０は、全分割領域中のすべての点の分解係数に基づいて分解係数の変化が所定閾値より小さくなるまで、前記領域分割ステップと、前記分解係数決定ステップと、等価原子番号決定ステップとをループするように構成される。
任意的に、サイクルの始めに必要なスキャン対象のすべての点の分解係数と等価原子番号とは、所定の（固定の）主成分材料の組み合わせにより決定される。
任意的に、領域分割モジュール４１０は、所定の原子番号関係が同じ副区間にあるすべての点を、原子番号関係チェーン中のスキャン対象の各点の等価原子番号の位置により１つの領域に分割するように構成され、副区間の両端の原子番号に対応する材料が、その領域の主成分材料の組み合わせを構成する主成分材料と決定されることを含む。
任意的に、分解係数決定モジュール４２０は、反復的方法を用いて実装されるように構成される。前記反復的方法は、ｔ＝１から開始して、領域ｔの主成分材料の組み合わせに基づき、領域ｔの各点の分解係数を計算し、次いで、領域ｔの各点の分解係数を次の計算の既知パラメータとして、領域ｔ＋１の主成分材料の組み合わせに基づいて領域ｔ＋１の各点の分解係数を計算し、ｔがＮまで増加したら終了し、各領域の各点の分解係数を取得する。Ｎは領域数を示す。
任意的に、領域ｔの主成分材料の組み合わせに基づき、ある領域ｔの各点の分解係数の計算方法は、様々なエネルギーの構成の投影データ（例えば、高エネルギーの投影データと、低エネルギーの投影データ）と、前記領域のエネルギースペクトルと、前記領域の主成分材料の組み合わせに対応する減衰係数関数とにしたがって、各光線が通過する各経路に沿って、前記領域の各点の分解係数の積分を計算することと、各光線が通過する各経路に沿って、前記領域の各点の分解係数の積分にしたがって、所定の再構成アルゴリズムで、前記領域の各点の分解係数を決定することとを含む。再構成アルゴリズムは、例えば、フィルタードバックプロジェクション（ｆｉｌｔｅｒｅｄ ｂａｃｋ ｐｒｏｊｅｃｔｉｏｎ（ＦＢＰ））であってもよい。
任意的に、等価原子番号決定モジュール４３０は、前記領域の各点の分解係数と、前記領域の主成分材料の組み合わせのすべての主成分材料の原子番号と原子量とにしたがって、分割領域の各点の等価電子密度を決定し、前記領域の各点の分解係数及び等価電子密度と、前記領域の主成分材料の組み合わせのすべての主成分材料の原子番号及び原子量とにしたがって、前記領域の各点の等価原子番号を決定するように構成されてもよい。
図５は、本開示の主成分材料の組み合わせに基づく分解装置の他の一実施形態を示す概略的構造図である。
図５に示すように、主成分材料の組み合わせに基づく分解装置５０は、メモリ５１０と、そのメモリ５１０に結合したプロセッサ５２０とを含む。プロセッサ５２０は、メモリ５１０に記憶された命令に基づき、前記実施形態のいずれかの分解方法を実行するように構成されている。
メモリ５１０は例えば、システムメモリ、固定不揮発性記憶媒体などを含む。システムメモリは、例えば、オペレーティングシステム、アプリケーション、ブートローダ、及びその他のプログラムを記憶する。
装置５０は、入出力インターフェース５３０、ネットワークインターフェース５４０、ストレージインターフェース５５０なども含んでもよい。これらのインターフェース５３０、５４０、５５０と、メモリ５１０と、プロセッサ５２０とは、例えば、バス５６０を介して接続されていてもよい。入出力インターフェース５３０は、ディスプレイ、マウス、キーボード、及びタッチスクリーンなどの入出力デバイスに対して接続インターフェースを提供する。ネットワークインターフェース５４０は、ネットワークされた様々なデバイスに接続インターフェースを提供する。ストレージインターフェース５５０は、ＳＤカードやＵＳＢフラッシュドライブなどの外部ストレージデバイスに対して接続インターフェースを提供する。
本開示は、コンピュータプログラムを記憶するコンピュータ可読記憶媒体であって、前記コンピュータプログラムがプロセッサにより実行され、上記の実施形態のいずれかの分解方法を実施する、コンピュータ可読記憶媒体も提供する。
当業者には言うまでもなく、本開示の実施形態は、方法、システム、又はコンピュータプログラム製品として提供し得る。したがって、本開示は、完全にハードウェアの実施形態、完全にソフトウェアの実施形態、またはハードウェアとソフトウェアを両方組み合わせた実施形態の形式を取ることができる。さらに、本開示は、コンピュータ利用可能プログラムコードを含む一以上のコンピュータ利用可能非一時的記憶媒体（ディスクストレージ、ＣＤ−ＲＯＭ、光学的ストレージなどを含む）の形式を取ることもできる。
上記の説明は、単に、本開示の好ましい実施形態であり、本開示を限定するものではない。本開示の精神と原理に含まれるいかなる修正、等価物、改良なども、本開示の保護範囲に含まれる。
Claims (9)
- 主成分材料の組み合わせに基づく分解方法であって、
マルチエネルギーＣＴによりスキャンされるスキャン対象の各点の分解係数と等価原子番号とを取得する取得ステップであって、前記点はピクセル点またはボクセル点である取得ステップと、
前記スキャン対象の各点の等価原子番号にしたがって、前記スキャン対象を複数の領域に分割する領域分割ステップであって、分割領域の各々で用いられる主成分材料の組み合わせは互いに異なるステップと、
各分割領域について、前記領域で用いられる前記主成分材料の組み合わせに基づいて、前記領域中の各点の分解係数を決定する分解係数決定ステップと、
すべての分割領域のすべての点の分解係数に基づく分解係数の変化が所定閾値より小さくない場合、前記分割領域中のすべての点の決定された分解係数と、前記分割領域中で用いられる前記主成分材料の組み合わせのすべての主成分材料の原子情報とに基づいて、前記分割領域の各々について、前記分割領域中のすべての点の等価原子番号を決定する等価原子番号決定ステップと、
全分割領域中のすべての点の分解係数に基づいて分解係数の変化が所定閾値より小さくなるまで、前記領域分割ステップと、前記分解係数決定ステップと、等価原子番号決定ステップとをループするステップと、
を含む方法。 - 前記領域分割ステップは、
所定の原子番号関係チェーンにおける前記スキャン対象の各点の等価原子番号の位置を決定するステップであって、前記原子番号関係チェーンの同じ副区間にある各点に対応する前記スキャン対象のすべての点を、１つの領域とし、副区間の両端の原子番号に対応する材料が、その領域の主成分材料の組み合わせを構成する主成分材料と決定されるステップを含み、
前記原子番号関係チェーンは複数の原子番号よりなる関係チェーンであって、前記原子番号関係チェーンは境界として複数の原子番号を有する複数の副区間に分割される、
請求項１に記載の方法。 - 前記分解係数決定ステップは、各分割領域について、
前記マルチエネルギーＣＴの光源により放射される様々なエネルギーの光線の投影データと、前記領域の光線のエネルギースペクトルと、前記領域の主成分材料の組み合わせに対応する減衰係数関数とにしたがって、各光線が通過する各経路に沿って、前記領域の各点の分解係数の積分を計算することと、
各光線が通過する各経路に沿って、前記領域の各点の分解係数の積分にしたがって、前記領域の各点の分解係数を決定することとを含む、
請求項１に記載の方法。 - 前記領域の光線のエネルギースペクトルは、前記マルチエネルギーＣＴの光源のエネルギースペクトルと、前記領域の外の他の領域のすべての点の分解係数と、前記他の領域の主成分材料の組み合わせに対応する減衰係数関数とにしたがって決定される、
請求項３に記載の方法。 - 前記等価原子番号決定ステップは、
前記領域の各点の分解係数と、前記領域の主成分材料の組み合わせのすべての主成分材料の原子番号と原子量とにしたがって、分割領域の各点の等価電子密度を決定することと、
前記領域の各点の分解係数及び等価電子密度と、前記領域の主成分材料の組み合わせのすべての主成分材料の原子番号及び原子量とにしたがって、前記領域の各点の等価原子番号を決定することとを含む、
請求項１に記載の方法。 - 前記取得ステップは、所定の主成分材料組み合わせに応じて、スキャン対象の全ての点の分解係数と等価原子番号とを取得する、
請求項１に記載の方法。 - 分解係数の変化が所定閾値より小さい場合、前記領域の各点の等価原子番号を用いて、各分割領域の物質識別を実行する、
請求項１に記載の方法。 - 主成分材料の組み合わせに基づく分解装置であって、
メモリと、
前記メモリに結合したプロセッサであって、前記メモリに記憶された命令に基づいて、請求項１ないし７いずれか一項に記載の方法を実行するように構成されたプロセッサとを有する、
分解装置。 - コンピュータプログラムを記憶するコンピュータ可読記憶媒体であって、前記コンピュータプログラムがプロセッサにより実行され、請求項１ないし７いずれか一項に記載の方法を実施する、コンピュータ可読記憶媒体。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711432840.5 | 2017-12-26 | ||
| CN201711432840.5A CN109959669B (zh) | 2017-12-26 | 2017-12-26 | 基材料分解方法和装置及计算机可读存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2019117193A JP2019117193A (ja) | 2019-07-18 |
| JP6917356B2 true JP6917356B2 (ja) | 2021-08-11 |
Family
ID=65364536
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018240727A Active JP6917356B2 (ja) | 2017-12-26 | 2018-12-25 | 主成分材料組み合わせによる分解方法および装置 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US10448904B2 (ja) |
| JP (1) | JP6917356B2 (ja) |
| CN (1) | CN109959669B (ja) |
| DE (1) | DE102018133435A1 (ja) |
| GB (1) | GB2571192B (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111134709B (zh) * | 2020-01-17 | 2021-09-14 | 清华大学 | 一种多能量ct基材料分解方法 |
| CN113100803A (zh) * | 2021-04-20 | 2021-07-13 | 西门子数字医疗科技(上海)有限公司 | 用于显示静脉血栓的方法、装置、计算机设备和介质 |
| JP2023009382A (ja) * | 2021-07-07 | 2023-01-20 | 浜松ホトニクス株式会社 | 元素分布測定装置および元素分布測定方法 |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7298812B2 (en) * | 2006-03-31 | 2007-11-20 | General Electric Company | Image-based material decomposition |
| US20090052621A1 (en) * | 2007-08-23 | 2009-02-26 | Deborah Joy Walter | Method and apparatus for basis material decomposition with k-edge materials |
| US8260023B2 (en) * | 2008-11-26 | 2012-09-04 | General Electric Company | Forward projection for the generation of computed tomography images at arbitrary spectra |
| US20100310036A1 (en) * | 2009-06-04 | 2010-12-09 | General Electric Company | Computed tomography method and apparatus |
| CN102313753B (zh) * | 2010-06-30 | 2014-07-16 | 清华大学 | 物品检测设备及其检测方法 |
| RU2587556C2 (ru) * | 2010-12-13 | 2016-06-20 | Конинклейке Филипс Электроникс Н.В. | Оценка плотности молочной железы |
| US9435900B2 (en) * | 2011-06-08 | 2016-09-06 | The Board Of Trustees Of The Leland Stanford Junior University | X-ray system utilizing alternating spectrum x-ray source in conjuction with a fixed spectrum separation filter approval |
| DE102012215515A1 (de) * | 2012-08-31 | 2014-03-06 | Klinikum Der Universität München | Verfahren zur Detektion von Schäden an Silikonimplantaten und Computertomographiegerät |
| US8811709B2 (en) * | 2012-11-14 | 2014-08-19 | General Electric Company | System and method for multi-material correction of image data |
| JP6310473B2 (ja) * | 2012-11-26 | 2018-04-11 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | プロジェクションデータからノイズ除去を行う方法、プロセッサ及びコンピュータプログラム |
| US8885910B2 (en) * | 2012-12-28 | 2014-11-11 | General Electric Company | Systems and methods for X-ray imaging |
| US20140236488A1 (en) * | 2013-02-19 | 2014-08-21 | Mindways Software, Inc. | Computed Tomography Calibration Systems and Methods |
| CN104700389B (zh) * | 2013-12-09 | 2019-08-13 | 通用电气公司 | 双能ct扫描图像中的对象识别方法 |
| US9977140B2 (en) * | 2015-01-09 | 2018-05-22 | Toshiba Medical Systems Corporation | More efficient method and apparatus for detector response correction and material decomposition of projection data obtained using photon-counting detectors |
| CN107427276B (zh) * | 2015-03-30 | 2021-01-26 | 株式会社日立制作所 | 图像生成装置、图像生成方法以及x射线ct装置 |
| DE102016221684A1 (de) * | 2015-11-27 | 2017-06-01 | Siemens Healthcare Gmbh | Verfahren und Bilddatenverarbeitungseinrichtung zum Verarbeiten eines Multi-Energie-Computertomographie-Bilddatensatzes |
| CN107356615B (zh) * | 2016-05-10 | 2020-01-21 | 清华大学 | 一种用于双能x射线ct的方法和系统 |
-
2017
- 2017-12-26 CN CN201711432840.5A patent/CN109959669B/zh active Active
-
2018
- 2018-12-21 DE DE102018133435.4A patent/DE102018133435A1/de active Pending
- 2018-12-21 GB GB1820957.7A patent/GB2571192B/en active Active
- 2018-12-25 JP JP2018240727A patent/JP6917356B2/ja active Active
- 2018-12-26 US US16/232,072 patent/US10448904B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109959669A (zh) | 2019-07-02 |
| US10448904B2 (en) | 2019-10-22 |
| GB2571192A (en) | 2019-08-21 |
| CN109959669B (zh) | 2021-02-23 |
| GB2571192B (en) | 2020-04-29 |
| DE102018133435A1 (de) | 2019-06-27 |
| US20190192090A1 (en) | 2019-06-27 |
| GB201820957D0 (en) | 2019-02-06 |
| JP2019117193A (ja) | 2019-07-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107356615B (zh) | 一种用于双能x射线ct的方法和系统 | |
| US8855395B2 (en) | Conditional likelihood material decomposition and methods of using the same | |
| US10139354B2 (en) | Spectral X-ray imaging | |
| JP5290145B2 (ja) | ダーティアイソトープのｐｅｔ再構成 | |
| JP6917356B2 (ja) | 主成分材料組み合わせによる分解方法および装置 | |
| Weidinger et al. | Polychromatic iterative statistical material image reconstruction for photon-counting computed tomography | |
| JP5853123B2 (ja) | 光子計数検出装置のための不完全性の動的モデリング | |
| CN110660111B (zh) | Pet散射校正及图像重建方法、装置及设备 | |
| US9651506B2 (en) | Spectral response effects (SRE) compensation method for photon counting detectors (PCDs) | |
| JP2010223963A (ja) | コンテナを検査する方法及びシステム | |
| CN106659452B (zh) | 在定量单光子发射计算机断层扫描中利用多个光电峰的重构 | |
| CN104050631A (zh) | 一种低剂量ct图像重建方法 | |
| KR102507596B1 (ko) | 저선량 다중 스펙트럼 x 선 단층 촬영을 위한 시스템 및 방법 | |
| Zhao et al. | Iterative beam hardening correction for multi-material objects | |
| Holbrook et al. | Deep learning based spectral distortion correction and decomposition for photon counting CT using calibration provided by an energy integrated detector | |
| JP2010204755A (ja) | 画像処理装置、画像再構成システム、画像処理方法およびプログラム | |
| Steuben et al. | X-ray marching for the computational modeling of tomographic systems applied to materials applications | |
| Hausladen et al. | 3D Tomography and Image Processing Using Fast Neutrons Final Report | |
| Blackston et al. | Progress Update on Iterative Reconstruction of Neutron Tomographic Images, 2019 | |
| Mason et al. | A Monte Carlo framework for low dose CT reconstruction testing | |
| Vijayalakshmi et al. | Comparison of algebraic reconstruction methods in computed tomography | |
| Iskender | X-ray CT scatter correction by a physics-motivated deep neural network | |
| Abir | Iterative CT reconstruction from few projections for the nondestructive post irradiation examination of nuclear fuel assemblies | |
| Bismark | On the application of the polychromatic statistical reconstruction technique to C-arm CT data | |
| Melli | Compressed Sensing Based Reconstruction Algorithm for X-ray Dose Reduction in Synchrotron Source Micro Computed Tomography |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200130 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20201014 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20201117 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210216 |
|
| 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: 20210629 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210719 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6917356 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |