(57)【要約】      （修正有）【課題】デノイズ精度を高めることのできる医用情報処理装置を提供する。【解決手段】医用情報処理装置２００は、導出部と、調整部と、除去部と、合成部とを備える。導出部は、複数の受信コイルの各々により収集される磁気共鳴信号に対応するデータに含まれるノイズに関する指標値を導出する。調整部は、導出部により導出された指標値に基づいて、磁気共鳴信号に対応するデータからノイズを除去する度合いを調整する。除去部は、調整部により調整された度合いに基づいて、磁気共鳴信号に対応するデータからノイズを除去する。合成部は、除去部によりノイズが除去された、複数の受信コイルの各々で収集された磁気共鳴信号に対応するデータを合成する。【選択図】図４
2020103890.tif 000002
【特許請求の範囲】
【請求項１】
複数の受信コイルの各々により収集される磁気共鳴信号に対応するデータに含まれるノイズに関する指標値を導出する導出部と、
前記導出部により導出された前記指標値に基づいて、前記磁気共鳴信号に対応するデータからノイズを除去する度合いを調整する調整部と、
前記調整部により調整された前記度合いに基づいて、前記磁気共鳴信号に対応するデータからノイズを除去する除去部と、
前記除去部によりノイズが除去された、前記複数の受信コイルの各々で収集された前記磁気共鳴信号に対応するデータを合成する合成部と
を備える医用情報処理装置。
【請求項２】
前記磁気共鳴信号に対応するデータは、前記磁気共鳴信号に基づく複数の画像のデータを含み、
前記導出部は、前記画像の信号対雑音比を前記指標値として導出し、
前記調整部は、前記信号対雑音比が小さくなるほど前記度合いを高くし、前記信号対雑音比が大きくなるほど前記度合いを低くする、
請求項１に記載の医用情報処理装置。
【請求項３】
前記除去部は、磁気共鳴信号に対応するデータが入力されたときに該データからノイズが除去されたデータを出力するように学習されたモデルに基づいて、前記磁気共鳴信号に対応するデータからノイズを除去する、
請求項１または２に記載の医用情報処理装置。
【請求項４】
前記調整部は、前記モデルの内部パラメータを調整することにより、前記度合いを調整する、
請求項３に記載の医用情報処理装置。
【請求項５】
前記内部パラメータは、前記モデルに含まれる少なくとも１つのノードの活性化関数のパラメータである、
請求項４に記載の医用情報処理装置。
【請求項６】
複数の受信コイルの各々により被検体から収集される磁気共鳴信号に基づく複数の画像を合成して合成画像を生成する合成部と、
前記合成部により生成された前記合成画像に含まれるノイズの強度分布を推定し、推定された前記ノイズの強度分布を均一化させる均一化分布を生成し、生成された前記均一化分布により示される信号を前記合成画像に付加することで、前記合成画像に含まれるノイズの強度分布を均一化する均一化処理部と、
前記均一化処理部によりノイズの強度分布が均一化された前記合成画像からノイズを除去する除去部と、
を備える医用情報処理装置。
【請求項７】
コンピュータが、
複数の受信コイルの各々により収集される磁気共鳴信号に対応するデータに含まれるノイズに関する指標値を導出し、
導出された前記指標値に基づいて、前記磁気共鳴信号に対応するデータからノイズを除去する度合いを調整し、
調整された前記度合いに基づいて、前記磁気共鳴信号に対応するデータからノイズを除去し、
前記ノイズが除去された、前記複数の受信コイルの各々で収集された前記磁気共鳴信号に対応するデータを合成する、
医用情報処理方法。
【請求項８】
コンピュータが、
複数の受信コイルの各々により被検体から収集される磁気共鳴信号に基づく複数の画像を合成して合成画像を生成し、
生成された前記合成画像に含まれるノイズの強度分布を推定し、
推定された前記ノイズの強度分布を均一化させる均一化分布を生成し、
生成された前記均一化分布により示される信号を前記合成画像に付加することで、前記合成画像に含まれるノイズの強度分布を均一化し、
前記ノイズの強度分布が均一化された合成画像からノイズを除去する、
医用情報処理方法。
【請求項９】
コンピュータに、
複数の受信コイルの各々により収集される磁気共鳴信号に対応するデータに含まれるノイズに関する指標値を導出させ、
導出された前記指標値に基づいて、前記磁気共鳴信号に対応するデータからノイズを除去する度合いを調整させ、
調整された前記度合いに基づいて、前記磁気共鳴信号に対応するデータからノイズを除去させ、
前記ノイズが除去された、前記複数の受信コイルの各々で収集された前記磁気共鳴信号に対応するデータを合成させる、
プログラム。
【請求項１０】
コンピュータに、
複数の受信コイルの各々により被検体から収集される磁気共鳴信号に基づく複数の画像を合成して合成画像を生成させ、
生成された前記合成画像に含まれるノイズの強度分布を推定させ、
推定された前記ノイズの強度分布を均一化させる均一化分布を生成させ、
生成された前記均一化分布により示される信号を前記合成画像に付加することで、前記合成画像に含まれるノイズの強度分布を均一化させ、
前記ノイズの強度分布が均一化された合成画像からノイズを除去させる、
プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、医用情報処理装置、医用情報処理方法、およびプログラムに関する。
【背景技術】
【０００２】
従来、ＭＲＩ（Magnetic Resonance Imaging）装置に設けられた複数の受信コイルにより得られる信号の平方二乗和（ＳｏＳ：Sum of Squares）をとって画像を生成し、さらに、生成された画像に対してデノイズ処理を行うことでノイズ成分が除去或いは抑制された画像を得る技術が知られている。デノイズ処理としては、例えば、ディープニューラルネットワークを利用することが可能である。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】米国特許出願公開第２０１８／００１８７５７号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明が解決しようとする課題は、デノイズ精度を高めることである。
【課題を解決するための手段】
【０００５】
実施形態の医用情報処理装置は、導出部と、調整部と、除去部と、合成部とを備える。導出部は、複数の受信コイルの各々により収集される磁気共鳴信号に対応するデータに含まれるノイズに関する指標値を導出する。調整部は、前記導出部により導出された前記指標値に基づいて、前記磁気共鳴信号に対応するデータからノイズを除去する度合いを調整する。除去部は、前記調整部により調整された前記度合いに基づいて、前記磁気共鳴信号に対応するデータからノイズを除去する。合成部は、前記除去部によりノイズが除去された、前記複数の受信コイルの各々で収集された前記磁気共鳴信号に対応するデータを合成する。
【図面の簡単な説明】
【０００６】
【図１】第１の実施形態に係る医用画像処理装置２００を含む医用画像処理システム１の構成の一例を示す図。
【図２】第１の実施形態に係る医用画像生成装置１００の一例を示す図。
【図３】第１の実施形態に係る医用画像生成装置１００に備えられるＲＦコイル１０８の配置の一例を示す図。
【図４】第１の実施形態に係る医用画像処理装置２００の一例を示す図。
【図５】第１の実施形態に係るデノイズモデルＭＤＬ１の構成の一例を示す図。
【図６】第１の実施形態に係るアクティベーション層３３０の活性化関数の一例を示す図。
【図７】第１の実施形態に係る原画像からデノイズ画像および合成画像を生成する構成の一例を示す図。
【図８】第１の実施形態における処理回路２１０の学習処理の一連の流れを示すフローチャート。
【図９】第１の実施形態における処理回路２１０の画像処理の一連の流れを示すフローチャート。
【図１０】第１の実施形態の変形例に係る医用画像生成装置１００の一例を示す図。
【図１１】第２の実施形態に係る医用画像処理装置２００の一例を示す図。
【図１２】第２の実施形態における処理回路２１０の画像処理の一連の流れを示すフローチャート。
【図１３】第２の実施形態の変形例に係る医用画像生成装置１００を示す図。
【発明を実施するための形態】
【０００７】
以下、図面を参照しながら、医用情報処理装置、医用情報処理方法、およびプログラムの実施形態について詳細に説明する。
【０００８】
（第１の実施形態）
図１は、第１の実施形態に係る医用画像処理装置２００を含む医用画像処理システム１の構成の一例を示す図である。図１に示すように、医用画像処理システム１は、例えば、医用画像生成装置１００と、医用画像処理装置２００とを備える。医用画像生成装置１００と医用画像処理装置２００とは、ネットワークＮＷを介して互いに接続される。ネットワークＮＷは、例えば、ＷＡＮ（Wide Area Network）やＬＡＮ（Local Area Network）、インターネット、専用回線、無線基地局、プロバイダなどを含む。医用画像処理装置２００は、「医用情報処理装置」の一例である。
【０００９】
医用画像生成装置１００は、例えば、ＭＲＩ装置などを含む。ＭＲＩ装置は、例えば、被検体（例えば人体）に磁場を与えて、核磁気共鳴現象によって被検体内の水素原子核から発生する電磁波を、コイルを利用して受信し、その受信した電磁波に基づく信号を再構成することで医用画像（ＭＲ画像）を生成する。以下の説明では、一例として医用画像生成装置１００がＭＲＩ装置であるものとして説明する。
【００１０】
医用画像処理装置２００は、一つまたは複数のプロセッサにより実現される。例えば、医用画像処理装置２００は、クラウドコンピューティングシステムに含まれるコンピュータであってもよいし、他の機器に依存せずに単独で動作するコンピュータ（スタンドアローンのコンピュータ）であってもよい。
【００１１】
［医用画像生成装置（ＭＲＩ装置）の構成例］
図２は、第１の実施形態に係る医用画像生成装置１００の一例を示す図である。図２に示すように、医用画像生成装置１００は、例えば、静磁場磁石１０１と、傾斜磁場コイル１０２と、傾斜磁場電源１０３と、寝台１０４と、寝台制御回路１０５と、送信回路１０７と、ＲＦコイル１０８と、受信回路１０９と、シーケンス制御回路１１０と、コンソール装置１２０とを備える。
【００１２】
静磁場磁石１０１は、中空の略円筒形状に形成された磁石である。静磁場磁石１０１は、内部の空間に一様な静磁場を発生させる。静磁場磁石１０１は、例えば、永久磁石や超伝導磁石などである。傾斜磁場コイル１０２は、中空の略円筒形状に形成されたコイルであり、静磁場磁石１０１の内側に配置される。傾斜磁場コイル１０２は、互いに直交するｘ，ｙ，ｚの各軸に対応する３つのコイルが組み合わされて形成される。ｚ軸方向は、寝台１０４の天板１０４ａの長手方向を表し、ｘ軸方向は、ｚ軸方向に直交し、医用画像生成装置１００が設置される部屋の床面に対して平行である軸方向を表し、ｙ軸方向は、床面に対して垂直方向である軸方向を表している。各軸方向に対応した３つのコイルは、傾斜磁場電源１０３から個別に電流を受けて、ｘ，ｙ，ｚの各軸に沿って磁場強度が変化する傾斜磁場を発生させる。なお、ｚ軸方向は、静磁場と同方向とする。
【００１３】
傾斜磁場電源１０３は、傾斜磁場コイル１０２に電流を供給する。ここで、傾斜磁場コイル１０２によって発生するｘ，ｙ，ｚの各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ、およびリードアウト用傾斜磁場Ｇｒにそれぞれ対応する。スライス選択用傾斜磁場Ｇｓは、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場Ｇｅは、空間的位置に応じて磁気共鳴信号の位相を変化させるために利用される。リードアウト用傾斜磁場Ｇｒは、空間的位置に応じて磁気共鳴信号の周波数を変化させるために利用される。
【００１４】
寝台１０４は、被検体ＯＢが載置される天板１０４ａを備える。寝台１０４は、寝台制御回路１０５による制御のもと、天板１０４ａを、被検体ＯＢが載置された状態で傾斜磁場コイル１０２の空洞（撮像口）内へ挿入する。通常、寝台１０４は、長手方向が静磁場磁石１０１の中心軸と平行になるように設置される。寝台制御回路１０５は、コンソール装置１２０による制御のもと、寝台１０４を駆動して天板１０４ａを長手方向および上下方向へ移動する。
【００１５】
ＲＦコイル１０８は、送信回路１０７からＲＦ（Radio Frequency）パルスの供給を受けて、高周波磁場を発生する。送信回路１０７は、対象とする原子核の種類および磁場の強度で決まるラーモア周波数に対応するＲＦパルスをＲＦコイル１０８に供給する。また、ＲＦコイル１０８は、高周波磁場の影響によって被検体ＯＢから発せられる磁気共鳴信号を受信する。磁気共鳴信号には、例えば、信号強度成分と位相成分が含まれる。ＲＦコイル１０８は、磁気共鳴信号を受信すると、受信した磁気共鳴信号を受信回路１０９へ出力する。
【００１６】
ＲＦコイル１０８は、例えば、ＭＲＩ装置の架台内に収められ、被検体ＯＢを取り囲むように構成された全身用コイルや、被検体ＯＢの撮像部位ごとに設けられる局所コイルである。以下では、ＲＦコイル１０８として局所コイルを例にとって説明するが、ＲＦコイル１０８の種別を限定することは意図しない。また、送信と受信をそれぞれ異なるＲＦコイルで行ってもよいし、ＲＦコイル１０８を送受信兼用として構成してもよい。なお、第１の実施形態において、ＲＦコイル１０８は、複数のコイルエレメントを有するコイルアレイである。
【００１７】
図３は、第１の実施形態に係る医用画像生成装置１００に備えられるＲＦコイル１０８の配置の一例を示す図である。図３では、ＲＦコイル１０８が、８つのコイルエレメント１０８ａ〜１０８ｈを備える例を示す。これらのコイルエレメント１０８ａ〜１０８ｈは、例えばｚ軸に沿って被検体ＯＢを取り囲むように配置される。各コイルエレメント１０８ａ〜１０８ｈは、被検体ＯＢから発せられた磁気共鳴信号を受信し、受信回路１０９へ出力する。
【００１８】
受信回路１０９は、コイルエレメント１０８ａ〜１０８ｈにより出力される磁気共鳴信号ごとに磁気共鳴データを生成する。例えば、受信回路１０９は、コイルエレメント１０８ａ〜１０８ｈにより出力されるアナログ信号である磁気共鳴信号の各々をデジタル変換することによって、デジタル信号である磁気共鳴データの組を生成する。また、受信回路１０９は、生成した磁気共鳴データの組をシーケンス制御回路１１０へ送信する。なお、受信回路１０９は、静磁場磁石１０１や傾斜磁場コイル１０２などを備える架台装置側に備えられていてもよい。
【００１９】
シーケンス制御回路１１０は、コンソール装置１２０により出力されるシーケンス情報に基づいて、傾斜磁場電源１０３、送信回路１０７および受信回路１０９を駆動することによって、被検体ＯＢを撮像する。シーケンス情報は、撮像処理を行うための手順を定義した情報である。シーケンス情報には、傾斜磁場電源１０３が傾斜磁場コイル１０２に供給する電流の大きさや電流を供給するタイミング、送信回路１０７がＲＦコイル１０８に送信するＲＦパルスの強さやＲＦパルスを印加するタイミング、受信回路１０９が磁気共鳴信号を検出するタイミングなどが定義された情報が含まれる。
【００２０】
なお、シーケンス制御回路１１０は、傾斜磁場電源１０３、送信回路１０７および受信回路１０９を駆動させ、受信回路１０９から磁気共鳴データを受信すると、受信した磁気共鳴データをコンソール装置１２０へ転送する。
【００２１】
コンソール装置１２０は、医用画像生成装置１００の全体を制御したり、磁気共鳴データを収集したりする。例えば、コンソール装置１２０は、通信インターフェース１２２と、入力インターフェース１２４と、ディスプレイ１２６と、処理回路１３０と、メモリ（ストレージ）１５０とを備える。
【００２２】
通信インターフェース１２２は、例えば、ＮＩＣ（Network Interface Card）などの通信インターフェースを含む。通信インターフェース１２２は、ネットワークＮＷを介して医用画像処理装置２００と通信し、医用画像処理装置２００との間で情報を送受信する。通信インターフェース１２２は、受信した情報を処理回路１３０に出力する。また、通信インターフェース１２２は、処理回路１３０による制御を受けて、ネットワークＮＷを介して接続された他の装置に情報を送信してもよい。
【００２３】
入力インターフェース１２４は、操作者から各種の入力操作を受け付けるインターフェースである。入力インターフェース１２４は、入力操作を受け付けると、その受け付けた入力操作を電気信号に変換して処理回路１３０に出力する。例えば、入力インターフェース１２４は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック、タッチパネルなどにより実現される。また、入力インターフェース１２４は、例えば、マイクなどの音声入力を受け付けるユーザインターフェースによって実現されてもよい。入力インターフェース１２４がタッチパネルである場合、後述するディスプレイ１２６は入力インターフェース１２４と一体として形成されてよい。
【００２４】
ディスプレイ１２６は、各種の情報を表示する。例えば、ディスプレイ１２６は、処理回路１３０によって生成された画像を表示したり、操作者からの各種の入力操作を受け付けるためのＧＵＩ（Graphical User Interface）などを表示したりする。例えば、ディスプレイ１２６は、ＬＣＤ（Liquid Crystal Display）や、ＣＲＴ（Cathode Ray Tube）ディスプレイ、有機ＥＬ（Electroluminescence）ディスプレイなどである。
【００２５】
処理回路１３０は、例えば、取得機能１３２と、生成機能１３４と、通信制御機能１３６と、表示制御機能１３８とを備える。これらの機能（構成要素）は、例えば、ＣＰＵ（Central Processing Unit）やＧＰＵ（Graphics Processing Unit）などのプロセッサ（あるいはプロセッサ回路）が、メモリ１５０に記憶されたプログラム（ソフトウェア）を実行することにより実現される。また、処理回路１３０の機能の一部または全部は、ＬＳＩ（Large Scale Integration）、ＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field-Programmable Gate Array）などのハードウェア（回路部：circuitry）によって実現されてもよいし、ソフトウェアとハードウェアの協働によって実現されてもよい。また、上記のプログラムは、予めメモリ１５０に格納されていてもよいし、ＤＶＤやＣＤ−ＲＯＭなどの着脱可能な記憶媒体に格納されており、記憶媒体がコンソール装置１２０のドライブ装置に装着されることで記憶媒体からメモリ１５０にインストールされてもよい。
【００２６】
メモリ１５０は、例えば、ＲＡＭ（Random Access Memory）、フラッシュメモリなどの半導体メモリ素子、ハードディスク、光ディスクなどによって実現される。これらの非一過性の記憶媒体は、ＮＡＳ（Network Attached Storage）や外部ストレージサーバ装置といったネットワークＮＷを介して接続される他の記憶装置によって実現されてもよい。また、これらの非一過性の記憶媒体は、ＲＯＭ（Read Only Memory）やレジスタなどの記憶装置によって実現されてもよい。
【００２７】
取得機能１３２は、シーケンス制御回路１１０から磁気共鳴データを取得する。磁気共鳴データは、上述したように、核磁気共鳴現象によって被検体ＯＢ内において発生した電磁波の信号（磁気共鳴信号）をデジタル化して得られる磁気共鳴データである。なお、前述したスライス選択用傾斜磁場Ｇｓ、位相エンコード用傾斜磁場Ｇｅ、及びリードアウト用傾斜磁場Ｇｒによって付与された情報に応じて例えば２次元や３次元に配置された磁気共鳴データは、ｋ空間データとも称される。
【００２８】
生成機能１３４は、取得機能１３２によって取得されたｋ空間データに対してフーリエ変換（例えば、逆フーリエ変換）などの処理を含む再構成処理を行うことで、ｋ空間データから再構成されたＭＲ画像を生成する。生成機能１３４は、コイルエレメント１０８ａ〜１０８ｈにより出力される磁気共鳴信号ごとの磁気共鳴データに対してフーリエ変換（例えば、逆フーリエ変換）などの処理を含む再構成処理を行うことで、コイルエレメント１０８ａ〜１０８ｈにより出力される磁気共鳴信号ごとのＭＲ画像（以下、原画像）の組を生成する。図３に示すように、生成機能１３４は、コイルエレメント１０８ａ〜１０８ｈにより出力される磁気共鳴信号の原画像の組として、原画像ａ〜ｈをそれぞれ生成する。
【００２９】
通信制御機能１３６は、生成機能１３４により原画像ａ〜ｈが再構成により生成されると、通信インターフェース１２２に医用画像処理装置２００と通信させ、その通信相手の医用画像処理装置２００に、再構成された原画像ａ〜ｈを送信する。また、通信制御機能１３６は、通信インターフェース１２２に医用画像処理装置２００と通信させ、その通信相手の医用画像処理装置２００から各種情報を受信させてよい。
【００３０】
表示制御機能１３８は、医用画像処理装置２００から受信した医用画像を、ディスプレイ１２６に表示させる。また、表示制御機能１３８は、生成機能１３４により生成された原画像ａ〜ｈをディスプレイ１２６に表示させてもよい。
【００３１】
［医用画像処理装置の構成例］
図４は、第１の実施形態に係る医用画像処理装置２００の一例を示す図である。医用画像処理装置２００は、医用画像生成装置１００から受信した原画像ａ〜ｈに対して、ノイズを除去或いは軽減するデノイズ処理を行い、デノイズ処理後の画像を合成することで最終的な画像（以下、合成画像）を生成する。以下においては、ノイズを除去或いは軽減することを、単にノイズを除去すると記載する。図４に示すように、医用画像処理装置２００は、例えば、通信インターフェース２０２と、入力インターフェース２０４と、ディスプレイ２０６と、処理回路２１０と、メモリ２３０とを備える。
【００３２】
通信インターフェース２０２は、例えば、ＮＩＣなどの通信インターフェースを含む。例えば、通信インターフェース２０２は、ネットワークＮＷを介して医用画像生成装置１００と通信し、医用画像生成装置１００から、再構成された医用画像を受信する。通信インターフェース２０２は、受信した医用画像を処理回路２１０に出力する。また、通信インターフェース２０２は、処理回路２１０による制御を受けて、ネットワークＮＷを介して接続された医用画像生成装置１００やその他の装置に情報を送信してもよい。他の装置とは、例えば、医師や看護師などの画像の読影者が利用可能な端末装置であってよい。
【００３３】
入力インターフェース２０４は、操作者からの各種の入力操作を受け付け、受け付けた入力操作を電気信号に変換して処理回路２１０に出力する。例えば、入力インターフェース２０４は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック、タッチパネルなどにより実現される。また、入力インターフェース２０４は、例えば、マイクなどの音声入力を受け付けるユーザインターフェースによって実現されてもよい。入力インターフェース２０４がタッチパネルである場合、後述するディスプレイ２０６は入力インターフェース２０４と一体として形成されてよい。
【００３４】
ディスプレイ２０６は、各種の情報を表示する。例えば、ディスプレイ２０６は、処理回路２１０によって生成された画像（後述するデノイズ画像や合成画像）を表示したり、操作者からの各種の入力操作を受け付けるためのＧＵＩなどを表示したりする。例えば、ディスプレイ２０６は、ＬＣＤや、ＣＲＴディスプレイ、有機ＥＬディスプレイなどである。
【００３５】
処理回路２１０は、例えば、取得機能２１２と、導出機能２１４と、パラメータ調整機能２１６と、デノイズ機能２１８と、画像合成機能２２０と、出力制御機能２２２と、学習機能２２４とを実行する。導出機能２１４は、「導出部」の一例であり、パラメータ調整機能２１６は、「調整部」の一例であり、デノイズ機能２１８は、「除去部」の一例であり、画像合成機能２２０は、「合成部」の一例である。
【００３６】
これらの機能（構成要素）は、例えば、ＣＰＵやＧＰＵなどのプロセッサ（あるいはプロセッサ回路）が、メモリ２３０に記憶されたプログラム（ソフトウェア）を実行することにより実現される。また、これらの複数の機能のうち一部または全部は、ＬＳＩ、ＡＳＩＣ、ＦＰＧＡなどのハードウェア（回路部：circuitry）によって実現されてもよいし、ソフトウェアとハードウェアの協働によって実現されてもよい。また、上記のプログラムは、予めメモリ２３０に格納されていてもよいし、ＤＶＤやＣＤ−ＲＯＭなどの着脱可能な記憶媒体に格納されており、記憶媒体が医用画像処理装置２００のドライブ装置に装着されることで記憶媒体からメモリ２３０にインストールされてもよい。
【００３７】
メモリ２３０は、例えば、ＲＡＭ、フラッシュメモリなどの半導体メモリ素子、ハードディスク、光ディスクなどにより実現される。これらの非一過性の記憶媒体は、ＮＡＳや外部ストレージサーバ装置といったネットワークＮＷを介して接続される他の記憶装置によって実現されてもよい。また、これらの非一過性の記憶媒体は、ＲＯＭやレジスタなどの記憶装置によって実現されてもよい。メモリ２３０には、例えば、デノイズモデル情報２３２、原画像の情報（以下、原画像情報２３４）、デノイズ画像の情報（以下、デノイズ画像情報２３６）、合成画像の情報（以下、合成画像情報２３８）などが格納される。
【００３８】
デノイズモデル情報２３２は、後述するデノイズモデルＭＤＬ１を定義した情報（プログラムまたはデータ構造）である。デノイズモデルＭＤＬ１は、ある画像が入力されたときに該画像からノイズが除去された画像を出力するように学習されたモデルである。デノイズモデルＭＤＬ１は、例えば、一つ以上のＤＮＮ（Deep Neural Network(s)）を含む。
【００３９】
デノイズモデル情報２３２には、例えば、デノイズモデルＭＤＬ１に含まれる各ＤＮＮを構成する入力層、一以上の隠れ層（中間層）、出力層の各々に含まれるニューロン（ユニットあるいはノード）が互いにどのように結合されるのかという結合情報や、結合されたニューロン間で入出力されるデータに付与される結合係数がいくつであるのかという重み情報などが含まれる。結合情報は、例えば、各層に含まれるニューロン数や、各ニューロンの結合先のニューロンの種類を指定する情報、各ニューロンを実現する活性化関数、隠れ層のニューロン間に設けられたゲートなどの情報を含む。ニューロンを実現する活性化関数は、例えば、入力符号に応じて動作を切り替える関数（ＲｅＬＵ（Rectified Linear Unit）関数やＥＬＵ（Exponential Linear Units）関数、クリッピング関数）であってもよいし、シグモイド関数や、ステップ関数、ハイパポリックタンジェント関数であってもよいし、恒等関数であってもよい。ゲートは、例えば、活性化関数によって返される値（例えば１または０）に応じて、ニューロン間で伝達されるデータを選択的に通過させたり、重み付けたりする。結合係数は、例えば、ニューラルネットワークの隠れ層において、ある層のニューロンから、より深い層のニューロンにデータが出力される際に、出力データに対して付与される重みを含む。また、結合係数は、各層の固有のバイアス成分などを含んでもよい。
【００４０】
例えば、デノイズモデルＭＤＬ１は、ＣＮＮ（Convolutional Neural Network）によって実現されてよい。ＣＮＮは、畳み込み層や、アクティベーション層などが多層に構成されていてよい。
【００４１】
図５は、デノイズモデルＭＤＬ１の構成の一例を示す図である。図５に示すように、デノイズモデルＭＤＬ１には、例えば、入力層３１０と、一以上の畳み込み層３２０と、一以上のアクティベーション層３３０と、出力層３４０とが含まれてよい。
【００４２】
例えば、原画像を、各画素に対応した要素をもつ行列とした場合、入力層３１０には、原画像に対応した行列が入力される。入力層３１０は、入力された行列に対して適時バイアス成分を加えるなどして、後段の畳み込み層３２０に出力する。
【００４３】
畳み込み層３２０は、入力された行列に対して、フィルタまたはカーネルと呼ばれる線形変換行列を、ある決められたストライド量でスライドさせながら積和演算を繰り返し、入力された行列から、線形変換行列との積和が要素値として対応付けられた複数の要素を含む行列を生成する。この際、畳み込み層３２０は、入力された行列の周囲に任意の値の要素を補間するパディング（例えばゼロパティング）を行って、畳み込み層３２０に入力された行列を、入力層３１０に入力される原画像の行列と行数および列数が同じ行列に変換してよい。そして、畳み込み層３２０は、生成した行列をアクティベーション層３３０に出力する。
【００４４】
アクティベーション層３３０は、畳み込み層３２０から入力された行列の各要素に対して活性化関数の計算処理を行い、その計算処理を行った行列を、後段のレイヤに出力する。
【００４５】
図６は、アクティベーション層３３０の活性化関数の一例を表す図である。図６に示すように、アクティベーション層３３０の活性化関数は、Ｓｏｆｔ‐Ｓｈｒｉｎｋａｇｅ関数であってよい。Ｓｏｆｔ‐Ｓｈｒｉｎｋａｇｅ関数は、例えば、以下の数式（１）によって表される。なお、アクティベーション層３３０の活性化関数は、Ｓｏｆｔ‐Ｓｈｒｉｎｋａｇｅ関数に代えて、Ｈａｒｄ‐Ｓｈｒｉｎｋａｇｅ関数であってもよい。
【００４６】
【数１】
2020103890.tif 000003
【００４７】
Ｓｏｆｔ‐Ｓｈｒｉｎｋａｇｅ関数あるいはＨａｒｄ‐Ｓｈｒｉｎｋａｇｅ関数は、入力値である要素値ｘが、ゼロを中心にした所定の正負の閾値±Ｔの範囲内である場合、ゼロを出力し、入力値である要素値ｘが、正の閾値＋Ｔを超える、あるいは負の閾値−Ｔ未満である場合、その要素値ｘに比例する値を出力する関数である。アクティベーション層３３０の活性化関数に、Ｓｏｆｔ‐Ｓｈｒｉｎｋａｇｅ関数やＨａｒｄ‐Ｓｈｒｉｎｋａｇｅ関数を適用することで、振幅が閾値Ｔよりも小さな画像信号、すなわち、ノイズである蓋然性が高い微弱な画像信号を、活性化関数の出力においてゼロとすることができる。
【００４８】
閾値Ｔは、入力画像に含まれるノイズのレベル（信号強度あるいは信号電力）に応じて変動するパラメータであり、例えば、以下の数式（２）によって表される。
【００４９】
【数２】
2020103890.tif 000004
【００５０】
数式（２）のＧは、入力画像に含まれるノイズのレベルを表すものであり、閾値Ｔの値を制御する信号（以下、制御信号）である。入力画像に含まれるノイズのレベルが大きいほど、制御信号Ｇは大きくなり、入力画像に含まれるノイズのレベルが小さいほど、制御信号Ｇは小さくなる。入力画像に含まれるノイズのレベルは、入力画像の信号対雑音比（ＳＮＲ：Signal-Noise Ratio）の大きさによって決定されてよい。
【００５１】
数式（２）のαは、ノイズレベルに乗算される重み係数である。重み係数αは、機械学習によって決定される。数式（２）に示すように閾値Ｔを変動させることで、アクティベーション層３３０の活性化関数を、ある一定の強度以上の信号に対して反応する活性化関数とすることができる。この結果、原画像に含まれるノイズの信号強度がばらついていても原画像からノイズを精度よく除去することができる。
【００５２】
出力層３４０は、前段の畳み込み層３２０およびアクティベーション層３３０によって処理された行列を出力する。
【００５３】
なお、図５に例示したデノイズモデルＭＤＬ１はあくまでも一例であり、例えば、プーリング層などが含まれていてもよい。プーリング層は、入力された行列の要素値を、その行列に含まれる全要素値の平均値や最大値といった代表値に置き換えることで、行列の次元数を圧縮する（減らす）。プーリング層は、次元数を圧縮した行列を後段のレイヤに出力する。
【００５４】
取得機能２１２は、通信インターフェース２０２に医用画像生成装置１００と通信させ、その通信相手の医用画像生成装置１００から再構成された原画像ａ〜ｈを取得する。取得機能２１２は、取得した原画像ａ〜ｈを、原画像情報２３４としてメモリ２３０に記憶させる。
【００５５】
導出機能２１４は、取得機能２１２によって取得された原画像ａ〜ｈに基づいて、原画像ａ〜ｈの各々に含まれるノイズに関する指標値を導出する。すなわち、導出機能２１４は、複数の受信コイルの各々により被検体ＯＢから収集される磁気共鳴信号に基づく複数の画像に含まれるノイズに関する指標値を、画像ごとに導出する。例えば、導出機能２１４は、ノイズに関する指標値として、ＳＮＲを導出する。ＳＮＲは、画像の信号強度を、ノイズの信号強度で除算した指標値である。ＳＮＲの導出方法の詳細について後述する。導出機能２１４は、原画像ａ〜ｈの信号強度、ノイズの大きさ、正規化処理時の利得など、ＳＮＲに関連するデータを導出してもよい。導出機能２１４は、原画像ａ〜ｈの付帯情報に基づいて、原画像ａ〜ｈの各々に含まれるノイズに関する指標値を導出してもよい。
【００５６】
パラメータ調整機能２１６は、導出機能２１４によって導出されたＳＮＲに基づいて、デノイズモデルＭＤＬ１の内部パラメータである制御関数Ｇを調整する。すなわち、パラメータ調整機能２１６は、導出機能２１４により導出された指標値であるＳＮＲに基づいて、複数の原画像ａ〜ｈの各々からノイズを除去する度合いを調整する。パラメータ調整機能２１６は、信号対雑音比が小さくなるほど度合いを高くし、信号対雑音比が大きくなるほど度合いを低くする。パラメータ調整機能２１６は、アクティベーション層３３０の少なくとも１つのノードに対して調整済みの制御関数Ｇを入力する。
【００５７】
なお、パラメータ調整機能２１６は、複数の原画像ａ〜ｈの画像に対するノイズを除去する度合い（デノイズ強度）をまとめて１つに決定してもよい。すなわち、決定された１つの度合いは、複数の原画像ａ〜ｈの画像の全てに適用可能である。この場合、導出機能２１４は、被検体ＯＢが天板１０４ａに載置され、且つ、ＲＦパルスを印加しない条件下で撮像を行うことにより得られたデータに含まれるノイズを求める。例えば、導出機能２１４は、ＲＦパルスを印加しない条件下で収集したｋ空間データを再構成した画像空間におけるガウシアンノイズの分布を数値化する。パラメータ調整機能２１６は、このようにして数値化したノイズに応じてノイズを除去する度合いを決定してよい。なお、ノイズの強度や分布は、撮像条件に応じて変化する場合があるため、撮像条件が変更されるごとに、ノイズを除去する度合いが調整されてよい。
【００５８】
デノイズ機能２１８は、デノイズモデル情報２３２が示すデノイズモデルＭＤＬ１を利用して、取得機能２１２によって取得された原画像ａ〜ｈの各々からノイズを除去した複数の画像（以下、デノイズ画像）を生成する。図７は、第１の実施形態に係る原画像からデノイズ画像および合成画像を生成する構成の一例を示す図である。図７に示すように、デノイズ機能２１８は、取得機能２１２によって取得された原画像ａ〜ｈからノイズを除去し、デノイズ画像ａ〜ｈをそれぞれ生成する。すなわち、デノイズ機能２１８は、デノイズモデルＭＤＬ１の出力層３４０により出力された行列を、原画像ａ〜ｈからノイズを除去したデノイズ画像ａ〜ｈとして取得する。デノイズ機能２１８は、生成したデノイズ画像ａ〜ｈを、デノイズ画像情報２３６としてメモリ２３０に記憶させる。
【００５９】
デノイズ画像ａ〜ｈを生成する場合、デノイズ機能２１８は、パラメータ調整機能２１６により調整されたパラメータに基づいて、デノイズモデルＭＤＬ１に設定される内部パラメータの値を調整することができる。すなわち、デノイズ機能２１８は、原画像ａ〜ｈの各々のＳＮＲに基づいて調整された制御関数Ｇにより、各原画像のデノイズの度合いを可変とすることができる。デノイズモデルＭＤＬ１は、例えば、プロセッサがデノイズモデルＭＤＬ１を実行することによって、デノイズモデルＭＤＬ１がデノイズ機能２１８の一部として実現されてよい。デノイズ機能２１８は、磁気共鳴信号に対応するデータが入力されたときに該データからノイズが除去されたデータを出力するように学習されたモデルに基づいて、磁気共鳴信号に対応するデータからノイズを除去する。
【００６０】
なお、デノイズ機能２１８は、ニューラルネットワークを用いたものに限られない。デノイズ機能２１８は、例えば、ロジスティック回帰分析、決定木分析、サポートベクターマシンに基づく技術などの任意の機械学習により生成されたモデルを用いてノイズを除去または軽減するものであってよい。また、デノイズ機能２１８は、原画像に対するウェーブレット変換によりノイズを除去または軽減するものであってよい。
【００６１】
また、デノイズ機能２１８は、機械学習の技術を用いたものに限られない。デノイズ機能２１８は、例えば、画像信号を統計的に分析し、画像信号をノイズ除去空間に射影したときの信号のスパース性に基づいてノイズを含む画像信号からノイズが除去された信号を復元することで、デノイズ処理を行うようにしてもよい。また、デノイズ機能２１８は、任意のフィルタリング技術により、デノイズ処理を行うようにしてもよい。また、デノイズ機能２１８は、画像信号を平均化することや、画像信号における特徴量を抽出することにより、デノイズ処理を行うようにしてもよい。
【００６２】
画像合成機能２２０は、デノイズ機能２１８により生成されたデノイズ画像ａ〜ｈを合成した合成画像を生成する。画像合成機能２２０は、例えば、平方二乗和を用いた合成アルゴリズムなどを用いて、合成画像を生成する。図７に示すように、画像合成機能２２０は、複数のデノイズ画像ａ〜ｈを合成して１つの合成画像を生成する。この合成画像が、医用画像生成装置１００の最終的な撮像画像となる。なお、医用画像生成装置１００がパラレルイメージングによりＭＲ画像の再構成を行う装置である場合、画像合成機能２２０は、デノイズ画像ａ〜ｈに基づく再構成処理において画像の折り返しを展開する処理などを行ってよい。画像合成機能２２０は、生成した合成画像を、合成画像情報２３８としてメモリ２３０に記憶させる。
【００６３】
出力制御機能２２２は、画像合成機能２２０により合成画像が生成されると、通信インターフェース２０２に医用画像生成装置１００と通信させ、その通信相手の医用画像生成装置１００に、合成画像を送信する。また、出力制御機能２２２は、合成画像をディスプレイ２０６に表示させてもよい。
【００６４】
学習機能２２４は、デノイズモデルＭＤＬ１に対して、ある学習データとする画像（以下、学習画像）を入力し、デノイズモデルＭＤＬ１によって出力されたデノイズ画像が、教師データとする画像（以下、教師画像）に近づくようにデノイズモデルＭＤＬ１を学習する。例えば、教師画像は、サンプリング周波数を大きくするなどしてＳＮＲを高くした画像であってよい。学習画像は、教師画像に既知のノイズを加えた画像であってよい。既知のノイズは、例えば、ガウシアンノイズであってよい。
【００６５】
例えば、学習機能２２４は、デノイズモデルＭＤＬ１によって出力されたデノイズ画像と、教師画像との差分が小さくなるように、畳み込み層３２０の線形変換行列の要素値や、アクティベーション層３３０の各ノードの活性化関数の重み係数αといった種々のパラメータを、ＳＧＤ（Stochastic Gradient Descent）、Ｍｏｍｅｎｔｕｍ  ＳＧＤ、AｄａＧｒａｄ、ＲＭＳｐｒｏｐ、ＡｄａＤｅｌｔａ、Ａｄａｍ(Adaptive moment estimation)などの勾配法を利用して学習する。
【００６６】
なお、上記においては、処理回路２１０の各機能が、医用画像生成装置１００から受信した原画像ａ〜ｈに対してデノイズ処理を行い、デノイズ処理後の画像を合成することで合成画像を生成する場合を例に挙げて説明したがこれに限られない。例えば、処理回路２１０の各機能は、医用画像生成装置１００から受信したｋ空間データ等の再構成前のデータ（コイルエレメント１０８ａ〜１０８ｈにより出力される磁気共鳴信号に対応するデータ）に対してデノイズ処理を行い、デノイズ処理後のデータを合成することで合成画像を生成してもよい。
【００６７】
すなわち、導出機能２１４は、複数の受信コイルの各々により収集される磁気共鳴信号に対応するデータに含まれるノイズに関する指標値を導出する。パラメータ調整機能２１６は、導出機能２１４により導出された指標値に基づいて、磁気共鳴信号に対応するデータからノイズを除去する度合いを調整する。デノイズ機能２１８は、パラメータ調整機能２１６により調整された度合いに基づいて、磁気共鳴信号に対応するデータからノイズを除去する。画像合成機能２２０は、デノイズ機能２１８によりノイズが除去された、複数の受信コイルの各々で収集された磁気共鳴信号に対応するデータを合成する。磁気共鳴信号に対応するデータは、磁気共鳴信号に基づく複数の画像のデータを含む。
【００６８】
［処理フロー（学習処理）］
以下、第１の実施形態における処理回路２１０の処理フローについて説明する。処理回路２１０の処理には、デノイズモデルＭＤＬ１を学習する学習処理と、学習済みのデノイズモデルＭＤＬ１を使用してノイズを除去する処理を行う画像処理とが含まれる。以下においては、まず、処理回路２１０の学習処理について説明する。図８は、第１の実施形態における処理回路２１０の学習処理の一連の流れを示すフローチャートである。本フローチャートの処理は、例えば、医用画像処理装置２００の操作者が、入力インターフェース２０４を操作して学習処理の開始を指示した場合に行われる。
【００６９】
まず、学習機能２２４は、デノイズモデルＭＤＬ１に対する制御信号Ｇ、アクティベーション層３３０の各ノードの活性化関数の重み係数αなどの各種パラメータを初期値に設定する（ステップＳ１００）。例えば、制御信号Ｇおよび重み係数αを１に設定する。
【００７０】
次に、学習機能２２４は、学習画像をデノイズモデルＭＤＬ１に入力し、その処理結果を得る（ステップＳ１０２）。学習画像は、例えば、ノイズを含まない画像またはＳＮＲが高い画像である教師画像に対して、既知のノイズを付加した画像である。学習画像は、例えば、このノイズと教師画像の画素値とを加算することで得られる。
【００７１】
例えば、学習機能２２４は、デノイズモデルＭＤＬ１の入力層の各入力端に、学習画像の各画素の画素値を入力する。この画素値は、デノイズモデルＭＤＬ１の中間層の各ノードを、重み付き加算、バイアス加算、及び活性化関数処理の演算により、入力層から出力層に向かって、値を変化させながら伝搬していく。そして、出力層の出力端に、デノイズモデルＭＤＬ１による処理を受けた画像（以下、処理結果画像）の画素値が処理結果として出力される。
【００７２】
次に、学習機能２２４は、処理結果画像と教師画像との間の誤差である訓練誤差を算出する（ステップＳ１０４）。訓練誤差は、例えば、処理結果画像と教師画像の各画素の平均二乗誤差、各画素の二乗誤差の総和などである。
【００７３】
次に、学習機能２２４は、例えば、誤差逆伝播法を用いて、算出した訓練誤差が小さくなるようにデノイズモデルＭＤＬ１の内部パラメータを更新する（ステップＳ１０６）。次に、学習機能２２４は、例えば、教師画像とは別途準備されたノイズを含まない画像またはＳＮＲが高い画像である正解画像と、この正解画像にノイズを加えた検証画像とを用いて、検証画像をデノイズモデルＭＤＬ１に入力した場合の処理結果画像と、正解画像との間の誤差である汎化誤差を算出する（ステップＳ１０８）。次に、学習機能２２４は、汎化誤差が極小値に達したか否かを判定する（ステップＳ１１０）。
【００７４】
学習機能２２４は、汎化誤差が極小値に達していないと判定した場合、再度、ステップＳ１０２以降の処理を繰り返す。一方、学習機能２２４は、汎化誤差が極小値に達していると判定した場合、学習後のデノイズモデルＭＤＬ１の内部パラメータをメモリ２３０に保存する（ステップＳ１１２）。以上により、本フローチャートの処理が終了する。
【００７５】
［処理フロー（画像処理）］
次に、処理回路２１０の画像処理について説明する。図９は、第１の実施形態における処理回路２１０の画像処理の一連の流れを示すフローチャートである。本フローチャートの処理は、例えば、医用画像生成装置１００により送信された原画像ａ〜ｈが取得機能２１２によって取得された場合に行われる。
【００７６】
まず、導出機能２１４は、取得機能２１２によって取得された原画像ａ〜ｈの各々のＳＮＲを導出する（ステップＳ２００）。例えば、導出機能２１４は、取得機能２１２によって取得された原画像ａ〜ｈの各々に関して、傾斜磁場コイル１０２が傾斜磁場を発生させた状態で、ＲＦコイル１０８にＲＦパルスを供給し、高周波磁場を出力させた結果得られる原画像（以下、ＲＦパルスあり原画像と称する）と、傾斜磁場コイル１０２が傾斜磁場を発生させた状態で、高周波磁場を出力せずに得られる原画像（以下、ＲＦパルスなし原画像と称する）との画素値の差分に基づいて、ＳＮＲを導出する。
【００７７】
また、導出機能２１４は、取得機能２１２によって取得された原画像ａ〜ｈの各々に関して、同じ被検体ＯＢに向けて高周波磁場が出力された結果得られる２つ以上のＲＦパルスあり原画像の画素値の差分に基づいて、原画像のＳＮＲを導出してもよい。この際、導出機能２１４は、中心スライスに近い２つのＲＦパルスあり原画像の画素値の差分に基づいて、原画像のＳＮＲを導出してよい。例えば、本スキャンがスキャンを１０回繰り返すシーケンスである場合、中心スライスとは、５回目または６回目のスキャンによって得られる原画像である。
【００７８】
このように、被検体ＯＢが同じであるという条件の下、複数回にわたってスキャンを行って得られた２つの原画像の差分をとることによって、磁気共鳴信号成分を互いにキャンセルし、受信系の熱雑音に起因したランダムノイズ成分を基にしてＳＮＲを求めることができる。
【００７９】
次に、パラメータ調整機能２１６は、原画像ａ〜ｈの各々に関して、導出機能２１４により導出されたＳＮＲに基づいて制御関数Ｇを決定し、デノイズモデルＭＤＬ１のアクティベーション層３３０の少なくとも１つのノードに対して調整済みの制御関数Ｇを入力することで、デノイズモデルＭＤＬ１の内部パラメータである閾値Ｔの値を調整する（ステップＳ２０２）。
【００８０】
次に、デノイズ機能２１８は、原画像ａ〜ｈの各々に関して、パラメータ調整機能２１６により調整された閾値Ｔが内部パラメータとして設定されたデノイズモデルＭＤＬ１を利用してノイズを除去し、デノイズ画像ａ〜ｈを生成する（ステップＳ２０４）。
【００８１】
次に、画像合成機能２２０は、デノイズ機能２１８により生成されたデノイズ画像ａ〜ｈを合成して、最終的な撮像画像である合成画像を生成する（ステップＳ２０６）。
【００８２】
次に、出力制御機能２２２は、通信インターフェース２０２を制御して、画像合成機能２２０により生成された合成画像を医用画像生成装置１００に送信する（ステップＳ２０８）。この際、出力制御機能２２２は、ディスプレイ２０６に合成画像を表示させてもよい。医用画像生成装置１００の通信インターフェース１２２が医用画像処理装置２００から合成画像を受信した場合、医用画像生成装置１００の表示制御機能１３８は、合成画像をディスプレイ１２６に表示させてよい。また、出力制御機能２２２は、通信インターフェース２０２を制御して、画像の読影者が利用可能な端末装置に合成画像を送信してもよい。
【００８３】
以上説明した第１の実施形態によれば、コイルエレメント１０８ａ〜１０８ｈにより出力される磁気共鳴信号に基づく原画像ａ〜ｈの各々に対してデノイズ処理を行い、ノイズが除去されたデノイズ画像を合成して合成画像を生成することで、原画像ａ〜ｈのノイズ分布の違いに起因したデノイズ処理の精度低下を防ぐことができる。また、本実施形態のデノイズ処理では、原画像ａ〜ｈの各々に含まれるノイズの量に応じてノイズを除去する度合いを可変とすることができる。このデノイズ処理では、ノイズ以外の信号成分には影響が無い。このため、デノイズ処理の精度を高めるともに、最終的に得られる画像の品質を維持することができる。
【００８４】
（第１の実施形態の変形例）
以下、第１の実施形態の変形例について説明する。上述した第１の実施形態では、医用画像生成装置１００と医用画像処理装置２００とが互いに異なる装置であるものとして説明したがこれに限られない。例えば、医用画像処理装置２００は、医用画像生成装置１００のコンソール装置１２０の一機能によって実現されてもよい。すなわち、医用画像処理装置２００は、医用画像生成装置１００のコンソール装置１２０によって仮想的に実現される仮想マシンであってもよい。この場合、医用画像生成装置１００は、「医用情報処理装置」の一例である。
【００８５】
図１０は、第１の実施形態の変形例に係る医用画像生成装置１００を示す図である。図１０に示すように、コンソール装置１２０の処理回路１３０は、上述した取得機能１３２と、生成機能１３４と、通信制御機能１３６と、表示制御機能１３８と、に加えて、導出機能２１４と、パラメータ調整機能２１６と、デノイズ機能２１８と、画像合成機能２２０と、学習機能２２４とを実行してよい。また、コンソール装置１２０のメモリ１５０には、デノイズモデル情報２３２と、原画像情報２３４と、デノイズ画像情報２３６と、合成画像情報２３８と、が格納されてよい。
【００８６】
以上説明した第１の実施形態の変形例によれば、医用画像生成装置１００単体で、デノイズ処理の精度を高めることができる。
【００８７】
（第２の実施形態）
以下、第２の実施形態について説明する。上述した第１の実施形態では、複数のコイルエレメント１０８ａ〜１０８ｈにより出力される磁気共鳴信号に基づく原画像ａ〜ｈの各々に対してデノイズ処理を行い、ノイズが除去されたデノイズ画像を合成する構成を説明した。本実施形態では、医用画像処理装置２００の処理回路２１０が、デノイズ処理を行うことなく原画像を合成することにより生成された合成画像に対してノイズの強度分布を均一化させるためのノイズ成分を付加する処理を行い、処理後の合成画像に対してデノイズ処理を行う。このため、構成などについては第１の実施形態で説明した図および関連する記載を援用し、詳細な説明を省略する。
【００８８】
［医用画像処理装置の構成例］
図１１は、第２の実施形態に係る医用画像処理装置２００の一例を示す図である。医用画像処理装置２００の処理回路２１０は、取得機能２１２と、導出機能２１４と、パラメータ調整機能２１６と、デノイズ機能２１８と、画像合成機能２２０と、出力制御機能２２２と、学習機能２２４とに加えて、ノイズ均一化機能２２６を実行する。ノイズ均一化機能２２６は、「均一化処理部」の一例である。
【００８９】
ノイズ均一化機能２２６は、画像合成機能２２０によって原画像を合成することにより生成された合成画像からノイズの強度分布を推定し、推定されたノイズ強度分布から合成画像のノイズ強度分布を均一化させるためのノイズ分布（以下、均一化分布）を生成し、均一化分布により示される信号を合成画像に対して付加することで、合成画像に含まれるノイズの強度分布を均一化する。そして、デノイズ機能２１８は、ノイズ均一化機能２２６によりノイズの強度分布が均一化された合成画像に対してデノイズ処理を行う。
【００９０】
［処理フロー（画像処理）］
次に、処理回路２１０の画像処理について説明する。図１２は、第２の実施形態における処理回路２１０の画像処理の一連の流れを示すフローチャートである。本フローチャートの処理は、例えば、医用画像生成装置１００により送信された原画像ａ〜ｈが取得機能２１２によって取得された場合に行われる。
【００９１】
まず、画像合成機能２２０は、取得機能２１２によって取得された原画像ａ〜ｈを合成して合成画像を生成する（ステップＳ３００）。すなわち、画像合成機能２２０は、複数の受信コイル（複数のコイルエレメント）の各々により被検体ＯＢから収集される磁気共鳴信号に基づく複数の原画像ａ〜ｈを合成して合成画像を生成する。なお、医用画像生成装置１００のコンソール装置１２０の処理回路１３０により実行される生成機能１３４が原画像ａ〜ｈを合成して合成画像を生成し、該合成画像を画像合成機能２２０に送信している場合、ステップＳ３００の処理は行われなくてよい。
【００９２】
次に、ノイズ均一化機能２２６は、画像合成機能２２０により生成された合成画像のノイズ強度分布を推定する（ステップＳ３０２）。例えば、ノイズ均一化機能２２６は、
取得機能２１２によって取得された原画像a〜hそれぞれの画像と、原画像a〜hを合成した合成画像との比較に基づいて、合成画像のノイズ強度分布を推定する。
【００９３】
次に、ノイズ均一化機能２２６は、推定したノイズ強度分布に基づいて、画像合成機能２２０により生成された合成画像のノイズ強度分布を均一化させるための均一化分布を生成する（ステップＳ３０４）。ノイズ均一化機能２２６は、生成した均一化分布をノイズ均一化情報としてメモリ２３０に記憶させてよい。次に、ノイズ均一化機能２２６は、生成した均一化分布を、画像合成機能２２０により生成された合成画像に付加することで、合成画像のノイズ強度分布を均一化する（ステップＳ３０６）。
【００９４】
次に、デノイズ機能２１８は、デノイズモデル情報２３２が示すデノイズモデルＭＤＬ１を利用して、ノイズ均一化機能２２６によりノイズ強度分布が均一化された合成画像からノイズを除去するデノイズ処理を行う（ステップＳ３０８）。
【００９５】
次に、出力制御機能２２２は、通信インターフェース２０２を制御して、デノイズ機能２１８によりデノイズ処理が行われた合成画像を医用画像生成装置１００に送信する（ステップＳ３１０）。この際、出力制御機能２２２は、ディスプレイ２０６に合成画像を表示させてもよい。医用画像生成装置１００の通信インターフェース１２２が医用画像処理装置２００から合成画像を受信した場合、医用画像生成装置１００の表示制御機能１３８は、合成画像をディスプレイ１２６に表示させてよい。また、出力制御機能２２２は、通信インターフェース２０２を制御して、画像の読影者が利用可能な端末装置に合成画像を送信してもよい。
【００９６】
以上説明した第２の実施形態によれば、原画像を合成することにより生成された合成画像に対してノイズの強度分布を均一化させるためのノイズ成分を付加する処理を行い、処理後の合成画像に対してデノイズ処理を行うことで、デノイズ処理の精度を高めることができる。
【００９７】
（第２の実施形態の変形例）
以下、第２の実施形態の変形例について説明する。上述した第２の実施形態では、医用画像生成装置１００と医用画像処理装置２００とが互いに異なる装置であるものとして説明したがこれに限られない。例えば、医用画像処理装置２００は、医用画像生成装置１００のコンソール装置１２０の一機能によって実現されてもよい。すなわち、医用画像処理装置２００は、医用画像生成装置１００のコンソール装置１２０によって仮想的に実現される仮想マシンであってもよい。この場合、医用画像生成装置１００は、「医用情報処理装置」の一例である。
【００９８】
図１３は、第２の実施形態の変形例に係る医用画像生成装置１００を示す図である。図１３に示すように、コンソール装置１２０の処理回路１３０は、上述した取得機能１３２と、生成機能１３４と、通信制御機能１３６と、表示制御機能１３８と、に加えて、導出機能２１４と、パラメータ調整機能２１６と、デノイズ機能２１８と、画像合成機能２２０と、学習機能２２４と、ノイズ均一化機能２２６とを実行してよい。また、コンソール装置１２０のメモリ１５０には、デノイズモデル情報２３２と、原画像情報２３４と、デノイズ画像情報２３６と、合成画像情報２３８と、ノイズ均一化情報２４０とが格納されてよい。
【００９９】
以上説明した第２の実施形態の変形例によれば、医用画像生成装置１００単体で、デノイズ処理の精度を高めることができる。
【０１００】
上記説明したいずれかの実施形態は、以下のように表現することができる。
プログラムを格納するストレージと、
プロセッサと、を備え、
前記プロセッサは、前記プログラムを実行することにより、
複数の受信コイルの各々により収集される磁気共鳴信号に対応するデータに含まれるノイズに関する指標値を導出し、
導出された前記指標値に基づいて、前記磁気共鳴信号に対応するデータからノイズを除去する度合いを調整し、
調整された前記度合いに基づいて、前記磁気共鳴信号に対応するデータからノイズを除去し、
前記ノイズが除去された、前記複数の受信コイルの各々で収集された前記磁気共鳴信号に対応するデータを合成する、
医用情報処理装置。
【０１０１】
上記説明したいずれかの実施形態は、以下のように表現することができる。
プログラムを格納するストレージと、
プロセッサと、を備え、
前記プロセッサは、前記プログラムを実行することにより、
複数の受信コイルの各々により被検体から収集される磁気共鳴信号に基づく複数の画像を合成して合成画像を生成し、
生成された前記合成画像に含まれるノイズの強度分布を推定し、
推定された前記ノイズの強度分布を均一化させる均一化分布を生成し、
生成された前記均一化分布により示される信号を前記合成画像に付加することで、前記合成画像に含まれるノイズの強度分布を均一化し、
前記ノイズの強度分布が均一化された合成画像からノイズを除去する、
医用情報処理装置。
【０１０２】
以上説明した少なくともひとつの実施形態によれば、処理回路２１０が、複数のＲＦコイル１０８の各々により被検体から収集される磁気共鳴信号に基づく複数の画像に含まれるノイズに関する指標値を、画像ごとに導出する導出機能２１４と、導出機能２１４により導出された指標値に基づいて、複数の画像の各々からノイズを除去する度合いを調整するパラメータ調整機能２１６と、パラメータ調整機能２１６により調整された度合いに基づいて、複数の前記画像の各々からノイズを除去するデノイズ機能２１８と、を備えることで、デノイズ精度を高めることができる。
【０１０３】
以上説明した少なくともひとつの実施形態によれば、処理回路２１０が、複数のＲＦコイル１０８の各々により被検体から収集される磁気共鳴信号に基づく複数の画像を合成して合成画像を生成する画像合成機能２２０と、画像合成機能２２０により生成された合成画像に含まれるノイズの強度分布を推定し、推定されたノイズの強度分布を均一化させる均一化分布を生成し、生成された均一化分布により示される信号を合成画像に付加することで、合成画像に含まれるノイズの強度分布を均一化するノイズ均一化機能２２６と、ノイズ均一化機能２２６によりノイズの強度分布が均一化された合成画像からノイズを除去するデノイズ機能２１８と、を備えることで、デノイズ精度を高めることができる。
【０１０４】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１０５】
１…医用画像処理システム、１００…医用画像生成装置、１０１…静磁場磁石、１０２…傾斜磁場コイル、１０４…寝台、１０５…寝台制御回路、１０７…送信回路、１０８…ＲＦコイル、１０８ａ〜ｈ…コイルエレメント、１０９…受信回路、１１０…シーケンス制御回路、１２０…コンソール装置、２００…医用画像処理装置、２０２…通信インターフェース、２０４…入力インターフェース、２０６…ディスプレイ、２１０…処理回路、２１２…取得機能、２１４…導出機能、２１６…パラメータ調整機能、２１８…デノイズ機能、２２０…画像合成機能、２２２…出力制御機能、２２４…学習機能、２２６…ノイズ均一化機能、２３０…メモリ
【図１】
2020103890.tif 000005
【図２】
2020103890.tif 000006
【図３】
2020103890.tif 000007
【図４】
2020103890.tif 000008
【図５】
2020103890.tif 000009
【図６】
2020103890.tif 000010
【図７】
2020103890.tif 000011
【図８】
2020103890.tif 000012
【図９】
2020103890.tif 000013
【図１０】
2020103890.tif 000014
【図１１】
2020103890.tif 000015
【図１２】
2020103890.tif 000016
【図１３】
2020103890.tif 000017