JP7320421B2 - 磁気共鳴イメージング装置、及び、自動撮像位置設定方法 - Google Patents

磁気共鳴イメージング装置、及び、自動撮像位置設定方法 Download PDF

Info

Publication number
JP7320421B2
JP7320421B2 JP2019185131A JP2019185131A JP7320421B2 JP 7320421 B2 JP7320421 B2 JP 7320421B2 JP 2019185131 A JP2019185131 A JP 2019185131A JP 2019185131 A JP2019185131 A JP 2019185131A JP 7320421 B2 JP7320421 B2 JP 7320421B2
Authority
JP
Japan
Prior art keywords
imaging
measurement
cross
magnetic resonance
dimensional
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
Application number
JP2019185131A
Other languages
English (en)
Other versions
JP2021058447A (ja
Inventor
俊 横沢
亨 白猪
尚子 永尾
周子 森分
彬栄 呉
Original Assignee
富士フイルムヘルスケア株式会社
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 富士フイルムヘルスケア株式会社 filed Critical 富士フイルムヘルスケア株式会社
Priority to JP2019185131A priority Critical patent/JP7320421B2/ja
Priority to CN202010472767.XA priority patent/CN112637478B/zh
Priority to US16/889,987 priority patent/US11255939B2/en
Publication of JP2021058447A publication Critical patent/JP2021058447A/ja
Application granted granted Critical
Publication of JP7320421B2 publication Critical patent/JP7320421B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/60—Control of cameras or camera modules
    • H04N23/695—Control of camera direction for changing a field of view, e.g. pan, tilt or based on tracking of objects
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/546—Interface between the MR system and the user, e.g. for controlling the operation of the MR system or for the design of pulse sequences
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • 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/10088—Magnetic resonance imaging [MRI]

Description

本発明は、磁気共鳴イメージング装置（以下、ＭＲＩ装置という）に係り、特にその撮像位置を自動設定する技術に関する。
ＭＲＩ装置は、被検体の任意の方向の断面を撮像することができるため、診断に適した断面画像を得ることができる。一般に、ＭＲＩにおいて撮像断面を設定する場合、まず位置決めのためにアキシャル、コロナル、サジタルの３つの直交断面の画像を位置決め画像（スカウト画像ともいう）を取得し、スカウト画像上で解剖組織を目印として撮像位置を決定する。撮像位置の決定は、従来はユーザーインターフェイスを介して手動操作で行われていたが、スカウト画像の画像処理により自動的に撮像位置を提示する機能（自動位置決め機能という）が開発されている（例えば、特許文献１）。
自動位置決め機能による撮像断面の決定は、スカウト画像に現われる解剖組織の特徴が画像処理によって比較的容易に抽出できるものである場合には、精度よく行うことができるが、複雑な三次元構造を持つ肩や膝などの関節では、一回取得したスカウト画像だけでは撮像位置の目印となる解剖組織の同定が難しい。このため、一回目のスカウト画像から二回目のスカウト画像の位置を決定し、二回目のスカウト画像で診断用の画像の撮像位置を決定する。
決定した撮像位置は、それをユーザーが確認するために表示装置等により提示される。この際、ユーザーが提示された撮像位置が許容できるか否かを判別する必要があるため、撮像位置は、それを決定した際の目印となる解剖組織が同定可能な断面上に提示される。例えば、肩検査であれば、撮像位置を提示する画像として、肩甲骨と上腕骨が同定できるアキシャル面画像、上腕骨の傾きが同定できるオブリークのサジタル面画像、及び肩甲骨と上腕骨の関節面が同定できるオブリークのコロナル面画像が必要となる。
このような要請に対し、肩関節のＭＲ検査において、３次元のスカウト画像を取得し、３Ｄ Ｄｅｆｏｒｍａｂｌｅ ｍｅｓｈ ａｄａｐｔｉｏｎを用いたモデルベースのアルゴリズムで目印を検知し、撮像断面を自動位置決めする技術（非特許文献１）や、膝関節の検査において、３次元のスカウト画像を取得し、複数個所の異なる解剖構造の位置関係を検出することで誤検出を低減するとともに、解剖学的構造の不変性を抽出する階層学習を実施して、個々の膝の向きや曲がりなどの違いに起因する誤差を低減する技術（非特許文献２）が提案されている。
特開２０１４－１２１５９８号公報
Ｃ．Ｊ． ｄｅｎ Ｈａｒｄｅｒら，"Ｃｏｎｓｉｓｔｅｎｔ ａｕｔｏｍａｔｅｄ ｓｃａｎ ｐｌａｎｎｉｎｇ ｏｆ ｓｈｏｕｌｄｅｒ" Ｐｒｏｃ．Ｉｎｔｌ．Ｓｏｃ．Ｍａｇ．Ｒｅｓｏｎ．Ｍｅｄ． １９（２００８），ｐ３６６５ Ｙｉｑｉａｎｇ Ｚｈａｎ，"Ｒｏｂｕｓｔ Ａｕｔｏｍａｔｉｃ Ｋｎｅｅ ＭＲ Ｓｌｉｃｅ Ｐｏｓｉｔｉｏｎｎｉｎｇ Ｔｈｒｏｕｇｈ Ｒｅｄｕｎｄａｎｔ ａｎｄ Ｈｉｅｒａｒｃｈｉｃａｌ Ａｎａｔｏｍｙ Ｄｅｔｅｃｔｉｏｎ"
上述したように従来の技術では、解剖学的構造が複雑な組織について自動位置決めの精度を高めるために３次元のスカウト画像を取得しているが、３次元撮像は撮像時間が長く、しかも撮像断面特定に必要ではない情報を多く含む。即ち過剰な情報を取得している。また決定した撮像断面を提示した後で、それが不適当と判断された場合、修正に手間がかかり、時間の長い３次元撮像を再度行わなければならない場合も生じる。
そこで、本発明は解剖学的構造が複雑な組織について、短時間で精度の高い自動位置決めが可能な技術を提供することを課題とする。
本発明は、複雑な組織の断面を自動化する際に必要となる処理には、左右を判定する処理と解剖学的特徴を抽出する処理が含まれ、これら処理に必要な情報、情報量が異なることに着目し、少ない情報で行う処理については短時間で取得できる計測を行うことより、計測する情報の取得効率を高め、精度を保ちながら自動位置決めの時間を短縮するものである。
すなわち、本発明のＭＲＩ装置は、核磁気共鳴を用いて被検体の断面の画像を取得する撮像部と、撮像位置を決定するスカウト画像を取得するスカウト撮像、及び、診断用の画像を取得する本撮像を前記撮像部が行うよう制御する撮像制御部と、前記スカウト撮像で取得したデータを用いて、前記本撮像の撮像断面位置を決定する撮像断面決定部と、を備え、前記撮像制御部は、前記スカウト撮像として、１次元又は２次元の計測データを取得する第１計測、及び、２次元の計測データを取得する第２計測を実行する制御を行い、前記撮像断面決定部は、前記第１計測で取得した計測データを用いて、前記被検体の左右を判定する左右判定部と、前記左右判定における判定結果と前記第１計測で取得した計測データと用いて前記第２計測における断面位置を算出するとともに、前記第２計測で取得した計測データを用いて前記本撮像における断面位置を算出する断面位置算出部と、を備える。
また本発明のＭＲＩ装置における自動撮像位置設定方法（自動位置決め方法）は、次のステップを含む。即ち、診断用画像を取得する本撮像に先立って、スカウト撮像の第１計測を実行し、１次元又は２次元の計測データを取得し、前記第１計測で取得した計測データを用いて、前記被検体の左右を判定し、前記左右判定における判定結果と前記第１計測で取得した計測データと用いて、スカウト撮像の第２計測における断面位置を算出し、当該断面位置にて前記第２計測を実行し、２次元の計測データを取得し、前記第２計測で取得した２次元の計測データを用いて前記本撮像における断面位置を算出する。
本発明によれば、複雑な組織の特徴抽出において前提となる左右判定を、１次元或いは２次元の計測データで行うことにより、冗長な計測データの取得を不要とし、全体としての撮像時間を短縮することができる。また左右判定に基づいて、それに続くスカウト撮像の断面位置を算出することにより、２次元のスカウト画像から撮像断面位置設定に必要な精度のよい情報を得ることができる。結果として、３次元データを取得する場合と同等の精度を保ちながら短時間で自動位置決めを行うことができる。
本発明のＭＲＩ装置の一実施形態の全体構成を示すブロック図。 本発明が適用されるＭＲＩ装置の一例を示すブロック図。 自動位置決め処理手順の概要を示すフロー図。 ＭＰＲ処理を説明する図。 実施形態１の自動位置決め機能を実行する計算機の機能ブロック図。 実施形態１の自動位置決めの手順の一例を示すフロー図で、２回目のスカウト撮像の断面位置決定までのフローを示す。 実施形態１の自動位置決めの手順の一例を示すフロー図で、本撮像の断面位置決定までのフローを示す。 １回目のスカウト撮像で得た計測データを示す図で、（ａ）は２次元データである肩部のアキシャル画像を示す図、（ｂ）は（ａ）の画像の一次元投影を示す図。 物体抽出部の一例として、Ａｄａｐｔｉｖｅ ｂｏｏｓｔｉｎｇ法のアルゴリズムを示す図。 １回目のスカウト画像に対する特徴抽出部の処理を説明する図で、（ａ）はマッチング処理、（ｂ）はその要部拡大図、（ｃ）は特徴抽出後に決定される撮像断面位置を示す図。 対称線検出部の処理を説明する図。 ２回目のスカウト画像及びＭＰＲ画像に対する特徴抽出部の処理を説明する図で、（ａ）はスカウトＡＸ画像、（ｂ）はＭＰＲ処理後のＳＡＧ画像、（ｃ）はＭＰＲ処理後のＣＯＲ画像を示す。 本撮像の撮像断面位置の提示例を示す図で、（ａ）はＡＸ面、（ｂ）はＳＡＧ面、（ｃ）はＣＯＲ面である。 実施形態１と従来技術との効果の差を説明する図。 実施形態２の自動位置決めの手順の一例を示すフロー図。 ２回目のスカウト画像及びＭＰＲ画像に対する特徴抽出部の処理を説明する図で、（ａ）はスカウトＡＸ画像、（ｂ）はＭＰＲ処理後のＳＡＧ画像、（ｃ）はＭＰＲ処理後のＣＯＲ画像を示す。 本撮像の撮像断面位置の提示例を示す図で、（ａ）はＡＸ面、（ｂ）はＳＡＧ面、（ｃ）はＣＯＲ面である。
以下、本発明のＭＲＩ装置の実施形態を、図面を参照して説明する。
図１に示すように、ＭＲＩ装置１００は、核磁気共鳴を用いて被検体の断面の画像を取得する撮像部１０と、撮像部１０が取得した画像に対し種々の演算を施す演算部２０と、撮像部１０や演算部２０の動作を制御する制御部３０とを備えている。
撮像部１０の構成は、公知のＭＲＩ装置の撮像部の構成と同様である。具体的には、図２に示すように、ＭＲＩ装置１００は、撮像部として、被検体（その検査部位）１０２が配置される空間に静磁場を発生させる静磁場磁石１０１、被検体１０２に対し高周波磁場を印加する送信用ＲＦコイル１０３、静磁場の均一度を高めるためのシムコイル１０４、静磁場空間に３軸方向の傾斜磁場を発生させる傾斜磁場コイル１０５、及び被検体１０２から発生する核磁気共鳴信号（高周波信号）を検出する受信用ＲＦコイル１０６などを備えている。また被検体１０２を静磁場空間に配置するための寝台１０７が設置されている。
送信用ＲＦコイル１０３は、高周波信号を発生する高周波発生部１１３に、シムコイル１０４はシム電源１１４に、傾斜磁場コイル１０５は傾斜磁場電源１１５に、それぞれ接続されている。また受信用ＲＦコイル１０６は、核磁気共鳴信号の検波及びデジタルへの変換等を行う受信器１１６に接続されている。高周波発生部１１３、シム電源１１４、傾斜磁場電源１１５、受信器１１６は、計算機１１０の指令のもとで動作するシーケンサ１１７からの制御信号により動作する。計算機１１０には、付属装置としてディスプレイ（表示装置）１１１やマウスやキーボード等の入力装置１１２及び記憶装置１１８などが備えられており、傾斜磁場や高周波パルスのパルス形状や印加タイミングなどを予め定めパルスシーケンスを記憶装置１１８から読み出すとともに、入力装置１１２を介してユーザーが設定した撮像条件（撮像断面、エコー時間、繰り返し時間、繰り返し回数）などをシーケンサ１１７に設定する。これらより撮像部１０によって所定の撮像或いは計測が実行される。撮像断面の方向は、３軸方向の傾斜磁場の組み合わせによって、位相エンコード方向、スライス方向、及び読み出し方向が決まり、それらを軸方向とする互いに直交する３つの断面（例えばアキシャル面、コロナル面、サジタル面）が決まる。
撮像部１０が実行する撮像には、診断用の画像を取得するための撮像（本撮像）のほかに、撮像時の断面位置（撮像位置）を決定のための位置決め撮像（スカウト撮像と呼ぶ）が含まれ、スカウト撮像には２つの計測、即ち、検査対象の左右を判定するために必要なデータを収集する第１計測（１回目のスカウト撮像）と、検査対象の組織や部位（まとめて組織という）の特徴を抽出するために必要なデータを収集する第２計測（２回目のスカウト撮像）とが含まれる。
演算部２０が行う演算には種々の演算が含まれるが、本実施形態のＭＲＩ装置においては、位置決め撮像で取得した計測データを用いて、画像取得のための撮像断面を決定する撮像断面決定部２００が含まれ、撮像断面決定部２００には、１回目のスカウト撮像で取得したデータ（第１計測データ）を用いて検査対象の左右を判定する左右判定部２１０と、スカウト撮像で取得したデータを用いて検査対象部位の特徴抽出を行う特徴抽出部２３０と、２回目のスカウト撮像及びその後に行われる本撮像における撮像断面の位置を算出する断面位置算出部２５０と、を備える。
制御部３０は、予め定められた検査プロトコル（検査対象に応じて撮像シーケンスの種類やその実行順序等を定めたもの）やユーザーの設定や指令に従って、撮像部１０を制御する撮像制御部３１０と、演算部２０の演算結果の提示やユーザーインターフェイスを介したユーザー設定の受付などを行う表示制御部３３０とを備える。撮像制御部３１０は、検査プロトコルの検査対象が膝関節や肩関節などの複雑な構造を持つ組織である場合、位置決め撮像として、撮像部１０に２回のスカウト撮像を実行させる。また２回目のスカウト撮像の結果を用いて撮像断面決定部２００が提示した撮像断面がユーザーに承認された場合、撮像部１０にその断面で本撮像を実行させる。
演算部２０及び制御部３０の機能は、予め定めたプログラムを、同一の或いは別個の計算機の中央処理装置（ＣＰＵ）で実行することにより実現することが可能である。また演算部２０の機能の一部は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）やＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）などのハードウェアで実現してもよい。図２に示す構成では、計算機１１０がその機能を実現する。
本実施形態のＭＲＩ装置は、検査対象が特定されると、予め定めた手順に従って検査対象の撮像断面を決定するための位置決め撮像を自動的に行い、本撮像の撮像断面を決定し、ユーザーに提示する。
以下、本実施形態による撮像断面決定の処理の流れを、図３を参照して説明する。
まず被検体１０２の検査部位の周囲に受信用ＲＦコイル１０６を取り付け、被検体１０２を静磁場空間内に配置し、撮像制御部３１０は、撮像部１０を制御し、１回目のスカウト撮像を実行し、第１計測データを取得する（Ｓ３１）。１回目のスカウト撮像は、例えば、１回のＲＦパルスを照射後に読み出し傾斜磁場を印加しながらエコー信号を収集する計測であり、このエコー信号を１次元フーリエ変換することによって読み出し方向の信号強度（軸方向への投影値）を示す１次元のデータが得られる。この際、読み出し方向（例えばＸ方向とする）は、被検体１０２の左右方向と概ね一致させることが好ましい。第２の計測は、例えば、少ない位相エンコード数のエコー信号を計測する。このエコー信号を２次元フーリエ変換することにより、第２計測データとして２次元画像データが得られる。この計測では読み出し方向は被検体の左右方向と一致させなくともよい。
撮像断面決定部２００は、撮像部１０から第１計測データを受け取り、まず左右判定部２１０が第１計測データから得られる読み出し方向の信号強度をもとに、画像の左右を判定する（Ｓ３２）。ついで断面位置算出部２５０が左右判定部２１０の判定結果をもとに２回目のスカウト撮像における撮像断面位置を算出する（Ｓ３３）。２回目のスカウト撮像では、アキシャル面画像（ＡＸ画像）、コロナル面画像（ＣＯＲ画像）及びサジタル面画像（ＳＡＧ画像）を取得するので、これら三断面について、それぞれ撮像位置を算出する。
具体的には、第１計測データが１次元データである場合には、左右判定結果を元に、２回目のスカウト撮像の一つの断面、例えばＡＸ面の断面位置を決定する。そして、この断面位置で取得した２次元データ（例えばＡＸ画像）から検査対象である部位の特徴（位置、傾き、形状など）を抽出し、その特徴に基づいて他の断面（ＳＡＧ面、ＣＯＲ面）の位置や傾き（オブリーク面）を決定する。また第１計測データが例えば一つの断面（ＡＸ面）の２次元データである場合には、この２次元データを用いて左右判定を行うとともに、この２次元の画像データから抽出した検査対象部位の特徴を用いて、２回目のスカウト撮像における当該一つの断面を含む複数の断面の断面位置を決定する。
２回目のスカウト撮像における撮像断面位置が決定されると、撮像制御部３１０は決定された位置に対応する傾斜磁場印加条件で２回目のスカウト撮像を行うように撮像部１０を制御する（Ｓ３４）。２回目のスカウト撮像では、所定位置のＡＸ画像と、組織の構造に沿ったオブリークのＳＡＧ画像及びオブリークのＣＯＲ画像を撮像断面として、２次元の画像を得る撮像を行う。通常各断面の撮像では複数の平行な面の画像（複数のスライスデータ）を取得する。特に、サジタル面及びコロナル面については、その後にＭＰＲ（Ｍｕｌｔｉ－Ｐｌａｎａｒ Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ）処理を行うために、それぞれ、複数のスライスデータを取得する。
撮像断面決定部２００は、撮像部１０から２回目のスカウト撮像の画像データ（第２計測データ）を受け取り、本撮像の撮像断面を決定する（Ｓ３５）。そのため、まず特徴抽出部２３０が、第２計測データに含まれる複数の画像から検査対象の組織の特徴を抽出し、画像の中心とすべき組織の位置や、複数のサジタル面画像及び複数のコロナル面画像から画像を切り出すための組織の傾きを算出する。この情報をもとに、複数のサジタル面画像及び複数のコロナル面画像に対し、それぞれＭＰＲ処理を行い特徴が含まれる面（切り出し面）を決定する。この面は、図４に示すように、例えば、複数のサジタル面画像から切り出した面であり、サジタル面に対し所定の傾きを持つオブリーク面である。ついで、切り出し面において検査対象の組織の特徴を用いて、その傾きを算出する。このコロナル面についても、検査対象の組織の特徴を用いて、その傾きを算出する。以上の処理によって求められた検査対象となる組織の傾きや位置をもとに、例えば所定の位置が画像中心となり且つ検査対象組織の傾きと平行となるように、本撮像の撮像断面を決定する。
撮像断面決定部２００が３つの断面を決定すると、表示制御部３３０はその情報を受け取り、ユーザーに提示する（Ｓ３６）。ユーザーへの提示方法は、例えば、ディスプレイに三断面の画像上に、推奨される撮像断面位置を重畳した画面を表示する。同じ画面に、例えば、その推奨撮像断面にて本撮像を開始するか否かのユーザー指示やユーザー調整を入力するユーザーインターフェイスを表示してもよい。ユーザーが本撮像の開始を入力した場合には（Ｓ３７）、本撮像を開始する（Ｓ３８）。ユーザー調整の場合には、提示した撮像断面位置の変更を受付け（Ｓ３９）、ユーザー調整後に本撮像を開始する。これらの表示或いはユーザーインターフェイスの表示態様は種々の態様を取ることができ、後述の実施形態において、いくつかの例を詳述する。
本実施形態のＭＲＩ装置によれば、自動位置決めの際に前提となる左右判定を一次元計測あるいは二次元計測で行い、その結果を利用して、本撮像断面位置決定に必要なデータを二次元計測で取得するので、２回目のスカウト撮像を含むスカウト撮像時間を大幅に低減することができる。また左右判定結果に基づいて２回目のスカウト撮像を行うことにより、二次元計測のスカウト撮像で得た画像における特徴抽出とそれに基づく撮像位置算出の精度を向上することができる。
以下、具体的な検査部位について、本発明の自動位置きめ技術を適用した実施形態について説明する。
［実施形態１］
本実施形態は、検査部位が肩関節である場合の自動位置決め処理の実施形態である。ここでは、撮像断面として、肩関節を中心として、上腕骨を含む三断面を自動設定する場合を説明する。
本実施形態においても装置の構成は、図１及び図２に示す構成と同様であり、また基本的な処理の流れは図３と同様である。本実施形態における計算機１１０（演算部２０及び制御部３０）の自動位置決めに関する詳細な構成を図５に示す。図５において、図１に示す要素と同じ要素は同じ符号で示し、重複する説明は省略する。
図５に示すように、撮像断面決定部２００の特徴抽出部２３０は、特定の組織を同定する処理を行うマッチング処理部２３１、特定の組織を抽出する物体抽出部２３３、抽出された物体について左右対称線を検出する対称線検出部２３５、複数のスライスデータからＭＰＲ処理により１枚の画像（ＭＰＲ画像）を切り出す画像切り出し部２３７を備える。また制御部３０の機能として、表示装置１１１を介して、撮像断面位置に関するユーザー変更を受付けるとともに変更内容を実行する撮像位置調整／実行部３２０を備える。
以下、撮像断面決定部２００の各部の処理を中心に、本実施形態における撮像断面決定処理を、図６を参照して説明する。図６は、処理の流れを示す図で、図６Ａは、１回目のスカウト撮像から２回目のスカウト撮像用の撮像位置決定まで、図６Ｂは、２回目のスカウト撮像から本撮像用の撮像位置決定までを示している。
＜Ｓ６１：１回目のスカウト撮像＞
撮像制御部３１０の制御のもと、撮像部１０は１回目のスカウト撮像を実行し、スカウト画像を取得する。この実施形態では、対象領域を選択して二次元撮像を行い、１ないし複数のスライスのＡＸ画像を取得するものとする。
＜Ｓ６２：左右判定処理＞
＜＜Ｓ６２－１：一次元投影画像作成＞＞
左右判定部２１０は、スカウト撮像で取得したＡＸ画像を入力とし、まず、それを体軸の左右軸に対応するｘ軸に投影した一次元投影画像を作成する。ＡＸ画像が複数のスライスデータからなる場合には、複数のスライスデータをそれぞれｘ軸に投影して得られる一次元投影画像を加算して、左右判定のための一次元投影画像とする。
＜＜Ｓ６２－２：左右判定＞＞
次いで左右判定部２１０は、一次元投影画像から検査部位が右肩か左肩かを判定する。例えば、右肩関節を検査する場合、受信用ＲＦコイルは右肩に装着されており、その状態で取得した画像は、図７（ａ）に示すように、右肩からの信号値が高く、それ以外の領域からの信号値は低い。即ち画像の左右方向（ｘ軸方向）において信号値に偏りがある。これを利用して、一次元投影像（図７（ｂ））のｘ方向の信号の偏りから左右判定を行う。具体的には、信号値に閾値ＳＴを設定し、信号値Ｓ（ｘ）が閾値ＳＴより高くなる位置（ｘ）の平均値ｍＡ（位置の平均値）を算出する。閾値ＳＴは、例えば、判別分析法等の手法により算出する。或いは単純に最大値と最小値との中央値を閾値としてもよい。
次に信号強度Ｓ（ｘ）が閾値ＳＴより低くなる位置（ｘ）の信号値の重み付き平均値ｍＢ（位置の平均値）を算出する。重みｗｉは各位置ｘにおける信号強度Ｓ（ｘ）の規格値とする。
ｍＢ＝Σ（ｗｉ×ｘｉ）／Ｎ
Ｎ：信号強度Ｓ（ｘ）が閾値ＳＴより低くなる位置（ｘ）の総数
ｗｉ＝Ｓ（ｘｉ）／｛ΣＳ（ｘｉ）／Ｎ｝
最後に平均値ｍＡとｍＢとを比較し、ｍＡ＞ｍＢのときは左肩、ｍＡ＜ｍＢのときは右肩と判定する。なお平均値ｍＡ、ｍＢの代わりに、信号値Ｓ（ｘ）が平均値を取る位置を用いる、或いはｍＡについては信号値がピーク値となる位置を用いるなど変更が可能である。撮像断面決定部２００は、判定結果をディスプレイに表示してもよい。これにより、ユーザーは、画像からは直ちに左右が判別できない場合にも、この表示によって左右のいずれが撮像されているか、検査プロトコルの対象が正しく撮像対象になっているかを確認することができる。
＜Ｓ６３：特徴抽出＞
左右判定後、特徴抽出部２３０は、１回目のスカウト画像（ＡＸ画像）を用いて、対象となる部位、ここでは肩関節、の特徴的な構造の探索と特徴抽出を行う。この探索の範囲は、例えば、左右判定により信号値が閾値より高いと判断された領域とする。
＜＜Ｓ６３－１：物体抽出＞＞
このステップでは、まず、物体抽出部２３３が、Ａｄａｐｔｉｖｅ ｂｏｏｓｔｉｎｇ法など公知の物体抽出方法を用いて、上腕骨の関節面を探索する。Ａｄａｐｔｉｖｅ ｂｏｏｓｔｉｎｇ法は、図８に示すように、比較的簡易な機械学習アルゴリズムを用いた物体抽出法の一つであり、入力データｘの特徴量から、各識別器ｈ１（ｘ）～ｈＮ（ｘ）にて正解か否かを判定し、最終的な判別を弱識別器の重み付き多数決で決定する。学習の過程で、入力データｘと予め正解と不正解の分かっている出力データｙを用いて、正解と不正解をもっとも精度よく分離する各各識別器の重みα１～αＮを算出する。上腕骨の関節面画像を抽出したい場合には、複数の画像から関節面画像を抽出して正解データとし、他断面の画像や上腕骨面の画像を黒塗りするなどして上腕骨の特徴を失わせた画像を不正解データとして用いる。さらに学習用のデータとして、元画像を－４５度～４５度まで５度間隔で回転させた画像など角度を異ならせた画像も加えてもよい。
このような機械学習アルゴリズムの学習過程で用いた入力データｘと同じサイズの画像データを、ＡＸ画像の探索範囲から抽出して入力データとして学習済み機械学習アルゴリズムに入力し、上腕骨の関節面画像と判定された画像データを抽出する。これにより、上腕骨関節面の候補の画像が複数得られる。
＜＜Ｓ６３－２：マッチング＞＞
上腕骨関節面の複数の候補から、マッチング処理部２３１は、上腕関節面の位置及び肩甲骨と上腕骨とを結ぶラインの傾きを算出する。
具体的には図９（ａ）、（ｂ）に示すように、マッチング処理によって肩甲骨の両端部（関節唇）の位置を特定し、これを囲む領域を関節面の位置８０１とする。また肩甲骨と上腕骨とを結ぶラインＬ１のテンプレートに対する傾きを求める。
＜＜Ｓ６３－３：対称線検出＞＞
さらに、対称線検出部２３５が、肩甲骨と上腕骨とを結ぶラインＬ１に沿って、ラインＬ１に垂直な対称線Ｌ３を検出し、上腕骨の中心位置を算出する。
対称線の検出は、以下のような手法を用いることができる。まず対象となるアキシャル面画像において、図１０に示すように、画像上に与えた中心線Ｌｃの左右両側に幅Ｗ、高さＨの領域ＲＯＩＡ、ＲＯＩＢを抽出し、ＲＯＩＡとＲＯＩＢの左右反転画像との一致度（或いはＲＯＩＢとＲＯＩＡの左右反転画像との一致度）を算出する。一致度の算出は、例えばＳＳＤ（Ｓｕｍ ｏｆ Ｓｑｕａｒｅｄ Ｄｉｆｆｅｒｅｎｃｅ）やＳＡＤ（Ｓｕｍ ｏｆ Ａｂｏｓｏｌｕｔｅ Ｄｉｆｆｅｒｅｎｃｅ）など公知の手法を用いることができる。ＲＯＩの大きさ（Ｗ、Ｈ）は、検査対象に応じて所定の予めプリセットしておくものとする。
このような一致度の算出を、ＲＯＩＡ、ＲＯＩＢを設定する位置を所定の探索範囲で変化させながら、一致度がもっとも高くなる中心線の位置を対称線とする。探索範囲は、例えば、図９（ｂ）の対称線Ｌ３を抽出する場合、ラインＬ１に沿った所定の範囲（ＰｓからＰｅまでの範囲）とし、各位置でＰｓとＰｅとを結ぶ線に対しオフセット距離Ｏを与えた位置を中心として中心線Ｌｃの角度φを変化させて角度毎の一致度を求める。つまりＲＯＩの中心の位置と中心線の角度の両方を変化させながら一致度を算出し、一致度が最大となる位置及び角度を求める。所定の範囲及びオフセット距離は、対象とする部位に応じて予めプリセットしておくものとする。このような探索によって一致度が最大となる位置及び角度で決まるラインを対称線Ｌ３とする。図９（ｂ）の例では、こうして決定した対称線Ｌ３とラインＬ１との交点を上腕骨の中心位置８０２とする。
＜Ｓ６４：断面位置算出＞
断面位置算出部２５０は、以上のステップＳ６３で得られた上腕骨関節面の位置８０１や上腕骨の中心位置８０２をもとに、２回目のスカウト撮像における推奨撮像位置を算出する。具体的には、上腕骨の中心位置８０２が中心となるように推奨アキシャル面（ＡＸ面）を算出する。また肩甲骨と上腕骨とを結ぶラインＬ１に平行なコロナル面（オブリークＣＯＲ面）と、それに垂直なサジタル面（オブリークＳＡＧ面）をそれぞれ２回目スカウト撮像の推奨撮像位置とする。オブリークＣＯＲ面とオブリークＳＡＧ面については、上腕骨の中心位置８０２を含む面を中心スライスとして、それぞれ複数のスライス面を設定する。なお各断面のＦＯＶやスライス数は、予め設定された撮像条件により設定されている。
＜Ｓ６５：推奨撮像位置提示＞
表示制御部３３０は、撮像断面決定部２００が決定した推奨断面をディスプレイ１１１に表示する。推奨断面の表示方法は特に限定されないが、例えば、図９（ｃ）に示すように、１回目のスカウトＡＸ画像の上に、２回目のスカウト撮像位置として、肩甲骨と上腕骨を結ぶ線Ｌ１に平行な複数のＣＯＲ面の位置、肩甲骨と上腕骨を結ぶ線Ｌ１に垂直な複数のＳＡＧ面の位置、上腕骨の中心位置８０２を中心とするＡＸ面の位置をそれぞれ提示する。ユーザーは、提示された推奨撮像断面で問題がないかを確認し、２回目のスカウト撮像に進むことができる。また表示制御部３３０は、推奨撮像位置の表示とともにユーザー調整を受付けるＧＵＩ等のユーザーインターフェイスを表示してもよく、その場合、ユーザーはユーザーインターフェイスを介して推奨撮像位置の調整を行うことができる。以上のステップＳ６１～Ｓ６５により（図６Ａ）、２回目のスカウト撮像の位置の自動設定が完了する。
＜Ｓ６６：２回目のスカウト撮像＞
撮像制御部３１０は、ユーザーからの２回目のスカウト撮像開始の指示を受けて、スカウト撮像を開始する。撮像部１０は、２回目のスカウト撮像では、上記ステップＳ６４で決定された各撮像断面（ＡＸ面、オブリークＣＯＲ面、オブリークＳＡＧ面）、或いはユーザーによる位置調整後の撮像断面の位置で撮像を行い、スカウト画像を取得する。
＜Ｓ６７、Ｓ６８：スカウト画像の特徴抽出及び撮像断面位置算出＞
撮像断面決定部２００は、２回目のスカウト撮像で取得した画像（ＡＸ画像、オブリークＳＡＧ画像、及び、オブリークＣＯＲ画像）を用いて、本撮像における撮像断面位置の算出を行う。本実施形態では、本撮像の撮像断面位置は、上腕骨の関節面を画像の中心として、ＡＸ面は上腕骨の方向に対し垂直な面となり、ＳＡＧ面及びＣＯＲ面は上腕骨の方向が画像の一つの軸（縦軸）となるように三断面を決定する。
＜＜Ｓ６７－１：ＡＸ画像を用いた処理＞＞
このため、まず２回目のスカウト撮像のＡＸ画像を入力として、ステップＳ６３－２と同様に、マッチング処理により、図１１（ａ）に示すように、上腕骨関節面の位置８１１を特定するとともに、肩甲骨と上腕骨とを結ぶラインＬ１１の傾きαを算出する。さらに、対称線抽出処理により、Ｌ１１に垂直な対称線Ｌ３３を抽出し、Ｌ１１と対称線Ｌ３３との交点である上腕骨の中心８１２を算出する。ラインＬ１１に垂直な面と上腕骨の中心位置８１２の情報は、次のステップＳ６７－２においてＳＡＧ面のＭＰＲ画像を作成する際の位置情報として用いられる。
＜＜Ｓ６７－２：ＳＡＧ画像を用いた処理＞＞
次に複数のスライスデータからなるオブリークＳＡＧ画像から図４に示したようにＭＰＲ画像を作成する。ＳＡＧ面のＭＰＲ画像は、ステップＳ６７－１で得られた情報を用い、上腕骨の中心位置を通り、肩甲骨と上腕骨を結ぶラインに垂直な平面とする。これにより、図１１（ｂ）に示すようなＳＡＧのＭＰＲ画像が作成される。このＭＰＲ画像において、さらに上腕骨の傾き（ＭＰＲ画像上の上腕骨の傾き）βを求める。上腕骨の傾きβを求める処理は、図１０を用いて説明した対称線の抽出処理と同様である。即ち、上腕骨の形状に合わせてサイズが予め定められたＲＯＩ（ＲＯＩＡ、ＲＯＩＢ）を画像上に設定し、所定の検索範囲でＲＯＩを移動しながら且つ各位置でのＲＯＩの角度を変更しながら、それぞれの位置及び角度で、ＲＯＩ内の画像とその左右反転画像との一致度を算出し、一致度が最大となる位置及び角度を求める。この位置と角度のときの左右二つのＲＯＩＡ及びＲＯＩＢの中心線を対称線とし、この対称線の角度を上腕骨の傾きβとする。
＜＜Ｓ６７－３：ＣＯＲ画像を用いた処理＞＞
次に、ステップＳ６７－２で求めた対称線の角度βと位置及びステップＳ６７－１で算出した肩甲骨と上腕骨とを結ぶラインＬ１１の傾きαを用いて、２回目のスカウト撮像で得られたオブリークＣＯＲ画像からＭＰＲ画像を作成する。このＭＰＲ画像もＳＡＧ面のＭＰＲ画像と同様に、図１１（ｃ）に示すように、上腕骨の関節面中心を通り、肩甲骨と上腕骨を結ぶラインＬ１１及び上腕骨の傾きβに平行な平面のＭＰＲ画像を作成する。上腕骨の傾きβの情報は、上述したオブリークＳＡＧ面画像の対称線抽出処理において求められている。ついで作成したＣＯＲ面のＭＰＲ画像から、このＣＯＲ面における上腕骨関節面の中心位置８１３と傾きγを算出する。この処理は、上述した対称線抽出処理と同様に行うことができる。即ち上腕骨について対称線を抽出し、抽出した対称線と直交する方向の傾きを、上腕骨関節面の傾きγとする。あるいは，マッチング処理にて上腕骨関節面の位置および傾きγを算出する方法もある。
＜Ｓ６８：撮像断面位置決定＞
断面位置算出部２５０は、ステップＳ６７－１～ステップＳ６７－３で求めた上腕骨関節面の位置８１３及び上腕骨の傾き（β、γ）をもとにアキシャル面、サジタル面、コロナル面の撮像位置を算出する。具体的には、アキシャル面の撮像位置は、図１２（ａ）に示すように、上腕骨関節面を画像の中心とし、肩甲骨と上腕骨を結ぶラインＬ１１と垂直な方向に上腕骨の傾きβに相当する角度だけ傾き、上腕骨中心を中心としてラインＬ１１の方向に沿って、傾きγに相当する角度だけ傾いたオブリークＡＸ面となる。
同様にサジタル面及びコロナル面の撮像位置は、それぞれ、上記オブリークＡＸ面に直交し、且つ互いに直交するオブリークＳＡＧ面、オブリークＣＯＲ面であり、図１２（ｂ）（ｃ）に示すように、上腕骨関節面を中心とし、上腕骨の軸に沿って傾いたオブリーク面となる。
＜Ｓ６９：推奨撮像位置提示＞
表示制御部３３０は、以上の処理によって決定した本撮像の撮像位置を表示装置に表示させる。撮像位置の表示例は、例えば図１２（ａ）～（ｃ）に示したように、２回目のＡＸ面スカウト画像、ＳＡＧ面のＭＰＲ画像、及びＣＯＲ面のＭＰＲ画像の上に、撮像位置を四角い枠線で囲って示すことができる。ユーザーは、本撮像で取得す画像の位置や範囲を確認することができる。また枠線をユーザーが変更可能なＧＵＩとしてもよく、これによりユーザーは枠線の位置や大きさを変えることで撮像位置の調整を行うことができる。
ユーザーによる変更を受付けた場合は（Ｓ７０）、断面位置算出部２５０がユーザー変更を反映した撮像断面位置を再度算出し、表示制御部３３０により表示装置１１１に表示してもよい。或いはユーザーによる変更を撮像位置調整／実行部３２０が受付けて、撮像位置を決定する撮像パラメータを変更してもよい。撮像制御部３１０は、最終的に決定した撮像断面位置にて撮像を開始する。
以上説明したように、本実施形態によれば、１回目のスカウト撮像を開始した後、自動的に本撮像の撮像断面位置を決定することができる。その際、１回目のスカウト撮像及び２回目のスカウト撮像に三次元撮像を採用することなく撮像断面位置決定に必要かつ十分な情報を収集することができるので、従来技術に比べ短時間で自動位置決め機能を達成することができる。
スカウト撮像から自動撮像断面位置設定までの時間を、従来技術と本実施形態を適用した場合とで比較したタイムチャートを図１３に示す。図示するように、スカウト撮像として三次元撮像を採用した従来技術では、撮像に４０～６０秒、自動位置設定に５～秒、最大で８０秒程度の時間がかかるのに対し、本実施形態によれば、撮像時間及び演算時間の両方を低減することができ、最大でも３４秒程度で自動位置設定を完了することができる。
上記の実施例では，２回目のスカウト撮像でＡＸ画像、ＳＡＧ画像、ＣＯＲ画像を取得する例について説明したが、ＡＸ画像は１回目に撮像しているため、２回目の撮像においては省略することもできる。これによって、より短い計測時間で撮像位置の自動設定が完了できる。
［実施形態２］
本実施形態は、検査部位が膝関節である場合の自動位置決め処理の実施形態である。ここでは、関節面の中心位置を画像の中心とし、ＡＸ面は大腿骨及び脛骨とほぼ垂直な面、ＣＯＲ面及びＳＡＧ面は大腿骨及び脛骨の方向とほぼ平行な面、となるように本撮像の断面位置を決定する。
本実施形態においても装置の構成は、実施形態１と同様であり、以下、実施形態１と異なる点を中心に、図１４を参照して説明する。図１４は、本実施形態の自動位置決め処理を示すフローである。
本実施形態では、１回目のスカウト撮像が一次元計測である場合を説明する。
＜Ｓ８１、Ｓ８２：左右判定＞
撮像制御部３１０の制御のもと、一次元計測のスカウト撮像を行った後（Ｓ８１）、左右判定部２１０は、一次元投影像（例えば図７（ｂ）に示したようなグラフ）を入力として、検査部位が右膝か左膝かの判定を行う（Ｓ８２）。本実施形態では、スカウト撮像で一次元データが得られるので、投影処理は不要である。左右の判定方法は実施形態１と同様に左右の軸に相当するｘ方向の信号値の強度の偏りから判定する。
＜Ｓ８３：一断面のスカウト撮像位置決定＞
断面位置算出部２５０は、左右判定部２１０の判定結果に基づき、左右どちらかに中心を合わせた２回目のスカウト撮像における推奨撮像位置（ＡＸ面の位置）を決定する（Ｓ８３）。２回目のスカウト撮像の撮像位置としては、例えば、一次元プロファイルから算出した信号強度Ｓ（ｘ）が閾値より高くなる位置ｘの平均値ｍ Ａ の位置をＦＯＶの中心とするＡＸ面を設定する。
この際、左右判定部２１０による判定結果をユーザーの確認のために表示装置１１１に表示してもよい。提示方法は、左右のいずれであるかを示してもよいし、１回目のスカウト撮像で得られた一次元プロファイルと判定結果を表示してもよい。ユーザーは、この表示により検査プロトコルの検査対象に対し正しく受信ＲＦコイルが装着されているかなどを確認することができる。
＜Ｓ８４～Ｓ８６：特徴抽出と他断面のスカウト撮像位置決定＞
撮像制御部３１０は設定された位置で、まずＡＸ画像（複数のスライスのＡＸ画像）、それぞれ直交するＣＯＲ画像（複数のスライスのＣＯＲ画像）、ＳＡＧ画像（複数のスライスのＳＡＧ画像）を取得する（Ｓ８４）。本撮像用のＣＯＲ面の設定のため、ステップＳ８２のスカウト撮像で取得したＡＸ画像を入力として、マッチング処理により、図１５（ａ）に示すように、内側顆と外側顆とを結ぶラインＬ４の傾き（θ４）を算出する（マッチング処理部２３１の処理）。
次にＣＯＲ面のスカウト撮像で取得した複数のスライスデータを用いて特徴抽出を行う（Ｓ８５）。まず、ＣＯＲ面の複数のスライスデータに対し、画像切り出し部２３７によるＭＰＲ処理を行い、図１５（ｂ）に示すように、内側顆と外側顆とを結ぶラインＬ４を通るＣＯＲ面のＭＰＲ画像を作成する。このＣＯＲ面のＭＰＲ画像において、マッチング処理部２３１により関節面の中心位置を求める。この中心位置を通り、ＣＯＲ面に直交する面をＳＡＧ面とする。
さらに、ＳＡＧ面の複数のスライスデータに対し関節面中心位置でＭＰＲ処理を行い、図１５（ｃ）に示すようなＳＡＧ面のＭＰＲ画像を作成し、本撮像の撮像断面位置決定に必要な特徴抽出を行う（Ｓ８６）。まず、作成したＳＡＧ面画像のＭＰＲ画像に対し、関節面上部及び下部のそれぞれについて、対称線検出部２３５が対称線抽出処理を行い、対称線Ｌ５、Ｌ６を求め、その傾きを算出する。関節面より上部の対称線の傾きは大腿骨の傾き（θ５）とみなすことができ、関節面より下部の対称線の傾きは脛骨の傾き（θ６）とみなすことができる。
＜Ｓ８７：本撮像の断面位置算出、Ｓ８８：表示＞
断面位置算出部２５０は、ステップＳ８４で得た内側顆と外側顆とを結ぶラインＬ４の傾き（θ４）、及びステップＳ８６で得た大腿骨及び脛骨の傾き（θ５、θ６）を用いて本撮像の撮像位置を算出し、提示する（Ｓ８８）。例えば、本撮像の各断面は、関節面の中心位置が画像の中心となるように設定され、ＡＸ面については、内側顆と外側顆とを結ぶラインＬ４と画像の横軸とが平行になり、且つ大腿骨の傾きと脛骨の傾きとを平均化した傾きに対し直交する面となるように設定される。ＳＡＧ面については、画像の縦軸が大腿骨の傾きと脛骨の傾きとの平均傾きと平行になり、且つ内側顆と外側顆とを結ぶラインＬ４の方向と直交する面となるように設定される。ＣＯＲ面については、大腿骨の傾きと脛骨の傾きとの平均と平行であって、且つ内側顆と外側顆とを結ぶラインＬ４と平行な面に設定される。
このような三断面の表示例を図１６（ａ）～（ｃ）に示す。図示する例では、三断面のそれぞれについて、スカウト撮像で得たＡＸ画像、ＭＰＲ処理したＣＯＲ画像、ＭＰＲ処理したＳＡＧ画像を表示し、それぞれに重畳して、各断面の位置（太線の四角）が表示される。本実施形態においてもユーザー変更を受け付けることができる（Ｓ８９）。
本実施形態によれば、自動位置決めの際に前提となる左右判定を一次元計測のみで行うので、２回目のスカウト撮像を含むスカウト撮像時間を大幅に低減することができる。またも左右判定結果に基づいて２回目のスカウト撮像を行うことにより、２回目のスカウト撮像で得た画像の特徴抽出とそれに基づく撮像位置算出の精度を向上することができる。
本実施形態では、一次元撮像にて左右を判定するため、左右の位置を確認するためのスカウト撮像が不要となる。一次元の撮像時間は１秒程度、その後の三断面のスカウト撮像時間は１０～１５秒、処理時間が２秒であるので、最大でも１８秒以内に自動位置決め設定を行うことができ、三次元撮像を行う従来技術に比べ時間の短縮効果が得られる。
以上、本発明の自動位置決め技術を、肩関節及び膝関節に適用した実施形態を説明したが、本発明はこれらの部位に限定されることなく、左右判定を必要とする部位、主として関節のＭＲＩ検査について同様に適用することができ、同様の効果を得ることができる。
１０：撮像部、２０：演算部、３０：制御部、１００：ＭＲＩ装置、１０１：静磁場磁石、１０３：送信用ＲＦコイル、１０５：傾斜磁場コイル、１０６：受信用ＲＦコイル、１１０：計算機、１１１：表示装置、１１２：入力装置、１１７：シーケンサ、１１８：記憶装置、２００：撮像断面決定部、２１０：左右判定部、２３０：特徴抽出部、２３１：マッチング処理部、２３３：物体抽出部、２３５：対称線検出部、２３７：画像切り出し部、２５０：断面位置算出部、３１０：撮像制御部、３３０：表示制御部。

Claims (17)

  1. 核磁気共鳴を用いて被検体の断面の画像を取得する撮像部と、
    撮像位置を決定するためのスカウト撮像、及び、診断用の画像を取得する本撮像を前記撮像部が行うよう制御する撮像制御部と、
    前記スカウト撮像で取得したデータを用いて、前記本撮像の撮像断面位置を決定する撮像断面決定部と、を備え、
    前記撮像制御部は、前記スカウト撮像として、１次元又は２次元の計測データを取得する第１計測、及び、２次元の計測データを取得する第２計測を実行する制御を行い、
    前記撮像断面決定部は、前記第１計測で取得した計測データを用いて、前記被検体の左右を判定する左右判定部と、
    前記左右判定部における判定結果と前記第１計測で取得した計測データとを用いて前記第２計測における断面位置を算出するとともに、前記第２計測で取得した計測データを用いて前記本撮像における断面位置を算出する断面位置算出部と、を備えたことを特徴とする磁気共鳴イメージング装置。
  2. 請求項１に記載の磁気共鳴イメージング装置であって、
    前記左右判定部は、前記第１計測で取得した計測データの、前記被検体の左右方向に沿った一次元方向の信号値をもとに、左右の判定を行うことを特徴とする磁気共鳴イメージング装置。
  3. 請求項１に記載の磁気共鳴イメージング装置であって、
    前記撮像制御部は、前記第１計測において、２次元の計測データを取得する制御を行い、
    前記左右判定部は、前記第１計測で取得した前記２次元の計測データを前記被検体の左右方向に沿った一次元方向に投影した一次元データを作成し、当該一次元データを用いて左右の判定を行うことを特徴とする磁気共鳴イメージング装置。
  4. 請求項３に記載の磁気共鳴イメージング装置であって、
    前記断面位置算出部は、前記第１計測で取得した一断面の計測データを用いて、前記一断面を含む複数断面について、前記第２計測における断面位置を算出することを特徴とする磁気共鳴イメージング装置。
  5. 請求項１に記載の磁気共鳴イメージング装置であって、
    前記撮像制御部は、前記第１計測において、前記被検体の左右方向に沿った一次元方向をエコー信号の読み出し方向とするパルスシーケンスを実行して、１次元の計測データを取得する制御を行うことを特徴とする磁気共鳴イメージング装置。
  6. 請求項５に記載の磁気共鳴イメージング装置であって、
    前記断面位置算出部は、前記左右判定部が前記第１計測で取得した一次元の計測データを用いて左右判定した結果に基づき、一断面について第２計測における断面位置を算出し、前記第２計測で取得した前記一断面の計測データを用いて、前記一断面以外の断面について前記第２計測における断面位置を算出することを特徴とする磁気共鳴イメージング装置。
  7. 請求項１に記載の磁気共鳴イメージング装置であって、
    前記撮像断面決定部は、前記第１計測または前記第２計測で取得した計測データから検査対象組織の特徴を抽出する特徴抽出部をさらに備えることを特徴とする磁気共鳴イメージング装置。
  8. 請求項７に記載の磁気共鳴イメージング装置であって、
    前記撮像断面決定部は、前記第２計測で取得した複数スライス分の計測データから、検査対象組織を含む１枚のスライスデータを切り出し、ＭＰＲ画像を作成する画像切り出し部をさらに備え、
    前記特徴抽出部は、前記画像切り出し部が作成したＭＰＲ画像から検査対象組織の特徴を抽出することを特徴とする磁気共鳴イメージング装置。
  9. 請求項７に記載の磁気共鳴イメージング装置であって、
    前記特徴抽出部は、前記検査対象組織の特徴として、前記検査対象組織の位置及び傾きの少なくとも一方をマッチング処理により算出するマッチング処理部を備えることを特徴とする磁気共鳴イメージング装置。
  10. 請求項７に記載の磁気共鳴イメージング装置であって、
    前記特徴抽出部は、前記検査対象組織の特徴として、機械学習アルゴリズムを用いて画像データから前記検査対象組織の位置を抽出する物体抽出部を備えることを特徴とする磁気共鳴イメージング装置。
  11. 請求項７に記載の磁気共鳴イメージング装置であって、
    前記特徴抽出部は、前記検査対象組織について左右対称となる中心線とその傾きを算出する対称線検出部を備えることを特徴とする磁気共鳴イメージング装置。
  12. 請求項１に記載の磁気共鳴イメージング装置であって、
    前記左右判定部が判定した結果及び／又は前記断面位置算出部が算出した断面位置を表示装置に表示させる表示制御部をさらに備えることを特徴とする磁気共鳴イメージング装置。
  13. 請求項１２に記載の磁気共鳴イメージング装置であって、
    前記表示制御部は、断面位置についてのユーザー変更を受付けるユーザーインターフェイスを備えることを特徴とする磁気共鳴イメージング装置。
  14. 磁気共鳴イメージング装置を用いて撮像する際の撮像断面位置を自動設定する方法であって、
    診断用画像を取得する本撮像に先立って、スカウト撮像の第１計測を実行し、１次元又は２次元の計測データを取得し、
    前記第１計測で取得した計測データを用いて、被検体の左右を判定し、
    前記左右判定における判定結果と前記第１計測で取得した計測データと用いて、スカウト撮像の第２計測における断面位置を算出し、当該断面位置にて前記第２計測を実行し、２次元の計測データを取得し、
    前記第２計測で取得した２次元の計測データを用いて前記本撮像における断面位置を算出することを特徴とする自動撮像位置設定方法。
  15. 請求項１４に記載の自動撮像位置設定方法であって、
    前記第１計測で２次元の計測データを取得し、当該２次元の計測データを、前記被検体の左右方向に沿った一次元方向に投影した１次元データを作成し、
    当該１次元データの信号強度をもとに前記被検体の左右を判定することを特徴とする自動撮像位置設定方法。
  16. 請求項１５に記載の自動撮像位置設定方法であって、
    前記第１計測で取得した一断面の２次元の計測データを用いて、前記第２計測における前記一断面を含む複数断面の断面位置を算出することを特徴とする自動撮像位置設定方法。
  17. 請求項１４に記載の自動撮像位置設定方法であって、
    前記第１計測で１次元の計測データを取得し、当該１次元の計測データの信号強度をもとに前記被検体の左右を判定することを特徴とする自動撮像位置設定方法。
JP2019185131A 2019-10-08 2019-10-08 磁気共鳴イメージング装置、及び、自動撮像位置設定方法 Active JP7320421B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2019185131A JP7320421B2 (ja) 2019-10-08 2019-10-08 磁気共鳴イメージング装置、及び、自動撮像位置設定方法
CN202010472767.XA CN112637478B (zh) 2019-10-08 2020-05-29 磁共振成像装置和自动拍摄位置设定方法
US16/889,987 US11255939B2 (en) 2019-10-08 2020-06-02 Magnetic resonance imaging apparatus, and automatic imaging position setting method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019185131A JP7320421B2 (ja) 2019-10-08 2019-10-08 磁気共鳴イメージング装置、及び、自動撮像位置設定方法

Publications (2)

Publication Number Publication Date
JP2021058447A JP2021058447A (ja) 2021-04-15
JP7320421B2 true JP7320421B2 (ja) 2023-08-03

Family

ID=75275011

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019185131A Active JP7320421B2 (ja) 2019-10-08 2019-10-08 磁気共鳴イメージング装置、及び、自動撮像位置設定方法

Country Status (3)

Country Link
US (1) US11255939B2 (ja)
JP (1) JP7320421B2 (ja)
CN (1) CN112637478B (ja)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090138020A1 (en) 2007-11-27 2009-05-28 Otismed Corporation Generating mri images usable for the creation of 3d bone models employed to make customized arthroplasty jigs
WO2010150783A1 (ja) 2009-06-25 2010-12-29 株式会社 日立メディコ 医用画像撮影装置
JP2014000393A (ja) 2012-05-22 2014-01-09 Toshiba Corp 医用画像診断装置及び画像表示装置
JP2014008245A (ja) 2012-06-29 2014-01-20 Ge Medical Systems Global Technology Co Llc 検出装置、医用装置、およびプログラム
JP2015142723A (ja) 2013-12-27 2015-08-06 株式会社東芝 磁気共鳴イメージング装置
JP2016119976A (ja) 2014-12-24 2016-07-07 株式会社東芝 医用画像診断装置
US20180329004A1 (en) 2017-05-12 2018-11-15 Canon Medical Systems Corporation Magnetic resonance imaging apparatus

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5944645B2 (ja) * 2010-11-02 2016-07-05 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
JP5971682B2 (ja) * 2011-03-02 2016-08-17 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
JP6218569B2 (ja) 2012-11-22 2017-10-25 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
JP6440980B2 (ja) * 2013-07-10 2018-12-19 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
JP6517031B2 (ja) * 2015-02-05 2019-05-22 キヤノンメディカルシステムズ株式会社 医用画像処理装置および磁気共鳴イメージング装置
JP6378426B2 (ja) * 2015-04-22 2018-08-22 株式会社日立製作所 磁気共鳴イメージング装置および磁気共鳴イメージング方法
EP3373024A1 (en) * 2017-03-08 2018-09-12 Esaote S.p.A. Method and device for automatic scan planning in mri
KR101971736B1 (ko) * 2017-08-18 2019-04-24 엘지전자 주식회사 이동단말기 및 그 제어방법

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090138020A1 (en) 2007-11-27 2009-05-28 Otismed Corporation Generating mri images usable for the creation of 3d bone models employed to make customized arthroplasty jigs
WO2010150783A1 (ja) 2009-06-25 2010-12-29 株式会社 日立メディコ 医用画像撮影装置
JP2014000393A (ja) 2012-05-22 2014-01-09 Toshiba Corp 医用画像診断装置及び画像表示装置
US20150084959A1 (en) 2012-05-22 2015-03-26 Kabushiki Kaisha Toshiba Medical image diagnosis apparatus and image display apparatus
JP2014008245A (ja) 2012-06-29 2014-01-20 Ge Medical Systems Global Technology Co Llc 検出装置、医用装置、およびプログラム
JP2015142723A (ja) 2013-12-27 2015-08-06 株式会社東芝 磁気共鳴イメージング装置
JP2016119976A (ja) 2014-12-24 2016-07-07 株式会社東芝 医用画像診断装置
US20180329004A1 (en) 2017-05-12 2018-11-15 Canon Medical Systems Corporation Magnetic resonance imaging apparatus
JP2018191697A (ja) 2017-05-12 2018-12-06 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置

Also Published As

Publication number Publication date
CN112637478A (zh) 2021-04-09
US20210103021A1 (en) 2021-04-08
US11255939B2 (en) 2022-02-22
CN112637478B (zh) 2022-04-19
JP2021058447A (ja) 2021-04-15

Similar Documents

Publication Publication Date Title
JP5230809B2 (ja) 医用画像撮影装置
JP5736459B2 (ja) 医用画像撮影装置
Harris et al. Correlations between the alpha angle and femoral head asphericity: implications and recommendations for the diagnosis of cam femoroacetabular impingement
JP5667066B2 (ja) Ｍｒスキャンの自動シーケンシャルプランニング
US8923584B2 (en) Method and system of automatic determination of geometric elements characterizing a bone deformation from 3D image
Konan et al. Is the frog lateral plain radiograph a reliable predictor of the alpha angle in femoroacetabular impingement?
Ahedi et al. Hip shape as a predictor of osteoarthritis progression in a prospective population cohort
JP2004524942A (ja) 断層撮像パラメータの自動指示
US9763636B2 (en) Method and system for spine position detection
US7561910B2 (en) Method for determination of the supported position of a patient in a magnetic resonance apparatus
Ghafurian et al. Computerized 3D morphological analysis of glenoid orientation
CN104717920B (zh) 磁共振成像装置以及磁共振成像方法
CN112946547A (zh) 用于从校准图像生成诊断扫描参数的系统和方法
Van Ijsseldijk et al. Three dimensional measurement of minimum joint space width in the knee from stereo radiographs using statistical shape models
JP7320421B2 (ja) 磁気共鳴イメージング装置、及び、自動撮像位置設定方法
JP2020005889A (ja) 自動位置決め機能付き磁気共鳴イメージング装置、磁気共鳴イメージング装置の制御方法、及び、撮像断面自動設定用プログラム
Neubert et al. Automatic segmentation of the glenohumeral cartilages from magnetic resonance images
US10820837B2 (en) System and method for angle measurements using magnetic resonance images
Kang et al. Preliminary evaluation of an MRI-based technique for displaying and quantifying bony deformities in cam-type femoroacetabular impingement
Blumenfeld et al. Automatic prospective registration of high-resolution trabecular bone images of the tibia
Horenstein et al. Measuring 3D growth plate shape: Methodology and application to cam morphology
US20240005480A1 (en) Methods and systems for automated saturation band placement
US20220133217A1 (en) Method and system for visualizing a property of cartilage
Mahendrakar et al. A Survey on Morphological Assessment of Knee Articular Cartilage from MR Images
Reisman et al. Robust local intervertebral disc alignment for spinal MRI

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20211013

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20221003

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230426

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230530

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230612

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: 20230627

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230724

R150 Certificate of patent or registration of utility model

Ref document number: 7320421

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150