JP5230809B2 - 医用画像撮影装置 - Google Patents

医用画像撮影装置 Download PDF

Info

Publication number
JP5230809B2
JP5230809B2 JP2011519903A JP2011519903A JP5230809B2 JP 5230809 B2 JP5230809 B2 JP 5230809B2 JP 2011519903 A JP2011519903 A JP 2011519903A JP 2011519903 A JP2011519903 A JP 2011519903A JP 5230809 B2 JP5230809 B2 JP 5230809B2
Authority
JP
Japan
Prior art keywords
image
recommended
section
anatomical feature
photographing
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
JP2011519903A
Other languages
English (en)
Other versions
JPWO2010150783A1 (ja
Inventor
俊 横沢
陽 谷口
良孝 尾藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
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 Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP2011519903A priority Critical patent/JP5230809B2/ja
Publication of JPWO2010150783A1 publication Critical patent/JPWO2010150783A1/ja
Application granted granted Critical
Publication of JP5230809B2 publication Critical patent/JP5230809B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/0037Performing a preliminary scan, e.g. a prescan for identifying a region of interest
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0042Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/483NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/543Control 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • A61B2576/02Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4528Joints
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4533Ligaments
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/546Interface between the MR system and the user, e.g. for controlling the operation of the MR system or for the design of pulse sequences
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Neurology (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明は、磁気共鳴イメージング装置等の医用画像撮影装置を用いた検査における撮影断面の位置決め技術および診断画像切り出し技術に関する。
3次元空間上で所望の位置および傾きの断面を撮影し、診断に用いる医用画像診断装置がある。中でも、磁気共鳴イメージング(Magnetic Resonance Imaging)装置は、主にプロトンの核磁気共鳴現象を利用した医用画像診断装置であり、静磁場に置かれた被検体に高周波磁場を印加して核磁化を励起し、計測した磁気共鳴信号を画像化する。MRI装置は、撮影部位に制限が無く、非侵襲に任意の断面の撮影が可能である。
MRI装置では、一般に、撮影する断面(撮影断面)を特定するスライス傾斜磁場を印加すると同時に、撮影断面内の磁化を励起させる励起パルスを与え、励起された磁化が発生する核磁気共鳴(エコー)信号を得る。磁化に位置情報を与えるため、励起からエコー信号を得るまでの間に、撮影断面内で互いに垂直な方向の位相エンコード傾斜磁場とリードアウト傾斜磁場とを印加する。
MRI装置のような任意の断面の撮影が可能な医用画像診断装置では、撮影する断面を予め設定しなければならない。このため、検査では、撮影断面を設定するスカウト画像を取得するスカウト撮影を、診断用画像を取得する本撮影の前に行う。
撮影断面は、撮影部位や対象とする疾患によって概ね決められているため、スカウト画像上で操作者がユーザインタフェースを介してマニュアル操作で行っている。ところが、検査は、異なる撮影断面の撮影を組み合わせて行ったり、定期的に同じ撮影断面での撮影を繰り返したりすることがある。異なる撮影断面の撮影を行う場合は、その度に操作者が設定する必要がある。また、同じ撮影断面の撮影を繰り返す場合、マニュアル操作では、毎回同じ撮影断面を設定することは難しい。いずれにしても、撮影断面の設定は撮影のたびに同じ作業を繰り返す必要があり、単純で煩雑な作業の一つに挙げられている。
このような撮影断面決定の作業の操作性を向上するため、例えば、視覚的に位置決め操作しやすい工夫を行うもの(例えば、特許文献1参照)、画像認識により自動的に位置決めを行うもの(例えば、非特許文献1および2参照)などが提案されている。自動化による効果としては、操作性の向上だけでなく,フォローアップ検査時の撮影断面の再現性向上が期待されている。
特開平05−269113号公報
『Automated planning of MRI neuro scans』(Young S, Bystrov D, Netsch T, Bergmans R, van Muiswinkel A, Visser F, Springorum R, Gieseke J著、Medical Imaging: Image Processing, Reinhardt JM, ed. Proc. of SPIE Vol. 6144, 2006 『Automated Scan Prescription for Brain MRI』(Itti L, Chang L, Ernst T著、Magnetic Resonance in Medicine 45;486−494, 2001
特許文献1に記載の手法は、操作者のマニュアル操作による部分がまだ残っているため、依然として再現性の低さと煩雑さは残る。また、非特許文献1および2に記載の手法では、スカウト画像として3次元のボリュームデータを計測する必要があるためスカウト撮影にかかる時間が従来よりも長くなる。さらに、スカウト撮影後に撮影断面を算出して決定するため、さらに自動化のための処理時間が必要となる。このため、全体としての検査時間が延長するだけでなく、自動化のための処理による待ち時間が操作者にストレスを与え、操作者の要求に十分に応えていない。
一方、近年、パラレルイメージングなどの高速撮影技術の向上に伴い、高速に3次元画像を取得する3次元撮影が行われている。3次元撮影では撮影範囲として3次元空間を設定すればよいため、断面を設定する2次元撮影に比べて容易である。しかし、撮影終了後に得られた3次元ボリュームデータから診断したい2次元画像を切り出す必要がある。この処理はMPR(Multi Planar Reconstruction)と呼ばれ、2次元画像を簡単かつ高速に抽出可能な画像処理システムが求められている。
本発明は、上記の問題を鑑みたもので、上記MRI装置のような3次元空間の任意の面を撮影断面として設定可能な医用画像撮影装置において、従来の検査の流れを変えることなく、また、検査時間を延長すること無く、撮影断面を自動的に算出するとともにMPRにおいて自動的に切り出し断面を算出可能な技術を提供することを目的とする。
本発明は、手動で撮影断面を設定する場合と同様の二次元スカウト撮影を行い、得られたスカウト画像を処理し、推奨撮影断面を算出する。処理のアルゴリズムおよび処理に用いる各種画像処理手法は、予め撮影対象部位および検査の種類毎に保持される。
具体的には、3次元空間の任意の面を撮影可能な医用画像撮影装置であって、互いに交差する2つのスライス面の内の一方のスライス面である第一のスライス面に平行な1枚以上の二次元画像からなる第一の画像群と、他方のスライス面である第二のスライス面に平行な1枚以上の二次元画像からなる第二の画像群とを取得する画像取得手段と、二次元画像から、予め定められた解剖学的特徴構造の情報を抽出する解剖学的特徴構造抽出手段と、撮影断面として推奨する推奨撮影断面を算出する推奨撮影断面算出手段と、撮影対象部位に応じて、前記推奨撮影断面を算出するために必要な推奨撮影断面算出情報を保持する推奨撮影断面算出情報記憶手段と、を備え、前記解剖学的特徴構造抽出手段は、前記第一の画像群から第一の解剖学的特徴構造の情報を抽出するとともに、前記第二の画像群から第二の解剖学的特徴構造の情報を抽出し、前記推奨撮影断面算出手段は、前記推奨撮影断面算出情報に従い、前記第一の解剖学的特徴構造の情報および前記第二の解剖学的特徴構造の情報を用いて前記推奨撮影断面を算出することを特徴とする医用画像撮影装置を提供する。
本発明によれば、上記MRI装置のような3次元空間の任意の面を撮影断面として設定可能な医用画像撮影装置において、従来の検査の流れを変えることなく、また、検査時間を延長すること無く撮影断面を自動的に算出するとともにMPRにおいて自動的に切り出し断面を算出できる。
第一の実施形態のMRI装置の典型的な構成を示すブロック図である。 第一の実施形態の計算機の機能ブロック図である。 (a)〜(c)は第一の実施形態のスカウト画像の一例を示す図である。 第一の実施形態の推奨撮影断面算出処理のフローである。 (a)は、第一の実施形態のコロナル正中線算出処理を説明するための図であり、(b)は、第一の実施形態のアキシャル画像決定処理を説明するための図である。 (a)〜(c)は、第一の実施形態のアキシャル正中線算出処理を説明するための図である。 (a)〜(c)は、第一の実施形態の正中面決定処理を説明するための図である。 (a)〜(c)は、第一の実施形態のテンプレートモデルの作成手順を説明するための図である。 (a)、(b)は、第一の実施形態のテンプレートモデルのフィッティングを説明するための図である。 (a)、(b)は、第一の実施形態のテンプレートモデルから推奨撮影断面算出手法を説明するための図である。 第一の実施形態の計測準備処理を説明するための図である。 (a)、(b)は、第一の実施形態の計測準備処理を説明するための図である。 (a)、(b)は、第一の実施形態の計測準備処理を説明するための図である。 第二の実施形態の計算機の機能ブロック図である。 第三の実施形態の計算機の機能ブロック図である。
<<第一の実施形態>>
以下、本発明を適用する第一の実施形態について説明する。本発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付し、その繰り返しの説明は省略する。
まず、本実施形態の磁気共鳴イメージング(MRI)装置について説明する。本実施形態のMRI装置100は、上述のように、静磁場に置かれた被検体103に高周波磁場を印加して、被検体103内の核磁化を励起し、発生する核磁気共鳴信号(エコー信号)を計測する。このとき、傾斜磁場を印加して計測する磁気共鳴信号に位置情報を与え、画像化(撮影)する。図1は、これを実現する、本実施形態のMRI装置100の典型的な構成を示すブロック図である。本実施形態のMRI装置100は、静磁場を発生するマグネット101と、傾斜磁場を発生する傾斜磁場コイル102と、被検体(生体)103に高周波磁場パルス(以下、RFパルス)を照射するRFコイル107と、被検体103から発生するエコー信号を検出するRFプローブ108と、マグネット101の発生する静磁場空間内で被検体(例えば、生体)103を載置する寝台(テーブル)115と、を備える。
さらに、本実施形態のMRI装置100は、傾斜磁場コイル102を駆動する傾斜磁場電源105と、RFコイル107を駆動する高周波磁場発生器106と、RFプローブ108で検出したエコー信号を受信する受信器109と、傾斜磁場電源105と高周波磁場発生器106とに命令を送り、それぞれ傾斜磁場および高周波磁場を発生させるとともに、検波の基準とする核磁気共鳴周波数を受信器109にセットするシーケンサ104と、検波された信号に対して信号処理を施す計算機110と、計算機110での処理結果を表示する表示装置111と、同処理結果を保持する記憶装置112と、操作者からの指示を受け付ける入力装置116と、を備える。
以上の構成を有するMRI装置100では、シーケンサ104の制御により、RFパルスがRFコイル107を通じて被検体103に印加されるとともに、スライス選択や位相エンコードなどの位置情報をエコー信号に与えるための傾斜磁場パルスが傾斜磁場コイル102によって印加される。また、被検体103から発生した信号はRFプローブ108によって受波され、検波された信号は計算機110に送られ、ここで画像再構成などの信号処理が行われる。なお、記憶装置112には、信号処理の結果だけでなく、必要に応じて、検波された信号自体、撮影条件等を記憶させてもよい。
また、MRI装置100は、静磁場均一度を調節する必要があるときには、シムコイル113と、シムコイル113を駆動するシム電源114をさらに備えてもよい。シムコイル113は、複数のチャネルからなり、シム電源114から供給される電流によりにより静磁場不均一を補正する付加的な磁場を発生する。静磁場均一度調整時にシムコイル113を構成する各チャネルに流す電流は、シーケンサ104により制御される。
なお、シーケンサ104は、上述のようにMRI装置100を構成する各部の動作を制御し、計測を実行させるもので、予めプログラムされたタイミング、強度で各部が動作するように制御を行う。上記プログラムのうち、特に、高周波磁場、傾斜磁場、信号受信のタイミングや強度を記述したものはパルスシーケンスと呼ばれる。計測は、パルスシーケンスとこれを制御するために必要な撮影パラメータとに従って行われる。パルスシーケンスは、予め作成され、記憶装置112に保持され、撮影パラメータは、操作者からユーザインタフェースを介して入力される。
また、計算機110は、受信した信号を処理する信号処理だけでなく、MRI装置100全体の動作の制御等を行う。さらに、本実施形態の計算機110は、記憶装置112とともに情報処理装置を構成し、撮影断面設定用の画像(スカウト画像)から、推奨する撮影断面(推奨撮影断面)を算出する。
これを実現する、本実施形態の計算機110と記憶装置112とにより構成される情報処理装置の機能ブロック図を図2に示す。本図に示すように、本実施形態の情報処理装置は、制御部200と記憶部300とを備える。制御部200は、ユーザインタフェース(UI)制御部210と、信号処理部220と、計測制御部230と、撮影断面決定部240と、を備える。また、記憶部300は、制御部200が実現する各種処理に必要な情報を記憶するもので、撮影情報記憶部310と、最適スカウト撮影情報記憶部320と、推奨撮影断面算出情報記憶部330とを備える。
撮影情報記憶部310には、MRI装置100により検査を実行するために必要な各種の情報、例えば、プロトコルが登録される。一般に、MRI装置100による検査では、スカウト画像を取得するスカウト撮影、静磁場不均一の調整やコイル感度補正のための準備撮影、決定された撮影断面の診断用の画像を取得する本撮影が順次実行される。各撮影は1以上の計測で構成され、各計測は、パルスシーケンスおよび撮影パラメータに従って行われる。プロトコルは、各検査における、撮影順、各撮影内での計測順、計測種別を規定したもので、検査の実行に先立ち、操作者により作成され、記憶される。プロトコルは、例えば、頭部、腰椎、膝、肩といった検査部位や対象とする疾患に応じて作成され、各撮影で実行するパルスシーケンスとその撮影パラメータを含む。
また、撮影情報記憶部310には、プロトコル生成に用いるパルスシーケンスや、操作者から入力される撮影パラメータ等が登録される。登録されるパルスシーケンスには、例えば、FSE(Fast Spin Echo)、GrE(GradientEcho)、EPI(Echo Planar Imaging)、などがあり、撮影パラメータには、TR(繰り返し時間)、TE(エコー時間)、FOV(撮影視野)、スライス厚、スライス枚数、複数の断面を撮影する場合は、その撮影順序などがある。
なお、プロトコルは、検査毎に操作者がユーザインタフェースを介して作成し、撮影情報記憶部310に登録するよう構成してもよいし、予め、検査部位および/または疾患毎に作成され、撮影情報記憶部310に保存されていてもよい。この場合は、操作者が、検査毎に保存されているプロトコルの中から抽出し、決定する。なお、検査部位毎に最適なプロトコルを対応付けて撮影情報記憶部310に保存し、操作者が撮影パラメータの設定時に検査部位を指定した際、当該検査部位に対応づけて記憶されているプロトコルが、最適なプロトコルの初期値として抽出されるよう構成してもよい。
最適スカウト撮影情報記憶部320には、検査部位毎の、また、検査部位毎に複数種の撮影断面の異なる検査種別がある場合は、検査部位および検査種別毎の、最適なスカウト撮影の手順を規定した最適スカウト撮影情報が登録される。具体的には、スカウト撮影で使用するパルスシーケンス、撮影パラメータ、撮影する断面、複数の断面を撮影する場合は、その撮影順が、登録される。なお、最適スカウト撮影は、操作者によってプロトコル作成時に選択および設定ができる。あるいは、上記の最適なプロトコルの初期値に含まれるよう構成してもよい。
推奨撮影断面算出情報記憶部330には、スカウト画像から推奨撮影断面を算出するために用いる情報が登録される。例えば、検査部位に応じて、スカウト画像から推奨撮影断面を生成するために実行するアルゴリズムが登録される。また、アルゴリズム内で実行される画像処理の種類、推奨撮影断面算出に用いる解剖学的な特徴構造等、検査部位に応じた解剖学的な特徴構造と推奨撮影断面との位置関係等が登録される。なお、解剖学的な特徴構造と推奨撮影断面との位置関係は、操作者が設定できる構成でもよい。
なお、複数の推奨撮影断面の設定が可能な部位の場合、目印の名称(OMライン、椎間板ライン等)、検査種別(ルーチン検査、てんかん検査等)ごとに解剖学的な特徴構造と推奨撮影断面との位置関係を撮影断面リストとして保持し、プロトコル作成時に操作者が、本撮影のパルスシーケンスごとに選択可能なように構成してもよい。これにより、操作者の選択に応じた撮影断面の設定が可能となる。なお、撮影断面リストからの選択時に、撮影断面と、解剖学的目印位置との関係を操作者が目視で確認できるよう、表示装置111に表示させるグラフィカルインタフェース(GI)を備えるよう構成してもよい。さらに、このGI上で位置の調整入力を受け付けるよう構成してもよい。
解剖学的特徴構造は、例えば、検査部位が頭部であれば、正中線、頭部輪郭、脳輪郭、脳梁、橋、脳幹、下垂体、斜台、腰椎であれば脊髄神経の傾き、椎間板の位置、膝であれば、大腿骨内側顆、大腿骨外側顆、大腿骨、頸骨の位置、大腿骨内側顆と大腿骨外側顆とを結んだ線、大腿骨と頸骨と間の関節面、肩であれば、棘上筋、骨頭、肩甲骨、肩峰、鎖骨の位置、棘上筋に平行な線、上腕骨に沿ったライン、骨頭と肩甲骨間の関節面の接線、骨頭と肩甲骨間を結ぶラインなどである。これらの解剖学的特徴構造は、検査部位に応じて予め定められ、登録される。
UI制御部210は、入力装置116および表示装置111などのユーザインタフェースを制御し、操作者に処理結果を提示する、操作者からの入力を受け付ける、といったユーザインタフェース処理を行う。
信号処理部220は、MRI装置100で取得したエコー信号を処理し、画像を再構成する。また、エコー信号を処理し、中心周波数やRF照射強度等の撮影に必要な制御値を算出し、装置にシーケンサ104に送信する。
計測制御部230は、操作者によって作成され、撮影情報記憶部310に登録されたプロトコルに従って、実行するパルスシーケンス、用いる撮影パラメータを読み出し、シーケンサ104に命令を与え、計測を実行する。
撮影断面決定部240は、スカウト画像に基づき、本撮影での撮影断面を決定する。また、本実施形態の撮影断面決定部240は、スカウト画像から推奨撮影断面を算出する推奨撮影断面算出部241を備える。また、推奨撮影断面算出部241は、二次元画像上で予め定められた解剖学的特徴構造を抽出する解剖学的特徴構造抽出部242と、2つの二次元画像上それぞれの解剖学的特徴構造を用い、新たなスライス面を決定するスライス面決定部243と、を備える。
推奨撮影断面算出部241は、最適スカウト撮影情報に従って実行されたスカウト撮影で得たスカウト画像から、推奨撮影断面算出情報記憶部330に登録されるアルゴリズムに従って、推奨撮影断面を算出する。
解剖学的特徴構造抽出部242は、推奨撮影断面算出情報記憶部330に記憶されたアルゴリズムに従って、スカウト画像などの二次元画像から推奨撮影断面算出情報記憶部330に記憶された画像処理により、推奨撮影断面算出情報記憶部330に記憶された解剖学的特徴構造を抽出する。
スライス面決定部243は、推奨撮影断面算出情報記憶部330に記憶されたアルゴリズムに従って、2つの互いに交差する2つのスライス面に平行な画像群から抽出した解剖学的特徴構造を用い、これらのスライス面いずれにも交差する新たなスライス面を決定する。
なお、計算機110は、CPUとメモリとを備え、計算機110が実現する制御部200の各機能は、CPUが記憶装置112に格納されたプログラムをメモリにロードして実行することにより実現される。また、記憶部300は、記憶装置112上に実現される。なお、各機能の全てまたはその一部は、MRI装置100とは独立して設けられる汎用の情報処理装置であって、MRI装置100とデータの送受信が可能な情報処理装置により実現されていてもよい。同様に、記憶部300の全てまたはその一部は、MRI装置100とは独立して設けられ、MRI装置100とデータの送受信が可能な外部記憶装置により実現されてもよい。
以下、推奨撮影断面算出情報記憶部330に登録されるアルゴリズムに従って推奨撮影断面算出部241により実行される推奨撮影断面算出処理を、具体例を用いて説明する。ここでは、頭部のルーチン検査の場合を例にあげて説明する。
頭部のルーチン検査では、本撮影は、T1強調画像、T2強調画像、FLAIR画像、拡散強調画像を取得するとする。撮影断面は、正中面に垂直で、OMラインに平行な傾きを持ち、全脳を覆う領域に設定される。OMラインとは、正中面の画像上の鼻根部と橋端下部とを結ぶラインであり、正中面に近い平面上の画像であるほど認識しやすいラインである。このため、撮影断面は、正中面画像上で決定する。
従って、最適スカウト撮影情報記憶部320には、頭部のルーチン検査の最適スカウト撮影情報として、正中面を決定し、正中面画像を得ることができるものが登録される。例えば、T1強調のGrE(グラディエントエコー)系のパルスシーケンスを使用し、コロナル面、アキシャル面、サジタル面の順番で各面に平行な複数のスライスの撮影を行うよう登録される。ここでのスライス枚数は限定されないが、以下、5枚のスライスの撮影を行う場合を例にあげて説明する。図3に、このパルスシーケンス使用して取得したスカウト画像の一例を示す。図3(a)がコロナル面に平行なスライスの画像(コロナル画像)301で、図3(b)がアキシャル面に平行なスライスの画像(アキシャル画像)302で、図3(c)がサジタル面に平行なスライスの画像(サジタル画像)303である。このとき用いた撮影パラメータは、TR/TE=30/1.5、FOV=270、FA=45、画素数128×128、スライス厚10mmである。それぞれの画像の左下には、座標軸304を示す。なお、以下、本明細書では、体軸方向をz軸方向、被検体103が水平にテーブル115上に載置されるものとし、z軸に垂直な2方向のうち、水平面に平行な方向をx軸方向、水平面に垂直な方向をy軸方向とする。
また、推奨撮影断面算出情報記憶部330には、コロナル画像およびアキシャル画像から正中面を決定し、サジタル画像から正中面の画像を生成し、その正中面の画像上で推奨撮影断面を算出するアルゴリズムが登録される。そして、正中線および正中面画像上の頭部輪郭、脳輪郭、脳梁、橋、脳幹、下垂体、斜台が解剖学的特徴構造として登録され、画像処理としては、各スカウト画像から正中線を抽出するもの、正中面の画像上で推奨撮影断面を特定するものが登録される。
頭部ルーチン検査における、推奨撮影断面算出情報記憶部330に登録されたアルゴリズムに従って推奨撮影断面算出部241により実行される推奨撮影断面算出処理の流れを、図4に示す。解剖学的特徴構造抽出部242は、コロナル画像上で正中線の傾きを計算する(ステップS401)。次に、推奨撮影断面算出部241は、コロナル画像上で、画像処理に適するアキシャル画像の位置を決定する(ステップS402)。そして、解剖学的特徴構造抽出部242は、ステップS402で決定された位置のアキシャル画像上で、正中線の位置を計算する(ステップS403)。そして、スライス面決定部243は、ステップS401で得た正中線の傾きと、ステップS402で得た正中線とから正中面を決定する(ステップS404)。
次に、推奨撮影断面算出部241は、サジタル画像を参照し、ステップS404で得た正中面の妥当性を解析し(ステップS405)、妥当であると判別した場合(ステップS406)、サジタル画像から正中面画像を作成する(ステップS407)。そして、推奨撮影断面算出部241は、正中面画像上で、推奨撮影断面を特定し(ステップS408)、処理を終了する。一方、ステップS406において、妥当でないと判別した場合は、ステップS407以降の処理を行わず、処理を終了する。
次に、推奨撮影断面算出処理における、各画像処理の詳細について説明する。これらの画像処理の手法は、予め推奨撮影断面算出情報記憶部330に登録される。まず、ステップS401における、解剖学的特徴構造抽出部242によるコロナル画像上で正中線の傾きを計算する処理(以下、コロナル正中線算出処理)について説明する。図5(a)は、本実施形態のコロナル正中線算出処理を説明するための図である。
コロナル正中線算出処理では、正中線の大まかな傾きを算出する。まず、取得した5枚のコロナル画像を全て加算し、加算画像501を作成する。そして、閾値により二値化処理を行い、加算画像501上で、頭部領域502とバックグラウンド領域503とを分離する。ここで、閾値には、例えば、加算画像上の全画素の画素値の平均値等を用いる。閾値以上の値を有する画素領域を頭部領域502、閾値未満の値を有する領域をバックグラウンド領域503とする。
そして、頭部領域502内のピクセル座標を抽出し、この座標値をプロットし(504)、一次関数で最小二乗フィッティングを行う。このフィッティングで得られた一次関数の傾き505を、コロナル画像上における正中線の大まかな傾きとする。これは、コロナル面において、頭部形状がほぼ左右対称である性質を利用した手法であり、画像のコントラストに左右されない。
ただし、コロナル正中線算出処理はこれに限られない。画像処理の各種の手法を用いることができる。例えば、画素値、画素の微分値の組み合わせによる評価関数を作成し、正中線を算出し、傾きを決定してもよい。
次に、図5(b)を用い、ステップS402の、推奨撮影断面算出部241による、画像処理に適するアキシャル画像の位置を決定する処理(アキシャル画像決定処理)を説明する。ここでは、図5(a)で作成した加算画像501を用いる。
まず、加算画像501の、x軸への一次元投影像511およびz軸への一次元投影像521を作成する。そして、x軸への一次元投影像511から、頭部を決定し、その両端のx座標とその間の幅512とを決定する。x座標と幅512とは、閾値処理により算出する。すなわち、一次元投影像511においてその画素値が閾値513以上の範囲の画素群を頭部の幅512と判別する。閾値513には、加算画像501上の全画素の画素値の平均値等、予め定められた手法で決定される値を用いる。
z軸への一次元投影像521については、その逆数の微分値522を計算し、その最大値となる点を求め、頭頂点531のz座標を決定する。そして、頭頂点531のx座標は、頭部の両端のx座標の中点とする。この加算画像501上で、頭頂点531から、上記コロナル正中線算出処理で決定した正中線の傾きに沿って幅512だけ離れた点532の座標を算出する。頭頂点531と点532との間の、所定の範囲内のスライス面のアキシャル画像を、画像処理に適するアキシャル画像として選択する。所定の範囲としては、例えば、頭頂点531と点532との距離(頭部幅512)をLとすると、頭頂点531から、L/4の位置と3L/4の位置との範囲533内で撮影されるアキシャル画像を1枚以上選択する。アキシャル画像の選択においては、撮影パラメータとして先に登録されている、5枚のアキシャル画像を取得するアキシャル面に平行なスライス面の各スライス位置534を用いる。
次に、ステップS403における、解剖学的特徴構造抽出部242による、アキシャル画像上で正中線の位置を計算する処理(アキシャル正中線算出処理)について、説明する。図6(a)〜(c)は、アキシャル正中線算出処理を説明するための図である。図6(a)に示すように、アキシャル画像選択処理において、選択されたアキシャル画像において、上記コロナル正中線算出処理と同様の手法で、閾値により二値化処理を行い、頭部領域とバックグラウンド領域とを分離する。さらに、頭部領域とされた領域を、その画素値の大きさに応じて2つ以上のグループに分ける。グループごとに画素値として一定の値を割り振った階調調整画像601をグループ数作成する。階調調整画像601を生成する処理により、血管や脂肪の領域といった大きな画素値を平滑化できる。その後、階調調整画像601をそれぞれ加算し、加算画像602を作成する。
次に、加算画像602において、その画像の重心603を計算する。重心603は、一般的に用いられている画像処理の手法で算出する。例えば、画素値で重みをつけて算出する、加算画像602において頭部領域として抽出した領域の画素値は一様として座標のみで計算する、といった手法を用いる。
次に、加算画像602のx軸への一次元投影像604を作成し、アキシャル画像選択処理と同様に閾値を用い、頭部の幅605、および両端のx座標を決定する。
次に、図6(b)に示すように、加算画像602を、重心603を中心に、予め定めた基準ラインの両側−90度から+90度の範囲で予め定めた刻みで回転させ、基準ラインとの成す角度606毎にx軸への一次元投影像604を作成する。基準ラインは例えばy軸とする。各投影像604において、先に決定した頭部の両端のx座標の中心を中心として、頭部の幅605の1/2の幅607の範囲で、輝度が最小の値となる点608のx座標およびその点の輝度を抽出する。そして、抽出した輝度を、その投影像604の角度606に対応づけてプロット609する。
プロット609の結果から、各投影像604の中で、点608の輝度が最小となるものの角度606を特定する。それを、アキシャル画像における正中線の傾きα610とする。そのときの一次元投影像604上の、輝度が最小となる点608のx座標Xminと、重心の座標(X0,Y0)とを用い、以下の式(1)に従って、正中面が通過する点P611の座標(XP,YP)をそれぞれ算出する。
XP=X0−(Xmin−X0)×Cosα
YP=Y0−(Xmin−X0)×Sinα (1)
点P611を通り、傾きが角度α610である直線を、アキシャル画像における正中線612とする。
以上のアキシャル正中線算出処理は、アキシャル画像上の正中線に対応する領域の画素の画素値(以後、簡単に正中線の画素値と呼ぶ。他も同様)が回りの脳実質の画素値に対して低いことを利用した手法である。逆に、正中線の画素値が脳実質の画素値より高くなる場合は、一次元投影像において、輝度の高い値を算出する処理を行う。
また、上記アキシャル正中線算出処理では、画素値の大きさを評価関数として用いている。しかし、これに限られない。例えば、一次元投影像604を微分したエッジ情報を用いるよう構成してもよい。また、一次元投影像604の標準偏差の情報を用いてもよい。さらに、これらを複合して処理を行ってもよい。これらを複合することにより、より正確な正中線抽出が可能となる。
アキシャル画像が脳低部に近い画像の場合は、正中線抽出に異なる手法を用いてもよい。すなわち、図6(c)に示すように、重心を通るx軸に平行な線を用い、アキシャル画像を前壁側の領域621と後壁側の領域622とに分ける。次に、前壁側の領域621の画素値の総和を算出し、重心を中心として画像を、予め定めた基準ラインの両側−90度から+90度の範囲で予め定めた刻みで回転させ、回転刻み毎に画素値の総和を算出する。計算した総和が最小になる回転角を正中線の傾きとし、重心を正中面の通る点とする。この手法は、脳低に近いアキシャル画像では、鼻腔領域の画素値が低く、眼球および小脳領域の画素値が高いことを利用している。また、脳低に近いアキシャル画像では、頭部の傾きに依存せずに正中面がほぼ重心を通ることも利用している。なお、アキシャル画像が脳低部に近い画像であるか否かは、アキシャル画像の位置で判断できる。例えば、選択されたアキシャル画像が頭頂側から数えて三枚目以下である場合は、脳低部に近いと判断できる。
なお、ここでは、処理の一例を示しており、正中線の位置を抽出する方法はこの限りではない。他の方法としては、画素値や画素の微分値の組み合わせによる評価関数を作成することで正中線を算出する方法などがある。
次に、ステップS404における、スライス面決定部243による、コロナル画像上の正中線の傾きと、アキシャル画像上の正中線の位置とから、正中面を決定する正中面決定処理について説明する。図7(a)〜(c)は、この正中面決定処理を説明するための図である。
図7(a)に示すように、コロナル正中線算出処理で決定したコロナル画像501上の正中線の傾き505と、アキシャル正中線算出処理で決定したアキシャル画像602上の正中線612とにより、正中面の初期位置を決定する。すなわち、コロナル画像501上の正中線の傾き505とアキシャル画像602上の正中線612との両方を満たす面を正中面の初期位置(初期正中面)とする。そして、図7(b)および図7(c)に示すように、初期正中面と、スカウト撮影で取得した5枚のコロナル画像501a、501b、501c、501d、501eおよび5枚のアキシャル画像601a、601b、601c、601d、601eそれぞれの交線701を求める。
次に、交線701上の画素値を評価関数とし、評価関数が最小となるよう、初期正中面の傾きや位置を変更し、調整する。そして、評価関数が最小となる傾きおよび位置を有する面を正中面とする。なお、調整には、例えば、最小二乗法などを用いる。また、評価関数に、画素値の1階微分、2階微分などを組み合わせてもよい。これにより、抽出精度を向上させることができる。さらに、評価の精度を向上させるため、交線701は、ほぼ脳領域のみにおける交線であることが望ましい。なお、脳領域は、例えば、コロナル画像においては、頭頂点531と点532の間、アキシャル画像では、重心を中心に、直径が頭部の幅605である円の範囲とする。また、コロナル画像およびアキシャル画像を、それぞれ画素値の平均値を閾値として二値化処理し、信号領域とバックグラウンド領域とに分け、信号領域を脳領域としてもよい。
次に、上記ステップS407における、推奨撮影断面算出部241によるサジタル画像から正中面画像を作成する、正中面画像作成処理について説明する。まず、スカウト撮影で得られた5枚のサジタル画像の画素値を用い、ステップS404において決定された正中面の、各画素値を補間により求める。この際、5枚のサジタル画像の位置は、撮影パラメータにおいて、スカウト撮影時のサジタル面に平行な5枚のスライスのスライス位置として設定されたものを用いる。
なお、5枚のサジタル画像の中で、その頭部領域が正中面に70%以上含まれている画像があれば、その画像のスライス位置を正中面とするよう構成してもよい。このように構成することにより、処理時間を短縮することができる。
次に、上記ステップS408における、推奨撮影断面算出部241により作成した正中面画像上で、推奨撮影断面を特定する断面特定処理について説明する。ここでは、正中面画像上の解剖学的特徴点の位置を自動認識し、予め定められた解剖学的特徴点の位置と、撮影断面との関係の情報に従って、撮影断面を決定する。
まず、正中面画像上の解剖学的特徴点を自動認識する手法について説明する。自動認識は、固体毎の形状の違いに対応して変形可能なテンプレートモデルを用い、フィッティング処理により行う。使用するテンプレートモデルは、複数の点で構成され、各点は、正中面画像上の解剖学的特徴点に対応する。このテンプレートモデルを正中面画像にフィッティングし、テンプレートモデルを構成する各点の位置を算出し、正中面画像上の解剖学的特徴点の位置を抽出する。
まず、使用するテンプレートモデル800の作成手順を図8(a)〜(c)を用いて説明する。まず、複数の被検体の正中面画像810を取得する。これらの正中面画像810において、頭部輪郭、脳梁、大脳下縁、下垂体、斜台、橋、延髄、小脳下縁等の境界など、複数の解剖学的特徴点(目印点)801の座標を抽出する。ここでの被検体は、健常者であることが好ましく、また、被検体の数は、抽出する目印点801の数に比例することが好ましい。抽出した各被検体の目印点801の分布の大きさや傾きは、回転、スケーリング、平行移動により標準化する。標準化した目印点801の位置に関し、複数の被検体間での平均値を算出し、これを図8(b)に示すように、標準モデル820とする。
次に、各被検体および標準モデル820の全目印点801の座標を一次元ベクトル化し、分散共分散行列を作成後、行列の固有値および固有ベクトルを算出する。ここで算出される固有ベクトルは、被検体の目印点801の分布において個体差が見られる方向を示すベクトルとなる。固有ベクトルをPi(i=1,…,N。Nは分散共分散行列のマトリックス数に依存)、標準モデル820の一次元ベクトルをXVCと表すと、個体差の変化傾向をパラメータ化したテンプレートモデルの一次元ベクトルXは以下の式(2)で表される。
Figure 0005230809
ここで、bは変化の大きさを決めるパラメータである。このように、この標準モデル820は、bの値を調整することにより、形状の個体差に対応できる。図8(c)に示すように、テンプレートモデル800は、bの値を調節し、標準モデル820の形状を変化させて得る。
なお、bの値は任意に設定できる。例えば、固有値をλとしたとき、bの範囲を以下の式(3)で表される範囲に設定してもよい。
−3√λ<b<3√λ (3)
このように、bの値に制限を設けることにより、標準モデル820を変形させてテンプレートモデル800を作成する際、初期形状から大きく逸脱することを防止できる。また、固有値が大きい固有ベクトルほど個体差が大きい傾向を示すため、固有値の大きな固有ベクトルのみ標準モデル820に組み込むよう構成しても良い。このように構成することにより、計算時間を短縮できる。
また、使用するテンプレートモデル800は、個体差を無視した固定形状のテンプレートモデルを用いてもよい。これにより、大まかな位置を高速に計算できる。
なお、テンプレートモデル800は、予め推奨撮影断面算出情報記憶部330に登録される標準モデル820を用い、操作者が作成する。また、予め作成し、推奨撮影断面算出情報記憶部330に保存しておいてもよい。テンプレートモデル800は、人種、性別、年齢等毎に作成し、保持してもよい。これにより、より精度の高い特徴抽出が可能となる。テンプレートモデル800を複数用意する場合は、被検体の情報として検査時に登録する情報のいずれかと対応付けて管理し、検査時に自動的に設定するよう構成してもよいし、使用するテンプレートモデル800を操作者が選択可能なインタフェースを設けてもよい。使用するテンプレートモデル800は、いつでも自由に入れ替え可能であるよう構成してもよい。
次に、上記テンプレートモデル800によるフィッティング行う際の初期配置について、図9(a)、(b)を参照して説明する。まず、図9(a)に示すように、正中面画像901の重心902を計算する。次に、z軸への一次元投影像903を作成し、得られた一次投影像903において頭部の長さ904を計算する。重心902を中心に頭部の長さ904で定まる範囲において、閾値処理により前壁側の輪郭座標905を抽出する。抽出した輪郭座標905を一次関数906で最少二乗フィッティングする。次に、得られた一次関数906に重心902より垂線を下ろし、垂線との交点907を中心に頭部の長さ904で定まる範囲において、閾値処理により再び前壁側の輪郭座標908を抽出する。さらに、抽出した輪郭座標908を再び一次関数910で最少二乗フィッティングし、その傾き911を求める。なお、一次関数910の傾き911は、閾値処理で頭部領域の座標を抽出し、一次関数でフィッティングする方法などで求めてもよい。次に、y軸への一次元投影像912を作成し、閾値処理にて頭部の幅913を算出する。図9(b)に示すように、テンプレートモデル800を正中面画像上に初期位置920に配置するときは、このようにして得られた重心902、傾き910、頭部の幅904に基づき、正中面画像上の同じ目印点801に近い位置で配置させるようにする。
次に、上記初期位置920に配置されたテンプレートモデル800を、正中面画像にフィッティングする手法について説明する。フィッティングは、テンプレートモデル800の各構成点である目印点801が、正中面画像上の対応する組織の解剖学的特徴点に収束するよう、初期位置920からその位置を少しずつ更新する。具体的には、それぞれの目印点801について、当該目印点801が移動すべき更新点を検索し、最小二乗法によりテンプレートモデル800全体を変形させる処理を繰り返し行う。更新点は、検索方向の直線上の画素値を評価関数として、評価関数の特異点の検索により決定する。検査方向や検索する幅は特徴点を有する組織に応じて設定する。頭部の場合、フィッティングは、まず、頭部輪郭上の点について、更新点の検索および回転、スケーリング、平行移動によりテンプレートモデル変形の処理を繰り返し、その後、テンプレートモデル800上の全ての点について、更新点の検索および回転、スケーリング、平行移動によるモデル変形の処理を繰り返す。最後に、テンプレートモデル800上の全ての点において、固体差を示すパラメータによる変形も加え、モデル変形の処理を繰り返す。
これらの処理により、テンプレートモデル800の各目印点801は、図10(a)に示すように、処理を実行している正中面画像901上に対応する解剖学的特徴点に収束する。そして、このときのテンプレートモデル800の各目印点801の座標が、撮影断面を決定するために必要な、解剖学的特徴点の座標となる。このようにして、正中面画像上の解剖学的特徴点を自動認識する。その後、推奨撮影断面算出情報記憶部330に登録される解剖学的な特徴構造と推奨撮影断面との位置関係に従って、推奨撮影断面を決定するために必要な解剖学的特徴点を抽出し、推奨撮影断面を決定する。例えば、頭部では、鼻根部、橋、頭部輪郭、脳輪郭、脳梁、脳幹を決定する点などが解剖学的特徴構造として登録され、鼻根部と橋端下部とを結ぶOMライン1002の傾きが推奨撮影断面の傾き、全脳の範囲1003が推奨撮影断面の範囲と登録される。これに基づき、図10(b)に示すように、正中面画像上におけるフィッティング後のテンプレートモデル800の各目印点801の位置より推奨撮影断面の傾き1004と範囲1005とを決定する。
なお、図10では上下方向の撮影範囲のみを示したが、前後方向の撮影範囲に関しても同様に位置情報を取得できる。また、OMラインに平行な面を推奨撮影断面の傾きとして決定しているが、これに限られない。例えば、頭部検査では、正中面画像上の前交連を結ぶAC−PCライン、脳幹に平行なラインなどが推奨撮影断面の傾きとして登録されていてもよい。
以上が、推奨撮影断面算出情報記憶部330に登録されるアルゴリズムに従って推奨撮影断面算出部241により実行される推奨撮影断面算出処理である。次に、この推奨撮影断面算出処理を組み込んだ、本実施形態のスカウト撮影開始操作から、本撮影開始操作までの計測準備処理の流れを、具体例を用いて説明する。先の推奨撮影断面算出処理と同様、頭部のルーチン検査の場合を例にあげて説明する。図11から図13は、本実施形態の計測準備処理を説明するための図であり、図11は、スカウト撮影の、図12および図13は、それ以降の処理の流れを説明するための図である。ここでは、操作者の動作、計算機110における処理、シーケンサ104の動作に分けて説明する。なお、頭部ルーチン検査のプロトコルは以下のとおりである。スカウト撮影を行い、その後、シミングと感度マップ補正のためのデータを取得するパルスシーケンスで構成される準備撮影を行う。そして、スカウト撮影の結果により決定された撮影断面で本撮影を行う。
操作者がスカウト撮影開始の指示を入力装置116を介して行う(ステップS1101)と、計算機110では、UI制御部210が操作者からの指示を受け付ける。そして、計測制御部230が、最適スカウト撮影情報記憶部320に記憶されている最適スカウト撮影情報のパルスシーケンスおよび撮影パラメータに従って、シーケンサ104に指令を出力する(ステップS1201)。シーケンサ104は、計算機110からの指示に従って、各部を動作させ、スカウト撮影を開始し(ステップS1301)コロナル面、アキシャル面、サジタル面の順番で、順次計測を行い、それぞれ得られたエコー信号を順に計算機110に送信する(ステップS1302、S1303、S1304)。
計算機110では、コロナル面の計測によるエコー信号を受け取ると、信号処理部220が、画像を再構成し、コロナル画像を得る(ステップS1202)。そして、得られたコロナル画像上で、推奨撮影断面算出部241は、正中線の傾きを計算する(ステップS1203)。また、推奨撮影断面算出部241は、コロナル画像上で、画像処理に適するアキシャル画像の位置を決定する(ステップS1204)。なお、本実施形態の推奨撮影断面算出部241は、計測制御部230によるアキシャル面の計測が完了する前に、ステップS1204までの処理を行う。
次に、計算機110では、アキシャル面の計測によるエコー信号を受け取ると、信号処理部220が、画像を再構成し、アキシャル画像を得る(ステップS1205)。そして、推奨撮影断面算出部241は、得られたアキシャル画像の中の、ステップS1204で決定した位置にあるアキシャル画像上で正中線を特定する(ステップS1206)。そして、推奨撮影断面算出部241は、ステップS1204で算出したコロナル画像上の正中線の傾きと、ステップS1206で算出したアキシャル画像上の正中線の位置とから、正中面を特定する(ステップS1207)。なお、本実施形態の推奨撮影断面算出部241は、計測制御部230によるサジタル面の計測が完了する前に、ステップS1207までの処理を行う。
計算機110では、サジタル面の計測によるエコー信号を受け取ると、信号処理部220が、画像を再構成し、サジタル画像を得る(ステップS1208)。ここで、UI制御部210は、再構成されたスカウト画像(コロナル画像、アキシャル画像、サジタル画像)とともに、スカウト撮影を再度行うか否かの指示を受け付ける表示を、表示装置111に表示させる(ステップS1209)。操作者は、表示装置111に表示されたスカウト画像を確認し(ステップS1102)、スカウト画像の取り直しの有無を指示する(ステップS1103)。一方、計算機110では、推奨撮影断面算出部241が、サジタル画像を参照しながら、ステップS1207で得た正中面の妥当性を解析し、判別する(ステップS1210)。取り直しの有無の指示は、例えば、UI制御部210が、表示装置111上に取り直しおよび処理続行の指示を受け付ける操作ボタンをスカウト画像とともに表示させ、その選択等によって受け付ける。
ここで、ステップS1210における、正中面の妥当性は、ステップS1208において再構成したサジタル画像から、正中面画像が作成可能であるか否か、ステップS1207で特定された正中面が解剖学的に妥当か否かにより決定する。正中面画像が作成可能であり、かつ、解剖学的に妥当である場合、正中面は妥当であると判別する。なお、作成可能であるか否かは、例えば、ステップS1207で特定された正中面が、ステップS1208で取得した5枚のサジタル画像の撮影範囲に含まれる割合で判別する。すなわち、計算により得た割合が、所定以上(例えば、75%以上)である場合、作成可能と判別し、それ以外は否と判別する。また、解剖学的に妥当であるか否かは、正中面位置が、コロナル画像およびアキシャル画像から決定される頭部中心から大きく外れた位置(例えば、頭部中心から正中面までの距離がコロナル画像から算出される頭部の幅の1/4以上)にあるか否かで判別する。すなわち、大きく外れた位置に該当する場合、妥当でないと判別する。
次に、ステップS1103で、スカウト撮影の取り直しは無し、すなわち、そのまま計測準備処理を進める指示を受け付けた場合の処理を、図12(a)および(b)を用いて説明する。図12(a)は、さらに、ステップS1210で、正中面が妥当であると判別された場合の処理の流れである。
本図に示すように、ステップS1103でスカウト撮影の取り直しは無しとの指示を受け付け、かつ、ステップS1210で正中面が妥当であると判別された場合、計算機110では、推奨撮影断面算出部241が、サジタル画像から、正中面の位置の画像(正中面画像)を生成し(ステップS1222)、生成した画像上で推奨撮影断面を算出する(ステップS1223)。その間、操作者は、準備撮影開始の指示を行う(ステップS1121)。計算機110では、UI制御部210が操作者からの指示を受け付け、計測制御部230が、準備撮影用として予め保持されているパルスシーケンスおよび撮影パラメータに従って、シーケンサ104に指令を出力する(ステップS1221)。シーケンサ104は、計算機110からの指示に従って、各部を動作させ、準備撮影を行う(ステップS1321)。なお、ここでは、順撮影開始後、計算機110がステップS1222およびS1223の処理を行うよう記載しているが、計算機110では、ステップS1222およびS1223の処理を、操作者の準備撮影開始の指示を待たずに進める。
操作者が本撮影の撮影断面設定画面を開く操作を行うと(ステップS1122)、計算機110では、ステップS1223で算出した推奨撮影断面を特定する情報、例えば、位置、傾き、本撮影で設定されている撮影パラメータ(スライス枚数、スライス厚、FOV等)に基づく推奨撮影断面範囲、スカウト画像等をUI制御部210が表示装置110に表示させる(ステップS1224)。操作者は、必要に応じて、表示された推奨撮影断面を調整する入力を行う(ステップS1123)。計算機110では、UI制御部210が調整の入力を受け付け、撮影断面決定部240が受け付けた調整量に従って、撮影断面位置を調整し、本撮影の最終的な撮影断面を決定する(ステップS1225)。このとき、本撮影の撮影断面に関連する撮影パラメータを決定した撮影断面に従って更新する構成でもよい。操作者から本撮影開始の指示を(ステップS1124)受け付けると、計算機110では、本撮影用として予め保持されているパルスシーケンスおよび撮影パラメータに従って、計測制御部230がシーケンサに指令を出し(ステップS1226)、シーケンサ104は、計算機110からの指示に従って、各部を動作させ、本撮影を行う(ステップS1322)。
なお、ステップS1225で操作者からの調整の入力を受け付けた場合、計算機110では、撮影断面決定部240が調整後の撮影断面を推奨撮影断面として表示装置111に表示させ、再度、調整の入力を受け付けるよう構成してもよい。この場合、操作者から本撮影開始の指示を受け付けるまで、調整の入力が可能なように表示装置111に調整後の撮影断面を推奨撮影断面として表示することを繰り返す。
次に、ステップS1103でスカウト撮影の取り直しは無しとの指示を受け付け、かつ、ステップS1210で妥当でないと判別された場合の処理を、図12(b)を用いて説明する。計算機110では、推奨撮影断面算出処理をそこで終了する。一方、操作者は、準備撮影開始の指示を行う(ステップS1131)。計算機110では、UI制御部210が操作者からの指示を受け付け、計測制御部230が、準備撮影用として予め保持されているパルスシーケンスおよび撮影パラメータに従って、シーケンサ104に指令を出力する(ステップS1231)。シーケンサ104は、計算機110からの指示に従って、各部を動作させ、準備撮影を行う(ステップS1331)。
操作者が本撮影の撮影断面設定画面を開く操作を行うと(ステップS1132)、計算機110では、UI制御部210が、従来同様、スカウト画像と設定された撮影パラメータ(スライス枚数、スライス厚、FOV等)を元に決定される推奨撮影断面位置を軸断で表示装置110に表示させる(ステップS1232)。操作者は、表示された軸断上で撮影断面の入力を行う(ステップS1133)。計算機110では、UI制御部210が調整の入力を受け付け、撮影断面決定部240が受け付けた入力に従って、撮影断面を決定し、本撮影の撮影断面に関連する撮影パラメータを決定した撮影断面に従って更新する(ステップS1233)。ここでも、上記図12(a)の場合と同様、何度でも調整の入力が可能なよう構成してもよい。操作者が本撮影開始の指示を行うと(ステップS1134)、計算機110では、本撮影用として予め保持されているパルスシーケンスおよび撮影パラメータに従って、計測制御部230がシーケンサに指令を出し(ステップS1234)、シーケンサ104は、計算機110からの指示に従って、各部を動作させ、本撮影を行う(ステップS1332)。
次に、ステップS1103において、スカウト撮影を取り直す旨の指示を受け付けた場合の処理を図13を用いて説明する。図13(a)は、さらに、ステップS1210の解析において、正中面が解剖学的に妥当であると判別された場合の計測準備処理の流れを説明するための図である。このとき、計算機110では、推奨撮影断面算出部241が、ステップS1207で特定した正中面を取得可能なスカウト撮影用の各断面を推奨スカウト撮影断面として算出する(ステップS1241)。そして、操作者がスカウト撮影断面設定画面を開く操作を行うと(ステップS1142)、計算機110では、UI制御部210がそれを受け付け、表示装置111にスカウト画像と算出した推奨スカウト撮影断面を特定する情報を表示する(ステップS1242)。ここで、操作者がスカウト撮影設定画面を開く操作は必ずしも行う必要はなく、撮り直す旨の指示を受けた際に、自動でスカウト撮影設定画面に切り替え、スカウト画像と推奨スカウト撮影位置を表示させる構成でもよい。操作者は、表示された推奨スカウト撮影断面を調整する入力を行う(ステップS1141)。計算機110では、UI制御部210が、調整の入力を受け付け、撮影断面決定部240が、受け付けた調整量に従って、推奨スカウト撮影断面を調整し、最終的なスカウト撮影断面を決定する(ステップS1243)。このとき、撮影断面に関連する撮影パラメータを決定した撮影断面に従って更新する構成でもよい。なお、ここでも、上記図12(a)の場合と同様、何度でも調整の入力が可能なよう構成してもよい。その後、それぞれ、ステップS1101、S1201、S1301以降の処理を繰り返す。
次に、ステップS1103において、スカウト撮影を取り直す旨の指示を受け付け、さらに、ステップS1210の解析において、正中面が解剖学的に妥当でないと判別された場合の計測準備処理の流れを図13(b)に従って説明する。このとき、操作者がスカウト撮影断面設定画面を開く操作を行うと(ステップS1151)、計算機110では、UI制御部210がそれを受け、表示装置111に3軸直交断面を表示する(ステップS1251)。操作者は、スカウト撮影断面を決定するため、表示された3軸直交断面を調整する入力を行う(ステップS1152)。計算機110では、UI制御部210が、調整の入力を受け付け、撮影断面決定部240が、受け付けた調整量に従って、スカウト撮影断面を決定する(ステップS1252)。ここで、撮影断面に関連する撮影パラメータを決定した撮影断面に従って更新する構成でもよい。なお、ここでも、上記図12(a)の場合と同様、何度でも調整の入力が可能なよう構成してもよい。その後、それぞれ、ステップS1101、S1201、S1301以降の処理を繰り返す。
以上が、頭部ルーチン検査の場合の、スカウト撮影開始操作から本撮影開始操作までの処理の流れである。
すなわち、頭部ルーチン検査では、計測制御部230と信号処理部220とは、互いに交差する2つのスライス面の内の一方のスライス面である第一のスライス面をコロナル面とし、このコロナル面に平行な1枚以上の二次元画像からなる第一の画像群と、他方のスライス面である第二のスライス面をアキシャル面として、このアキシャル面に平行な1枚以上の二次元画像からなる第二の画像群とを取得する。また、第一のスライス面および第二のスライス面の双方に交差する第四のスライス面をサジタル面とし、このサジタル面に並行な1枚以上の二次元画像からなる第四の画像群を取得する。
そして、解剖学的特徴構造抽出部242は、第一の画像群から、第一の解剖学的特徴構造として正中線の傾きを、第二の画像群から正中線を決定する。そして、スライス面決定部243は、解剖学的特徴構造抽出部242が決定した正中線の傾きおよび正中線から、第三のスライス面として正中面を決定する。
推奨撮影断面算出部241は、第四の画像群から決定した正中面の二次元画像を第三の画像として生成する。そして、推奨撮影断面が含むべき複数の解剖学的特徴点から構成されるテンプレートモデル800を用い、第三の画像上でこれらの解剖学的特徴点の位置を特定し、推奨撮影断面算出情報記憶部330に予め登録される解剖学的な特徴構造と推奨撮影断面との位置関係に基づき、推奨撮影断面を決定する。
以上説明したように、本実施形態によれば、スカウト撮影に平行して画像処理を実行し、推奨撮影断面を算出する。また、実行するスカウト撮影は、撮影断面をスカウト画像から手動で設定する場合と同様のものである。従って、従来の検査の流れを変更することなく、また、新たな処理を追加することによる処理時間の延長無しに、操作者に推奨撮影断面を提示することができる。このため、自動的に推奨撮影断面を得るにあたり、作業性が悪化することがない。
また、スカウト撮影の取り直しを行う場合であっても、図13(a)に示すように、より適切なスカウト撮影断面を設定後に取り直しを行うため、効率的にスカウト撮影の取り直しを行うことができる。さらに、上記図12(b)および図13(b)のように、推奨撮影断面または推奨スカウト撮影断面を算出しない場合であっても、処理時間の延長がなく、従来の手順と同じであるため、検査効率が悪化することはない。
なお、上記実施形態では、頭部のように、一般に、スカウト画像を取得する直交する3つのスライス面とは大きく異なるオブリーク面を撮影断面に設定する場合の推奨撮影断面算出処理を例にあげて説明した。ここで、撮影断面が、スカウト画像を取得する3つのスライス面のいずれかに略合致する面となる場合の、推奨撮影断面算出処理を、具体例をあげて説明する。ここでも、推奨撮影断面算出処理は、推奨撮影断面算出情報記憶部330に登録されるアルゴリズム、画像処理、解剖学的特徴構造に従って推奨撮影断面算出部241が実行する。以下、それぞれ説明する。
まず、腰椎検査の場合について説明する。腰椎検査では、サジタル面(腰椎サジタル面)のT1強調画像およびT2強調画像と、椎間板ラインのT1強調画像およびT2強調画像とを取得する。すなわち、撮影断面は、腰椎サジタル面と、椎間板ラインを含む面(椎間板ライン面)の2つである。腰椎検査サジタル面は、脊髄神経に平行なスライス面とし、椎間板ライン面は、第1から第5腰椎それぞれの下に位置する各椎間板の傾きに合わせた面とする。従って、腰椎検査サジタル面は、腰椎位置および脊髄神経の傾きにより決定される。また、椎間板ライン面は、椎間板の位置および脊髄神経の傾きにより決定される。
従って、最適スカウト撮影情報記憶部320には、腰椎検査の最適スカウト撮影情報として、例えば、プロトン強調のGrE系パルスシーケンスを用い、アキシャル面、コロナル面、サジタル面の順に、各面に平行な複数のスライス(例えば、5スライス)の撮影を行うよう登録される。また、推奨撮影断面算出情報記憶部330には、撮影断面処理では、アキシャル画像上で腰椎位置を決定し、コロナル画像上で脊髄神経の傾きを決定し、これらから、腰椎サジタル面を決定するアルゴリズムと画像処理と解剖学的特徴構造、および、腰椎検査サジタル面画像上で撮影対象の椎間板の位置を決定し、先にコロナル画像上で決定した脊髄神経の傾きを用い、椎間板ライン面を決定するアルゴリズムと画像処理と解剖学的特徴構造が、登録される。
以下、腰椎検査における推奨撮影断面算出部241による推奨撮影断面算出処理の詳細について説明する。
まず、スカウト撮影の最初に取得するアキシャル画像において、画像の重心の計算などにより腰椎位置、すなわち、脊髄神経の位置を決定する。次に、コロナル画像から画素値の評価関数を用いて脊髄神経の傾きを求める。そして、アキシャル画像で求めた脊髄神経の位置を通り、コロナル画像で求めた脊髄神経の傾きを有する面を、腰椎サジタル面の推奨撮影断面とする。なお、ここまでの処理を、本撮影前までに完了させる。
次に、スカウト撮影で得たサジタル画像から、腰椎検査サジタル面の画像を補間等により作成する。なお、腰椎検査サジタル面の画像は、本撮影により取得したものを用いてもよい。そして、腰椎検査サジタル面画像上でエッジ強調や画素値の評価関数に基づき椎間板の位置を特定する。求めた椎間板の位置を通り、脊髄神経の傾きに直交する面を、椎間板ライン面の推奨撮影断面とする。
次に、膝検査の場合について説明する。膝検査では、コロナル面(膝コロナル面)のT2*強調画像、T1強調画像と、サジタル面(膝サジタル面)のT2*強調画像、T1強調画像と、前十字靭帯の診断面のT2強調画像とを取得する。すなわち、撮影断面は、膝コロナル面と、膝サジタル面と、前十字靭帯の診断面と、の3面である。
最適スカウト撮影情報記憶部320には、膝検査の最適スカウト撮影情報として、GrE系パルスシーケンスを用い、アキシャル面、サジタル面、コロナル面の順番で、それぞれ、複数のスライス(例えば、5スライス)の撮影を行うよう登録される。
以下、推奨撮影断面算出情報記憶部330に膝に対応づけて登録されるアルゴリズムおよび画像処理に従って、推奨撮影断面算出部241により実行される、膝検査における推奨撮影断面算出処理の詳細について説明する。
まず、スカウト撮影の最初に取得するアキシャル画像において、大腿骨内側顆と大腿骨外側顆とを結ぶ線を特定する。次に、サジタル画像において、大腿骨と脛骨の関節面の接線に垂直なラインを特定する。そして、大腿骨内側顆と大腿骨外側顆とを結ぶ線に平行で、かつ、大腿骨と脛骨の関節面の接線に垂直な平面を膝コロナル面の推奨撮影断面とする。なお、推奨撮影断面の決定にあたっては、頭部ルーチン検査における正中面決定処理と同様の手法等を用いる。これは、以下においても同様である。
次に、スカウト撮影で得たアキシャル画像から大腿骨内側顆と大腿骨外側顆とを結ぶ線に垂直な線を特定する。また、スカウト撮影で得たコロナル画像から、膝コロナル面の画像(膝コロナル画像)を補間等により作成する。なお、膝コロナル画像は、本撮影により取得したものを用いてもよい。そして、膝コロナル面画像上で大腿骨と脛骨の関節面の接線を特定する。そして、大腿骨と脛骨の関節面の接線に垂直で、かつ、大腿骨内側顆と大腿骨外側顆とを結ぶ線に垂直な平面を、膝サジタル面の推奨撮影断面とする。
次に、スカウト撮影で得たサジタル画像から、膝サジタル面の画像(膝サジタル画像)を補間等により作成する。なお、膝サジタル画像は、本撮影により取得したものを用いてもよい。膝サジタル画像から前十字靭帯が描写されている画像を選択し、膝サジタル画像上の前十字靭帯に沿ったラインを特定する。次に、膝コロナル画像から前十字靭帯が描写されている画像を選択し、膝コロナル画像上の前十字靭帯に沿ったラインを特定する。そして、膝サジタル画像上の前十字靭帯に平行で、かつ、膝コロナル画像上の前十字靭帯に平行な平面を、前十字靭帯の診断面の推奨撮影断面とする。
次に、肩検査の場合について説明する。肩検査では、コロナル面(肩コロナル面)のT1強調画像、T2強調画像と、サジタル面(肩サジタル面)のT1強調画像、T2強調画像と、アキシャル面(肩アキシャル面)のT1強調画像、T2強調画像とを取得する。すなわち、肩検査の撮影断面は、肩コロナル面と肩サジタル面と肩アキシャル面の3面である。
最適スカウト撮影情報記憶部320には、肩検査の最適スカウト撮影情報として、GrE系パルスシーケンスを用い、アキシャル面、サジタル面、コロナル面の順番で、それぞれ複数のスライス(例えば、5スライス)の撮影を行うよう登録される。
以下、推奨撮影断面算出情報記憶部330に肩に対応づけて登録されるアルゴリズムおよび画像処理に従って、推奨撮影断面算出部241により実行される、肩検査における推奨撮影断面算出処理の詳細について説明する。
まず、スカウト撮影の最初に取得するアキシャル画像から、棘上筋が描写されている画像を選択し、当該画像上で棘上筋に平行な線を特定する。次に、スカウト撮影で得たサジタル画像上で上腕骨に沿ったラインを特定する。そして、棘上筋に平行で、かつ、上腕骨に平行な平面を、肩コロナル面の推奨撮影断面とする。
次に、スカウト撮影で得たアキシャル画像から骨頭と肩甲骨とが描写されている画像を選択し、当該画像上で関節面の接線に垂直な線を特定する。また、スカウト撮影で得たコロナル画像から、肩コロナル面の画像(肩コロナル画像)を補間等により作成する。なお、肩コロナル画像は、本撮影により取得したものを用いてもよい。そして、肩コロナル画像上で骨頭と肩甲骨とを結ぶラインに垂直なラインを特定する。そして、骨頭と肩甲骨とを結ぶ線に垂直で、かつ、骨頭と肩甲骨の関節面の接線に垂直な平面を、肩サジタル面の推奨撮影断面とする。
次に、スカウト撮影で得たサジタル画像から、肩サジタル面の画像(肩サジタル画像)を補間等により作成する。なお、肩サジタル画像は、本撮影により取得したものを用いてもよい。そして、肩サジタル画像上で上腕骨に沿ったラインに垂直な直線を特定する。次に、肩コロナル画像上で肩峰と鎖骨とを結ぶ直線を特定する。そして、上腕骨に垂直で、かつ、肩峰と鎖骨とを結ぶ直線に平行で、かつ、関節面の接線に平行な平面を、肩アキシャル面の推奨撮影断面とする。
以上説明した各部位の推奨撮影断面算出処理を組み込んだ、検査全体の流れは上記頭部ルーチン検査の場合と同様である。なお、これらの検査部位において、正中面に相当する面は、腰椎の場合、腰椎サジタル面、膝の場合、膝コロナル面、肩の場合、肩コロナル面である。従って、ステップS1210では、これらの面の妥当性を解析し、判別する。
すなわち、腰椎検査の場合、計測制御部230と信号処理部220とは、第一のスライス面をアキシャル面とし、第二のスライス面をコロナル面とし、第四のスライス面をサジタル面とし、これらに平行な1枚以上の二次元画像群を、それぞれ、第一の画像群、第二の画像群、第四の画像群として取得する。そして、解剖学的特徴構造抽出部242は、第一の解剖学的特徴構造としてアキシャル画像上の脊髄神経の位置、第二の解剖学的特徴構造としてコロナル画像上の脊髄神経の傾きを、それぞれ、抽出する。スライス面決定部243は、脊髄神経の位置を通り脊髄神経の傾きを有する面を、腰椎サジタル面、すなわち、第三のスライス面と決定する。推奨撮影断面算出部241は、第三の画像群である第三のスライス面の画像を補間により作成する。解剖学的特徴構造抽出部242は、腰椎サジタル面上で椎間板の位置を第三の解剖学的特徴構造として抽出する。そして、スライス面決定部243は、第一の解剖学的特徴構造または第一の解剖学的特徴構造と、第三の解剖的特徴構造とから、第五のスライス面を決定し、推奨撮影断面算出部241は、第五のスライス面を推奨撮影断面とする。
また、膝検査の場合、計測制御部230と信号処理部220とは、第一のスライス面をアキシャル面とし、第二のスライス面をサジタル面とし、第四のスライス面をコロナル面とし、これらに平行な1枚以上の二次元画像群を、それぞれ、第一の画像群、第二の画像群、第四の画像群として取得する。そして、解剖学的特徴構造抽出部242は、第一の解剖学的特徴構造としてアキシャル画像上の大腿骨内側顆と大腿骨外側顆とを結んだ線を、第二の解剖学的特徴構造としてサジタル画像上の大腿骨と脛骨の間接面の接線の垂線を、第四の解剖学的特徴構造としてアキシャル画像上の大腿骨内側顆と大腿骨外側顆とを結んだ線の垂線を、それぞれ、抽出する。スライス面決定部243は、第一の解剖学的特徴構造と第二の解剖学的特徴構造とにより特定される膝コロナル面を第三のスライス面と決定する。なお、推奨撮影断面算出部241は、第三のスライス面の画像として第三の画像を補間等により作成する。そして、解剖学的特徴構造抽出部242は、第三の画像上で大腿骨と脛骨の間接面の接線を第三の解剖学的特徴構造として特定する。
さらに、スライス面決定部243は、第三の解剖学的特徴構造と第四の解剖学的構造とにより特定される膝サジタル面を第五のスライス面と決定し、推奨撮影断面算出部241は、第五のスライス面の画像を第五の画像として補間等により作成する。解剖学的特徴構造抽出部242は、膝サジタル画像上の前十字靱帯に沿ったラインを第五の解剖学的特徴構造として抽出し、膝コロナル画像上で前十字靱帯に沿ったラインを第六の解剖学的特徴構造として抽出する。スライス面決定部243は、第五の解剖学的特徴構造に平行で、かつ、第六の解剖学的特徴構造に平行な平面を第六のスライス面と決定し、推奨撮影断面算出部241は、第六のスライス面を推奨撮影断面とする。
肩検査の場合、計測制御部230と信号処理部220とは、第一のスライス面をアキシャル面とし、第二のスライス面をサジタル面とし、第四のスライス面をコロナル面とし、これらに平行な1枚以上の二次元画像群を、それぞれ、第一の画像群、第二の画像群、第四の画像群として取得する。そして、解剖学的特徴構造抽出部242は、第一の解剖学的特徴構造としてアキシャル画像上の棘上筋に平行な線を、第二の解剖学的特徴構造としてサジタル画像上の上腕骨に沿ったラインを、第四の解剖学的特徴構造としてアキシャル画像上の関節面の接線に垂直な線を、それぞれ、抽出する。スライス面決定部243は、第一の解剖学的特徴構造と第二の解剖学的特徴構造とにより特定される肩コロナル面を第三のスライス面と決定する。なお、推奨撮影断面算出部241は、第三のスライス面の画像として第三の画像を補間等により作成する。そして、解剖学的特徴構造抽出部242は、第三の画像上で骨頭と肩甲骨とを結ぶラインに垂直なラインを第三の解剖学的特徴構造として特定する。
さらに、スライス面決定部243は、第三の解剖学的特徴構造と第四の解剖学的構造とにより特定される肩サジタル面を第五のスライス面と決定し、推奨撮影断面算出部241は、第五のスライス面の画像を第五の画像として補間等により作成する。解剖学的特徴構造抽出部242は、肩サジタル画像上の上腕骨に沿ったラインに垂直な直線を第五の解剖学的特徴構造として抽出し、肩コロナル画像上で肩峰と鎖骨とを結ぶ直線を第六の解剖学的特徴構造として抽出する。スライス面決定部243は、第五の解剖学的特徴構造に平行で、かつ、第六の解剖学的特徴構造に平行な平面を第六のスライス面と決定し、推奨撮影断面算出部241は、第六のスライス面を推奨撮影断面とする。
なお、本実施形態の頭部ルーチン検査では、頭部における組織間のコントラストが明瞭となるT1強調画像を取得するパルスシーケンスによりスカウト画像を取得する場合を例にあげて説明しているが、用いるパルスシーケンスはこれに限られない。例えば、T2強調画像を取得するパルスシーケンスを用いてもよい。ただし、この場合、画像の画素値のパターンが本実施形態とは異なるため、画像処理における評価関数を調節する必要がある。
また、本実施形態では、スカウト撮影と並行して推奨撮影断面算出部241による推奨撮影断面算出処理を行う場合を例にあげて説明しているが、これに限られない。例えば、全スカウト撮影終了後に推奨撮影断面算出処理を行うよう構成してもよい。
また、スライス面毎に、次のスライス面の撮影前に、当該スライス面の画像に対する推奨撮影断面算出処理を行い、その結果を次のスライス面の撮影に反映するよう構成してもよい。この場合、例えば、被検体103の体位を認識し、次の面の計測の位置を調整する機能を備える。頭部ルーチン検査の例では、コロナル面の計測後、得られた再構成画像から、上記範囲にアキシャル画像が収まるようにアキシャル撮影のz軸方向のスライス位置を調整し、アキシャル面の計測を行う。さらに、サジタル面の計測は、コロナル画像、アキシャル画像から計算された正中面の位置を含むようにスライス位置を調整後、行う。このように構成することにより、頭部が傾いてセッティングされた場合であっても、スカウト撮影を撮り直す必要がない。
さらに、スカウト撮影の前に、x軸、y軸、z軸それぞれのゼロ位相エンコードの核磁気共鳴信号を取得し、被検体103の大まかな位置および、大きさを特定後、スカウト撮影の位置を決定するよう構成してもよい。
また、本実施形態では、全スカウト撮影終了後、操作者に、処理を続行するか、または、スカウト撮影を再度実行するか選択させる場合を例にあげて説明した。しかし、これに限られない。推奨撮影断面算出部241は、選択を受け付ける構成を省略し、そのまま処理を進め、推奨撮影断面を算出し、例えば、頭部ルーチン検査であれば、正中面の解析結果が妥当でない場合は、推奨撮影断面とともに精度に信頼性がないことを伝えるメッセージを表示するよう構成してもよい。あるいは、正中面の解析結果が妥当である場合は処理を進め、妥当でない場合は、スカウト撮影を再度行うよう推奨するメッセージを表示装置111に表示し、処理を終えるよう構成してもよい。さらに、正中面の解析結果が妥当でない場合、推奨撮影断面算出部241が正中面を作成可能な位置を算出し、この位置で自動的にスカウト撮影を再度行うよう構成してもよい。さらに、正中面位置の解析時に、推奨撮影断面算出部241が正中面が大きく傾いていると判断した場合、この時点でスカウト撮影の再実行を推奨するメッセージを表示するよう構成してもよい。これらの構成によれば、操作者が画像を確認しなくとも、スカウト画像の善し悪しを判断可能であるため、検査効率が向上する。
また、本実施形態では、本撮影の前に操作者が撮影断面設定画面を開く操作を行い、表示される推奨撮影断面上で操作者からの調整の入力を受け付ける場合を例にあげて説明した。しかし、これに限られない。撮影断面決定部240は、推奨撮影断面が妥当である場合は、そのまま推奨撮影断面を撮影断面と決定し、本撮影を実行するよう構成してもよい。例えば、頭部ルーチン検査の場合、正中面が妥当と解析された場合、位置決め設定の入力を受け付ける画面を表示させることなく、本撮影開始の指示を待つ。妥当でないと解析された場合のみ、本撮影開始の指示受け付けの前に、位置決め設定用の入力を受け付ける画面を表示するか、または、上述のように精度が低いこと、あるいは、スカウト撮影再実行を推奨するメッセージを表示する。さらに、推奨撮影断面上で調整の入力を受け付ける機能の使用の有無を選択可能なよう構成してもよい。このように構成することで、算出した正中面が妥当である場合、操作者の処理量および検査全体の処理時間が低減する。
また、本実施形態では、推奨撮影断面の算出まで自動的に行う場合を例にあげて説明しているが、これに限られない。例えば、撮影断面の設定は、操作者が行うよう構成してもよい。この場合、例えば、頭部ルーチン検査であれば、推奨撮影断面算出部241は、正中面画像の作成まで行う。そして、操作者が本撮影前に位置決め設定用の画面を開く操作をした場合、撮影断面決定部240は、正中面画像を表示させ、操作者からの入力を受け付ける。操作者は、推奨撮影断面算出部241が算出した正中面画像上で、撮影断面を設定する。このとき、算出した正中面画像とともに、取得したスカウト画像(コロナル画像、アキシャル画像、サジタル画像)を表示装置111に表示するよう構成してもよい。このように構成することにより、操作者が手動で撮影断面を設定する場合の操作性が向上する。また、頭部が傾いてセッティングされていた場合に、従来は複数枚にまたがって取得された正中面画像が一枚の画像で確認可能となり、位置決め作業が容易にとなる。
また、本実施形態では、推奨撮影断面算出情報記憶部330に撮影断面リストを予め保持する場合を例にあげて説明しているが、これに限られない。操作者が新たな撮影断面リスト項目を追加可能な機能を有してもよい。例えば、操作者が表示装置111および入力装置116を介して、目印となる解剖学的特徴構造の位置が既知の画像上で撮影断面を設定する。そして、その関係を任意の名称で撮影断面リストに登録する。設定時には、過去に撮影したスカウト画像や現在の患者を撮影したスカウト画像、あるいは標準的なヒトのスカウト画像等を用いる。いずれの画像を用いる場合でも,解剖学的特徴構造を手動または画像処理にて自動で抽出する。使用するインタフェースとして、専用のインタフェースを持つよう構成してもよいし、撮影断面決定用のインタフェースを用いるよう構成してもよい。この場合、入力した情報を撮影断面リストに保存する機能を追加する。また、過去に行った検査で設定した撮影断面を撮影断面リストに保存するよう構成してもよい。
また、撮影断面リストは、検査毎に操作者が登録するよう構成してもよい。まず、本撮影のパルスシーケンスの撮影パラメータ設定時に、設定されたパラメータ(スライス枚数、スライス厚、FOV等)をもとに撮影断面を軸断で表示する。このとき、解剖学的特徴構造の位置が既知のスカウト画像を一緒に表示する。画像は過去の撮影等で取得したスカウト画像でもよいし、標準的なヒトのスカウト画像でもよい。次に、操作者は一緒に表示された画像の解剖学的特徴構造を目印に手動で撮影断面を設定する。このとき、計算機110は手動設定した位置を本撮影の撮影断面とし、自動認識する解剖学的特徴構造の位置との位置関係を保存する。以上の動作は、操作者が従来の検査感覚で撮影断面位置を設定できるため、事前準備の煩雑さを感じることがない。また、特別なユーザインタフェースが必要でない。さらに、検査の過程で撮影断面位置設定を含めたプロトコルの変更が可能となる。
<<第二の実施形態>>
次に、本発明を適用する第二の実施形態について説明する。本実施形態のMRI装置は基本的に第一の実施形態と同様の構成を有する。ただし、本実施形態では、推奨撮影断面または推奨スカウト撮影断面を表示させた後に操作者から受け付けた調整量を学習データとして記憶し、その後の処理に反映する機能を有する。以下、第一の実施形態と異なる構成に主眼をおいて説明する。
図14は、本実施形態の計算機110Aと記憶装置112Aとにより構成される情報処理装置の機能ブロック図である。本実施形態の情報処理装置は、第一の実施形態と同様の構成に加え、それぞれ、制御部200Aの撮影断面決定部240が学習機能部244を、記憶部300Bが学習データ記憶部340を備える。また、推奨撮影断面算出部241の代わりに推奨撮影断面算出部241Aを備える。
本実施形態においても、計算機110Aは、CPUとメモリとを備え、計算機110Aが実現する制御部200Aの各機能は、CPUが記憶装置112Aに格納されたプログラムをメモリにロードして実行することにより実現される。また、記憶部300Aは記憶装置112A上に実現される。各機能の全てまたはその一部は、MRI装置100とは独立して設けられる汎用の情報処理装置であって、MRI装置100とデータの送受信が可能な情報処理装置により実現されていてもよい。同様に、記憶部300Aの全てまたはその一部は、MRI装置100とは独立して設けられ、MRI装置100とデータの送受信が可能な外部記憶装置により実現されてもよい。
学習機能部244は、上記ステップS1225、S1243において、推奨撮影断面算出部241Aにより算出された推奨撮影断面および推奨スカウト撮影断面に操作者が加えた調整量を抽出し、学習データ記憶部340に学習データとして保持する。このとき、撮影対象部位に対応づけて学習データを登録する。
学習データ記憶部240に登録される学習データは、操作者が推奨撮影断面に対して加えた角度および位置の調整量である。学習機能部244は、調整がなされる毎に調整量を収集し、撮影対象部位に対応づけて学習データ記憶部340に登録されている学習データを最新のものに更新する。
なお、学習データとして、複数検査にわたって変更を加えた量(調整量)の平均値を登録するよう構成してもよい。また、複数検査にわたって収集された調整量の中の頻度の高い値のみで算出された平均値を登録するよう構成してもよい。頻度の高い値のみで算出された平均値を登録することにより、例外的に位置調整をした検査の調整量を学習しないようにすることができる。これにより、通常の検査における補正値としての精度を高めることができる。さらに、調整量に閾値を設け、閾値内の値のみで算出した平均値を登録するよう構成してもよい。また、被験者の年齢,性別,検査部位の大きさ,セッティングの向き,解剖学的構造の特徴量などで学習データをクラスタリングし、同じクラスタ内で調整値の平均値を算出し、登録するよう構成してもよい。このように構成することにより、例えば、個体差に応じた補正値を学習データとして得ることができる。
推奨撮影断面算出部241Aは、登録された学習データを次回の算出時の補正値とする。すなわち、上記検査準備処理において、ステップS1223、ステップS1241でそれぞれ推奨撮影断面および推奨スカウト撮影断面を算出後、学習データ記憶部340を参照し、撮影対象部位について学習データが登録されている場合、当該学習データを用いて算出結果を補正し、補正後の断面をそれぞれ、推奨撮影断面および推奨スカウト撮影断面とする。
以上説明したように、本実施形態によれば、第一の実施形態と同様の効果に加え、さらに高い精度で推奨撮影断面を得ることができる。例えば、施設毎、操作者毎に、設定する撮影断面が異なる場合、それぞれに合わせた撮影断面を自動的に出力できる。
なお、調整量を学習するか否かを操作者に選択させるインタフェースを設けるよう構成してもよい。学習機能部244は、操作者が調整を学習することを選択した場合のみ、上記処理を行う。さらに、学習データを推奨撮影断面および推奨スカウト撮影断面に反映させるか否かを操作者が選択可能なように構成してもよい。推奨撮影断面算出部241は、操作者が学習データを反映させることを選択した場合のみ、上記処理を行う。なお、選択は、計測毎、撮影毎、検査毎のいずれで可能なよう構成してもよい。
また、推奨撮影断面または推奨スカウト撮影断面からの調整量だけでなく、例えば、頭部ルーチン検査であれば、最終的なテンプレートモデルの形状やテンプレートモデルの初期配置、対象とする被検体なども学習データとして保存するよう構成してもよい。
さらに、パルスシーケンスおよび撮影パラメータが同じ本撮影を異なるプロトコルで規定される検査で行う場合、操作者の操作傾向である学習データを共有するか否かを選択可能なよう構成してもよい。このように構成することによって、異なるプロトコル間で同じ本撮影を行う場合、学習機能の適用を操作者の好みにあわせることができる。
また、学習機能部244は、推奨撮影断面位置の設定に用いることもできる。上記第一の実施形態で説明したように、初期設定として軸断で表示された撮影断面について操作者が加えた変更を記憶し、次回からの処理に反映させる。この場合、操作者は予め推奨撮影断面位置の設定もしなくとも良い利点がある。
<<第三の実施形態>>
次に、本発明を適用する第三の実施形態を説明する。本実施形態のMRI装置は、基本的に上記各実施形態のいずれかと同様の構成を有する。ただし、本実施形態では、第一の実施形態および第二の実施形態の推奨撮影断面算出部をMPR処理に用いる。以下、上記各実施形態のいずれかと異なる構成に主眼をおいて、本実施形態を説明する。
図15は、本実施形態の計算機110Bと記憶装置112Bとにより構成される情報処理装置の機能ブロック図である。本実施形態の情報処理装置は、基本的に上記各実施形態の情報処理装置のいずれかと同様であるが、上記各実施形態の情報処理装置の構成に加え、制御部200Bに、MPR処理部260を備える。なお、学習機能部244および学習データ340は備えなくてもよい。
MPR処理部260は、MPR処理用のインタフェース画面(MPRIF画面)を生成し、UI制御部210を介して表示装置111に表示させる。MPRIF画面には、撮影した3次元データの表示領域と、操作者からの入力を受け付ける指定受付領域とを備える。ここでは、操作者から、診断したい断面(切出断面)を特定する情報と、FOV、スライス枚数、スライス厚さ、スライス間隔などの切出断面を特定するパラメータ(断面パラメータ)と、の指定を受け付ける。
なお、切出断面を特定する情報は、例えば、頭部の場合、OMラインに平行な面、AP−ACラインに平行な面などの情報である。また、断面パラメータの指定は、例えば、リストから選択する、表示領域に表示された画像上で操作者が手動で設定する、などいずれの方式であってもよい。リストから選択する方式の場合、グラフィカルなインタフェースによる視覚的なサポートがあってもよい。また、過去に切り出し切出断面を表示し、単に同意の有無を受け付けるインタフェースを有するよう構成してもよい。この場合、同意しない意思を受け付けた場合のみ、断面パラメータの指定を受け付ける画面に移行する。本構成にすることで、過去に設定した切出断面と同じ断面を容易に指定できる。
MPR処理部260は、操作者からMPR処理開始の指示を受け付けると、まず、3次元データから、少なくとも2面の、互いに交差する面の画像を作成する。例えば、頭部の場合、コロナル面に平行な画像、アキシャル面に平行な画像、サジタル面に平行な画像を作成する。推奨撮影断面算出部241は、作成された各画像に対し、上記各実施系形態と同様の手順で同様の処理を行い、推奨撮影断面として、切出断面を算出する。そして、MPR処理部260は、算出された切出断面の画像を、3次元データから生成し、UI制御部210は、生成された画像を表示装置に表示する。
なお、本実施形態では、例えば、頭部の場合、正中面の位置が解剖学的に妥当でない位置と判別された場合、エラーメッセージを表示するよう構成する。この場合、例えば、切出断面を算出せずに、直交三軸断面を表示し、操作者に手動で切出断面を設定させる。
また、本実施形態においても、処理のアルゴリズム、画像処理の種類、切出断面と解剖学的特徴構造との位置関係などの情報は、予め推奨撮影断面算出情報記憶部330に格納される。
また、MPR処理部260が、各断面の画像を、最適スカウト撮影情報記憶部320に記憶されている、スカウト撮影と同様の順に行うよう構成し、必要な画像の生成が終わると、推奨撮影断面算出部241は、MPR処理部260による画像の生成処理と並行して推奨撮影断面作成処理を行うよう構成してもよい。
以上説明したように、本実施形態によれば、3次元データから、短時間で容易に所望の診断画像を生成することができる。
以上説明したように、上記各実施形態によれば、撮影断面を設定するにあたり、手動で設定する場合と同じ二次元スカウト撮影と同時に画像処理を実行し、推奨撮影断面を算出し、提示するため,従来の検査の流れを変更することなく、また、推奨撮影断面算出処理による時間の延長を発生させることなく、撮影断面の位置決めの操作性を向上させることができる。さらに、推奨撮影断面を算出する処理は、MPRにおいて切出断面を自動で算出する処理にも応用でき、後処理における検査効率も向上する。
なお、上記各実施形態では、MRI装置に本発明を適用する場合を例にあげて説明しているが、装置はこれに限られない。上記各実施形態にかかる推奨撮影断面算出部は、3次元空間の任意の面を撮影断面として設定可能な各種の医用画像撮影装置、画像撮影装置の推奨撮影断面算出に用いることができる。
100:MRI装置、101:マグネット、102:傾斜磁場コイル、103:被検体(生体)、104:シーケンサ、105:傾斜磁場電源、106:高周波磁場発生器、107:RFコイル、108:RFプローブ、109:受信器、110:計算機、111:表示装置、112:記憶装置、113:シムコイル、114:シム電源、115:寝台(テーブル)、116:入力装置、200:制御部、200A:制御部、200B:制御部、210:UI制御部、220:信号処理部、230:計測制御部、240:撮影断面決定部、241:推奨撮影断面算出部、241A:推奨撮影断面算出部、242:解剖学的特徴構造抽出部、243:スライス面決定部、244:学習機能部、260:MPR処理部、300:記憶部、300A:記憶部、301:コロナル画像、302:アキシャル画像、303:サジタル画像、304:座標軸、310:撮影情報記憶部、320:最適スカウト撮影情報記憶部、330:推奨撮影断面情報記憶部、340:学習データ記憶部、501:加算画像、502:頭部領域、503:バックグラウンド領域、504:プロット、505:一次関数の傾き、511:x軸への一次元投影像、512:頭部の幅、513:閾値、521:z軸への一次元投影像、522:逆数の微分値、531:頭頂点、532:点、533:範囲、534:スライス位置、601:階調調整画像、602:加算画像、603:重心、604:x軸への一次元投影像、605:頭部の幅、606:角度、607:1/2の幅、608:輝度最小点、609:プロット、610:傾きα、611:正中面通過点P、612:正中線、621:前壁側の領域、622:後壁側の領域、800:テンプレートモデル、801:交線、801:目印点、810:正中面画像、820:標準モデル、901:正中面画像、902:重心、903:z軸への一次元投影像、904:頭部の長さ、905:輪郭座標、906:一次関数、907:垂線との交点、908:輪郭座標、910:一次関数、911:一次関数の傾き、913:頭部の幅、920:初期位置、1002:OMライン、1003:全脳の範囲、1004:推奨撮影断面の傾き、1005:推奨撮影断面の範囲

Claims (21)

  1. 3次元空間の任意の面を撮影可能な医用画像撮影装置であって、
    互いに交差する2つのスライス面の内の一方のスライス面である第一のスライス面に平行な1枚以上の二次元画像からなる第一の画像群と、他方のスライス面である第二のスライス面に平行な1枚以上の二次元画像からなる第二の画像群とを取得する画像取得手段と、
    二次元画像から、予め定められた解剖学的特徴構造の情報を抽出する解剖学的特徴構造抽出手段と、
    撮影断面として推奨する推奨撮影断面を算出する推奨撮影断面算出手段と、
    撮影対象部位に応じて、前記推奨撮影断面を算出するために必要な推奨撮影断面算出情報を保持する推奨撮影断面算出情報記憶手段と、を備え、
    前記解剖学的特徴構造抽出手段は、前記第一の画像群から第一の解剖学的特徴構造の情報を抽出するとともに、前記第二の画像群から第二の解剖学的特徴構造の情報を抽出し、
    前記推奨撮影断面算出手段は、前記推奨撮影断面算出情報に従い、前記第一の解剖学的特徴構造の情報および前記第二の解剖学的特徴構造の情報を用いて前記推奨撮影断面を算出すること
    を特徴とする医用画像撮影装置。
  2. 請求項1記載の医用画像撮影装置であって、
    前記推奨撮影断面算出手段は、
    スライス面を決定するスライス面決定手段と、
    再構成画像から前記スライス面決定手段が決定したスライス面の二次元画像を生成する画像生成手段と、
    画像上の解剖学的な特徴点の位置を特定する特徴点決定手段と、をさらに備え、
    前記スライス面決定手段は、前記第一の解剖学的特徴構造および前記第二の解剖学的特徴構造を用い、前記推奨撮影断面算出情報に従って第三のスライス面を決定し、
    前記画像取得手段は、前記第一のスライス面および前記第二のスライス面の双方に交差する第四のスライス面に平行な1枚以上の二次元画像からなる第四の画像群をさらに取得し、
    前記画像生成手段は、前記第四の画像群から、前記第三のスライス面の二次元画像を第三の画像として生成し、
    前記特徴点決定手段は、前記推奨撮影断面算出情報に従い、前記第三の画像上で前記推奨撮影断面が含むべき複数の解剖学的な特徴点の位置を特定し、
    前記推奨撮影断面算出手段は、第三のスライス面に交差し、前記複数の解剖学的特徴点を含む面を前記推奨撮影断面とすること
    を特徴とする医用画像撮影装置。
  3. 請求項1記載の医用画像撮影装置であって
    前記推奨撮影断面算出手段は、スライス面を決定するスライス面決定手段をさらに備え、
    前記スライス面決定手段は、前記第一の解剖学的特徴構造および前記第二の解剖学的特徴構造を用い、前記推奨撮影断面算出情報に従って第三のスライス面を決定し、
    前記推奨撮影断面算出手段は、前記第三のスライス面を前記推奨撮影断面とすること
    を特徴とする医用画像撮影装置。
  4. 請求項3記載の医用画像撮影装置であって、
    前記推奨撮影断面算出手段は、再構成画像から前記スライス面決定手段が決定したスライス面の二次元画像を生成する画像生成手段をさらに備え、
    前記画像取得手段は、前記第一のスライス面および前記第二のスライス面の双方に交差する第四のスライス面に平行な1枚以上の二次元画像からなる第四の画像群をさらに取得し、
    前記画像生成手段は、前記第四の画像群から、前記第三のスライス面の二次元画像を第三の画像として生成し、
    前記解剖学的特徴構造抽出手段は、前記推奨撮影断面算出情報に従い、前記第三の画像から第三の解剖学的特徴構造の情報を抽出し、
    前記スライス面決定手段は、前記第一の解剖学的特徴構造の情報および前記第二の解剖学的特徴構造の情報のいずれか一方と、前記第三の解剖学的特徴構造の情報とを用い、前記推奨撮影断面算出情報に従って、第五のスライス面を決定し、
    前記推奨撮影断面算出手段は、前記第五のスライス面を第二の前記推奨撮影断面とすること
    を特徴とする医用画像撮影装置。
  5. 請求項3記載の医用画像撮影装置であって、
    前記推奨撮影断面算出手段は、再構成画像から前記スライス面決定手段が決定したスライス面の二次元画像を生成する画像生成手段をさらに備え、
    前記画像取得手段は、前記第一のスライス面および前記第二のスライス面の双方に交差する第四のスライス面に平行な1枚以上の二次元画像からなる第四の画像群をさらに取得し、
    前記画像生成手段は、前記第四の画像群から、前記第三のスライス面の二次元画像を第三の画像として生成し、
    前記解剖学的特徴構造抽出手段は、前記推奨撮影断面算出情報に従い、前記第三の画像から第三の解剖学的特徴構造の情報を抽出するとともに、前記第一の画像群から第四の解剖学的特徴構造の情報を抽出し、
    前記スライス面決定手段は、前記第三の解剖学的特徴構造の情報および前記第四の解剖学的特徴構造の情報を用い、前記推奨撮影断面算出情報に従って、第五のスライス面を決定し、
    前記推奨撮影断面算出手段は、前記第五のスライス面を第二の前記推奨撮影断面とすること
    を特徴とする医用画像撮影装置。
  6. 請求項5記載の医用画像撮影装置であって、
    前記画像生成手段は、前記第二の画像群から前記第五のスライス面の二次元画像を第五の画像として生成し、
    前記解剖学的特徴構造抽出手段は、前記推奨撮影断面算出情報に従い、前記第五の画像から第五の解剖学的特徴構造の情報を抽出するとともに、前記第三の画像から第六の解剖学的特徴構造の情報を抽出し、
    前記スライス面決定手段は、前記第五の解剖学的特徴構造の情報および前記第六の解剖学的特徴構造の情報を用い、前記推奨撮影断面算出情報に従って、第六のスライス面を決定し、
    前記推奨撮影断面算出手段は、前記第六のスライス面を第三の前記推奨撮影断面とすること
    を特徴とする医用画像撮影装置。
  7. 請求項1から6いずれか1項記載の医用画像撮影装置であって、
    前記推奨撮影断面算出情報記憶手段は、前記解剖学的特徴構造の情報の抽出に用いる画像処理の種別を記憶し、
    前記解剖学的特徴構造抽出手段は、前記推奨撮影断面算出情報記憶手段に記憶される画像処理により、前記解剖学的特徴構造の情報を抽出すること
    を特徴とする医用画像撮影装置。
  8. 請求項2記載の医用画像撮影装置であって、
    前記推奨撮影断面算出情報記憶手段は、前記撮影断面が含むべき複数の解剖学的な特徴点に対応する複数の座標点で構成されるテンプレートモデルを記憶し、
    前記特徴点決定手段は、前記第三の画像上で、前記テンプレートモデルをフィッティングすることにより、前記複数の解剖学的特徴点の位置を特定すること
    を特徴とする医用画像撮影装置。
  9. 請求項1から8いずれか1項記載の医用画像撮影装置であって、
    前記推奨撮影断面算出手段が算出した前記推奨撮影断面を表示する表示手段と、
    前記表示手段に表示された推奨撮影断面を変更するための調整量の入力を受け付ける入力手段と、
    前記入力手段を介して調整量を受け付けた場合、前記推奨撮影断面を、前記受け付けた調整量で調整し、調整後の断面を撮影断面と決定する撮影断面決定手段と、をさらに備えること
    を特徴とする医用画像撮影装置。
  10. 請求項9記載の医用画像撮影装置であって、
    撮影対象部位に対応付けて前記入力手段を介して受け付けた調整量を記憶する学習データ記憶手段をさらに備え、
    前記推奨撮影断面算出手段は、
    前記学習データ記憶手段に撮影を行う対象の撮影部位に対応付けて前記調整量が記憶されている場合、前記表示手段に表示する前に、前記算出した推奨撮影断面を前記調整量で調整すること
    を特徴とする医用画像撮影装置。
  11. 請求項1から10いずれか1項記載の医用画像撮影装置であって、
    予め定められたシーケンスに従って、前記画像取得手段による前記二次元画像の取得を制御するスカウト撮影制御手段をさらに備え、
    前記スカウト撮影制御手段は、前記第一の画像群の取得後に前記第二の画像群を取得し、
    前記解剖学的特徴構造抽出手段は、前記第一の解剖学的特徴構造の情報の抽出を、前記第二の画像群の取得の前に終えること
    を特徴とする医用画像撮影装置。
  12. 請求項1から10いずれか1項記載の医用画像撮影装置であって、
    予め定められたシーケンスに従って、前記画像取得手段による前記二次元画像の取得を制御するスカウト撮影制御手段をさらに備え、
    前記スカウト撮影制御手段は、前記第一の画像群の取得後に前記第二の画像群を取得し、
    前記解剖学的特徴構造抽出手段は、前記第一の解剖学的特徴構造の情報の抽出を、前記第二の画像群を取得するシーケンスの開始前に終え、
    前記スカウト撮影制御手段は、前記第一の解剖学的特徴構造の情報に基づき、第二の画像群を取得する位置を調整すること
    を特徴とする医用画像撮影装置。
  13. 請求項3項記載の医用画像撮影装置であって、
    予め定められたシーケンスに従って、前記画像取得手段による前記二次元画像の取得を制御するスカウト撮影制御手段と、
    前記推奨撮影断面算出手段が算出した推奨撮影断面から撮影断面を決定する撮影断面決定手段と、
    前記撮影断面決定手段が決定した撮影断面の診断用画像を取得する本撮影を行う本撮影制御手段と、をさらに備え、
    前記解剖学的特徴構造抽出手段は、前記推奨撮影断面算出情報に従い、前記診断用画像から第三の解剖学的特徴構造の情報を抽出し、
    前記スライス面決定手段は、前記第一の解剖学的特徴構造の情報および前記第二の解剖学的特徴構造の情報のいずれか一方と、前記第三の解剖学的特徴構造の情報とを用い、前記推奨撮影断面算出情報に従って、第五のスライス面を決定し、
    前記推奨撮影断面算出手段は、前記第五のスライス面を第二の前記推奨撮影断面と決定すること
    を特徴とする医用画像撮影装置。
  14. 請求項11から13いずれか1項記載の医用画像撮影装置であって、
    撮影対象部位に対応づけて、前記第一のスライス面、前記第二のスライス面をそれぞれ特定する情報を保持するスカウト撮影情報保持手段をさらに備え、
    前記スカウト撮影制御手段は、撮影対象部位に応じて前記スカウト撮影情報保持手段に保持されている情報に従って、前記画像取得手段を制御すること
    を特徴とする医用画像撮影装置。
  15. 請求項11から14いずれか1項記載の医用画像撮影装置であって、
    前記撮影断面決定手段が決定した撮影断面の適否を判別する適否判別手段と、
    前記適否判別手段で否と判別された場合、前記推奨撮影断面に基づき、前記画像取得手段が二次元画像を取得するスライス面を決定し、当該画像取得手段に当該スライス面の画像を取得させる再スカウト撮影指示手段と、をさらに備えること
    を特徴とする医用画像撮影装置。
  16. 請求項1から15記載の医用画像撮影装置であって、
    当該医用画像撮影装置は、静磁場におかれた被検体に高周波磁場を印加して核磁化を励起し、計測した核磁気共鳴信号を画像化する磁気共鳴イメージング装置であること
    を特徴とする医用画像撮影装置。
  17. 請求項16記載の医用画像撮影装置であって、
    3次元ボリュームデータを取得する3次元撮影制御手段をさらに備え、
    前記第一の画像群および前記第二の画像群は、前記3次元撮影制御手段で取得した3次元ボリュームデータから、切り出した断面の画像であること
    を特徴とする医用画像撮影装置。
  18. 3次元空間の任意の面を撮影可能な医用画像撮影装置において、当該医用画像撮影装置が取得した画像であって、互いに交差する2つのスライス面の内の一方のスライス面である第一のスライス面に平行な1枚以上の二次元画像からなる第一の画像群と、他方のスライス面である第二のスライス面に平行な1枚以上の二次元画像からなる第二の画像群とから、撮影断面として推奨する推奨撮影断面を算出する画像処理装置であって、
    撮影対象部位に応じて、前記推奨撮影断面を算出するために必要な推奨撮影断面算出情報を保持する推奨撮影断面算出情報記憶手段と
    前記推奨撮影断面情報に従って、前記第一の画像群から第一の解剖学的特徴構造の情報を抽出するとともに、前記第二の画像群から第二の解剖学的特徴構造の情報を抽出する解剖学的特徴構造抽出手段と、
    前記第一の解剖学的特徴構造の情報と前記第二の解剖学的特徴構造の情報とを用い、前記推奨撮影断面情報に従って、前記推奨撮影断面を算出する推奨撮影断面算出手段と、を備えること
    を特徴とする画像処理装置。
  19. 3次元空間の任意の面を撮影可能な医用画像撮影装置が備えるコンピュータに、
    互いに交差する2つのスライス面の内の一方のスライス面である第一のスライス面に平行な1枚以上の二次元画像からなる第一の画像群を再構成するとともに、前記第一の画像群の再構成後に他方のスライス面である第二のスライス面に平行な1枚以上の二次元画像からなる第二の画像群を再構成する画像再構成ステップと、
    前記画像再構成ステップにおいて、前記第一の画像群の再構成を終えると、当該第一の画像群から第一の解剖学的特徴構造の情報を抽出し、前記第二の画像群の再構成を終えると、当該第二の画像群から第二の解剖学的特徴構造の情報を抽出する解剖学的特徴構造抽出ステップと、
    前記第一の解剖学的特徴構造の情報と前記第二の解剖学的特徴構造の情報とから、前記推奨撮影断面を算出する推奨撮影断面算出ステップと、を実行させるためのプログラム。
  20. 3次元空間の任意の面を撮影可能な医用画像撮影装置において、撮影断面として推奨する推奨撮影断面算出方法であって、
    互いに交差する2つのスライス面の内の一方のスライス面である第一のスライス面に平行な1枚以上の二次元画像からなる第一の画像群を再構成するとともに、前記第一の画像群の再構成後に他方のスライス面ある第二のスライス面に平行な1枚以上の二次元画像からなる第二の画像群を再構成する画像再構成ステップと、
    前記画像再構成ステップにおいて、前記第一の画像群の再構成を終えると、当該第一の画像群から第一の解剖学的特徴構造の情報を抽出し、前記第二の画像群の再構成を終えると、当該第二の画像群から第二の解剖学的特徴構造の情報を抽出する解剖学的特徴構造抽出ステップと、
    前記第一の解剖学的特徴構造の情報と前記第二の解剖学的特徴構造の情報とから、前記推奨撮影断面を算出する推奨撮影断面算出ステップと、を備えること
    を特徴とする推奨撮影断面算出方法。
  21. 請求項20記載の推奨撮影断面算出方法であって、
    前記画像再構成ステップは、前記第二の画像群の再構成後に前記第一のスライス面および前記第二のスライス面の双方に交差する第四のスライス面に平行な1枚以上の二次元画像からなる第四の画像群を再構成するステップを備え、
    前記推奨撮影断面算出ステップは、
    前記第一の解剖学的特徴構造の情報と前記第二の解剖学的特徴構造の情報とから第三のスライス面を決定する第三のスライス決定ステップと、
    前記第三のスライス面の二次元画像を、前記第四の画像群から第五の画像として生成する画像生成ステップと、
    前記第五の画像上で前記撮影断面が含むべき複数の解剖学的な特徴点の位置を特定する特徴点特定ステップと、
    前記第三のスライス面に交差し、前記複数の解剖学的特徴点を含む面を前記撮影断面として決定する決定ステップと、を備えること
    を特徴とする推奨撮影断面算出方法。
JP2011519903A 2009-06-25 2010-06-22 医用画像撮影装置 Active JP5230809B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011519903A JP5230809B2 (ja) 2009-06-25 2010-06-22 医用画像撮影装置

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2009151338 2009-06-25
JP2009151338 2009-06-25
JP2011519903A JP5230809B2 (ja) 2009-06-25 2010-06-22 医用画像撮影装置
PCT/JP2010/060557 WO2010150783A1 (ja) 2009-06-25 2010-06-22 医用画像撮影装置

Publications (2)

Publication Number Publication Date
JPWO2010150783A1 JPWO2010150783A1 (ja) 2012-12-10
JP5230809B2 true JP5230809B2 (ja) 2013-07-10

Family

ID=43386552

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011519903A Active JP5230809B2 (ja) 2009-06-25 2010-06-22 医用画像撮影装置

Country Status (4)

Country Link
US (1) US8693760B2 (ja)
JP (1) JP5230809B2 (ja)
CN (1) CN102802519B (ja)
WO (1) WO2010150783A1 (ja)

Families Citing this family (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5415772B2 (ja) * 2009-01-07 2014-02-12 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 正中面決定装置および磁気共鳴イメージング装置
JP5575454B2 (ja) * 2009-10-29 2014-08-20 株式会社東芝 磁気共鳴イメージング装置
JP2011212148A (ja) * 2010-03-31 2011-10-27 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置、スライス位置設定方法、およびプログラム
JP5925998B2 (ja) * 2011-04-28 2016-05-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置およびプログラム
WO2012165485A1 (ja) 2011-05-31 2012-12-06 株式会社東芝 磁気共鳴イメージング装置
CN103764024B (zh) * 2011-08-25 2016-01-06 株式会社日立医疗器械 医用图像摄影装置
JP5766070B2 (ja) * 2011-08-31 2015-08-19 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 画像データ処理装置、磁気共鳴装置、画像データ処理方法、およびプログラム
JP6054024B2 (ja) * 2011-11-29 2016-12-27 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー スライス位置設定装置、磁気共鳴装置、スライス位置設定方法、およびプログラム
US11286463B2 (en) 2012-03-08 2022-03-29 Advanced ReGen Medical Technologies, LLC Reprogramming of aged adult stem cells
JP5806175B2 (ja) * 2012-07-06 2015-11-10 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー ラベリング領域決定装置および磁気共鳴装置並びにプログラム
JP2014030621A (ja) 2012-08-03 2014-02-20 Sony Corp 情報処理装置、プログラム及び生体信号測定セット
US9076246B2 (en) * 2012-08-09 2015-07-07 Hologic, Inc. System and method of overlaying images of different modalities
US9646376B2 (en) 2013-03-15 2017-05-09 Hologic, Inc. System and method for reviewing and analyzing cytological specimens
RU2523929C1 (ru) * 2013-03-18 2014-07-27 Корпорация "САМСУНГ ЭЛЕКТРОНИКС Ко., Лтд." Система и способ для автоматического планирования видов в объемных изображениях мозга
CN104146766B (zh) * 2013-05-13 2016-09-28 东芝医疗系统株式会社 扫描装置、扫描方法和医学图像设备
JP6404029B2 (ja) 2013-10-01 2018-10-10 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置及び撮像計画プログラム
JP6041781B2 (ja) * 2013-10-11 2016-12-14 富士フイルム株式会社 医用画像処理装置およびその作動方法並びに医用画像処理プログラム
WO2015148534A1 (en) * 2014-03-24 2015-10-01 Advanced ReGen Medical Technologies, LLC Cell free compositions for cellular restoration and methods of making and using same
JP6353073B2 (ja) 2013-12-20 2018-07-04 アドヴァンスド リジェン メディカル テクノロジーズ,エルエルシー 細胞回復のための組成物並びにその作製及び使用方法
JP6482842B2 (ja) 2013-12-27 2019-03-13 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
JP6358816B2 (ja) * 2014-03-04 2018-07-18 株式会社日立製作所 医用画像診断装置およびそれに用いる画像処理方法
US10061979B2 (en) 2014-09-28 2018-08-28 Koninklijke Philips N.V. Image processing apparatus and method
CN105701438B (zh) * 2014-11-26 2020-06-23 东芝医疗系统株式会社 医学图像处理装置和方法以及医学成像设备
WO2016203295A1 (en) * 2015-06-19 2016-12-22 Synaptive Medical (Barbados) Inc. A medical imaging system for determining a scan orientation
JP6975131B2 (ja) * 2015-07-23 2021-12-01 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ボリュームイメージの対話型平面スライスからの内視鏡誘導
JP6999575B2 (ja) 2016-04-29 2022-01-18 アドヴァンスド リジェン メディカル テクノロジーズ,エルエルシー マイクロrna組成物並びにその作製及び使用方法
JP6765894B2 (ja) * 2016-08-02 2020-10-07 キヤノンメディカルシステムズ株式会社 医用画像処理装置及び医用画像診断装置
CN106725570B (zh) * 2016-12-30 2019-12-20 上海联影医疗科技有限公司 成像方法及系统
DE102017202601A1 (de) * 2017-02-17 2018-08-23 Siemens Healthcare Gmbh Bestimmung von Markierungen auf quantitativen Bilddaten
US11903691B2 (en) * 2017-07-18 2024-02-20 Siemens Healthineers Ag Combined steering engine and landmarking engine for elbow auto align
CN110728274A (zh) * 2018-06-29 2020-01-24 通用电气公司 医疗设备计算机辅助扫描方法、医疗设备及可读存储介质
JP7211725B2 (ja) 2018-07-09 2023-01-24 富士フイルムヘルスケア株式会社 自動位置決め機能付き磁気共鳴イメージング装置、磁気共鳴イメージング装置の制御方法、及び、撮像断面自動設定用プログラム
CN110855963B (zh) * 2018-08-21 2022-03-08 视联动力信息技术股份有限公司 一种视频数据的投影方法和装置
WO2020262683A1 (ja) * 2019-06-28 2020-12-30 富士フイルム株式会社 医用画像処理装置、方法およびプログラム
JP7320421B2 (ja) * 2019-10-08 2023-08-03 富士フイルムヘルスケア株式会社 磁気共鳴イメージング装置、及び、自動撮像位置設定方法
CN112716476A (zh) * 2019-10-28 2021-04-30 通用电气精准医疗有限责任公司 磁共振成像方法及系统、计算机可读存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01155837A (ja) * 1987-12-14 1989-06-19 Mitsubishi Electric Corp Mri断層位置指定方法
JPH05212013A (ja) * 1992-02-07 1993-08-24 Hitachi Ltd 部位の決定方法および正中線の決定方法
JPH06189935A (ja) * 1992-09-16 1994-07-12 Toshiba Medical Eng Co Ltd 位置決めスキャン方法及び磁気共鳴イメージング装置
JPH0966046A (ja) * 1995-08-31 1997-03-11 Shimadzu Corp 医用画像診断装置
JP2007319689A (ja) * 2006-06-02 2007-12-13 Siemens Ag 心臓の表示方法および磁気共鳴装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0775627A (ja) * 1993-06-11 1995-03-20 Hitachi Ltd 磁気共鳴診断装置における体動追従計測方法
JP2003245360A (ja) * 2002-02-26 2003-09-02 Piolax Medical Device:Kk ステント設計支援装置、ステント設計支援方法、ステント設計支援プログラム、及びステント設計支援プログラムを記録した記録媒体
JP2004283483A (ja) * 2003-03-25 2004-10-14 Hitachi Medical Corp 磁気共鳴イメージング装置
JP2006000127A (ja) * 2004-06-15 2006-01-05 Fuji Photo Film Co Ltd 画像処理方法および装置並びにプログラム
JP4512824B2 (ja) * 2005-02-04 2010-07-28 国立大学法人 千葉大学 画像処理装置及びそれに用いられるプログラム
JP2006325747A (ja) * 2005-05-24 2006-12-07 Toshiba Corp 画像処理装置、画像処理システム、x線ct装置、及び画像処理プログラム
US7590440B2 (en) * 2005-11-14 2009-09-15 General Electric Company System and method for anatomy labeling on a PACS
JP4863893B2 (ja) * 2007-01-26 2012-01-25 株式会社日立メディコ 磁気共鳴イメージング装置
US8543186B2 (en) * 2007-03-01 2013-09-24 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP5127305B2 (ja) * 2007-05-31 2013-01-23 株式会社東芝 磁気共鳴イメージング装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01155837A (ja) * 1987-12-14 1989-06-19 Mitsubishi Electric Corp Mri断層位置指定方法
JPH05212013A (ja) * 1992-02-07 1993-08-24 Hitachi Ltd 部位の決定方法および正中線の決定方法
JPH06189935A (ja) * 1992-09-16 1994-07-12 Toshiba Medical Eng Co Ltd 位置決めスキャン方法及び磁気共鳴イメージング装置
JPH0966046A (ja) * 1995-08-31 1997-03-11 Shimadzu Corp 医用画像診断装置
JP2007319689A (ja) * 2006-06-02 2007-12-13 Siemens Ag 心臓の表示方法および磁気共鳴装置

Also Published As

Publication number Publication date
US8693760B2 (en) 2014-04-08
WO2010150783A1 (ja) 2010-12-29
CN102802519A (zh) 2012-11-28
CN102802519B (zh) 2014-12-31
US20120093385A1 (en) 2012-04-19
JPWO2010150783A1 (ja) 2012-12-10

Similar Documents

Publication Publication Date Title
JP5230809B2 (ja) 医用画像撮影装置
JP5736459B2 (ja) 医用画像撮影装置
JP5667066B2 (ja) Mrスキャンの自動シーケンシャルプランニング
JP5788318B2 (ja) 磁気共鳴イメージング装置および撮像スライス決定方法
JP5178521B2 (ja) 磁気共鳴イメージング(mri)データを取得するシステム及び方法
US9746535B2 (en) Magnetic resonance imaging apparatus
KR101659142B1 (ko) 영상 데이터를 생성하는 방법
KR101652387B1 (ko) 영상 데이터를 생성하는 방법
US20110228998A1 (en) System and method for automatic computation of mr imaging scan parameters
US9070212B2 (en) Medical imaging apparatus and imaging slice determination method
US20070038070A1 (en) Determining patient-related information on the position and orientation of MR images by the individualisation of a body model
JP2018198682A (ja) 磁気共鳴イメージング装置及び磁気共鳴画像処理方法
JP2020039507A (ja) 磁気共鳴撮像装置、画像処理装置、及び、画像処理方法
US11810227B2 (en) MR image reconstruction based on a-priori information
JP2007164592A (ja) モデリング装置、領域抽出装置およびプログラム
JP6188764B2 (ja) 磁気共鳴イメージング装置
JP7320421B2 (ja) 磁気共鳴イメージング装置、及び、自動撮像位置設定方法
EP4339879A1 (en) Anatomy masking for mri
JP2022164567A (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
Lee Patient-Specific Rib Cage and Thoracic Volume Reconstruction in Spine Patients Using Biplanar Radiographs

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20121127

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130123

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130319

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160329

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 5230809

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250