特開2023-133213磁気共鳴イメージング方法、磁気共鳴イメージング装置及びプログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023133213
(43)【公開日】2023-09-22
(54)【発明の名称】磁気共鳴イメージング方法、磁気共鳴イメージング装置及びプログラム
(51)【国際特許分類】
A61B 5/055 20060101AFI20230914BHJP
【ＦＩ】
A61B5/055 376
【審査請求】未請求
【請求項の数】9
【出願形態】ＯＬ
(21)【出願番号】P 2023034352
(22)【出願日】2023-03-07
(31)【優先権主張番号】17/690,834
(32)【優先日】2022-03-09
(33)【優先権主張国・地域又は機関】US
(71)【出願人】
【識別番号】594164542
【氏名又は名称】キヤノンメディカルシステムズ株式会社
(74)【代理人】
【識別番号】110001771
【氏名又は名称】弁理士法人虎ノ門知的財産事務所
(72)【発明者】
【氏名】サミア デーヴ シャルマ
【テーマコード（参考）】
4C096
【Ｆターム（参考）】
4C096AB07
4C096AB11
4C096AD06
4C096AD12
4C096AD24
4C096BA15
4C096DA11
4C096DA30
4C096DC21
(57)【要約】
【課題】アーチファクトの低減を基にしたＭＲＩにおいて改善されたノイズ低減を実現すること。
【解決手段】実施形態に係る磁気共鳴イメージング方法は、アーチファクトの低減を基にした磁気共鳴イメージングを実行する磁気共鳴イメージング方法であって、スキャン領域のスライス方向に沿って複数のスライスを収集するステップと、収集された複数のスライスにおける信号分散の空間的範囲を推定するステップと、推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成するステップとを含む。
【選択図】図２
【解決手段】実施形態に係る磁気共鳴イメージング方法は、アーチファクトの低減を基にした磁気共鳴イメージングを実行する磁気共鳴イメージング方法であって、スキャン領域のスライス方向に沿って複数のスライスを収集するステップと、収集された複数のスライスにおける信号分散の空間的範囲を推定するステップと、推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成するステップとを含む。
【選択図】図２
【特許請求の範囲】
【請求項1】
アーチファクトの低減を基にした磁気共鳴イメージングを実行する磁気共鳴イメージング方法であって、
スキャン領域のスライス方向に沿って複数のスライスを収集するステップと、
収集された複数のスライスにおける信号分散の空間的範囲を推定するステップと、
推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成するステップと
を含む、磁気共鳴イメージング方法。
スキャン領域のスライス方向に沿って複数のスライスを収集するステップと、
収集された複数のスライスにおける信号分散の空間的範囲を推定するステップと、
推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成するステップと
を含む、磁気共鳴イメージング方法。
【請求項2】
前記収集された複数のスライスそれぞれを２次元スライスから３次元空間ボリュームスライスに空間的にエンコードすることで、一連の３次元空間ボリュームスライスを生成するステップをさらに含む、
請求項１に記載の磁気共鳴イメージング方法。
請求項１に記載の磁気共鳴イメージング方法。
【請求項3】
前記再構成画像を生成するステップは、二乗和平方根又は複素数値和を計算することによって、前記収集された複数のスライスの信号を結合することを含む、
請求項１に記載の磁気共鳴イメージング方法。
請求項１に記載の磁気共鳴イメージング方法。
【請求項4】
前記信号分散の空間的範囲を推定するステップは、双極子応答関数を用いて、前記信号分散の空間的範囲を推定することを含む、
請求項１に記載の磁気共鳴イメージング方法。
請求項１に記載の磁気共鳴イメージング方法。
【請求項5】
前記収集された複数のスライスから、最高の信号強度を有するスライスを特定するステップと、
前記最高の信号強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定するステップと
をさらに含み、
前記再構成画像を生成するステップは、前記最高の信号強度を有するスライスの信号と、前記少なくとも１つの近傍スライスの信号とを結合することを含む、
請求項１に記載の磁気共鳴イメージング方法。
前記最高の信号強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定するステップと
をさらに含み、
前記再構成画像を生成するステップは、前記最高の信号強度を有するスライスの信号と、前記少なくとも１つの近傍スライスの信号とを結合することを含む、
請求項１に記載の磁気共鳴イメージング方法。
【請求項6】
前記少なくとも１つの近傍スライスを特定するステップは、双極子応答関数を用いて、前記最高の信号強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定することを含む、
請求項５に記載の磁気共鳴イメージング方法。
請求項５に記載の磁気共鳴イメージング方法。
【請求項7】
前記アーチファクトは、メタルアーチファクトである、
請求項１～６のいずれか１つに記載の磁気共鳴イメージング方法。
請求項１～６のいずれか１つに記載の磁気共鳴イメージング方法。
【請求項8】
アーチファクトの低減を基にした磁気共鳴イメージングを実行する磁気共鳴イメージング装置であって、
スキャン領域のスライス方向に沿って複数のスライスを収集し、
収集された複数のスライスにおける信号分散の空間的範囲を推定し、
推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成する
処理装置を備える、磁気共鳴イメージング装置。
スキャン領域のスライス方向に沿って複数のスライスを収集し、
収集された複数のスライスにおける信号分散の空間的範囲を推定し、
推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成する
処理装置を備える、磁気共鳴イメージング装置。
【請求項9】
アーチファクトの低減を基にした磁気共鳴イメージングを実行するプログラムであって、
スキャン領域のスライス方向に沿って複数のスライスを収集する手順と、
収集された複数のスライスにおける信号分散の空間的範囲を推定する手順と、
推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成する手順と
をコンピュータに実行させる、プログラム。
スキャン領域のスライス方向に沿って複数のスライスを収集する手順と、
収集された複数のスライスにおける信号分散の空間的範囲を推定する手順と、
推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成する手順と
をコンピュータに実行させる、プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本明細書及び図面に開示の実施形態は、磁気共鳴イメージング方法、磁気共鳴イメージング装置及びプログラムに関する。
【0002】
具体的には、本開示は、メタルアーチファクトの低減を基にした磁気共鳴イメージングにおいて改善されたノイズ低減を実現するための方法、装置及び非一時的コンピュータ可読記憶媒体に関する。
【背景技術】
【0003】
従来、医用イメージング装置は、患者の体内の臓器及び組織の画像を生成する。例えば、磁気共鳴イメージング（Magnetic Resonance Imaging：ＭＲＩ）は、電波、磁場及び傾斜磁場を使用して、内部臓器及び組織の画像を生成する。これらの画像が生成されると、医師は、当該画像を患者のけが又は病気を診断するために使用することができる。
【0004】
ここで、患者の体内に存在する金属インプラント又は金属の物体は、ＭＲＩスキャン中にいくつかの課題を生じさせる。１つの大きな課題は、ＭＲＩスキャンによって収集されるデータからメタルアーチファクトを補正、又は、少なくとも低減することである。ＭＲＩスキャンによって収集されるデータからメタルアーチファクトを低減するために利用されるメタルアーチファクト補正方法は、金属インプラントの強い磁化率によって発生する信号の歪みを解消するために、異なる軸のデータを収集する必要がある。このようなメタルアーチファクト補正方法は、異なる軸に沿った複数の領域を考慮して、信号の歪みを解消する。しかしながら、そのような領域の多くは、ノイズだけを含んでいることが多い。したがって、既知のメタルアーチファクト補正方法は、メタルアーチファクトを低減又は補正するための効果的な手法を提供するものではない。
【0005】
図１は、金属インプラントの強い磁化率によって生じ得る信号の歪みにより良好に対処するために、ＭＲＩのプロセスの一部として、異なる軸に沿ってデータが収集されることを示している。メタルアーチファクト補正方法の一例であるＳＥＭＡＣ（Slice Encoding for Metal Artifact Correction）法では、データ収集のプロセスにおいて、２次元（２Ｄ）スライスが励起され、その後、当該２Ｄスライスが３次元（３Ｄ）ボリュームにエンコードされる。ＳＥＭＡＣ法におけるデータ収集の結果は、（ｘ，ｙ，ｚ，スライス）次元の４次元（４Ｄ）ボリューム１００である。図１は、Ｚ軸１０４及びスライス軸１０２を示している。なお、図１にはｘ次元及びｙ次元は図示されていないが、当業者には、４Ｄボリュームの軸系ではｘ軸及びｙ軸が図１に入ることが理解される。加えて、簡略化のために、コイルの次元が省略されている。
【0006】
さらに、ＳＥＭＡＣ法では、データ結合によって、再構成の処理が実行される。具体的には、ＳＥＭＡＣ法では、データ結合のプロセスの一部として、（ｘ，ｙ）の位置ごとに、（ｚ，スライス）次元における線１０６に沿って全てのデータが結合される。このＳＥＭＡＣ法におけるデータ結合は、（ｚ，スライス）次元における線１０６に沿って二乗和平方根（Root-Sum-of-Squares：ＲＳＳ）を計算することによって、簡単に行われる。このような簡単なＲＳＳによるデータ結合は、線１０６上の多くのスライスが大抵はノイズを含んでいることから、ノイズを生じさせる。したがって、主な技術的課題は、ノイズを最小化しながらも信号を維持するように、画像データを結合する方法を開発することである。
【0007】
なお、上述した「背景技術」の説明は、本開示の背景の概要を示すためのものである。当該背景技術の部分に記載されている範囲内において、発明者の研究は、出願時に先行技術として認められない下記説明の態様と同様に、明示的にも暗示的にも、本発明に対する先行技術として認められるものではない。
【先行技術文献】
【非特許文献】
【0008】
【非特許文献1】Koch KM，Lorbiecki JE，Hinks RS，King KF，”A multispectral three-dimensional acquisition technique for imaging near metal implants”，Magnetic Resonance in Medicine，An Official Journal of the International Society for Magnetic Resonance in Medicine，2009 Feb，61(2)，381-90
【非特許文献2】Lu W，Pauly KB，Gold GE，Pauly JM，Hargreaves BA，”SEMAC: slice encoding for metal artifact correction in MRI”，Magnetic Resonance in Medicine，An Official Journal of the International Society for Magnetic Resonance in Medicine”2009 Jul，62(1)，66-76
【非特許文献3】Lu W，Pauly KB，Gold GE，Pauly JM，Hargreaves BA，”Slice encoding for metal artifact correction with noise reduction”，Magnetic resonance in medicine，2011 May，65(5)，1352-7
【発明の概要】
【発明が解決しようとする課題】
【0009】
本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、アーチファクトの低減を基にしたＭＲＩにおいて改善されたノイズ低減を実現することである。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置付けることもできる。
【課題を解決するための手段】
【0010】
実施形態に係るＭＲＩ方法は、アーチファクトの低減を基にしたＭＲＩを実行するＭＲＩ方法であって、スキャン領域のスライス方向に沿って複数のスライスを収集するステップと、収集された複数のスライスにおける信号分散の空間的範囲を推定するステップと、推定された信号分散の空間的範囲に基づいて、前記収集された複数のスライスの信号を結合することによって、アーチファクトが補正された前記スキャン領域の再構成画像を生成するステップとを含む。
【図面の簡単な説明】
【0011】
添付図面とともに以下の詳細な説明を検討し、当該説明を参照して本開示をより良く理解するにつれて、本開示の完全な理解及び本開示に付随する多くの利点を容易に把握することができるであろう。
【0012】
【発明を実施するための形態】
【0013】
本開示は、メタルアーチファクトに起因するノイズの低減を基にした画質改善のための装置、方法及び非一時的コンピュータ可読記憶媒体に関する。
【0014】
実施形態によれば、本開示は、磁気共鳴イメージング（ＭＲＩ）データにおけるアーチファクト（例えば、メタルアーチファクト）を補正する方法であって、方法は、身体部分に関連するスキャン領域のスライス方向に沿って複数のスライスを取得することと、スライス方向に沿って取得された複数のスライスの信号分散の空間的範囲を推定することと、信号分散の推定された空間的範囲に基づいて、スライス方向に沿って取得された複数のスライスの信号を結合して、スキャン領域の再構成された画像を生成することと、を含むが、これに限定されない、方法に関する。
【0015】
実施形態によれば、本開示は、磁気共鳴イメージング（ＭＲＩ）データにおけるアーチファクト（例えば、メタルアーチファクト）を補正する装置であって、装置は、身体部分に関連するスキャン領域のスライス方向に沿って複数のスライスを取得し、スライス方向に沿って取得された複数のスライスの信号分散の空間的範囲を推定し、信号分散の推定された空間的範囲に基づいて、スライス方向に沿って取得された複数のスライスの信号を結合して、スキャン領域の再構成された画像を生成する、ように構成された処理回路、を含むが、これに限定されない、装置に関する。
【0016】
実施形態によれば、本開示は、コンピュータによって実行されると、磁気共鳴イメージング（ＭＲＩ）データにおけるアーチファクト（例えば、メタルアーチファクト）を補正する方法をコンピュータに実行させる、コンピュータ可読命令を格納する非一時的コンピュータ可読記憶媒体であって、方法は、身体部分に関連するスキャン領域のスライス方向に沿って複数のスライスを取得することと、スライス方向に沿って取得された複数のスライスの信号分散の空間的範囲を推定することと、信号分散の推定された空間的範囲に基づいて、スライス方向に沿って取得された複数のスライスの信号を結合して、スキャン領域の再構成された画像を生成することと、を含むが、これに限定されない、非一時的コンピュータ可読記憶媒体に関する。
【0017】
なお、上記段落は、総論を示すものであり、添付の特許請求の範囲を限定することを意図したものではない。記載されている実施形態は、以下の詳細な説明を添付の図面とともに参照することによって、さらなる利点とともに最も良く理解されるであろう。
【0018】
また、以下で用いられている「１つの」の用語は、１つ又は１つ以上として定義される。また、以下で用いられている「複数の」の用語は、２つ又は２つ以上として定義される。また、以下で用いられている「他の」の用語は、少なくとも２つ目以上として定義される。また、以下で用いられている「含む」及び／又は「有する」の用語は、構成すること（つまり、オープンランゲージ）として定義される。また、本文書を通して、「一実施形態」、「特定の実施態様」、「実施形態」、「実装」、「例」又は同様の用語を用いた記載は、本実施形態に関連して記載されている特定の特性、構造又は特徴が、本開示の少なくとも１つの実施形態に含まれることを意味する。したがって、このような語彙表現又は本明細書を通して各所に記載されている語彙表現は、必ずしも全てが同一の実施形態に関するとは限らない。さらに、当該特定の特性、構造又は特徴は、１つ以上の実施形態において、限定されることなく任意の適切な方法で組み合わされてもよい。
【0019】
ここで、図２を参照すると、図２は、メタルアーチファクトの低減を基にしたＭＲＩにおいて改善されたノイズ低減を実現するための方法を示している。この方法では、図１と同様に、データ収集のプロセスにおいて、２次元（２Ｄ）スライスが励起され、その後、当該２Ｄスライスが３次元（３Ｄ）ボリュームにエンコードされる。データ収集の結果は、（ｘ，ｙ，ｚ，スライス）次元の４次元（４Ｄ）ボリューム１００である。図２は、Ｚ軸２０４及びスライス軸２０２を示している。図２に示すように、この方法では、収集されたスライスが解析されることによって、あるエリアの（ｚ軸方向における）特定の部分を共有する全てのスライスではなく、局所的な数のスライスに補正を集中させる。具体的には、この方法では、データ結合のプロセスの一部として、最初に、中心スライス２１２が特定される。例えば、中心スライス２１２は、最高の信号強度を有するスライスであってもよい（即ち、当該中心スライス２１２は、患者の体内における金属インプラントがありそうな位置を示す）。次に、双極子応答関数（dipole response function)の半値全幅(Full Width at Half Maximum：ＦＷＨＭ）を利用することによって、近傍サイズ２０８が決定される（後に図４を参照して説明する）。近傍サイズ２０８が決定されると、当該近傍サイズ２０８に含まれるスライスが特定される。この例では、近傍スライス２０６が特定される。したがって、（ｘ，ｙ）の位置ごとに、（ｚ，スライス）次元における線２１０に沿った全てのデータが結合される。このデータ結合は、（ｚ，スライス）次元における線２１０に沿ってＲＳＳを計算することによって、簡単に行われる。図示されているように、最高の信号強度を有する中心スライス２１２及びその近傍スライス２０６のみが、データが結合される部分として考慮される。このように、中心スライス２１２及び近傍スライス２０６のみを考慮することによって、スライス軸２０２に沿った、患者の体内における金属インプラントの位置を提供しない他のスライスは利用されないようになる。これにより、メタルアーチファクトの低減を行うためのより良好な推定結果が得られる。
【0020】
この近傍結合の技術は、金属誘起磁化率の影響を受けるボクセルの信号分散が限定的であることに基づいている。この技術の利点は、当該信号分散が（ランダムに分布するのではなく、）密接にクラスタ化されることである（後に図３及び図４を参照して説明する）。これは、双極子応答関数が有する局所的な性質の結果である。双極子応答関数は、高磁化率源（例えば、金属）がＢ０磁場にどのような影響を与えるかを数学的に表す。高磁化率物体は、信号分散を引き起こすが、その分散は、（ｚ，スライス）次元における線に沿って局所的である。また、双極子応答関数に関する情報を使用して、考慮されるべきスライスの範囲が決定されてもよい。例えば、双極子応答関数のＦＷＨＭが計算され、データ収集のスライス厚を使用して、当該ＦＷＨＭ内にあるスライスの点がデータ結合のために保持されてもよい。
【0021】
具体的には、図３は、水／組織を含む参照領域のボクセル値に基づいて既知の媒体における参照領域の信号分散を説明するために、制御された環境下で生成されたＭＲＩスキャン参照画像３０２を示している。当該既知の媒体は、水又は患者の臓器の組織であり得る。例として、当該参照領域は水の媒体３０４中の点３０８であり、当該点３０８の信号分散がグラフ表示３１０によって表されている。グラフ表示３１０は、ｘ軸及びｙ軸にわたる点３０８の信号分散を表している。ここで、ｘ軸は、スライス軸のパラメータを表し、ｙ軸は、信号強度のパラメータを表している。結果的に、グラフ表示３１０は、水という既知の媒体では、信号分散がほとんどクラスタ化され、局所化されていることを示している。
【0022】
また、図４は、金属の媒体を含む別の参照領域のボクセル値に基づいて既知の媒体における参照領域の信号分散を説明するために、制御された環境下で生成されたＭＲＩスキャン参照画像４０２を示している。当該既知の媒体は、患者の体内における金属インプラントであり得る。例として、当該参照領域は金属インプラント４０４中の点４０８であり、当該点４０８の信号分散がグラフ表現４１０によって表されている。結果的に、グラフ表現４１０は、双極子応答関数が、高磁化率媒体領域（例えば、金属インプラント領域）がＢ０磁場にどのような影響を与えるかを数学的に表すことを示している。このように、金属インプラント領域（高磁化率物体とも称される）は、信号分散を引き起こすが、その分散は、（ｚ，スライス）次元における線に沿って局所的である。このグラフ表現４１０の信号分散は、ランダムに分布しているわけではない。信号分散のクラスタ化は、既知の双極子応答関数を用いたＭＲ物理特性と一致している。また、双極子応答関数に関する情報を使用して、考慮されるべきスライスの範囲が決定されてもよい。例えば、双極子応答関数のＦＷＨＭが計算され、データ収集のスライス厚を使用して、当該ＦＷＨＭ内にあるスライスの点がデータ結合のために保持されてもよい（後に図５を参照して詳細に説明する）。図２に戻り、ＦＷＨＭ及びデータ収集のスライス厚を利用して、近傍サイズ２０８が計算される。
【0023】
【0024】
まず、ステップ５０２において、図２を参照して図示及び上述したように、スキャン領域のスライス方向（スライス軸２０２）に沿って複数のスライスが収集される。
【0025】
次に、ステップ５０４において、データ結合のプロセスの一部として、収集された複数のスライスにおける信号分散の空間的範囲が推定される。具体的には、データ結合のプロセスの一部として、最初に、中心スライス２１２が特定される。例えば、中心スライス２１２は、最高の信号強度を有するスライスであってもよい。即ち、中心スライス２１２は、患者の体内における金属インプラントの位置を示す。次に、（例えば、双極子応答関数のＦＷＨＭを利用することによって、）近傍サイズ２０８（信号分散の空間的範囲とも称される）が推定される。近傍サイズに含まれるスライスの数は、１つであってもよいし、スライスの全数より少ない数であってもよい。３ｍｍのスライス厚を使用する臨床スキャンのための好ましい実施形態では、中心スライスの両隣それぞれの方向に２つのスライスを含む近傍サイズが利用される。
【0026】
次に、ステップ５０６において、推定された空間的範囲に基づいて、スライス方向に沿って収集された複数のスライスの信号が結合されることによって、メタルアーチファクトが補正されたスキャン領域の再構成画像が生成される。具体的には、近傍サイズ２０８を決定すると、近傍サイズ２０８に含まれるスライスが特定される。この例では、近傍スライス２０６が特定される。したがって、各（ｘ，ｙ）位置について、（ｚ，スライス）次元における線２１０に沿った全てのデータが結合される。このデータ結合は、（ｚ，スライス）次元における線２１０に沿って二乗和平方根（ＲＳＳ）を計算することによって、簡単に行うことができる。または、このデータ結合は、（ｚ，スライス）次元における線２１０に沿って複素数値和を計算することによって行われてもよい。したがって、この技術では、最高の信号強度を有する中心スライス２１２及びその近傍スライス２０６のみが、データ結合のプロセスの一部として考慮される。このように、中心スライス２１２及び近傍スライス２０６のみを考慮することによって、スライス軸２０２に沿った、患者の体内における金属インプラントの位置を提供しない他のスライスは利用されないようになる。これにより、メタルアーチファクトの低減を行うためのより良好な推定結果が得られる。
【0027】
ここで、図１及び図２を参照しつつ、図６Ａについて説明する。具体的には、図６Ａは、（１）スキャン画像に対して線１０６に沿って図１のノイズ低減の方法を適用して生成された画像６０２、及び、（２）当該スキャン画像に対して線２１０に沿って図２のノイズ低減の方法を適用して生成された画像６０４を示している。また、画像６０６は、画像６０２と画像６０４との違いを説明する目的で、理想的な状態（ノイズの無い状態）で生成された参照画像である。
【0028】
画像６０２は、図１のメタルアーチファクト補正方法を適用して生成された画像である。また、領域６１０は、金属インプラントを示しており、領域６０８は、スキャン画像の水／組織領域を示している。画像６０２は、さらに、四角で囲まれた関心領域６１２を示しており、当該関心領域は、図６Ｂに詳細に示されている。
【0029】
画像６０４は、図２のメタルアーチファクト補正方法を適用して生成された画像である。また、領域６１４は、金属インプラントを示しており、領域６１６は、スキャン画像の水／組織領域を示している。画像６０４は、さらに、四角で囲まれた関心領域６１８を示しており、当該関心領域は、図６Ｂに詳細に示されている。
【0030】
画像６０６は、前述したように、参照画像である。また、領域６２０は、金属インプラントを示しており、領域６２２は、スキャン画像の水／組織領域を示している。画像６０６は、さらに、四角で囲まれた関心領域６２４を示しており、当該関心領域は、図６Ｂに詳細に示されている。
【0031】
ここで、図６Ｂは、図６Ａの関心領域６１２、６１８及び６２４の拡大画像を示している。図６Ｂから分かるように、関心領域６１２と関心領域６１８とを比較すると、関心領域６１８は、ノイズが大幅に少ない。また、関心領域６１８と参照画像の関心領域６２４とを比較すると、関心領域６１８と参照画像の関心領域６２４との間のノイズの違いが小さいことが分かる。したがって、図６Ｂは、図２のメタルアーチファクト補正方法を適用して生成された画像である画像６０４が、明らかなノイズ低減を提供するものであり、その結果として、メタルアーチファクトを低減するためのより効率的な方法を提供するものであることを示している。
【0032】
図６Ｃは、図６Ａの画像６０２、６０４及び６０６上の１つの位置で生成された信号を示している。図６Ｃは、画像６０２、６０４及び６０６を通る線６２６を示している。また、グラフ図６２８は、画像６０２、６０４及び６０６上の線６２６の位置に対応して生成された信号を示している。例として、画像６０２上の線６２６の位置に対応する信号６３０が生成され、画像６０４上の線６２６の位置に対応する信号６３２が生成され、画像６０６上の線６２６の位置に対応する信号６３４が生成されている。楕円６３６及び６３８は、ノイズ低減の方法による結果の違いが大きい位置を強調しており、図２の方法は、目標であるノイズの無い信号に近い結果を提供している。ここで、楕円６３６は、図２の方法が、図１の既知の方法と同様に高信号を維持することを示している。また、楕円６３８は、図２の方法が、低信号を有する領域（即ち、ノイズのある領域内）において、目標であるノイズの無い信号に近い結果を提供することを示している。概して、図２の方法は、高信号の領域で信号を維持しつつ、低信号（即ち、ノイズ）の領域でノイズを低減することを意図している。
【0033】
ここで、図１及び図２を参照しつつ、図７Ａについて説明する。具体的には、図７Ａは、（１）スキャン画像の関心領域７０８を通る線に沿って図１のノイズ低減の方法を適用して生成された画像７０２、及び、（２）当該スキャン画像の対応する関心領域７１２を通る線に沿って図２のノイズ低減の方法を適用して生成された画像７０４を示している。また、画像７０６は、画像７０２と画像７０４との違いを説明する目的で、理想的な状態（ノイズの無い状態）で生成された参照画像である。
【0034】
画像７０２は、図１のメタルアーチファクト補正方法を適用して生成された画像である。また、関心領域７０８は、スキャン画像における第１の金属インプラント関心領域を示しており、関心領域７１０は、第２の金属インプラント関心領域を示している。例として、金属インプラントは、金属ねじであるが、他の金属インプラントも含まれ得る。関心領域７０８及び７１０は、それぞれ図７Ｂ及び図７Ｃに詳細に示されている。
【0035】
画像７０４は、図２のメタルアーチファクト補正方法を適用して生成された画像である。また、関心領域７１２は、スキャン画像における第１の金属インプラント関心領域を示しており、関心領域７１４は、第２の金属インプラント関心領域を示している。関心領域７１２及び７１４は、それぞれ図７Ｂ及び図７Ｃに詳細に示されている。画像７０６は、前述したように、参照画像である。また、関心領域７１６は、スキャン画像における第１の金属インプラント関心領域を示しており、関心領域７１８は、第２の金属インプラント関心領域を示している。
【0036】
ここで、図７Ｂは、図７Ａの関心領域７０８、７１６及び７１２の拡大画像を示している。図７Ｂから分かるように、関心領域７０８と関心領域７１２と比較すると、関心領域７１２は、ノイズが明らかに少ない。また、関心領域７１２と参照画像の関心領域７１６とを比較すると、関心領域７１２と参照画像の関心領域７１６との間のノイズの違いが明らかに小さいことが分かる。一方、図７Ｃは、図７Ａの関心領域７１０、７１４及び７１８の拡大画像を示している。図７Ｃから分かるように、関心領域７１４と関心領域７１０と比較すると、関心領域７１４は、ノイズが明らかに少ない。また、関心領域７１４と参照画像の関心領域７１８と比較すると、関心領域７１４と参照画像の関心領域７１８との間のノイズの違いが明らかに小さいことが分かる。したがって、図７Ｂ及び図７Ｃは、図２のメタルアーチファクト補正方法を適用して生成された画像である画像７０４が、明らかなノイズ低減を提供するものであり、その結果として、メタルアーチファクトを低減するためのより効率的な方法を提供するものであることを示している。
【0037】
図７Ｄは、図７Ａの画像７０２、７０４及び７０６上の１つの位置で生成された信号を示している。図７Ｄは、画像７０２、７０４及び７０６を通る線７２０を示している。また、グラフ図７２２は、画像７０２、７０４及び７０６上の線７２０の位置に対応して生成された信号を示している。例として、画像７０２上の線７２０の位置に対応する信号７２４が生成され、画像７０４上の線７２０の位置に対応する信号７２６が生成され、画像７０６上の線７２０の位置に対応する信号７２８が生成されている。楕円７３０及び７３２は、ノイズ低減の方法による結果の違いが大きい位置を強調しており、図２の方法は、目標であるノイズの無い信号に近い結果を提供している。ここで、楕円７３０は、図２の方法が、図１の既知の方法と同様に高信号を維持することを示している。また、楕円７３２は、図２の方法が、低信号を有する領域（即ち、ノイズのある領域）において、目標であるノイズの無い信号に近い結果を提供することを示している。
【0038】
図８は、本開示の方法５００を実施可能な医用イメージング装置８６０の例示的な実施形態を示している。医用イメージング装置８６０は、少なくとも１つのスキャン装置８６２、それぞれが特別に構成されたコンピュータ装置（例えば、特別に構成されたデスクトップコンピュータ、特別に構成されたラップトップコンピュータ、特別に構成されたサーバ）である１つ以上の画像生成装置８６４と、表示装置８６６と、を含む。
【0039】
スキャン装置８６２は、被検体（例えば、患者）の領域（例えば、エリア、ボリューム、スライス）をスキャンすることによってスキャンデータを収集するように構成されている。スキャンの手段は、例えば、ＭＲＩ、コンピュータ断層撮影（Computed Tomography：ＣＴ）、ポジトロン放射断層撮影（Positron Emission Tomography：ＰＥＴ）、Ｘ線透過撮影及び超音波検査等である。
【0040】
１つ以上の画像生成装置８６４は、スキャン装置８６２からスキャンデータを取得し、当該スキャンデータに基づいて被検体の部位の画像を生成する。１つ以上の画像生成装置８６４は、画像を生成するために、例えば、中間画像生成中又は最終画像再構成中にスキャンデータに対して再構成処理を実行してもよい。再構成処理の例としては、圧縮センシング（Compressed Sensing：ＣＳ）、ＧＲＡＰＰＡ（GeneRalized Autocalibrating Partially Parallel Acquisitions）、ＣＧ－ＳＥＮＳＥ（Conjugate Gradient-SENSitivity Encoding）、ＳＥＮＳＥ（SENSitivity Encoding）、ＡＲＣ（Autocalibrating Reconstruction for Cartesian imaging）、ＳＰＩＲＩＴ（Iterative Self-consistent Parallel Imaging Reconstruction）及びＬＯＲＡＫＳ（LOw-RAnk modeling of local K-Space neighborhoods）等が挙げられる。
【0041】
一実施形態では、１つ以上の画像生成装置８６４は、画像を生成した後に、当該画像を表示装置８６６に送信し、表示装置８６６は、当該画像を表示する。
【0042】
別の実施形態では、上記に加えて、１つ以上の画像生成装置８６４は、同じスキャンデータから複数の画像を生成してもよい。１つ以上の画像生成装置８６４は、同じスキャンデータから異なる再構成プロセスを用いて２つの画像を生成してもよく、一方の画像の解像度が、他方の画像の解像度より低くてもよい。また、１つ以上の画像生成装置８６４は、１つの画像を生成してもよい。
【0043】
ここで、図９を参照すると、図９は、ＭＲＩ装置９７０の非限定的な実施例を示している。図９に示すＭＲＩ装置９７０は、架台９７１（概略断面図で示す）、及び、インターフェースを介して相互に接続された様々な関連するＭＲＩシステム構成要素９７２を含む。通常、少なくとも架台９７１は、シールドルーム内に配置される。図９に示すＭＲＩ装置９７０の幾何学的構造は、静磁場Ｂ０磁石９７３、Ｇｘ、Ｇｙ、Ｇｚ傾斜磁場コイルセット９７４、及び、全身用ＲＦコイル（Whole-Body RF Coil：ＷＢＣ）９７５を含む実質的に同軸の円筒状構造を含む。この円筒状構造に含まれる要素の水平軸に沿って、患者寝台９７８によって支持された患者９７７の頭部を実質的に取り囲むように、イメージングボリューム９７６が示されている。
【0044】
イメージングボリューム９７６において、患者の頭部（本明細書では、例えば、「スキャン被検体」又は「被検体」と称される）により近接するように、１つ以上の小型のアレイＲＦコイル９７９が接続される。当業者には理解されるように、ＷＢＣ（全身用ＲＦコイル）と比べて相対的に小さい、特定の身体部分（例えば、腕、肩、肘、手首、膝、脚、胸、脊椎等）に対してカスタマイズされた表面コイル等のＲＦコイル及び／又はアレイが用いられる。ここでは、このような小型のＲＦコイルをアレイコイル（Array Coil：ＡＣ）又はフェーズドアレイコイル（Phased-Array Coil：ＰＡＣ）と呼ぶ。これらは、イメージングボリューム９７６にＲＦ信号を送信するように構成された少なくとも１つのコイルと、イメージングボリューム９７６内の患者の頭部等の被検体からＲＦ信号を受信するように構成された複数の受信コイルとを含む。
【0045】
ＭＲＩ装置９７０は、ディスプレイ９８０、キーボード９８１及びプリンタ９８２に接続された入出力ポートを有するＭＲＩシステムコントローラ９８３を含む。当業者には理解されるように、ディスプレイ９８０は、制御入力も提供するようなタッチスクリーン型であってもよい。また、マウス又は他の入出力デバイスが用いられてもよい。
【0046】
ＭＲＩシステムコントローラ９８３は、ＭＲＩシーケンスコントローラ９８４とインターフェースを介して接続されており、Ｇｘ、Ｇｙ、Ｇｚ傾斜磁場コイルドライバ９８５に加えて、ＲＦ送信機９８６、及び、（同一のＲＦコイルが送信及び受信の両方に使用される場合は）送信／受信スイッチ９８７を順次に制御する。ＭＲＩシーケンスコントローラ９８４は、パラレルイメージングを含むＭＲＩ（核磁気共鳴イメージング、即ち、ＮＭＲイメージングとしても知られている）の技法を実行するための適切なプログラムコード構造９８８を含む。ＭＲＩシーケンスコントローラ９８４は、パラレルイメージングを伴うＭＲＩ用、又は、パラレルイメージングを伴わないＭＲＩ用に構成される。さらに、ＭＲＩシーケンスコントローラ９８４は、１つ以上の準備スキャン（プレスキャン）シーケンス、及び、メインスキャン磁気共鳴（Magnetic Resonance：ＭＲ）画像（診断画像と称される）を収集するためのスキャンシーケンスを実行させる。プレスキャンによって得られるＭＲデータは、例えば、全身ＲＦコイル９７５及び／又はアレイＲＦコイル９７９の感度マップ（コイル感度マップ又は空間感度マップと称されることもある）を計測するため、及び、パラレルイメージング用の展開マップを計測するために用いられる。
【0047】
ＭＲＩシステム構成要素９７２は、ディスプレイ９８０に送信される処理済みの画像データを生成するようにＭＲＩデータプロセッサ９９０に入力を提供するＲＦ受信機９８９を含む。また、ＭＲＩデータプロセッサ９９０は、以前に生成されたＭＲデータやＭＲ画像、マップ（例えば、コイル感度マップ、パラレルイメージング用の展開マップ、歪みマップ等）、システム構成パラメータ等が格納されたメモリ９９１、及び、ＭＲＩ画像再構成プログラムコード構造９９２及び９９３にアクセスするように構成されている。
【0048】
一実施形態では、ＭＲＩデータプロセッサ９９０は、処理回路によって実現される。処理回路は、例えば、特定用途向け集積回路（Application Specific Integrated Circuit ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple Programmable Logic Device：ＳＰＬＤ）、複雑プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）、及び、フィールドプログラマブルゲートアレイ（Field Programmable Gate Array：ＦＰＧＡ））等のデバイス、並びに、本開示に記載された機能を実行するように構成された他の回路構成要素を含む。
【0049】
ＭＲＩデータプロセッサ９９０は、本明細書に記載の方法５００等、ＭＲＩ画像再構成プログラムコード構造９９２及び９９３に含まれている１つ以上の命令の１つ以上のシーケンスを実行する。別の方法として、当該命令は、ハードディスク又はリムーバブルメディアドライブ等の別のコンピュータ可読媒体から読み取られてもよい。ＭＲＩ画像再構成プログラムコード構造９９２及び９９３に含まれる命令のシーケンスを実行するために、マルチプロセッシング構成における１つ以上のプロセッサが用いられてもよい。別の実施形態では、ソフトウェアの命令の代わりに、又は、ソフトウェアの命令と組み合わせて、ハードワイヤード回路が用いられてもよい。このように、本開示の実施形態は、ハードウェア回路及びソフトウェアの何らかの特定の組み合わせに限定されるものではない。
【0050】
さらに、本明細書で使用されている「コンピュータ可読媒体」の用語は、実行用のＭＲＩデータプロセッサ９９０に命令を供給する任意の非一時的な媒体を意味する。コンピュータ可読媒体は、不揮発性の媒体や揮発性の媒体を含む多くの形態のものであり得るが、これらに限定されるものではない。不揮発性の媒体としては、例えば、光学ディスク、磁気ディスク及び光磁気ディスク、又は、リムーバブルメディアのドライブが挙げられる。また、揮発性の媒体には、ダイナミックメモリが含まれる。
【0051】
また、図９には、上述したＭＲＩ画像再構成プログラムコード構造９９３が記憶された記憶装置（メモリ）が一般化されて示されている。ここで、ＭＲＩ画像再構成プログラムコード構造９９３は、ＭＲＩ装置９７０に含まれる種々のデータ処理コンポーネントによって呼び出し可能な非一時的コンピュータ可読記憶媒体に記憶されている。当業者であれば理解されるように、ＭＲＩ画像再構成プログラムコード構造９９３の記憶装置は、分割されて、少なくとも部分的に、ＭＲＩシステム構成要素９７２の中で、通常の動作中に最も即時にＭＲＩ画像再構成プログラムコード構造９９３を必要とする異なる処理コンピュータに直接接続される（即ち、共有的に記憶されて、ＭＲＩシステムコントローラ９８３に直接接続されるのではない）ように構成されていてもよい。
【0052】
さらに、図９に示されるＭＲＩ装置９７０は、以下に記載される例示的な実施形態を実行するために利用される。上述したＭＲＩシステム構成要素９７２は、異なる論理集合の「ボックス」に分割されていてもよく、典型的には、多数のデジタルシグナルプロセッサ（Digital Signal Processor：ＤＳＰ）、マイクロプロセッサ及び特殊目的処理回路（例えば、高速Ａ／Ｄ（Analog to Digital）変換、高速フーリエ変換、アレイ処理等用の処理回路）を含む。これらの各処理回路は、典型的には、物理的なデータ処理回路がクロックサイクル（又は、所定数のクロックサイクル）ごとに１つの物理的状態から他の物理的状態へ進行する、クロック制御された「ステートマシン」である。
【0053】
さらに、動作の過程において、処理回路（例えば、ＣＰＵ、レジスタ、バッファ、演算ユニット等）の物理的状態が１つのクロックサイクルから他のクロックサイクルへ漸進的に変化するだけでなく、データ記憶媒体（例えば、磁気記憶媒体内のビット記憶部）の関連する物理的状態も、そのようなシステムの動作中に１つの状態から他の状態へ変化する。例えば、画像再構成処理、及び／又は、時には画像再構成マップ（例えば、コイル感度マップ、展開マップ、ゴーストマップ、歪みマップ等）の生成処理の最後に、物理的な記憶媒体に含まれる、コンピュータ可読かつアクセス可能なデータ値記憶部のアレイが、何らかの先行の状態（例えば、全ての均一な「０」の値又は全て「１」の値）から新しい状態へ変化する。ここで、当該アレイの物理的部分の物理的状態は、最小値と最大値との間で変化することで、実世界の物理的イベント及び物理的状態（例えば、イメージングボリューム空間にわたる患者の内部の物理的構造）を表現する。当業者であれば理解されるように、コンピュータ制御プログラムコードの特定の構造がＭＲＩ装置９７０の１つ又は複数のＣＰＵによって命令レジスタに順次ロードされて実行された際にＭＲＩ装置９７０内で動作可能な状態の特定のシーケンスを発生させてその状態を変化させるように、データ値が記憶された当該アレイは、物理的構造を表現し、また、物理的構造を構成する。
【0054】
自明なように、上述した内容を踏まえて、多くの変更及び変形が可能である。したがって、本願が開示する技術は、添付の特許請求の範囲内で、本明細書で具体的に記載されている方法とは別の方法で実施されてもよいことが理解される。
【0055】
また、本開示の実施形態は、以下の付記に記載されるようであってもよい。
【0056】
（付記１）磁気共鳴イメージング（ＭＲＩ）データにおけるアーチファクトを補正する方法であって、方法は、身体部分に関連するスキャン領域のスライス方向に沿って複数のスライスを収集することと、スライス方向に沿って収集された複数のスライスの信号分散の空間的範囲を推定することと、信号分散の推定された空間的範囲に基づいて、スライス方向に沿って収集された複数のスライスの信号を結合して、スキャン領域の再構成された画像を生成することと、を含む、方法。
【0057】
（付記２）収集された複数のスライスのそれぞれを、２次元スライスから３次元（３Ｄ）空間ボリュームスライスに空間的にエンコードして、一連の３Ｄ空間ボリュームスライスを生成すること、をさらに含むが、これに限定されない、（付記１）に記載の方法。
【0058】
（付記３）取得された複数のスライスの信号を結合することが、二乗和計算及び複素数値和の少なくとも１つを使用して、収集された複数のスライスの信号を結合することを含むが、これに限定されない、（付記１）及び（付記２）のいずれか１つに記載の方法。
【0059】
（付記４）スライス方向に沿って収集された複数のスライスの信号分散の空間的範囲を推定することが、３Ｄ空間双極子応答関数を使用してスライス方向に沿って収集された複数のスライスの信号分散の空間的範囲を推定することを含むが、これに限定されない、（付記１）～（付記３）のいずれか１つに記載の方法。
【0060】
（付記５）スライス方向に沿って収集された複数のスライスから、最高の画素強度を有するスライスを特定することと、３Ｄ空間双極子応答関数に基づいて、最高の画素強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定することと、をさらに含み、信号分散の推定された空間的範囲に基づいて、スライス方向に沿って収集された複数のスライスの信号を結合することが、（ａ）最高の画素強度を有するスライスの信号と（ｂ）少なくとも１つの近傍スライスの信号を結合することを含むが、これに限定されない、（付記１）～（付記４）のいずれか１つに記載の方法。
【0061】
（付記６）最高の画素強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定することが、３Ｄ空間双極子応答関数に基づいて、最高の画素強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定することを含むが、これに限定されない、（付記５）に記載の方法。
【0062】
加えて、本発明は、アーチファクトがメタルアーチファクトである、（付記１）～（付記６）のいずれかに記載の方法も含む。
【0063】
（付記７）磁気共鳴イメージング（ＭＲＩ）データにおけるアーチファクトを補正する装置であって、装置は、身体部分に関連するスキャン領域のスライス方向に沿って複数のスライスを収集し、スライス方向に沿って収集された複数のスライスの信号分散の空間的範囲を推定し、信号分散の推定された空間的範囲に基づいて、スライス方向に沿って収集された複数のスライスの信号を結合して、スキャン領域の再構成された画像を生成する、ように構成された処理回路、を含むが、これに限定されない、装置。
【0064】
（付記８）処理回路が、収集された複数のスライスのそれぞれを、２次元スライスから３次元（３Ｄ）空間ボリュームスライスに空間的にエンコードして、一連の３Ｄ空間ボリュームスライスを生成するように構成された処理回路、をさらに含むが、これに限定されない、（付記７）に記載の装置。
【0065】
（付記９）収集された複数のスライスの信号を結合するように構成された処理回路が、二乗和計算及び複素数値和の少なくとも１つを使用して、収集された複数のスライスの信号を結合するように構成された処理回路を含むが、これに限定されない、（付記７）及び（付記８）のいずれか１つに記載の装置。
【0066】
（付記１０）スライス方向に沿って収集された複数のスライスの信号分散の空間的範囲を推定するように構成された処理回路が、３Ｄ空間双極子応答関数を使用してスライス方向に沿って収集された複数のスライスの信号分散の空間的範囲を推定するように構成された処理回路を含むが、これに限定されない、（付記７）～（付記９）のいずれか１つに記載の装置。
【0067】
（付記１１）スライス方向に沿って収集された複数のスライスから、最高の画素強度を有するスライスを特定するように構成された処理回路と、３Ｄ空間双極子応答関数に基づいて、最高の画素強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定するように構成された処理回路と、をさらに含み、信号分散の推定された空間的範囲に基づいて、スライス方向に沿って収集された複数のスライスの信号を結合するように構成された処理回路が、（ａ）最高の画素強度を有するスライスの信号と（ｂ）少なくとも１つの近傍スライスの信号を結合するように構成された処理回路を含むが、これに限定されない、（付記７）～（付記１０）のいずれか１つに記載の装置。
【0068】
（付記１２）最高の画素強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定するように構成された処理回路が、３Ｄ空間双極子応答関数に基づいて、最高の画素強度を有するスライスに近傍する少なくとも１つの近傍スライスを特定するように構成された処理回路を含む、（付記１１）に記載の装置。
【0069】
加えて、本発明は、アーチファクトがメタルアーチファクトである、（付記７）～（付記１２）のいずれかに記載の装置も含む。
【0070】
（付記１３）コンピュータによって実行されると、磁気共鳴イメージング（ＭＲＩ）データにおけるアーチファクトを補正する方法をコンピュータに実行させる、コンピュータ可読命令を格納する非一時的コンピュータ可読記憶媒体であって、方法は、（付記１）から（付記６）のいずれか１つに記載の方法を含むが、これに限定されない、非一時的コンピュータ可読記憶媒体。
【0071】
加えて、本発明は、アーチファクトがメタルアーチファクトである、（付記１３）に記載の非一時的コンピュータ可読記憶媒体も含む。
【0072】
したがって、前述の議論は、本発明の単なる例示的な実施形態を開示し、記載するにすぎない。当業者には理解されるように、本発明は、その精神又は本質的特徴から逸脱することなしに、他の特定の形態で具体化されてもよい。したがって、本発明の開示は、本発明の範囲並びに他の特許請求の範囲を限定するものではなく、例示することを意図している。本明細書での教示の容易に認識可能な変形を含む本開示は、発明の主題が公共の用に供されないように、前述の特許請求の範囲の用語の範囲を、部分的に定義する。
【0073】
以上説明した少なくとも一つの実施形態によれば、アーチファクトの低減を基にしたＭＲＩにおいて改善されたノイズ低減を実現することができる。
【0074】
いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【0075】
９７０ ＭＲＩ装置
９９０ ＭＲＩデータプロセッサ
９９０ ＭＲＩデータプロセッサ
【図1】
【図2】
【図3】
【図4】
【図5】
【図6A】
【図6B】
【図6C】
【図7A】
【図7B】
【図7C】
【図7D】
【図8】
【図9】