▶ 清華大学の特許一覧 ▶ 同方威視技術股▲分▼有限公司の特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-04-27
(45)【発行日】2022-05-11
(54)【発明の名称】ＣＴ画像を再構成する方法及びデバイス、並びに記憶媒体
(51)【国際特許分類】
G01N 23/046 20180101AFI20220428BHJP
A61B 6/03 20060101ALI20220428BHJP
G06T 1/00 20060101ALI20220428BHJP
G06T 7/00 20170101ALI20220428BHJP
【ＦＩ】
G01N23/046
A61B6/03 350F
G06T1/00 290B
G06T7/00 350C
【請求項の数】
17
(21)【出願番号】P 2018139744
(22)【出願日】2018-07-25
(65)【公開番号】P2019028072
(43)【公開日】2019-02-21
【審査請求日】2020-06-24
(31)【優先権主張番号】201710616651.7
(32)【優先日】2017-07-25
(33)【優先権主張国・地域又は機関】CN
(73)【特許権者】
【識別番号】502192546
【氏名又は名称】清華大学
【氏名又は名称原語表記】Ｔｓｉｎｇｈｕａ Ｕｎｉｖｅｒｓｉｔｙ
【住所又は居所原語表記】Ｔｓｉｎｇｈｕａ Ｕｎｉｖｅｒｓｉｔｙ，Ｈａｉｄｉａｎ Ｄｉｓｔｒｉｃｔ，Ｂｅｉｊｉｎｇ １０００８４，Ｐ．Ｒ．Ｃｈｉｎａ
(73)【特許権者】
【識別番号】503414751
【氏名又は名称】同方威視技術股▲分▼有限公司
(74)【代理人】
【識別番号】110001195
【氏名又は名称】特許業務法人深見特許事務所
(72)【発明者】
【氏名】▲シン▼ 宇 翔
(72)【発明者】
【氏名】張 麗
(72)【発明者】
【氏名】李 薦 民
(72)【発明者】
【氏名】陳 志 強
(72)【発明者】
【氏名】顧 建 平
【審査官】佐藤 仁美
(56)【参考文献】
【文献】国際公開第２０１６／１９２６１２（ＷＯ，Ａ１）
【文献】特開２０１６－１９８５０７（ＪＰ，Ａ）
【文献】特開２００７－１９０１８２（ＪＰ，Ａ）
【文献】Kerstin Hammernik et al.，A Deep Learning Architecture for Limited-Angle Computed Tomography Reconstruction，Bildverarbetung fur die Medizin 2017 (First Online March 2017)，2017年03月01日，Pages 92-97，https://doi.org/10.1007/978-3-662-54345-0_25
(58)【調査した分野】(Int.Cl.，ＤＢ名)
ＩＰＣ
Ａ６１Ｂ ６／００－Ａ６１Ｂ ６／１４、
Ｇ０１Ｎ ２３／００－Ｇ０１Ｎ ２３／２２７６、
Ｇ０６Ｔ １／００－Ｇ０６Ｔ １／４０、
Ｇ０６Ｔ ３／００－Ｇ０６Ｔ ９／４０
(57)【特許請求の範囲】
【請求項1】
被検査物にＣＴ走査を行って、投影データを得るステップと、複数の畳み込み層を含む第１畳み込みニューラルネットワークを利用して投影データを処理し、処理後の投影データを得るステップと、
処理後の投影データに逆投影操作を行って、再構成された画像を得るステップと、
第２畳み込みニューラルネットワークを利用して、再構成された画像を処理し、結果画像を得るステップとを含む、ＣＴ画像を再構成する方法。
【請求項2】
前記ＣＴ走査は、検知器アンダーサンプリング走査、スパース角度走査、内部再構成走査、有限角度走査、および直線軌跡走査のうちの１つであり、前記第１畳み込みニューラルネットワークは、プーリング層のない畳み込みニューラルネットワークである、請求項１に記載の方法。
【請求項3】
前記ＣＴ走査は、円周走査あるいは螺旋走査であり、前記第１畳み込みニューラルネットワークは、該当する畳み込み層の後に設けられた複数のプーリング層と全結合層とをさらに含む、請求項１に記載の方法。
【請求項4】
第１畳み込みニューラルネットワークを利用して投影データを処理する前に、ランプフィルターを用いて投影データにフィルタリングを行うステップをさらに含む、請求項１に記載の方法。
【請求項5】
第２畳み込みニューラルネットワークを利用して、再構成された画像に局所的な平滑処理を行って前記結果画像を得る、請求項１に記載の方法。
【請求項6】
第１畳み込みニューラルネットワークにおける畳み込み層の畳み込みカーネルは、１つの次元が検知器画素シーケンスであり、もう１つの次元が走査角度であり、且つ、第１畳み込みニューラルネットワークにおける畳み込み層の畳み込みカーネルは、検知器画素シーケンスの次元におけるスケールと走査角度の次元におけるスケールとが個別に設定される、請求項１に記載の方法。
【請求項7】
第１畳み込みニューラルネットワークにおける畳み込み層の畳み込みカーネルは、検知器画素シーケンスの次元におけるスケールが、走査角度の次元におけるスケールより大きい、請求項６に記載の方法。
【請求項8】
第１畳み込みニューラルネットワークは、少なくとも３つの畳み込み層を含み、畳み込み層のそれぞれは、畳み込み処理された投影データに非線形演算を行うための１つのアクティブ関数を有する、請求項１に記載の方法。
【請求項9】
第１畳み込みニューラルネットワークは、畳み込み層により処理された投影データに逆投影演算を行うための逆投影層をさらに含む、請求項１に記載の方法。
【請求項10】
第１畳み込みニューラルネットワークにおいて、逆投影層に最も近い畳み込み層の畳み込みカーネルの長さおよび幅のサイズパラーメータは、１＊１である、請求項９に記載の方法。
【請求項11】
第２畳み込みニューラルネットワークは、画像領域において再構成された画像を処理するための画像領域初期畳み込み層及び末端畳み込み層を含む、請求項１に記載の方法。
【請求項12】
画像領域初期畳み込み層に含まれる畳み込み層のそれぞれはアクティブ関数を有し、末端畳み込み層はアクティブ関数を有しない、請求項１１に記載の方法。
【請求項13】
被検査物をＣＴ走査して、投影データを得るＣＴ走査装置と、複数の畳み込み層を含む第１畳み込みニューラルネットワークを利用して投影データを処理し、処理後の投影データを得て、処理後の投影データに逆投影操作を行って再構成された画像を得るように配置されたプロセッサーとを含み、
前記プロセッサーは、さらに、第２畳み込みニューラルネットワークを利用して、再構成された画像を処理し、結果画像を得るように配置される、ＣＴ画像を再構成するデバイス。
【請求項14】
前記ＣＴ走査装置は、検知器アンダーサンプリング走査、スパース角度走査、内部再構成走査、有限角度走査、および直線軌跡走査のうちの１つを実行し、前記第１畳み込みニューラルネットワークは、プーリング層のない畳み込みニューラルネットワークである、請求項１３に記載のデバイス。
【請求項15】
前記ＣＴ走査装置は、円周走査あるいは螺旋走査を実行し、前記第１畳み込みニューラルネットワークは、該当する畳み込み層の後に設けられた複数のプーリング層と全結合層とをさらに含む、請求項１３に記載のデバイス。
【請求項16】
前記プロセッサーは、さらに、第２畳み込みニューラルネットワークを利用して、再構成された画像に局所的な平滑処理を行って前記結果画像を得るように配置される、請求項１３に記載のデバイス。
【請求項17】
プロセッサーにより実行されることで、複数の畳み込み層を含む第１畳み込みニューラルネットワークを利用して、被検査物にＣＴ走査を行って得られた投影データを処理し、処理後の投影データを得るステップと、
処理後の投影データに逆投影操作を行って、再構成された画像を得るステップと、
第２畳み込みニューラルネットワークを利用して、再構成された画像を処理し、結果画像を得るステップとを実現する、コンピュータープログラムを記憶した、コンピューター読み取り可能な媒体。
【発明の詳細な説明】
【技術分野】
【0001】
本開示の実施例は、放射結像に関し、具体的には、ＣＴ画像を再構成する方法及びデバイス、並びに記憶媒体に関する。
【背景技術】
【0002】
Ｘ線ＣＴ結像システムは、医療、保安検査、産業非破壊検出などの分野において広く用いられている。放射線源と検知器は、一定の軌道に沿って一連の投影データを採集し、画像再構成アルゴリズムによる復元により、被検査物の線形減衰係数の三次元空間分布を得ることができる。ＣＴ画像再構成過程は、検知器により採集されたデータから線形減衰係数分布を回復することであって、ＣＴ結像の核心的なステップである。現在、実際の応用には、フィルタ補正逆投影（Ｆｉｌｔｅｒｅｄ Ｂａｃｋ－Ｐｒｏｊｅｃｔｉｏｎ）、ＦＤＫ（Ｆｅｌｄｋｍａｐ－Ｄａｖｉｓ－Ｋｒｅｓｓ）類の解析再構成アルゴリズムと、ＡＲＴ（Ａｌｇｅｂｒａ Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ Ｔｅｃｈｎｉｑｕｅ）、ＭＡＰ（Ｍａｘｉｍｕｍ Ａ Ｐｏｓｔｅｒｉｏｒ）等の反復再構成方法とが主に用いられる。
【0003】
Ｘ線ＣＴ結像の要求がますます多様化されることにつれて、放射量の低減に対する要求もますます高くなり、従来の再構成方法が実現可能な画像の品質はもう限界に達している。新たなＣＴ画像再構成技術の開発が求められている。
【発明の概要】
【発明が解決しようとする課題】
【0004】
従来技術における１つまたは複数の問題に鑑みて、再構成画像の品質を高めることができるＣＴ画像を再構成する方法及びデバイス、並びに記憶媒体を提案する。
【課題を解決するための手段】
【0005】
本開示の一方面によると、被検査物にＣＴ走査を行って、投影データを得るステップと、複数の畳み込み層を含む第１畳み込みニューラルネットワークを利用して、投影データを処理し、処理後の投影データを得るステップと、処理後の投影データに逆投影操作を行って、再構成された画像を得るステップと、を含む、ＣＴ画像を再構成する方法を提案する。
【0006】
本開示のある実施例によると、前記ＣＴ走査は、検知器アンダーサンプリング走査、スパース角度走査、内部再構成走査、有限角度走査、および直線軌跡走査のうちの１つであり、前記第１畳み込みニューラルネットワークは、プーリング層のない畳み込みニューラルネットワークである。
【0007】
本開示のある実施例によると、前記ＣＴ走査は、円周走査あるいは螺旋走査であり、前記第１畳み込みニューラルネットワークは、該当する畳み込み層の後に設けられた複数のプーリング層と全結合層とをさらに含む。
【0008】
本開示のある実施例によると、前記方法は、第２畳み込みニューラルネットワークを利用して、再構成された画像を処理し、結果画像を得るステップをさらに含む。
【0009】
本開示のある実施例によると、前記方法は、第１畳み込みニューラルネットワークを利用して投影データを処理する前に、ランプフィルターを用いて投影データにフィルタリングを行うステップをさらに含む。
【0010】
本開示のある実施例によると、第２畳み込みニューラルネットワークを利用して、再構成された画像に局所的な平滑処理を行って前記結果画像を得る。
【0011】
本開示のある実施例によると、第１畳み込みニューラルネットワークにおける畳み込み層の畳み込みカーネルは、１つの次元が検知器画素シーケンスであり、もう１つの次元が走査角度であり、且つ、第１畳み込みニューラルネットワークにおける畳み込み層の畳み込みカーネルは、検知器画素シーケンスの次元におけるスケールと走査角度の次元におけるスケールとが個別に設定される。
【0012】
本開示のある実施例によると、第１畳み込みニューラルネットワークにおける畳み込み層（例えば、第１畳み込み層）の畳み込みカーネルは、検知器画素シーケンスの次元におけるスケールが、走査角度の次元におけるスケールより大きい。
【0013】
本開示のある実施例によると、第１畳み込みニューラルネットワークは、少なくとも３つの畳み込み層を含み、畳み込み層のそれぞれは畳み込み処理された投影データに非線形演算を行うための１つのアクティブ関数を有する。
【0014】
本開示のある実施例によると、第１畳み込みニューラルネットワークは、畳み込み層により処理された投影データに逆投影演算を行うための逆投影層をさらに含む。
【0015】
本開示のある実施例によると、第１畳み込みニューラルネットワークにおいて、逆投影層に最も近い畳み込み層の畳み込みカーネルの長さおよび幅のサイズパラーメータは、１＊１である。
【0016】
本開示のある実施例によると、第２畳み込みニューラルネットワークは、画像領域において再構成された画像を処理するための画像領域初期畳み込み層及び末端畳み込み層を含む。
【0017】
本開示のある実施例によると、画像領域初期畳み込み層に含まれる畳み込み層のそれぞれはアクティブ関数を有し、末端畳み込み層はアクティブ関数を有しない。
【0018】
本開示の別の一方面によると、被検査物をＣＴ走査して、投影データを得るＣＴ走査装置と、複数の畳み込み層を含む第１畳み込みニューラルネットワークを利用して投影データを処理し、処理後の投影データを得て、処理後の投影データに逆投影操作を行って再構成された画像を得るように配置されたプロセッサーと、を含む、ＣＴ画像を再構成するデバイスを提案する。
【0019】
本開示のある実施例によると、前記ＣＴ走査装置は、検知器アンダーサンプリング走査、スパース角度走査、内部再構成走査、有限角度走査、および直線軌跡走査のうちの１つを実行し、前記第１畳み込みニューラルネットワークは、プーリング層のない畳み込みニューラルネットワークである。
【0020】
本開示のある実施例によると、前記ＣＴ走査装置は、円周走査あるいは螺旋走査を実行し、前記第１畳み込みニューラルネットワークは、該当する畳み込み層の後に設けられた複数のプーリング層と全結合層とをさらに含む。
【0021】
本開示のある実施例によると、前記プロセッサーは、さらに、第２畳み込みニューラルネットワークを利用して、再構成された画像を処理し、結果画像を得るように配置される。
【0022】
本開示のある実施例によると、前記プロセッサーは、さらに、第２畳み込みニューラルネットワークを利用して、再構成された画像に局所的な平滑処理を行って前記結果画像を得るように配置される。
【0023】
本開示のさらに他の一方面によると、プロセッサーにより実行されることで、複数の畳み込み層を含む第１畳み込みニューラルネットワークを利用して、被検査物にＣＴ走査を行って得られた投影データを処理し、処理後の投影データを得るステップと、処理後の投影データに逆投影操作を行って、再構成された画像を得るステップと、を実現するコンピュータープログラムを記憶した、コンピューター読み取り可能な媒体を提案する。
【発明の効果】
【0024】
本開示の上記の実施例の態様を利用して、より高い品質のＣＴ画像を再構成することができる。
【0025】
本発明をよりよく理解するために、以下の図面に基づいて本発明を詳細に説明する。
【図面の簡単な説明】
【0026】
【図1】本開示の一実施例のＣＴ画像を再構成するデバイスの構成模式図を示す。
【図3】本開示の実施例によるデバイスが投影データを取得する正弦図の例を示す。
【図4】本開示の実施例によるデバイスにおいて用いられる畳み込みニューラルネットワークの構成模式図を示す。
【図5】本開示の実施例によるデバイスにおいて用いられる別の畳み込みニューラルネットワークの構成模式図を示す。
【図6A】本開示の実施例によるデバイスにおいて用いられるフィルターカーネルのサイズ模式図を示す。
【図6B】本開示の実施例によるデバイスにおいて用いられるフィルターカーネルのサイズ模式図を示す。
【図6C】本開示の実施例によるデバイスにおいて用いられるフィルターカーネルのサイズ模式図を示す。
【図7】本開示の実施例による方法を記述する模式的フローチャートである。
【図8】本開示の別の一実施例による有限角度ＣＴ走査を実現する走査装置の模式図である。
【図9】本開示のさらに別の一実施例によるスパース角度走査方式を実現する走査装置の模式図である。
【図10】本開示のさらに別の一実施例による内部再構成走査方式を実現する走査装置の模式図である。
【図11】本開示のさらに別の一実施例による検知器アンダーサンプリング走査方式を実現する走査装置の模式図である。
【図13】本開示のさらに別の一実施例による直線軌跡ＣＴ走査を実現する走査装置の模式図を示す。
【発明を実施するための形態】
【0027】
以下、本発明の具体的な実施例について詳細に説明する。注意すべきことは、以下の実施例は例示的なものであり、本発明を限定するものではない。以下の説明においては、本発明の理解を容易にするために、多くの特定的な細部について記述しているが、当業者にとって明らかに分かるように、必ずこれらの特定の細部を採用して本発明を実現することではない。その他の実例においては、本発明を混同させることを避けるために、周知の回路、材料または方法に対する具体的な記述を省略した。
【0028】
明細書の全般にかけて、「一実施例」、「実施例」、「一例」または「例」とは、当該実施例または例に記述された特定的な特徴、構成または特性を組み合わせたものは、本発明の少なくとも一実施例に含まれていることを意味する。従って、明細書の全般にかけて、記載された「一実施例において」、「実施例において」、「一例」または「例」は、必ずしも同一の実施例又は例示を指すものとは限らない。また、いずれの適宜な組合せ及び／又はサブ組合せによって、特定な特徴、構造又は特性を１つ又は複数の実施例又は例示に組み合わせることができる。また、当業者であれば理解できるように、ここで使用する「及び／又は」という専門用語は、一つまたは複数の関連の列記項目の何れかと全ての組み合わせを含む。
【0029】
従来技術の１つまたは複数の問題について、本開示の実施例はＣＴ画像を再構成する方法及びデバイス、並びに記憶媒体を提案する。本開示の実施例によると、１つの畳み込みニューラルネットワークを利用して、ＣＴ走査により得られた投影データ（正弦図として表示され、１つの次元は検知器画素シーケンスであり、もう１つの次元は走査角度である）を処理し、処理後の投影データ（例えば、投影領域における特徴が強化された正弦図）を得る。そして、このように処理して得られた投影データに逆投影を行って、ＣＴ画像を得る。このようにして、画像の品質を向上でき、特に、検知器アンダーサンプリング、スパース角度走査、内部再構成、有限角度走査または直線軌跡走査等のような他の投影データが不完全である場合に、再構成により品質の高いＣＴ画像が得られる。他の実施例においては、別の畳み込みニューラルネットワークを利用して、さらに、ＣＴ画像を処理して最終画像を得てもよい。
【0030】
図１は、本開示の一実施例のＣＴ画像を再構成するデバイスの構成模式図を示す。図１に示すように、本実施例によるＣＴデバイスは、Ｘ線源１０と、機械的運動装置５０と、検知器・データ採集システム２０と、制御・再構成装置６０とを含み、被検査物４０に対してＣＴ走査および再構成を行う。
【0031】
Ｘ線源１０は、例えば、Ｘ線機器であり、結像の分解能によって、適切なＸ線機器焦点サイズを選択する。他の実施例においては、Ｘ線機器を使用せず、直線型加速器等を使用してＸ線ビームを発生してもよい。
【0032】
機械的運動装置５０は、ステージ、フレーム及び制御システム等を含む。ステージは、回転中心の位置を調整するように並進可能である。フレームは、Ｘ線源（Ｘ線機器）１０、検知器及び回転中心の三者を合わせるように並進可能である。本実施例においては、回転ステージ、固定フレームの円周走査軌跡あるいは螺旋軌跡に従って記述する。ステージとフレームとの運動は相対的な運動であるので、ステージを静止させ、フレームを回転させる方式で本実施例を実現する方法を採用してもよい。
【0033】
検知器・データ採集システム２０は、Ｘ線検知器及びデータ採集回路等を含む。Ｘ線検知器には、固体検知器を使用してもよく、気体検知器あるいは他の検知器を使用してもよく、本開示の実施例はこれらに限られない。データ採集回路は、読み出し回路、採集トリガー回路及びデータ伝送回路等を含む。
【0034】
制御・再構成装置６０は、例えば、制御プログラム及び再構成プログラムがインストールされているコンピューターデバイスを含み、機械的回動、電気的制御、安全連鎖制御等を含むＣＴシステム運行過程の制御を担当するとともに、投影データからのＣＴ画像の再構成等を行う。
【0035】
図２は、図１に示すような制御・再構成デバイスの構成模式図を示す。図２に示すように、検知器・データ採集システム２０により採集されたデータは、インターフェースユニット６７及びバス６８を介して記憶デバイス６１に記憶される。ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）６２には、コンピューターデータプロセッサーの配置情報及びプログラムが記憶されている。ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｒｓｓ Ｍｅｍｏｒｙ）６３は、プロセッサー６５の動作中に各種のデータを一時的に記憶するために用いられる。なお、記憶デバイス６１には、データ処理のためのコンピュータープログラム、例えば、ＣＴ画像を再構成する計算プログラムなどが、さらに記憶されている。内部バス６８には、上記の記憶デバイス６１、ＲＯＭ６２、ＲＡＭ６３、入力装置６４、プロセッサー６５、表示デバイス６６及びインターフェースユニット６７が接続されている。
【0036】
ユーザーがキーボードー及びマウスのような入力装置６４により操作命令を入力した後、コンピュータープログラムの指令コードはＣＴ画像を再構成するアルゴリズムを実行するようにプロセッサー６５に命令し、再構成結果を得た後、それをＬＣＤ表示器のような表示デバイス６６に表示させたり、または、プリントのようなハードコピーの形で処理結果を直接的に出力させたりする。
【0037】
本開示の実施例によると、上記のデバイスを利用して被検査物に対してＣＴ走査を行い、投影データを得る。通常、このような投影データは二次元の画像の形で表示できる。図３は本開示の実施例により得られた投影データの例を示す。図３に示すような正弦図の横軸方向は検知器画素シーケンス（例えば、１から２５６まで）を示し、縦軸は角度（例えば、１度から３６０度まで）を示す。そして、デバイスにおけるプロセッサー６５に再構成プログラムを実行させるように制御し、第１畳み込みニューラルネットワークを利用して投影データを処理し、処理後の投影データ得て、さらには、処理後の投影データに逆投影操作を行って、再構成された画像を得る。
【0038】
上記のように、本開示の実施例においては、投影領域において畳み込みニューラルネットワークを利用して投影データを処理した後、逆投影操作を行ってＣＴ画像を再構成する。畳み込みニューラルネットワークは、畳み込み層、プーリング層、および全結合層を含んでもよい。畳み込み層は入力データ集合のキャラクタリゼーションを認識し、畳み込み層のそれぞれは、１つの非線形アクティブ関数演算を有する。プーリング層は、特徴に対する示しを精練し、典型的な操作は平均プーリング及び最大化プーリングを含む。１層または複数層の全結合層は、高次の信号の非線形総合演算を実現し、全結合層も非線形アクティブ関数を有する。よく使われる非線形アクティブ関数には、Ｓｉｇｍｏｉｄ、Ｔａｎｈ、ＲｅＬＵなどがある。
【0039】
図４は、本開示の実施例によるデバイスにおいて用いられる畳み込みニューラルネットワークの構成模式図を示す。図４に示すように、畳み込みニューラルネットワーク（ＣＮＮ）は、異なる層がスタックされてなってもよく、それらは入力データを出力データに変換できる。例えば、ＣＴ走査を行って得られた投影データ（例えば、図４に示す正弦図４１０）を処理して、処理後の投影データ（例えば、図４に示す処理後の正弦図４５０）を得る。
【0040】
図４に示すような畳み込みニューラルネットワークは、複数の畳み込み層、例えば、第１畳み込み層４２０、第２畳み込み層４３０、…、第ｎ＋１畳み込み層４４０を含み、ｎは自然数である。これらの畳み込み層はＣＮＮの核心的な構造セルである。各々の畳み込み層のパラーメータはラーニング可能な畳み込みカーネル（または、単に畳み込みカーネルと称する）の集合により構成され、畳み込みカーネルのそれぞれは、一定の受容野を持つとともに、入力データのデプスの全体にかけて延ばしている。フォーワード過程において、畳み込みカーネルのそれぞれを入力データの幅及び高さに従って畳み込みを行い、畳み込みカーネルの要素と入力データとの点乗積を算出し、当該畳み込みカーネルの二次元的なアクティブマップを生成する。結果として、ネットワークは、入力されたある空間的な位置においてある具体的な種類の特徴を見たときしかアクティブにすることができない畳み込みカーネルをラーニングできる。
【0041】
全ての畳み込みカーネルのアクティブマップをデプス方向においてスタックして、畳み込み層の全部の出力データを形成する。従って、出力データにおける要素のそれぞれは、入力中の小さな領域が見え、同じアクティブマップにおける他の畳み込みカーネルとパラーメータを共用する畳み込みカーネルの出力として解釈できる。
【0042】
【数1】
【0043】
第１畳み込み層に類似する複数の畳み込み層を設けてもよく、各々の畳み込み層の厚みはK（ｎ）とする。
【0044】
第２畳み込み層４３０は第１畳み込み層４２０と似て、第２畳み込み層の厚みはK(2)である。例えば、C(n)を利用して前の畳み込み層により出力された正弦図をさらに処理する。図４に示すように、第２畳み込み層４３０の厚みは３であり、ここで、全ての畳み込みカーネルは未定のネットワークパラーメータであり、これらのネットワークパラーメータは畳み込みニューラルネットワークに対するトレーニングにより得られる。以上で言及された３は例示に過ぎなく、同様に、１つの畳み込み層に対して２４個の畳み込みカーネルを設置してもよく、この場合、第２畳み込み層（ネットワークの第２層）も２４個の正弦図が得られる。ここで、アクティブ関数としてＲｅＬｕを選択する。別の一実施例において、第２畳み込み層において、形状／サイズパラーメータが１８＊１＊１である畳み込みカーネルを用いて、厚み方向において畳み込みを行って第２層の出力を得て、アクティブ関数としてＲｅＬｕを使用する。
【0045】
【数2】
【0046】
図５は、本開示の実施例によるデバイスにおいて用いられる別の畳み込みニューラルネットワークの構成模式図を示す。図５に示すように、入力された投影データ５１０、畳み込み層５２０、畳み込み層５３０及び畳み込み層５４０は、図４に示された入力された正弦図４１０、第１畳み込み層４２０、第２畳み込み層４３０及び第ｎ＋１畳み込み層４４０と似ているので、ここで説明を省略する。図５に示す構成は、図４の構成と比べて、画像領域初期畳み込み層（図５に示す例示は１つの畳み込み層５７０であるが、本発明はこれに限られない）と末端畳み込み層（畳み込み層５８０）とを含む第２畳み込みニューラルネットワークをさらに設け、画像領域において再構成されたＣＴ画像に対し例えば局所的な平滑処理を行って結果画像を出力する点で異なる。
【0047】
【数3】
【0048】
上記の構成においては、プーリング層を含んでいないが、プーリング層を設けてもよく、例えば、ＣＴ走査により完全なデータを得た場合に（例えば、全検知器３６０円周走査）、１つまたは複数の畳み込み層の後にプーリング層を設けてもよい。しかし、ＣＴ走査により非完全なデータが得られた場合には（例えば、検知器アンダーサンプリング、スパース角度サンプリング、有限角度、内部再構成、または、直線軌跡等の方式）、プーリング層を設けなくてもよい。
【0049】
【数4】
【0050】
具体的には、以上のように畳み込みニューラルネットワークを構築した後、ネットワークパラーメータのトレーニングを行う。例えば、基本的な走査対象の数学的モデルを構築して、実際のシステムに従ってモデリングを行って、ＣＴシミュレーションデータを生成する。また、本開示の実施例によると、ネットワークの入力として複数の走査対象のＣＴシミュレーションデータを使用し、フラグとして走査対象の画像真値を使用して、ネットワークパラーメータをトレーニングしてもよい。他の実施例においても、実際のシステムにおいて、対象に走査を行って、ＣＴ走査データを取得し、このネットワークに入力して再構成テストを行う。テスト結果について、例えば、既知の局所的な平滑領域について、局所的平滑を行うような指向性のある画像処理を行う。処理後の画像をフラグとして、さらに、ネットワークをトレーニングし、ネットワークパラーメータの微細調整を実現する。
【0051】
図７は、本開示の実施例による方法を記述する模式的フローチャートである。図７に示すように、ステップＳ７１０において、被検査物をＣＴ走査して、投影データを得る。ここで、ＣＴ走査は、シングルエネルギーのものであってもよく、マルチエネルギーのものであってもよく、本開示の実施例はこれらに限られない。
【0052】
ステップＳ７２０において、ランプフィルターを用いて投影データにフィルタリングを行う。例えば、ＲＬランプフィルターを用いて投影データにフィルタリングを行う。ここで、別のフィルターを用いてもよく、または、フィルタリングを行わなくてもよいことは、当業者にとって明らかである。
【0053】
ステップＳ７３０において、第１畳み込みニューラルネットワーク４２０／４３０／４４０を利用して、投影領域において投影データを処理し、処理後の投影データを得る。例えば、図４に示すように、トレーニングにより得られた畳み込みニューラルネットワークを利用して投影データを処理し、処理後の正弦図を得る。
【0054】
ステップＳ７４０において、処理後の投影データに逆投影操作を行い、例えば、図４に示された逆投影層４６０において逆投影操作を行って、再構成されたＣＴ画像を得る。
【0055】
上記のように、別の実施例として、ＣＴ画像を得た後、さらに後処理を行ってもよく、例えば、ステップＳ７５０において、第２畳み込みニューラルネットワーク５７０／５８０を利用して、再構成されたＣＴ画像を処理して、結果画像を得てもよい。例えば、ここで、再構成された画像に対して、局所的な平滑処理、または、分割、エッジ強化及び均等化等のような他の画像処理操作を行ってもよい。
【0056】
以上、主に３６０度円周走査により完全な投影データを得た場合について記述したが、上記の態様が非完全な投影データの場合、例えば、検知器アンダーサンプリング、スパース角度サンプリング、有限角度、内部再構成、或者直線軌跡走査等方式にも適用できることは、当業者にとって明らかである。
【0057】
図８は、本開示の別の一実施例による有限角度ＣＴ走査を実現する走査装置の模式図である。図８に示すように、放射線源１０から発されたＸ線は視野４５における被検査物４０を透過した後、検知器３０により受けられて電気信号に転換され、さらには、減衰値を表すデジタル信号に転換され、投影データとして、コンピューターにおいて再構成される。上記の態様によると、被検査物４０に対して有限角度のＣＴ走査（例えば、１３０度）を行っても、品質の高い画像を再構成することができる。
【0058】
図９は、本開示のさらに別の一実施例によるスパース角度走査方式を実現する走査装置の模式図である。図９に示すように、放射線源１０から発されたＸ線は視野４５における被検査物４０を透過した後、検知器３０により受けられて電気信号に転換され、さらには、減衰値を表すデジタル信号に転換され、投影データとして、コンピューターにおいて再構成される。上記の態様によると、被検査物４０に対して幾つかの回転位置におけるＣＴ走査（例えば、６個所）を行っても、品質の高い画像を再構成することができる。このようにして、被検査物に対してスパース角度ＣＴ走査を行っても、非完全な投影データから品質の高い画像を再構成することができる。
【0059】
図１０は、本開示のさらに別の一実施例による内部再構成走査方式を実現する走査装置の模式図である。図１０に示すように、放射線源１０から発されたＸ線は視野４５における被検査物４０の一部を透過した後、検知器３０により受けられて電気信号に転換され、さらには、減衰値を表すデジタル信号に転換され、投影データとして、コンピューターにおいて再構成される。上記の態様によると、被検査物４０に対して内部再構成ＣＴ走査を行っても、品質の高い画像を再構成することができる。
【0060】
図１１は、本開示のさらに別の一実施例による検知器アンダーサンプリング走査方式を実現する走査装置の模式図である。図１１に示すように、放射線源１０から発されたＸ線は視野４５における被検査物４０を透過した後、検知器３０により受けられて電気信号に転換され、さらには、減衰値を表すデジタル信号に転換され、投影データとして、コンピューターにおいて再構成される。この例において、検知器３０は、アンダーサンプリングの状態に設置されて、例えば、各々の検知器ユニットを予定の距離に離間してアンダーサンプリングを実現する。このようにして、上記の態様によると、被検査物４０に対して検知器アンダーサンプリングＣＴ走査を行っても、品質の高い画像を再構成することができる。
【0061】
図１２は、図８、図９、図１０および図１１に示すような走査方式に係る不完全な投影データを記述する模式図である。図１２に示すように、角度スパースサンプリングＣＴ走査、有限角度ＣＴ走査、検知器アンダーサンプリングＣＴ走査及び内部再構成ＣＴ走査を行って得られた投影データはいずれも不完全なものである。投影データが不完全にもかかわらず、上記の態様によると、これらの不完全な投影データから品質の高い画像を再構成することができる。
【0062】
以上、角度スパースサンプリング走査等の方式について記述したが、本開示の方法が直線軌跡ＣＴ走査システムにおいても同様に適用できることは、当業者にとって明らかである。図１３は、本開示のさらに別の一実施例による直線軌跡ＣＴ走査を実現する走査装置の模式図を示す。
【0063】
図１３に示すように、从射線源１０から発されたＸ線は視野における被検査物４０を透過した後、検知器３０により受けられて電気信号に転換され、さらには、減衰値を表すデジタル信号に転換され、投影データとして、コンピューターにおいて再構成される。この例において、被検査物４０は、検知器に平行なコンベアにおいて直線の軌跡に沿って運動する。検知器は水平方向において放射線源との開き角をできるだけ大きくし、鉛直方向において物体を覆う。例えば、検知器アレイは放射線源に対向する側に置かれて、放射線の水平の開き角θは９０度以上と求められ、直線軌跡ＣＴ走査の投影データを得る。上記の態様によると、被検査物４０に対して直線軌跡ＣＴ走査を行っても、品質の高い画像を再構成することができる。
【0064】
本開示の実施例は、畳み込みニューラルネットワークに基づくＸ線ＣＴ再構成方法を提案しており、データ情報を深く発掘して、畳み込みニューラルネットワークとシステム向けのパラーメータを形成し、効率的なＣＴ画像再構成方法を得られる。
【産業上の利用可能性】
【0065】
本開示の方法は、異なるＣＴ走査モード及びシステム構成に柔軟的に適用でき、医学診断、産業非破壊検出及び保安検査分野に適用できる。
【0066】
以上の詳細な記述においては、模式図、フローチャート及び／又は例によりＣＴ画像を再構成する方法とデバイスの多数の実施例を記述した。このような模式図、フローチャート及び／又は例が１つまたは複数の機能及び／又は操作を含む場合、当業者は、このような模式図、フローチャートまたは例のそれぞれの機能及び／又は操作が各種の構成、ハードウェア、ソフトウェア、ファームウェアまたは実質的にこれらの任意的な組み合わせにより単独及び／又はともに実現できることを理解すべきである。一実施例において、本発明の実施例に記載された主題の幾つかの部分はＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）、ＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）、または他の集積回路により実現できる。しかしながら、当業者は、ここで開示された実施例のある方面の全体または一部が等価的に集積回路において実現でき、一台または複数の台のコンピューターにおいて実行される１つまたは複数のコンピュータープログラム（例えば、一台または複数の台のコンピューターにおいて実行される１つまたは複数のプログラムとして実現でき）として実現でき、１つまたは複数のプロセッサーにおいて実行される１つまたは複数のプログラム（例えば、１つまたは複数のマイクロプロセッサーにおいて実行される１つまたは複数のプログラムとして実現でき）として実現でき、ファームウェアとして実現でき、あるいは実質的に上記形態の任意的な組み合わせとして実現できる、ことを理解すべきであり、且つ、当業者は、本開示から電路の設計及び／又はソフトウェアー及び／又はファームウェアコードの書き込み能力を持つ。また、当業者は、本開示に記載された主題のメカニズムが複数の形態のプログラム製品として配布でき、且つ、実際に配布を実行するための信号キャリア媒体がどのような種類であっても、本開示に記載された主題の例示的な実施例を適用できる、ことを理解すべきである。信号キャリア媒体の例としては、ＦＤ、ＨＤＤ、ＣＤ、ＤＶＤ、デジタルテープ、コンピューターメモリなどの記録可能な媒体、及びデジタル及び／又はアナログ通信媒体（例えば、光ファイバーまたは光ケーブル、導波路、有線通信回線、無線通信線路など）のような伝送型媒体を含むことができるが、これに限られない。
【0067】
幾つかの典型的な実施例を参照して本発明を記述したが、使用された用語は説明及び例示的なものであり、制限的なものではないことに注意すべきである。本発明は発明の趣旨又は実質を逸脱することなく、いろんな形態で具体的に実施できるので、上記の実施例は前述の細部に限られなく、添付の請求の範囲に限定された精神及び範囲において広義的に解釈すべきである。従って、請求の範囲又はそれと均等な範囲内の全ての変化及び改良は、添付の請求の範囲に含まれるべきである。
【符号の説明】
【0068】
１０ Ｘ線源、２０ 検知器・データ採集システム、４０ 被検査物、５０ 機械的運動装置、６０ 制御・再構成装置。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6A】
【図6B】
【図6C】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
