(57)【特許請求の範囲】
【請求項１】
水平磁場型の磁場発生装置に目標の磁場分布が与えられた領域があり、その領域の磁場分布の誤差磁場成分を低減し、目標の磁場分布に近づける磁場調整法において、
調整手段として、受動的に磁化する鉄片などの磁性体又は外部磁場に依存しない永久磁石を配置する曲面の磁場調整機構を持ち、
所定数の点において磁場計測を行い、目標磁場との差である誤差磁場を算出し、その誤差を近似的に補正できるリング状の磁場調整機構面上の電流ポテンシャル分布を求め、その電流ポテンシャル分布を磁気モーメントに換算し、その磁気モーメントに相当する磁性体片を配置する磁場調整作業を特徴とする磁場調整法。
【請求項２】
請求項１の磁場調整法において、特異値分解により得た基底である固有分布関数の中から分布関数を選択し、その組み合わせで近似的に誤差磁場を補正する電流ポテンシャルの分布を表現することを特徴とする磁場調整法。
【請求項３】
請求項２の磁場調整法で、近似的に補正する電流ポテンシャルから目標磁場を与えた領域のおよび磁場計測点の補正磁場量を計算し、磁場調整作業後の目標磁場からの残留誤差磁場を求め、固有分布関数の選択の妥当性を確認すると共に、目標の残留誤差磁場以内にとなる選択を行うことを特徴とする磁場調整法。
【請求項４】
請求項１から３のいずれか１項に記載の磁場調整法において、電流ポテンシャルを磁気モーメントに比例する量として、鉄片量密度に換算し、その換算した分布に従って鉄片を配置することを特徴とする磁場調整法。
【請求項５】
請求項１から４のいずれか１項に記載の磁場調整法において、誤差磁場の補正に必要な電流ポテンシャル分布を求める固有分布関数の選択を、特異値の大きさの順に並べて番号付けした番号（次数）と、誤差磁場に含まれる固有分布の強さの相関図（スペクトル図）上で選択することを特徴とする磁場調整法。
【請求項６】
請求項１から５のいずれか１項に記載の磁場調整法において、電流ポテンシャル、または、磁気モーメントもしくは鉄片や永久磁石の量について、密度分布を鉄片を配置する磁場調整機構面上に等高線を含む表示を行い、その表示に従って鉄片を作業者が配置することを特徴とする磁場調整法。
【請求項７】
請求項６の磁場調整法において、等高線とともに鉄片を配置する面を多角形で分割し、分割した領域毎に、磁気モーメントの大きさもしくは鉄片量や永久磁石量を、面積積分値で、等高線と共に、もしくは等高線無しで示すことを特徴とする磁場調整法。
【請求項８】
請求項６の磁場調整法において、等高線で示した山もしくは盆地部をまとめて積算し、
その量を山もしくは盆地部内の一カ所もしくは複数箇所に分散して配置することを特徴とする磁場調整法。
【請求項９】
請求項１から８のいずれか１項に記載の磁場調整方法で、磁場計測から磁気モーメントの大きさもしくは鉄片量や永久磁石量の配置までの計算と作業を繰り返し実行することを特徴とする磁場調整法。
【請求項１０】
請求項９の磁場調整法で、繰り返し計算と作業で、誤差磁場の大きさと共に、特異値分解で得た磁場分布を表す基底である固有分布関数について、個々の強度の大きさを調べて、磁場調整の進展を把握することを特徴とする磁場調整法。
【請求項１１】
請求項３の残留誤差磁場の代表値、たとえば最小最大値の差を目標もしくは計測磁場の平均磁場強度で割った値を、請求項５の相関図、または請求項６，７，８の磁場補正作業で配置する磁性体又は永久磁石の大きさ表示と共に、表示することを特徴とする磁場調整法。
【請求項１２】
磁場発生用のコイルや磁性体の起磁力源を含む電磁石の請求項３の磁場調整法で、目標の誤差磁場以下とする条件での、請求項６，７の磁気モーメントの大きさもしくは鉄片量や永久磁石量の配置が可能であることで、磁場調整が正常に実行できることを、磁場調整作業開始時に判断できることを特徴とする磁石品質把握法。
【請求項１３】
磁場発生用のコイルや磁性体の起磁力源を含む電磁石の設計において、
水平磁場型の磁場発生装置に目標の磁場分布が与えられた領域があり、その領域の磁場分布の誤差磁場成分を低減し、目標の磁場分布に近づける起磁力源配置調整法において、
調整手段として、受動的に磁化する鉄片などの磁性体又は外部磁場に依存しない永久磁石を配置する曲面の磁場調整機構を持ち、
仮定して配置した起磁力源での、所定数の点における磁場分布を計算し、目標磁場との差である誤差磁場を算出し、その誤差を近似的に補正できるリング状の磁場調整機構面上の電流ポテンシャル分布を求め、その電流ポテンシャル分布を磁気モーメントに換算し、その磁気モーメントに相当するように前記起磁力源の配置を補正することを特徴とする起磁力源配置調整法。
【請求項１４】
水平磁場型の磁場発生装置に目標の磁場分布が与えられた領域があり、その領域の磁場分布の誤差磁場成分を低減し、目標の磁場分布に近づける磁場調整のための表示装置において、
前記表示装置は、
受動的に磁化する鉄片などの磁性体又は外部磁場に依存しない永久磁石を配置する曲面の磁場調整機構を持ち、
所定数の点において計測された磁場と目標磁場との差である誤差磁場を算出し、その誤差を近似的に補正できるリング状の磁場調整機構面上の電流ポテンシャル分布を求め、その電流ポテンシャル分布を磁気モーメントに換算し、その磁気モーメントに相当する磁性体片の配置位置を表示することを特徴とする表示装置。
【請求項１５】
前記表示装置は、
特異値分解により得た基底である固有分布関数の中から分布関数を選択し、その組み合わせで近似的に誤差磁場を補正する電流ポテンシャルの分布を表現することを特徴とする請求項１４に記載の表示装置。
【請求項１６】
前記表示装置において、
近似的に補正する電流ポテンシャルから目標磁場を与えた領域のおよび磁場計測点の補正磁場量を計算し、磁場調整作業後の目標磁場からの残留誤差磁場を求め、固有分布関数の選択の妥当性を確認すると共に、目標の残留誤差磁場以内にとなる選択を行うことを特徴とする請求項１５に記載の表示装置。
【請求項１７】
前記表示装置は、
電流ポテンシャルを磁気モーメントに比例する量として、鉄片量密度に換算し、その換算した分布に従って鉄片の配置を表示することを特徴とする請求項１４から請求項１６のいずれか１項に記載の表示装置。
【請求項１８】
前記表示装置は、
誤差磁場の補正に必要な電流ポテンシャル分布を求める固有分布関数の選択を、特異値の大きさの順に並べて番号付けした番号（次数）と、誤差磁場に含まれる固有分布の強さの相関図（スペクトル図）上で選択することを特徴とする請求項１４から請求項１７のいずれか１項に記載の表示装置。
【請求項１９】
前記表示装置は、
電流ポテンシャル、または、磁気モーメントもしくは鉄片や永久磁石の量について、密度分布を鉄片を配置する磁場調整機構面上に等高線を含む表示をすることを特徴とする請求項１４から請求項１８のいずれか１項に記載の表示装置。
【請求項２０】
前記表示装置は、
等高線とともに鉄片を配置する面を多角形で分割し、分割した領域毎に、磁気モーメントの大きさもしくは鉄片量や永久磁石量を、面積積分値で、等高線と共に、もしくは等高線無しで表示することを特徴とする請求項１９に記載の表示装置。
【請求項２１】
前記表示装置は、
等高線で示した山もしくは盆地部をまとめて積算し、その量を山もしくは盆地部内の一カ所もしくは複数箇所に分散して配置することを特徴とする請求項１９に記載の表示装置。
【請求項２２】
前記表示装置は、
磁場計測から磁気モーメントの大きさもしくは鉄片量や永久磁石量の配置までの計算と
作業を繰り返し実行することを特徴とする請求項１４から請求項２１のいずれか１項に記
載の表示装置。
【請求項２３】
前記表示装置は、
繰り返し計算と作業で、誤差磁場の大きさと共に、特異値分解で得た磁場分布を表す基底である固有分布関数について、個々の強度の大きさを調べて、磁場調整の進展を把握することを特徴とする請求項２２に記載の表示装置。
【請求項２４】
前記表示装置は、
請求項１６の残留誤差磁場の代表値、たとえば最小最大値の差を目標もしくは計測磁場の平均磁場強度で割った値を、請求項１８の相関図、または請求項１９，２０，２１の磁場補正作業で配置する磁性体又は永久磁石の大きさ表示と共に表示することを特徴とする表示装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超電導磁石装置に係り、核磁気共鳴断層撮像装置（Magnetic Resonance Imaging）に関する。
【背景技術】
【０００２】
核磁気共鳴を利用した診断では、磁場強度と診断箇所が対応しているので、マグネットシステムが発生する磁場強度に要求される精度は磁場強度の百万分の１程度の変動が問題にされる精度である。ＭＲＩ装置における磁場には大別して３種類がある。それらは
（１）時間的に定常で空間的にも一定な磁場で、通常０.１から数テスラ以上の強さである、撮像を行う空間（通常直径で３０−４０cmの球もしくは楕円体の空間）ないで数ppm程度の変動範囲である。
（２）１秒程度以下の時定数で変化して、空間的に傾斜した磁場。
（３）核磁気共鳴に対応した周波数（数ＭＨｚ以上）の高周波の電磁波によるもの。
【０００３】
このうち（１）の磁場は時間的に一定で、かつ空間的にも人体の断層撮影を行う領域では極めて高精度に均一性が磁場強度に要求される。高精度とは、たとえば４０cm直径の撮像空間ＦＯＶ（Field of View）で、±１.５ppmのように百万分の１のオーダの精度が要求される。このようにきわめて高精度な均一性が要求される磁場分布は、磁石を製作・励磁の後に磁場を精度よく調整する必要がある。一般に製作誤差による誤差磁場は均一磁場に要求される許容誤差磁場に比べて１０００倍以上大きい。製作後の据え付け時に要求される磁場調整（シミング）は数１００ppmから数ppmへの誤差磁場の低減を行うことになり、きわめて高精度な磁場調整装置およびその手法が要求される。
【０００４】
従来の手法としては線形計画法を用いた方法によるシミングがある。たとえば、特開２００１−８７２４５号公報や特開２００３−１６７９４１号公報に記述され、実際の機器の調整に適用された手法がある。しかし、この線形計画法による手法には、以下の課題がある。
（１）詳細な磁場計算を行うために多くの計算時間が必要。
（２）ここの鉄片や電流の設置や変化に高精度の磁場に対応した精度を要求。
（３）誤ったシミング作業を行った場合に、誤った箇所の特定が難しく回復に手間がかかる。
【０００５】
さらに、図２に示すように球面調和関数を用いて磁場分布を調整していることで課題が発生する。図２は従来の磁場調整法の例を示す図であり、球面関数を用いている（特開２００１−８７２４５号公報）。
【０００６】
球面調和関数は、球面上では直交して基底を構成しているが、磁場調整機構上また非球面の磁場評価面上では相互の干渉があり、精度の良い球面調和関数分布の磁場を発生させようとすると、細かい磁場調整機構上の調整が要求される。たとえば、均一な磁場分布が球面調和関数のもっとも次数の低い分布であるが、実際にこの分布を正確に出すことは、完全に磁場調整領域を取り囲む磁場調整機構でなければ不可能であり、従来技術が対象としているＭＲＩではこのような磁場調整機構は存在しない。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００１−８７２４５号公報
【特許文献２】特開２００３−１６７９４１号公報
【特許文献３】特開２００１−３２７４７８号公報
【非特許文献】
【０００８】
【非特許文献１】M. ABE, T. NAKAYAMA, S. OKAMURA, K. MATSUOKA ,“A new technique to optimize coil winding path for the arbitrarily distributed magnetic field and application to a helical confinement system”, Phys. Plasmas. Vol.10 No.4 (2003)1022.
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明の課題は、磁場調整装置とその手法で、前記述の問題点の解決を含み、かつ調整の進展状況と最終的な誤差磁場がどの程度に低減できるかの見通しを調整作業中に確認しながら、確実に調整を完了できる手法と装置を提供することにある。素早く調整を完了するために誤った作業時にも容易に自動的に訂正できる機能を含む手法およびその手法を含み磁場調整作業の指針を表示する装置を提供する。
【課題を解決するための手段】
【００１０】
目標磁場に対して曲面もしくは平面などの任意面上の電流分布を求める手法には、論文に記載されている電流ポテンシャルを用いる方法（非特許文献１）がある。この計算手法は論文中でＤＵＣＡＳと名付けられている。この手法を、特にこの手法の中で用いられている電流ポテンシャルと特異値分解の考え方を応用して磁場調整を行う。
【００１１】
非特許文献１のＤＵＣＡＳでは、補正すべき誤差磁場として入力する磁場分布は、プラズマ閉じこめの理論で決めた目標磁場と仮定した電流ポテンシャルなどで計算した磁場分布との差、つまり数値計算で求めた値であるが、本発明では実装置を対象とするので、目標磁場と計測磁場の差を誤差磁場とし、多数点の計測磁場を扱い誤差磁場分布を把握する。
【００１２】
また、非特許文献１では、電流ポテンシャルＴの分布を求め、電流密度ベクトルｊが電流ポテンシャルＴと面の法線のベクトル積で、（▽Ｔ）×ｎで電流が与えられるので、Ｔの等高線を線電流もしくはコイル形状としていたが、本発明では磁気モーメント分布もしくは鉄片密度分布とする。
【発明の効果】
【００１３】
本発明によれば、高精度な磁場を生成するＭＲＩ装置を低コストで製造することができる。また、ＭＲＩでなくても、高精度な磁場を必要とする磁石の磁場調整法に適用できる。
【図面の簡単な説明】
【００１４】
【図１】本発明の好適な一実施例である磁場調整フロー図を示す。
【図２】従来方法のシミングフロー図を示す。
【図３】本発明の好適な一実施例による磁場補正に必要な電流ポテンシャルと磁場調整用の磁化鉄片量の換算の考え方を示す図である。
【図４】本発明の一実施例で利用する計算体系の一般的体系例の図である。
【図５】本発明の一実施例である磁場調整に利用するＭＲＩ用磁石の磁場調整機構の配置図を示す。
【図６】本発明を図５の磁場調整機構に適用する計算モデルの図を示す。
【図７】本発明による磁場分布のスペクトル図を、シミングによる到達可能均一度と共に示す図であり、（ａ）はシミング前、（ｂ）はシミング後のスペクトルを示す図を示す。
【図８】本発明によるシミトレイへの磁場補正用の鉄片量配置の表示例を、電流ポテンシャル等高線と共に表示する図である。
【図９】本発明による鉄片表示のために升目内磁気モーメント計算および鉄片量換算の考え方を示す図である。
【図１０】本発明を、起磁力配置設計手法に用いた場合のフロー図である。
【図１１】本発明の一実施例の対称とする水平磁場型ＭＲＩの全体図。
【図１２】本発明の一実施例の対称とする水平磁場型ＭＲＩの磁石装置の断面図でシムトレイの位置と撮像領域の関係を示す図。
【図１３】水平磁場型ＭＲＩ装置のシムトレイを模擬する計算体系の図、シミング用の鉄量を求めるための電流ポテンシャルを評価する面をリング状に配置している。
【図１４】シムトレイの模式図で、位置に依存して磁場の微調整用の鉄片量を変化させている。
【図１５】電流ポテンシャル評価面上の電流ポテンシャルの計算結果と鉄片物量の配置を示す図（下側）。升目で区切られた領域の磁気モーメントを鉄量に換算して上側に示すシムトレイ上の同じ位置に鉄片を計算結果の物量で配置する。
【図１６】升目内の磁気モーメントを計算する考え方を示す図。シムトレイの升目内に相当する領域で節点の磁気モーメントを面積積分する。
【図１７】升目内の鉄片配置を示す図。鉄片物量は計算結果の物量だけ升目内に配置するが、あらかじめ用意したいくつかの物量に分割して配置する
【図１８】磁気モーメントを電流ループに置き換えた場合の概念図。小コイルであるループが囲む面積と電流で磁気モーメントを調整する。
【発明を実施するための形態】
【００１５】
以下、本発明の実施形態を図面を用いて説明する。
【００１６】
図３は、電流ポテンシャルと小コイル３（電流ループ４ｃ），永久磁石片４ｐと、それに磁化鉄片の等価性を示す。図である。鉄片４として、ボルト状の形状と板状の形状を示している。また、鉄片でなく磁石であっても表面の磁化電流２を介して等価であるが、この場合には周囲の磁場に依らず永久磁石４ｐの磁化で磁化電流２の方向が決まっている。
【００１７】
図３（ａ）は、計算時の有限要素１２と接点１１および電流ポテンシャルＴによる電流２１、図３（ｂ）は、小コイル３を流れる電流１による磁気モーメント生成、図３（ｃ）は、磁化した鉄片４による磁化電流２による磁気モーメントを示す。上側にボルト状の形状、下側に板状の形状を示す。これらは、鉄などの受動的に磁化する磁性体でなく、自発的に磁化している永久磁石でも、その磁化電流を磁化の程度で調整して考えると、等価である。但し、永久磁石の場合図３（ｄ）にように周囲の磁場とは無関係に磁化の方向が存在して、磁化電流も存在する。図３（ａ）のように、電流ポテンシャルＴがある接点１１に値があるとすると、大きさＴの電流ポテンシャルによる電流２１がその周りの接点との間に旋回して流れていると解釈できる。つまり、図３（ｂ）の小コイル３の電流ループに電流１が流れている状況と等価である。また、これは磁化した鉄片４の表面に磁化電流２がｊm（Ａ／ｍ）が流れている右側の状況と同等である。つまり、ＤＵＣＡＳで電流分布を表現するために用いられている電流ポテンシャル値Ｔは単位として［Ａ］の次元を持つが、これは磁気モーメント［Ａｍ2］の密度［１／ｍ2］で、［Ａ］の次元を持つと考えることもできる。一方、充分磁化した鉄片４は、磁化電流が囲む面積と磁力線方向の長さの積に磁気モーメントが比例するので、体積に比例する磁気モーメントを持つ。つまり、磁場調整時には電流ポテンシャルＴは、鉄片４の密度［単位面積あたりの重量ｇ／ｍ2もしくは体積cc／cm2］に比例する量である。この性質を利用し、また従来手法の球面調和関数でなくＤＵＣＡＳ中で用いられている特異値分解で得た固有分布関数と特異値を利用する。
【００１８】
これにより、ＤＵＣＡＳを利用し、磁場発生装置を対象として、磁場を調整する支援計算を行い、その調整用に配置する鉄片もしくは磁気モーメントの配置を表示する装置を提供する。作業者は表示に従って調整を進めることで目標の磁場分布に調整することができる。
【００１９】
本発明は任意の磁場分布を目標磁場とすることができるが、以下では主に目標磁場が一様に均一な磁場であることとして議論する。しかし、目標磁場が分布を持つかどうかは以下の議論に影響は与えない。単に、議論の理解を容易化するためである。
【００２０】
誤差磁場Ｂerr（ｒ）は位置の関数であるが、本発明では固有分布関数の組み合わせとして考える。つまり、
Ｂerr（ｒ）＝ΣＣmψm（ｒ）                         …（１）
である。従来法ではルジャンドル多項式とか球面調和関数を用いている。本発明では、特異値分解による分布関数を用いる。加算する関数ψmとその係数Ｃmの決め方を具体的に説明していく。
【００２１】
本発明の議論は、一般的な体系として図４の体系を考える。図４は、本実施例の計算体系を示す図である。電流ポテンシャル評価面１３と磁場計測評価点の集合１４より構成されている。一般的には電流ポテンシャル評価面１３は複数枚の面もあり得るがここでは各１枚として議論する。また磁場評価点は必ずしも面を構成しないがここでは面上の点として示している。
【００２２】
計測点ｊには３次元の磁場成分Ｂxj，Ｂyj，Ｂzjがあるが、１点の計測は計測位置とその位置で定義されている単位ベクトルｐで計測磁場成分を示す。空間の点としては１点でも、本発明の中では３個のデータということもある。
【００２３】
またＭＲＩ装置のように均一磁場を得ようとする場合には、磁場の軸方向の主成分のみを一定にすることになる。ＭＲＩでは磁場強度が一定であることが重要であるが、主成分以外が非常に弱いために、主成分の磁場がほぼ磁場強度に等しいためである。
【００２４】
計測値と目標磁場との差である誤差磁場には計測データは複数で、全体として列ベクトルとなり、Ｂeと表す。誤差磁場Ｂeは、計測磁場Ｂmと均一磁場に調整する時に調整する磁場強度Ｂtgの差ある。
【００２５】
計測点ｊに対応した誤差磁場は、成分がＢejであるベクトルＢeであり、
Ｂej＝Ｂtg−Ｂmj                                    …（２）
である。
【００２６】
特異値分解を適用する一般的な体系は図４の様である。磁場評価点の領域があり、その位置で磁場を計測する。ＣＣＳ面に磁場調整用の鉄片を配置する。この面をＭＲＩではシムトレイ面と呼ぶ。
【００２７】
この面上の鉄片密度と誤差磁場補正の関係を説明する。面を三角要素で分割し、その接点に電流ポテンシャルを割りあてる。これは非特許文献１に記述の通りである。磁場評価点の計測データを要素に持つ磁場ベクトルとＣＣＳ面上の電流ポテンシャルを要素に持つ電流ポテンシャルベクトルの間の関係は、
Ｂ＝ＡＴ                                            …（３）
である。この式は、電流面上の接点の電流ポテンシャル値を要素に持つベクトルＴから、磁場評価点の磁場の応答を示す式で、行列Ａはｍ（磁場計測点数）行ｎ（接点数）列である。
【００２８】
この行列Ａに接点の拘束条件を加えて独立な接点の電流ポテンシャルから磁場評価点への応答行列Ａ′とした行列を特異値分解して磁場分布と電流ポテンシャルの固有分布関数のセットを得る。つまり
磁場分布の基底である固有分布
ｕ1，ｕ2，ｕ3                                       …（４）
と電流ポテンシャルの基底である固有分布
ｖ1，ｖ2，ｖ3                                       …（５）
であり、ｕjとｖjには
λjｕj＝Ａ・ｖj                                     …（６）
の関係があるが、ここでλjは特異値である。また添え字ｊは特異値の大きさの順に固有分布に番号を付けた順番の数値である。一つの番号に対して、電流ポテンシャル分布と磁場分布を示す基底ベクトルが各一個ずつ対応している。この一つの番号に関連した２つの基底ベクトルと一つの特異値をまとめて一つの固有モードと呼ぶことにする。また、順番の番号ｊは固有モードの次数である。
【００２９】
特異値の大きな低次の固有モードは、単位電流ポテンシャル分布ｖjあたりの磁場強度がλjｕjで、あることから理解できるように、大きな磁場を発生できるといえる。一方、特異値の小さな固有モードで電流ポテンシャルを変化させても磁場の分布はわずかであるということになる。この性質は後に説明するが本磁場調整法に重要な役目を果たす。式（１）の分布関数としては特異値分解で得た固有ベクトルの分布を使う。
【００３０】
本発明で次数ｊ番目の固有モードに対応する誤差磁場の補正法または調整法（シミングにより誤差磁場を低減）を説明する。基底電流ポテンシャル分布ｖjの何倍の補正するべきかの係数Ｄjは誤差磁場分布から求めることができる。その大きさは、
Ｃj＝Ｂe・ｕj                                       …（７）
から
Ｄj＝−Ｃj／λj                                     …（８）
である。つまりｊ番目固有分布の誤差磁場はＤjｖjの電流ポテンシャル分布を与えることで、完全に補正できる。
【００３１】
電流ポテンシャルと鉄片密度の関係を次に述べる。鉄片は表面の磁化電流を考慮して磁気モーメントに置き換えることができる。鉄片の表面の磁化電流ｊm（Ａ／ｍ）は、
ｊm＝Ｍ／μ0                                        …（９）
が流れている。ここでＭは磁化（Ｔ）である。鉄片が飽和状態に有ればほぼＭは２.１Ｔ程度である。従って、ｊmは約１.７×１０6Ａ／ｍである。従って１立方メートルの体積の鉄は約１.７×１０6Ａｍ2（１７０Ａcm2／１cc）の磁気モーメントを持つ。この値は磁
石の種類、特に磁場強度に依存するので、個々に検討する必要はある。しかし、磁石の磁場強度が１Ｔ程度を超える磁石では、鉄片はほぼ飽和に近い状態に磁化していると考えて差し支えない。この時には鉄の磁気モーメントは体積に比例する。
【００３２】
この換算でｊ番目固有モードの誤差磁場を打ち消すには
ｄj＝−ｖjＣj／（λjｊm）                         …（１０）
の成分ｄjk（ｊ番目固有分布関数に対応するｋ番目鉄片補正点の鉄片量）に相当する体積密度（ｍ＝ｍ3／ｍ2）で鉄片を配置する。また、電流で磁場補正を行う場合には電流密度ベクトルｊが電流ポテンシャルＴと面の法線のベクトル積で、（▽Ｔ）×ｎで与えられる電流を配置する。
【００３３】
以上が基本的な補正の方法である。１つの誤差磁場の固有分布（ｕj）による成分の補正である。本発明の補正では、磁場分布の磁場調整手段の分布関数（ｖj，ｕj）がそれぞれの基底であり、かつ一対一で対応していることが特徴で、一つの固有分布成分を補正するためには、ただ一つの調整手段の分布関数を調整する。
【００３４】
本発明による手法でも補正すべき誤差磁場の固有分布関数はいくつも存在する。多数の固有モードから補正すべき固有モードの選択方法と補正法に上記の方法を拡張する。これについての基本的な考え方は次の項目である。
（１）小さい電流ポテンシャル（つまり少ない鉄片量）で大きな磁場を補正できる固有モードから選択する。この指標が特異値λjである。特異値は本計算体系では固有分布毎の単位電流ポテンシャルあたりの磁場強度であるので、この特異値の小さい固有分布は選択しない。また言い換えると特異値は単位鉄片量あたりの磁場強度に比例する値であるとも言える。一般に、できるだけ小さな物量で均一磁場を生成したいために、特異値の大きな固有分布を、調整に用いる。
（２）計測磁場に含まれる固有磁場分布の成分強度が小さなものは無視できる。内積［式（７）］で計算する成分強度が、目標とする均一磁場レベルで許容される誤差磁場に比べて十分小さな強度であれば、あえて補正する必要はない。特異値が大きな固有分でも、成分強度が小さい場合にはシミングに用いる必要はなく、選択しない。
（３）特に補正を必要と作業者が判断する固有分布関数を個々に選択し、内積で求める強度または人為的に決めた強度で補正する。たとえば、誤差磁場分布のピークが重なり局所的に大きな誤差磁場が発生している場合には、人為的にピークを下げる補正を適切な固有分布関数の選択と大きさで補正する。
（４）選択した固有分布関数の電流ポテンシャル成分を補正した後の均一度（到達均一度）を求め固有分布関数の選択が妥当であるかどうかを判定する。到達均一度が十分でなければ固有分布関数の選択を再検討する。ここで均一度は、磁場評価領域の複数計測点のなかで、磁場強度の最大，最小の差、言い換えると、誤差磁場の最大ピークから最小ピークの差について、平均磁場に対する割合を示したもので、ＭＲＩでは通常１／百万（ppm）のオーダで議論する。
（５）目標磁場を変更すると、誤差磁場に含まれる各固有分布の強度と、残差として残る磁場の強度、つまり到達均一度も変化するので、固有分布選択には目標磁場も考慮が必要である。
（６）数回から数１０回の回数を繰り返して調整を行う。これは、調整機構の精度が、目標とする磁場精度に比べて通常荒いために繰り返し操作して磁場精度を向上させる。たとえば、ＭＲＩの磁場調整（シミング）では、１マイクロＴの精度で磁場調整を行う必要があるが、調整以前の誤差磁場は数ｍＴ程度のものである。これを一気に調整することにした場合、調整のために配置する鉄片は１／１０００より細かい精度での物量の管理が要求されるが、実作業でこの精度を管理することは容易でない。そこで、本発明によると一回の調整では１／１０程度以下の精度での物量管理で、回数を追う毎に誤差磁場を低下させ、最終的な磁場精度との相対比を「低下させ、１／１０以下の物量管理でも十分な最終磁場精度を得る。
【００３５】
次に、鉄片で磁場調整を行う場合を考え、鉄片配置とこれまでの説明の関係を述べる。選択した固有分布に対応する補正量Ｄはそれぞれの固有分布関数による補正量の和である。
Ｄ＝Σｄj＝Σ−ｖjＣj／（λjｊm）                 …（１１）
【００３６】
ここで、和Σは選択した固有分布関数に対して実行される。この補正を実行した後に撮像領域の磁場分布がどのようになるかは容易に計算予測できる。
【００３７】
一つの方法は、固有分布関数の磁場分布の関数から求める方法である。次式となる。
Ｂshim＝Ｂe−ΣＣjｕj                             …（１２）
【００３８】
ここで和Σは選択した固有分布関数に対して実行される。
【００３９】
他の一つは、再構成した電流ポテンシャルにより再構成される方法である。式（１２）で電流ポテンシャルの補正分ΔＴは
ΔＴ＝Σ−ｖjＣj／λj                             …（１３）
であり、補正後の誤差磁場分布は
Ｂshim＝Ｂe−ＡΔＴ                               …（１４）
で与えられる。この２つの方法は同じ計算結果を与える。ここでも和Σは選択した固有分布関数に対して実行される。この計算法で、磁場調整後の到達均一度を予測して目標の精度で磁場調整が進行しているかどうかを判定する。
【００４０】
判定には、到達均一度だけでなく、シミングに必要な鉄片物量も参考にして、過大な鉄片が必要ならば、固有分布関数の選択を再検討する。どのような選択を行っても、過大な鉄片量が必要と計算される場合には、磁石の設計もしくは製作不良と判断できる。この機能は
（ａ）製作完了時の品質管理
（ｂ）コイルや磁性体配置の設計時に、起磁力配置の設計が妥当性，配置の再検討要否の検討に利用できる。必要な鉄量は式（９）に関連して議論したようにおおよそ１７０Ａcm2／１ccで換算可能である。また、補正に必要な磁気モーメントは、鉄片を配置する領域で
∫Ｔｄｓ＝ΣＴｉΣＳｉｊ／３.０（Ａｍ2）          …（１５）
の面積積分であるが、積分領域については実例で議論する。また離散化した表示はメッシュで分割して計算しているために実際の計算内容を示している。前に書いている和Σは積分領域内の節点ｉで行われ、Ｓｉｊはｉ番目節点に属する要素ｊの面積である。三角要素であるために１／３がｉ番目節点に寄与すると考えている。後ろに書いている和Σは、ｉ番目節点が属する三角要素ｊについて実行される。以下では簡単に
Ｓｉ＝ΣＳｉｊ／３.０                             …（１６）
として議論する。
【００４１】
以上の議論では電流ポテンシャルを変数として議論してきたが、節点が関与する面積Ｓｉを最初から考え、
ＴｉＳｉ→Ｍｉ                                    …（１７）
と変換すれば、磁気モーメントの大きさＭｉを扱う議論になるが、これまでの議論は単に面積を大きさだけの倍率で変換するのみで、特異値分解を利用する議論としては全く同じとなる。この場合、式１６の積分は単純にその領域に属する磁気モーメントの和となる。
【００４２】
前記の本発明における改良項目の作用を順に記述する。
【００４３】
項目（１）の固有モードの選択は特に低次側の固有分布関数について誤差磁場を補正するために選択する。比較的少量の鉄片で磁場が補正できる範囲で低次の分布関数を選択する。低次側のみと言っても通常数十から数百の個数の固有分布関数を選択する。固有分布関数の鉄片（電流ポテンシャル）配置に従って磁場を補正することで、選択しなかった固有分布に対しては大きな影響、新たな誤差磁場を与えないように、補正できる。これは、特に、選択しなかった高次（番号付けした大きな番号の固有分布）を乱さないという点で利点がある。つまり磁場調整を行っているときに、選択しなかった高次固有分布が乱れて作業が繁雑になることはない。
【００４４】
特異値分解で選択した低次固有分布関数は少量の鉄片で補正できるが、高次を変化させるには多くの鉄量を必要とする。高次部分を乱さない理由には、分布が直交していること以外に、高次に多くの量の鉄片が必要という理由もある。つまり、少量の鉄片を配置する低次分布関数の補正では、誤差的に配置が乱れても高次成分の強度を変化させるには至らない。この意味でも低次側の固有分布から選択して補正していく。
【００４５】
また、低次は補正できる磁場が特異値に比例して大きいために、効率的につまり少量の鉄片で磁場調整つまりシミングが実行できる。
【００４６】
項目（２）は補正の必要でない固有分布関数成分は補正しない。しかし、選択して補正量に含めても、補正量が小さいため、前述のように高次成分を乱すこともないために、問題は発生しない。
【００４７】
項目（３）は、鉄片配置量と磁場分布を調整するために選択を調整する。鉄片のみで磁場を補正する場合には、鉄片を取り去る負の鉄片量による調整が実行しにくい場合もある。一方、高次分布は鉄片を配置しても小さな磁場しか発生しない。つまり、高次成分の鉄片を配置して、低次の補正時に鉄片を取り去る余地を作る。また、正のピークから負のピーク値の範囲で均一度を定義すると、特にピーク部分のみ集中して均一度の表示が悪くなる場合もある。この場合には人為的に適当な補正成分を加える。これにより目標の均一度に達しやすくなる。
【００４８】
項目（４）で、目標の精度で磁場を調整できるかどうかをチェックできる。選択した固有分布関数に対して磁場を補正したときに、繰り返し補正していけば最終的に目標の均一度に達することが必要である。本手法ではすでに述べた計算方法でどの程度までの均一度が得られるかを推定できる。この推定に従って、固有分布関数の選択を変更するかどうか、また均一度が非常に悪い値しか得ることができない場合は製作上の問題があり、品質に問題があると判断できる。品質上の問題は補正が実行し難い高次成分で発生しやすいが、本発明の手法では特異値分解で成分を分けているので、高次の固有モードで発生している問題を容易に発見できる。
【００４９】
項目（５）は均一とする磁場強度をいくらの強度にするかの選択である。目標磁場を変更しながら固有分布関数の選択を行い、到達均一度と鉄片量をチェックして、均一度がよく、鉄片の配置も容易な目標磁場を選ぶ。鉄片の配置が容易ということは、単に量が少ないと言うことではなく、必要な比較的低次の分布関数が十分補正可能で、負の鉄量の計算値で鉄片が配置できない領域が存在しない配置である。
【００５０】
項目（６）は計測から鉄片配置の作業を繰り返して、磁場調整を完了させる。磁場調整を実行する幾何学的配置にも依存するが、選択する特異値の大きさは低次と高次で大きい場合には４桁程度の差がある。つまり、磁場調整中に１００cc程度の物量を扱う調整から
０.０１ccを扱う調整まで、固有分布関数の選択で変化していく。一方、鉄片の大きさを
１／１０を下る精度で管理することは実用上手間をとる。そこで、繰り返し調整を実行して、磁場の残差を補正して一回ごとの調整は１／１０程度の精度でも、繰り返しで良い均一度を出せるようにする。繰り返し中に、初期には高次まで固有分布を選択して磁場調整を行って物量の大きな調整を行い、徐々に次数の上限を下げて行くようにする。次数を下げるときには、高次部分が十分補正できたことを確認する。また次数を下げるに従って、補正として計算結果が要求する物量も減少してくる。従って物量の１／１０程度である調整精度も向上していく。
【００５１】
繰り返しで磁場調整を行うが、項目（２）に関連して説明したように、選択しなかった高次の固有モードは乱されることはない。このため、予知する均一度が繰り返し中に変化することはない。
【００５２】
以上のように、シムトレイ面の電流ポテンシャルから撮像領域に置いた磁場評価点の磁場強度への応答行列の特異値分解で得る固有分布関数を基礎にする本発明の手法によれば、誤差磁場に対して少ない補正量で、調整後の磁場を予知しながら、目標の磁場均一度が得られる。繰り返しシミング作業を行う必要がある。また繰り返し作業途中では、特に高次分布のシミングを行っている場合には低次分布関数に対応した誤差磁場成分が増大し、均一度が見かけ上悪化する場合があるが、本手法では固有分布関数の成分強度を確認することで、選択した固有分布の高次側成分が補正され、弱くなり、磁場調整が進展していることが確認できる利点もある。これは作業者の作業妥当性確認にも有利である。また、調整終了時に到達できる均一度も把握できるために製作上の誤差で目標の精度までの磁場調整が不可能で無いことも確認しながら安心して調整を進めることが可能な利点もある。
【実施例１】
【００５３】
実施例を説明する。第一の実施例として垂直磁場を持つ開放型ＭＲＩ装置の磁場調整（シミング）への適用を説明する。図５はＭＲＩ装置の磁石が発生する磁場をシミング（磁場調整）の体系である。この図は磁場（磁力線）の方向が垂直方向を向いている開放型機を想定している。開放型ＭＲＩ磁石の概念形状も図５中に示した。上下に分割された磁石装置６２があり、その間を連結注が結んでいる。これらの中には、断熱用の真空を確保する真空容器６２ｃ，輻射シールド６２ｄ，極低温容器６２ｅ，磁気遮蔽用コイル６２ｂを含むコイル６２ａ群が存在する。被検診者用ベッド６１の上に被検診者が横たわり、核磁気共鳴断層撮影を行う。
【００５４】
撮像領域の空間（磁場計測評価領域）６があり、その表面もしくは内部の面も含む磁場評価点の磁場分布を均一に調整（シミング）する。図１は本実施例において、磁場分布を調整するシミングのフローを示す。ＭＲＩの撮像領域の磁場調整に適用する実施例である。地面に対して垂直方向の磁場成分の強さで、撮像領域の上下に磁場調整機構の面（シムトレイ５）があり、その面上に鉄片４を配置する。
【００５５】
ＭＲＩ装置のシミングに本実施例を適用した場合のメッシュ生成例を図６に示す。この実施例では球体面上が磁場計測評価点の集合１４で、数百点の磁場計測点が配置されている。球体の上下の円盤面がシミングを実行するときに鉄片４を配置する面の計算モデル、つまり電流ポテンシャル評価面１３である。図の右側に有限要素で概略を示したように、この面上に接点を持つ三角要素で構成される有限要素計算の体系を構成している。
【００５６】
図１の破線内の予備計算部分１Ｂをシミング作業に先だって計算しておく、特異値分解計算ステップ３２Ｓを含み、計算メッシュ生成ステップ３１Ｓ、特異値分解の結果である固有分布関数と特異値の保存ステップ３３Ｓを含む。この部分は数千点の電流ポテンシャル値に対応した接点から撮像領域の数百点程度の磁場計測点への応答行列Ａの特異値分解を含む予備計算部分１Ｂであり、比較的長時間の計算時間を必要とする。このため、磁石の体系に合わせた計算体系でシミング用の固有分布関数を計算しておき、シミング作業中の計算時間を短縮する。
【００５７】
あらかじめ計算したデータは計算機の記憶領域に保存ステップ３３Ｓで保存しておき、必要時に読み出し（特異値分解結果の読み出しステップ１６Ｓ）て利用する。つまり計算機中には、磁場分布の基底ベクトル群である固有分布関数が数個以上、また電流面の分布関数である基底ベクトル群が同数、そして両者の大きさの変換情報である特異値が同数の組み合わせで保存されている。
【００５８】
磁石の励磁後、しばらく時間が経過した後に磁場調整（シミング）開始１１Ｓとなる。図１のフローに従って作業が行われる。磁場計測ステップ１２Ｓし、磁場分布データ保存ステップ１３Ｓと磁場データ読み出しステップ１４Ｓの後、磁場均一度判断ステップ１５Ｓで均一度良好かどうかを判断する。十分な均一度であればシミングの必要はなく、磁場調整終了ステップ４０Ｓのステップとなる。これは、十分な均一度で利用されていた装置が保守時などで消磁した後に再励磁した場合にはあり得るが、新しい磁石では製作誤差で数百から千ppm程度の均一度であり、磁場調整（シミング）が必要と判断される。
【００５９】
そこで、固有モード選択と目標磁場決定ステップ１７Ｓに移る。次のステップ１８では、各固有モードの強度Ｃj，補正電流ポテンシャルΔＴ，補正鉄片配置と補正磁場分布、および到達可能均一度の計算等の式（１−１４）の計算を選択した固有モードに対して行う。
【００６０】
次にステップは固有モード選択妥当性判断のための表示ステップ１９Ｓである。ステップ１８の計算結果を表示して固有モード選択の妥当性を判断する。表示は主に２つがあるが、一つは図７に示し、他の一つは図８に示した。
【００６１】
図７は、式（７）に示すような式で求めた誤差磁場に含まれる磁場の固有分布の強さを縦軸に、固有モードの次数を横軸に示した図でスペクトルと呼ぶことにする。縦軸は対数スケールで示している。尚、図７には、固有モード選択の範囲，致達可能均一度も示す。また、図８は、シミング作業用の鉄片配置量の表示例を、電流ポテンシャル等高線と共に示す。
【００６２】
本例で示す計算内容は、図５，図６の体系であるが、磁場評価点を直径４０cmの面上と
し、この表面の誤差磁場を均一度２０ppm以下とする均一な磁場を目標としている。
【００６３】
この図７のスペクトル表示を参考にして、補正する固有分布関数を選択する。図中でｘが個々の固有モードに対応しているが、○で囲ったものが選択した固有モード１５である。○で囲ってないものは非選択の固有モードである。この選択は既に記述した方法で行う。補正をする分布関数を選択すると、その誤差磁場成分を計測誤差磁場から差し引くことで、到達できる均一度を計算して予測できる。図７では上部に楕円で囲ったところに、到達可能均一度１７として示している。
【００６４】
２つのスペクトルを示したが、図７（ａ）はシミング前で、図７（ｂ）はシミング後である。シミング前にはこの例では７２６ppmであり、スペクトル図では低次モードの誤差
磁場成分が大きいことがわかる。○は次数８０以下で、計測精度のほぼ下限以上の強さの誤差磁場成分として選択した固有モードである。この例では、選択した固有モードを補正すれば１５.２５ppmと予測されている。固有モード選択の次数上限を示す線２２と固有モード選択で強さ下限を示す線２３をスペクトルの図である図７上で表示して固有モードを選択する。
【００６５】
予測した到達可能な均一度が不十分であれば再度固有分布関数の選択を再検討する。固有分布関数の数を調整、つまり固有分布関数選択範囲の番号上下限や固有モード強度Ｃjの下限を調整する。また個々に選択した固有分布関数の補正割合を調整するなどの選択肢もある。
【００６６】
ステップ１９Ｓでの他のもう一つの表示は、図８の鉄片配置量の指示図で、シミング可能かどうかのチェックに利用する。図の円形は図５に示したシムトレイ５を示している。上下の２枚あるが、この図では下側のシムトレイである。図中の升目７はシムトレイ５に配置された区画であり、それぞれに番地が割りつけてある。図７では左右方向にＡ，Ｂ，Ｃ，・・・、上下方向に１，２，３，・・・で番地を指定する。升目７内の数値は、升目に配置する鉄体積１８を示している。図８では０.１cc単位である。升目には５cc程度の鉄片は十分配置できる構造であり、表示の量は十分少なく、配置可能である。繰り返しの調整中には取り扱う鉄片量が徐々に小さくなるので、当初の１／１０，１／１００，１／１０００と小さい単位で表示していく。
【００６７】
図８には、シムトレイ５を計算もモデル化した電流ポテンシャル評価面１３上で、升目７，升目の鉄片物量１０に加え、電流ポテンシャルの等高線１９を示した。電流ポテンシャル等高線１９をコイル形状と考えると、この形状のコイルで誤差磁場は補正できる。このことは先に述べた既公表論文に書かれている。本発明での等高線表示には別の利点がある。特異値分解で求めた分布関数は面上に広がりをもって鉄片もしくは磁気モーメントの配置を要求する。しかし等高線のピーク８や等高線の谷９付近に最も多くの配置（もしくは除去）を要求している。この２つ性質を利用して磁場調整用の鉄片配置位置を柔軟に考える。配置や除去の制限がなければ等高線ピークの周りの同符号で閉じた電流ポテンシャル等高線１９の鉄片量を加算して等高線のピーク８位置付近に配置（から除去）する。またその付近が、例えばシムトレイ固定用の支持具が有るなどして、配置できない場合には、同じ量の配置を、閉じた等高線領域内で他の部分に配置（から除去）してもかまわない。図７の等高線のピーク８はＬ−Ｍ間の線と７−８間線上にピークがあるが、その周りのＬ，Ｍと８，９の交差部分の升目には合計５＋３＋２＋１＝１１の量の配置が要求されている。そこで、本発明ではＬ−Ｍ間の線と７−８間線上に等高線のピーク８位置に１１の量の鉄片を配置して良い。このような配置を行うことで作業量低減と、配置位置精度が緩和されることで、作業が容易になる。
【００６８】
升目内の磁気モーメントもしくは鉄量の計算方法を図９参照で述べる。図９は、本実施例の磁場調整用計算で接点ポテンシャル値から磁気モーメント，鉄片量への換算概念を示す。既に式（９）に続く記述で、鉄の体積と磁気モーメントは比例する事を述べた。また電流ポテンシャルは単位面積当たりの磁気モーメントを示していると解釈できる。そこで、有る領域の鉄片量を求めるには、その領域で電流ポテンシャルＴを面積積分し、磁気モーメントその領域に必要な磁気モーメントとし、既に述べたように鉄量に換算する。図９は、図８に示した升目７と接点の関係を模式的に示した。×で示した点が接点である。連続な分布関数ではないので、例えば、図中の式のように接点とその接点に対応する面積の積を加算して升目７内の磁気モーメントとする。接点に対応する面積は例えば接点が属している要素の１／３（△要素の場合）とする等の方法がある。
【００６９】
升目７の大きさおよび計算時の要素サイズについてコメントする。升目７の大きさは図８に示す鉄量配置分布の分解能を持つ程度の細かさが要求される。均一度を得るために必要な固有モードの次数の上限付近で等高線分布を確認し、そのピークや谷部分のサイズの大きさより小さなものとする。一方、細かい升目にすると手間がかかることになる。図８では等高線のピーク８と等高線の谷９部分の最小の大きさに比べてほぼ同程度としている。同程度であるために、細かい等高線の部分では升目７のみでは十分な分解能が無いところもある。この場合には配置する位置を等高線のピークと谷の位置を参考にして位置を調整して鉄片を配置する。有限要素の大きさは接点の個数が升目内にいくつ存在するかで決める。既に述べたように鉄片配置量の精度は１／１０程度で良く、繰り返しで均一度を高めていく。接点の個数も５個程度以上有れば、対応する面積に若干の誤差が発生しても、十分な精度であると考える。図８の例では片面が１５００程度の接点。図６のように上下のシムトレイを電流評価面として、２枚合計で３０００接点程度以上と考える。
【００７０】
図７（ａ）と図８の鉄片量表示による確認で十分なシミングができると予測されたので、この例では鉄片配置作業ステップ２２Ｓに移った。この予測の根拠は、到達可能均一度１７が目標値に比べて十分よく、鉄片配置が可能な量であることである。
【００７１】
シミング可否の判断ステップ２０Ｓで到達可能均一度もしくは鉄片量が不適と判断されると再度固有モード選択と目標磁場決定ステップ１７Ｓに戻る。しかし、各種条件を変更しても目標の磁場調整可能かを判断する磁場調整可能判断ステップ２１Ｓで磁場調整不可と判断されると、磁石は不良であり、修理・調整のステップ４１Ｓに入る。
【００７２】
図７（ｂ）はシミング終了時のスペクトルである。１７ppmの均一度までシミングできている。これは当初予測した１５ppm台と比較して、大差ない均一度で、良い精度で均一度の予測ができていたと言える。図７（ａ）のスペクトルの均一度に至るまでには繰り返し作業を図１のフローのように行っている。この繰り返しの必要性については既に説明したが、後で、実例で説明する。
【００７３】
ステップ２１Ｓで磁場調整可能性を判断するが、この内容を説明する。ステップ２１を経由してステップ１７に戻り、固有モード選択を再検討しても十分な均一度が、適切な補正量（シム鉄片量）で得られない場合もある。つまりこれは磁石の製作精度などが不十分で、磁場が不良で、目標の均一度を得ようとすると大量の鉄片の配置が必要で、実際には不可能な場合である。この評価によって、磁場調整を行うことなく磁場の不良を見知できる。不良であれば、適切な手直しを行うが、補正量の分布から問題の箇所を推定することもできる。また手直し不可であれば、製品の不良と判断でき、磁場調整を繰り返して人手を使うことなく、判断ができる利点も本発明に含まれている。
【００７４】
到達可能な均一度の予測が十分な均一度を持っていれば、補正に必要な鉄片分布計算結果を、紙面印刷やプロジェクトによる拡大表示などで出力し、その分布に従って、シミング用鉄片を配置する作業を行う。シミング作業で配置する鉄片量には、その量や位置に誤差があり、また鉄片磁化の程度も鉄片の材料性質や磁石内の磁場分布に依存するため、電流ポテンシャルからの鉄片換算に誤差も存在する。このため、一度の作業では到達可能な均一度には達しない。そのため、図１のように繰り返し作業を行い、磁場を均一に近づけていく。
【実施例２】
【００７５】
第２の実施例を説明する。既に製作後の品質の検査に使えることを既に述べたが、この手法は同様の判定で磁石の設計に利用できる。この場合のフローを図１０に示す。本実施例では、計算上で磁場調整を行い、目標磁場精度を達成できることを確認することで起磁力配置設計に適用する。起磁力配置検討開始５１Ｓの後、起磁力配置を仮定するステップ５２Ｓする。この起磁力配置に基づいて磁場計算ステップ５３Ｓを行う。また起磁力配置からシミングトレイの配置に基づいて特異値分解を実行し、その結果を保存しておく。この予備計算部分１Ｂは第１の実施例と同じである。１Ｂの部分は起磁力配置仮定ステップ５２Ｓに基づき、起磁力配置改善判断ステップ５６Ｓでシムトレイ変更が必要であると判断された場合のみ、実行する。この予備計算部１Ｂは図１と同様である。５６Ｓで既存の特異値分解データが利用できるかどうかを、判断し、その結果によっては単に特異値分解結果のデータセットを読み出す１６Ｓのみとする。
【００７６】
磁場計算結果を読み出しステップ１４Ｓ、磁場均一度を判断するチェック１５Ｓし、既に良好な均一度であればＭＲＩ磁石用起磁力配置候補案ステップ５７Ｓとなる。一般には、磁場分布から既に説明した本発明の手法でシミングにより均一度が良くなるかを判断する。均一度がシミングにより十分には改善されない場合、またはシミングに必要な鉄片量が過大でシミングが不可と判断される場合は、再度、起磁力配置仮定ステップ５２Ｓに戻る。この判断部分は図１の１７〜２１Ｓのステップと同じである。３Ｂで示す。磁場調整計算部分３Ｂも、図１の３Ｂと同じである。
【００７７】
シミングが可能であれば、さらに５６Ｓのステップで起磁力配置補正して改善するかどうかを、シミングに必要な鉄片量を参考にして、また磁石全体の構造設計を参考にして検討する。起磁力配置の再検討を行わなければ、候補の起磁力配置案５７Ｓとなる。また起磁力配置を補正する場合には再度起磁力配置仮定の部分に戻る。起磁力配置を再検討する場合は、例えば超伝導コイルの経験磁界が過大であるとか、支持構造が難しい電磁力である場合等が考えられる。
【００７８】
このようにして、本発明のシミングを計算上で仮想的に行い、起磁力配置の候補を得る。起磁力配置の設計には、十分な均一度が得られると判定された起磁力配置については、起磁力量，電磁力さらに応力などの全体設計を行い、さらに磁石の成立性を判断する。成立性が難しいとわかれば、再度起磁力配置を仮定するところから始める。
【００７９】
図１では磁場調整手段として、磁化した鉄片４の磁気モーメントを利用する方法で記述して、鉄片配置作業ステップ２２Ｓとしている。しかし、図３で説明したように磁化した鉄片は小コイル３による電流と等価である。そこで、図８の升目状に小コイルを並べ、その電流１を、本手法で計算した磁気モーメント分布に従って調整することも鉄片配置作業ステップ２２Ｓの代替えとして可能である。
【００８０】
本発明の手法では負の量を要求する磁場調整が要求されることがある。小コイルによる電流調整で、磁場調整を行う場合には、電流極性を変えることで、また永久磁石を使う場合には方向を変えることで対応できるが、鉄片の磁化は周囲の磁場環境で決まり、極性を変更できない。この場合の負の量は次のように考える。高次部分まで選択している場合には負の量の部分は配置しない。正の物量部分のみ配置することで、誤差磁場の凹凸振動空間的波長が２倍程度になり、負の物量を配置しなかったことで打ち消されなかった誤差磁場成分の固有モード次数は約２倍の次数の高次側に移り、特異値の減少で、磁場強度も弱まり、磁場調整上無視できるようになる。しかし、比較的低次の固有モード選択となって磁場調整を行っている場合には、負の量の要求は通常既に高次までの磁場調整で配置した鉄片の量を削減する方法で調整できる。しかし時には既に升目内に鉄量が零となっている場合には、その近傍から取り去る。近傍とは等高線の閉じた線の領域からである。それでも採るべき鉄片がない場合には、人為的に特定の固有モードを補正に必要な大きさからはずし、負の量が無くなるようにする。特定の固有モードは高次のモードから選択することで、磁場に与える影響を小さくできる。
【００８１】
図１や図１０の実施例で、磁場調整計算部分３Ｂは特にまとまったソフトウエアとしておくことで、磁場調整を行う時には特異値分解結果の保存データと共にして機動性の良い磁場調整の支援ツールとなる。
【実施例３】
【００８２】
以上の計算方法，シミングの手順は垂直磁場型の磁石装置６２でも水平磁場型の磁石装置６２でも同じ考え方で実行できるので、水平磁場型のＭＲＩ用磁石装置６２にも、磁場調整（シミング），品質の管理、および起磁力配置設計に応用できる。但し、磁石装置６２の形状が異なるために、計算手順や鉄片４配置位置などの考え方が異なる点も存在する。そのため以下で実施例３として説明する。実施例３では、図１１に示す水平磁場型ＭＲＩ用の磁石装置６２に適用する。この場合図１１の磁石装置中心部を貫通するボア６２ｆ（筒状の中空穴）に、筒状にシミングに利用する領域（シムトレイ５）を配置する。シムトレイ５面は図１２に断面で示すように被健診者を取り囲むような筒状の領域となる。筒状のシムトレイ５を取り囲みやはり筒状の断熱用の真空容器６２Ｃで囲まれ、その内部に輻射シールド６２ｄと極低温容器６２ｅを持ち、極低温容器６２ｅの中にはコイル群６２ａが配置されている。以上の内容物を含む磁石装置６２がシムトレイ５を取り囲んでいる。撮像領域６はこの磁石の幾何学的な中心であり、３つの直交する対称軸の交点を中心とする点線で囲まれた領域である。その中に傾斜磁場１９を発生させながら撮像する。傾斜磁場コイルもシムトレイ５と同様の領域に配置されるが、ここでは省略している。シムトレイ５はこの傾斜磁場コイル集合体の内部に配置されることが多い。
【００８３】
このシムトレイ５の筒の断面は必ずしも円形断面でなくても良い。特異値分解を利用するので、シムトレイ上の磁性体から撮像領域の磁場への応答行列は任意の計算体系に応用できる。例えば特許文献３の傾斜磁場コイルの断面形状のコイル群の間にこのシムトレイを配置する場合には楕円断面の筒状となるが、本手法では既に記述した手法でシミングを実行できる。
【００８４】
この水平磁場型のＭＲｉ装置用磁石装置のシムトレイ５のメッシュは図６の垂直磁場でのメッシュと異なる。その理由は、筒状のシムトレイ５を構成する面が、磁場に平行となるためである。図６ではシムトレイ面が磁場に垂直の面であった。
【００８５】
ＤＵＣＡＳ計算を利用するには、磁化の方向が面に垂直であることが必要であり、鉄片は周囲の磁場の影響で磁化し、その磁化の方向はほぼ周囲の磁場の静磁場磁力線方向６５である。したがって、鉄片４の磁化の方向はシムトレイ面の面内方向になる。その解決のため、図１３のようにリング状の電流ポテンシャル評価面１３を軸方向（水平方向の磁場の磁力線方向）に多数配置する。面の端とならない節点１１（面内節点と呼ぶ事にする）を周回方向に多数確保しておく。面内の節点１１は図３の左部に示したように、面内節点の周りに周回電流を想定することに等価であり、磁気モーメントＭを配置することになる。
【００８６】
水平磁場機への実施で、シムトレイ面の様子をのべる。垂直磁場機では図５のように平面の対であった鉄片４を配置する領域（シムトレイ５）は、水平磁場機では図８のように磁石のボア６２ｆに配置する筒状の位置になり、図１４のように鉄片４を周回方向と軸方向に分布させて配置する。この配置は既に述べた垂直磁場機（開放型機）と共通の計算方法で求めた鉄の物量１０に従って配置する。鉄片４は図中のように、軸方向位置（磁場方向位置）と筒状の断面の周回方向位置で分布させて配置する。この様子を図１４では鉄片厚さを変化させて描いている。
【００８７】
垂直磁場型ＭＲＩ用の磁石装置６２で図８であった鉄片配置の分布と補正用の電流ポテンシャル分布を示す図は、図１５下部分に書くように周回方向を示す角度座標及び角度方向番地と軸方向の位置を占めす位置座標と番地で示す分布図となる。図１５上部分にはシムトレイ５を書いている。シムトレイ５にも電流ポテンシャル評価面１３と同じ番地を配置し、計算で求めた番地の物量を配置していく。周回方向に２４分割、軸方向（静磁場磁力線方向６５）にＡ〜Ｎの１４分割して示しているが、この番地の分割方法は必要な磁場精度（均一度）達成するために、補正を行う次数（固有モードの数）を検討し、その次数の分布を再現できる数の分割を必要とする。周回方向，軸方向共に１５〜３０程度の数である。升目７には図８のように鉄片４の量を示す数字も書かれるが、ここでは表示を省略して図１５を書いている。
【００８８】
垂直磁場型ＭＲＩ用の磁石装置６２垂直磁場機では、升目を配置する面と補正用の電流ポテンシャルを計算する面（従ってシムトレイ面）が同じ面上に存在したが、水平磁場機では、磁場方向との関連で、電流ポテンシャル評価面１３はシムトレイ５の筒状面とは異なっている。そのため、図１５の表示は、各リング状の電流評価面１３（図１４）で計算した電流ポテンシャル分布の中間節点１１ｃの値を、筒状面に磁石の幾何学的中心を通る線で投影して表示する。全ての中間節点１１ｃを同じ筒状面に配置しておくと、この表示には都合が良いことは考慮すべきである。
【００８９】
鉄片４の配置量を表示する鉄片物量１０の計算について、図１６を使って説明する。鉄片配置量は電流ポテンシャルの面積分で計算されるが、電流評価面（電流ポテンシャル評価面１３）がシムトレイ面と異なっているので、図１６のようにいくつかのリング状の面を跨って積算する。従って、升目７に相当する領域に１個以上のリング状電流面が配置されるように、リング状電流面数を多めに配置すると都合が良い。
【００９０】
実際のシミングでの升目７枠内の様子を説明する。図１７に示すように、実際の升目内は、実施例１と実施例２の場合でも同様で、いくつかの体積の異なる鉄片４を配置することを想定する。図１７には２種類のケースを示しているが、水平磁場用の磁石装置６２には主に下側（図（ｂ））に書いた板状の鉄片４を組み合わせる方法を採る。これは、図１４で議論したように鉄片配置のシムトレイ５の形状に関連した鉄片形状である。
【００９１】
その磁気モーメントの大きさはＭｆ１〜Ｍｆ３のように異なる。図９や図１６に示した必要な磁気モーメント
Ｍｆ＝ΣＴｉ×Ｓｉ（Ａｍ2）                       …（１８）
を発生できるように、枠内の鉄片物量１０（体積）を、体積の異なる鉄片４の組み合わせで調整する。この式で、和は升目７に含まれる位置の電流ポテンシャル評価面を跨いで実行される。また、Ｔｉは枠内の接点ｉの電流ポテンシャル値（Ａ）で、Ｓｉはその接点に付属する電流ポテンシャル評価面１３面上の要素面積である。接点は複数の要素に付属するので、ここで示す三角要素では各要素の１／３が個々の接点１１に属すると考えて問題はない。鉄片物量１０と磁気モーメントＭの換算方法は既に図３で議論した１７０Ａcm2／１cc程度である。升目７に必要とされた磁気モーメントを鉄片の体積に換算して必要な体積を升目内に配置する。磁場が弱く、鉄片が磁気飽和してない場合には、磁化Ｍは飽和磁化とは異なり、この換算係数も異なるが、この場合には材料の磁化曲線（Ｍ−Ｈ曲線，Ｍ＝磁化強度Ｔ，Ｈ＝磁界の強さＡ／ｍもしくはＴ）を参考にしてＭの大きさを決めておく。これらの考え方は垂直磁場（開放型）機用の磁場調整と同じ考え方である。
【００９２】
図１７の上側部（図１７（ｂ））に書いた升目７の様子には鉄片４や永久磁石４Ｐのボルト使った場合の配置を模式的に書いている。この場合でもいくつかの磁気モーメント（Ｍｆ１〜Ｍｆ３，Ｍｐ１〜Ｍｐ３）の大きさに相当する物量を組み合わせて、目標の量とする。
【実施例４】
【００９３】
実施例１と実施例２では鉄片７を必要な磁気モーメントを発生する物量で配置する例を示したが、既に言及したように、配置する鉄片物量１０として負の量が要求され、升目７の中に取り除ける鉄量が十分でないもしくは零となっている場合も発生する。前記のように対応しても十分な均一度にならない場合が発生する場合、鉄の代わりに永久磁石４Ｐもしくは電流ループ４Ｃを使う。これらは正の鉄物量に対しても適用しても問題は無いが、鉄の磁化で対応できる場合には、安価にシミングできる鉄片を用いることが望ましい。この様子は図１７の上の部分に示した。
【００９４】
図１８に示す電流ループでは、ループに流す電流による磁気モーメントＭｃ、原理式に基づき、ループ内の面積をＳｌとすると、
Ｍｃ＝電流×Ｓｌ
である。この磁気モーメントが必要とされる磁気モーメントと同じになるように電流を電源１０から、符号も考慮して調整する。
【００９５】
実際のシミングに用いる磁気モーメントの発生方法とその材料について、図１７に示した鉄や永久磁石以外に、ニッケルやコバルトなどの他の磁性体も考えられる。この場合には、永久磁石で述べたように磁化曲線を調べて、その換算磁気モーメントを求め、その上で必要な体積を求め、升目内に配置していく。
【００９６】
磁化曲線で磁化が求められない場合も存在する。例えばシミング用の磁性体が配置されると周囲の磁場（特に磁性体が周囲に有れば）を乱して、元の磁場とは異なる場合が発生するケースである。この場合には磁性体の磁化の程度を計測することが望ましい。例えば、実際に磁化が不明な磁性体片を配置する前後で、周囲の磁場を計測し、磁化が既知の場合の配置による磁場変化と比較する。もしくは計算上の磁場変化と比較することになる。また、磁化曲線が既知の最良のみで有れば、詳細な非線形磁場計算を行い、その計算結果から配置した鉄片の磁化４計算値として利用する。
【００９７】
実施例１〜４によれば、最終的な到達均一度を予測しながら、また磁石の品質を確認し、間違いを自動的に補正して、計測，補正鉄配置計算・配置を繰り返すことにより確実な磁場調整を行うことができる。また、高い磁場精度が要求される磁石の起磁力配置設計にも用いることができる。
【産業上の利用可能性】
【００９８】
本発明は、医療診断用に用いる核磁気共鳴断層写真装置（ＭＲＩ）などのように、コイルと鉄などの磁性体を配置して磁場を発生する磁石装置において、所望の磁場強度の分布に磁場を調整する方法と装置を提供する。特にＭＲＩ等のように核磁気共鳴応用装置では被計測領域において、極めて高精度に均一化する方法と装置を提供する。特に、誤差磁場を鉄片を配置することで補正し、磁場強度を均一化するシミングと言われる作業で、誤差磁場分布と鉄片配置分布をそれぞれの直行基底の組み合わせで均一磁場分布に補正する。
【符号の説明】
【００９９】
１  電流
１Ｂ  予備計算部分
２  磁化電流
２Ｂ  磁場計測部分
３  小コイル
３Ｂ  磁場調整計算部分
４  鉄片
４Ｃ  電流ループ
４Ｐ  永久磁石
５  シムトレイ
６  磁場計測評価領域（撮像領域）
７  升目
８  等高線のピーク
９  等高線の谷
１０  鉄片物量
１１  接点
１１Ｓ  磁場調整開始ステップ
１２  有限要素
１２Ｓ  磁場計測ステップ
１３  電流ポテンシャル評価面
１３Ｓ  計測磁場保存ステップ
１４  磁場計測評価点の集合
１４Ｓ  磁場データ読み出しステップ
１５  選択した固有モード
１５Ｓ  均一度判断ステップ
１６  非選択の固有モード
１６Ｓ  特異値分解結果読み出しステップ
１７  到達可能均一度
１７Ｓ  固有モード選択と目標磁場決定ステップ
１８  升目に配置する鉄体積
１８Ｓ  固有モード強度，補正電流ポテンシャル，鉄片量，補正磁場分布および到達可能均一度計算ステップ
１９  電流ポテンシャル等高線
１９Ｓ  スペクトル，到達可能均一度および鉄片配置量計算ステップ
２０Ｓ  シミング可否の判断ステップ
２１  電流ポテンシャルによる電流
２１Ｓ  品質良否判断ステップ
２２  固有モード選択の次数上限を示す線
２２Ｓ  鉄片配置作業ステップ
２３  固有モード選択で強さ下限を示す線
３１Ｓ  計算メッシュ生成ステップ
３２Ｓ  特異値分解計算ステップ
３３Ｓ  特異値分解結果保存ステップ
４０Ｓ  磁場調整終了ステップ
４１Ｓ  修理・調整ステップ
５１Ｓ  起磁力配置検討開始ステップ
５２Ｓ  起磁力配置仮定ステップ
５３Ｓ  磁場計算ステップ
５４Ｓ  磁場計算結果保存ステップ
５５Ｓ  シムトレイ変更要否判断ステップ
５６Ｓ  起磁力配置改善判断ステップ
５７Ｓ  ＭＲＩ磁石用起磁力配置候補案ステップ
６０  被検診者
６１  被検診者用ベッド
６２  磁石装置
６２ａ  コイル
６２ｂ  能動磁気シールド用コイル
６２ｃ  真空容器
６２ｄ  輻射シールド
６２ｅ  極低温容器
６２ｆ  磁石中空穴（ボア）
６３  連結柱
６４  電源
６５  静磁場磁力線方向
６６  傾斜磁場ベクトル
【図１】
0005427565.tif 000002
【図２】
0005427565.tif 000003
【図３】
0005427565.tif 000004
【図４】
0005427565.tif 000005
【図５】
0005427565.tif 000006
【図６】
0005427565.tif 000007
【図７】
0005427565.tif 000008
【図８】
0005427565.tif 000009
【図９】
0005427565.tif 000010
【図１０】
0005427565.tif 000011
【図１１】
0005427565.tif 000012
【図１２】
0005427565.tif 000013
【図１３】
0005427565.tif 000014
【図１４】
0005427565.tif 000015
【図１５】
0005427565.tif 000016
【図１６】
0005427565.tif 000017
【図１７】
0005427565.tif 000018
【図１８】
0005427565.tif 000019