磁場が陽子線に与える影響を正確に予測:MRiPTのための治療計画システムを検証

· 論文解説

🎧 この記事を音声で聴く(AI生成Podcast)

解説対象論文: Dosimetric validation of proton beam deflection and field shape deformation in a measured transverse magnetic field of an in-beam MR scanner using a novel dose calculation functionality in a research version of a commercial treatment planning system
Physics in Medicine and Biology|被引用数: 0(2026年9月15日時点)

磁気共鳴イメージング(MRI)と陽子線治療(PT)を組み合わせたMR統合陽子線治療(MRiPT)は、高精度ながん治療の未来を切り開く可能性を秘めています。しかし、MRIが発する磁場(MF)は、陽子線の軌道に影響を与え、線量分布を変化させるという課題があります。本研究は、この磁場下での陽子線挙動を商業治療計画システム(TPS)であるRayStationが正確に予測できるかを、実測された3D磁場マップを用いて包括的に検証しました。その結果、TPSが磁場による陽子線偏向と照射野変形を高精度で再現できることが示され、MRiPTの臨床実装に向けた重要な基盤が確立されました。

研究サマリー(FINER & PICO)
Feasible実施可能性 商業TPSの研究版と実測磁場マップを用いて、MRiPTの線量計算機能の実装と検証に成功した。
Ethical倫理面 臨床実装を目指すMRiPTの線量計算精度向上は、対象者の治療安全性に貢献するため倫理的である。
Interesting面白さ 磁場が陽子線に与える複雑な影響を、詳細な実測データと高度な計算モデルで解明している点が興味深い。
Novel新規性 実際のMRスキャナーで測定された3D磁場マップを商業TPSに実装し、線量計算を実験的に検証した初の試み。
Relevant切実さ MRiPTの臨床実装における不可欠な課題である、磁場下での陽子線治療計画の基盤を確立する。
Measurable測定可能性 複数のエネルギーと位置で線量分布を測定し、ガンマ解析により定量的に線量計算の精度を検証した。
Modifiable派生・改善 研究版TPSに新機能を実装しており、実測データに基づく検証結果を反映して機能改善が可能である。
Structured文章構造 空中および水等価ファントム内の単純な照射パターンから臨床的な治療野まで、段階的な検証戦略を採用。
PICO / PECO対象・介入・比較・結果 陽子線治療において、MRイメージングの磁場が存在する場合の陽子線輸送と線量分布の正確な計算能力を検証する。

MRiPT時代の幕開け:陽子線治療の精度を高める新たな挑戦

どうも、Beyond the Pixelです。がん治療において、陽子線治療はその高い精度で注目されています。しかし、その恩恵を最大限に引き出すためには、治療中の対象者の体の動きや位置の変化をリアルタイムで把握し、計画に反映させる「画像誘導」が不可欠です。これまではX線ベースの画像診断が主流でしたが、軟部組織のコントラストが低いことや、追加のX線被ばくがあるという課題がありました。そこで登場したのが、優れた軟部組織コントラストと非電離放射線という利点を持つ「磁気共鳴イメージング(MRI)」を陽子線治療装置に統合する「MR統合陽子線治療(MRiPT)」という画期的なコンセプトです。

Figure 1
Figure 1. Figure 1. Schematic of the beamline and the MRiPT prototype system. 3D view (top) and orthogonal解説: MRiPTプロトタイプシステムの全体像と陽子線ビームライン、MRアイソセンター(MRisoc)およびビームアイソセンター(PBSisoc)の位置関係を示しています。3Dビュー(上)、前面図(中)、上面図(下)で構成され、静磁場(B0)の向きも示されています。

しかし、MRIを陽子線治療に統合することは、磁場(MF)が陽子線(陽電荷を帯びた粒子)の軌道に影響を与えるという特有の課題を伴います。この研究では、この磁場による陽子線の偏向や照射野の変形といった影響を、治療計画システム(TPS)が正確にモデル化できるかを実験的に検証しました。これにより、MRiPTの臨床応用に向けた重要な一歩が確立されたと言えるでしょう。

精密な磁場マップの作成とTPSへの統合

本研究の最初のステップは、MRiPTプロトタイプシステムに組み込まれた0.32テスラ(T)のMRIスキャナーが生成する3D磁場を、高解像度で正確にマッピングすることでした。磁場はスキャナーの中心だけでなく、その周辺の「フリンジフィールド」と呼ばれる不均一な領域まで詳細に測定する必要があります。これは、建物の構造材や陽子線治療装置自体が磁場に影響を与える可能性があるため、メーカー提供のシミュレーションマップだけでは不十分だからです。研究チームは、ハンドヘルド型の3軸ホールプローブ磁力計と赤外線光学トラッカーを組み合わせたMR互換のマッピングフレームワークを使用し、約1時間で161 × 41 × 25 cm³という広範囲の3D磁場マップを取得しました。

Figure 2
Figure 2. Figure 2. Schematic of the 3D MF mapping setup within the in-beam MRiPT system. An optical解説: インビームMRiPTシステム内での3D磁場マッピング設定の模式図です。光学トラッキングシステムとホール磁力計が使用され、測定ボリュームが示されています。これは、陽子線PBSノズルのビーム出口窓からMRスキャナーの遠位端まで広がる範囲です。

この測定された磁場マップは、その後、商業用TPSであるRayStationの研究バージョンに実装されました。この際、TPSの座標系に正確に変換し、磁場の方向や強度が正しく反映されていることを確認するための入念な検証が行われました。特に、陽子線の軌道がローレンツ力によってどのように偏向するかという物理的な正確性も検証されました。

Figure 4
Figure 4. Figure 4. High-resolution 3D magnetic flux density map of the in-beam MR scanner. (a) Measured 3D magnetic vector field distribution with zoomed detail showing vector orientation. (b) 3D representation of the magnitude of the measured magnetic flux density map (|B0|) with coordinate axes indicated. (c) Representative central front and top view slices of the magnetic flux density magnitude, demonstrating homogeneous field distribution in the imaging volume and strong gradients toward the fringe fields. The colour scale is identical to that shown in (b).解説: インビームMRスキャナーの高解像度3D磁束密度マップを示しています。(a) 測定された3D磁気ベクトル場分布で、ベクトルの向きを拡大表示しています。(b) 測定された磁束密度マップの大きさを3Dで表現し、座標軸を示しています。(c) 磁束密度の大きさの代表的な中央前面および上面スライスで、画像化ボリューム内の均一な磁場分布と、フリンジフィールドに向かう強い勾配を示しています。
Figure 5
Figure 5. Figure 5. Comparison between measured and implemented 3D MF map components in RayStation. Front and top views of the three components (BX, BY, BZ) are shown. The implemented maps are named according to the DICOM convention of RayStation and reproduce the correct orientation relative to the treatment plan geometry.解説: RayStationに実装された3D磁場マップの各コンポーネントと測定された磁場マップとの比較を示しています。BX、BY、BZの3つのコンポーネントについて、前面図と上面図が示されており、実装されたマップがRayStationのDICOM規約に従って正しい向きと大きさで再現されていることが確認できます。

RayStationの陽子線モンテカルロ線量計算エンジンは、この磁場の影響を考慮するために拡張されました。陽子線が磁場を通過する際、そのエネルギーと磁場強度に応じて軌道が変更されます。シミュレーションでは、陽子線のステップ長を1mm以下に制限し、各ステップでの偏向角度が5.7度以下になるように設定することで、高精度な計算を実現しています。この新しい計算機能が、実際の磁場環境下でどれだけ正確な線量分布を予測できるかが、次の重要な検証ポイントとなりました。

Figure 6
Figure 6. Figure 6. Illustration of the effect of the MF map of the in-beam MR scanner on proton dose distribution. (a) Schematic top view of the MRiPT system with the mapped MF volume indicated. Green lines represent the expected beam divergence in the absence of the MF, while orange lines illustrate the deflected beam path in the presence of the MF (left). The trapezoidal deformation of the radiation field at the maximum dose depth is shown in beam’s eye view (BEV) (right). (b) 2D top- view heat map of the MF, visualised in RayStation and expressed in mT. The offset between PBSisoc and MRisoc is verified by visualising the central point of the implemented map. (c) Top view of calculated beam trajectory in RayStation for a 100 MeV square radiation field (shown in green for the case without MF), deflected by the transverse MF (left). Corresponding 2D dose distribution in BEV at the Bragg peak depth in a water phantom (right), demonstrating the trapezoidal deformation induced by the MF.解説: インビームMRスキャナーの磁場マップが陽子線線量分布に与える影響を示しています。(a) MRiPTシステムの上面模式図で、磁場がない場合の予想されるビーム発散(緑線)と、磁場が存在する場合の偏向したビーム経路(オレンジ線)を示しています。(b) RayStationで可視化された磁場マップの2D上面ヒートマップです。(c) 100MeVの正方形放射線野に対するRayStationでの計算されたビーム軌道(磁場がない場合は緑色で表示)と、ブラッグピーク深さでの2D線量分布が示され、磁場による台形変形が確認できます。

実験による線量計算の定量的な検証

線量計算の精度を定量的に評価するため、研究チームは段階的な検証戦略を採用しました。まず、単純な照射パターン(スポットパターンと正方形の放射線野)を用いて、空気中での陽子線輸送における磁場の影響を評価しました。陽子線のエネルギーは100、150、220 MeVの3種類で、2D電離箱アレイを用いて複数の位置で線量分布を測定しました。測定結果は、磁場が陽子線のエネルギーや位置に応じて横方向への偏向を引き起こし、放射線野の形状が変形することを示しました。特に、陽子線のエネルギーが低いほど、また磁場内を通過する距離が長いほど、偏向が顕著になることが確認されました。例えば、PBSアイソセンター(PBSisoc)から最も離れた位置Dでは、100 MeVの陽子線で182.7 ±3 mmという大きな横方向偏向が測定されました。

Table 1
Table 1. Table 1. Measured lateral displacement of the central proton dose spots from experiments in-air at four positions (A–D) relative to the PBS isocentre (PBSisoc), for beam energies of 100, 150, and 220 MeV.解説: 空気中での実験で測定された、PBSisocからの4つの位置(A-D)における中心陽子線スポットの横方向変位を、100、150、220 MeVのビームエネルギー別に示しています。

これらの測定値とRayStationによる計算結果を比較したところ、横方向のビームシフトは概ね1~3mmの範囲内で一致しており、実験誤差の範囲内であることが示されました。特に重要なのは、MRIのアイソセンター(MRisoc)周辺の臨床的に重要な位置においても、計算が測定値を高い精度で再現したことです。

Figure 7
Figure 7. Figure 7. Comparison of in-air measured and calculated horizontal beam deflections. Data points with error bars (±3 mm) represent OCTAVIUS 1500XDR measurements and positioning uncertainties, while solid lines show beam trajectories calculated by RayStation for proton energies of 100, 150, and 220 MeV. The upper panel depicts the full measurement range (positions A–D). The lower panel provides a magnified view around the MR isocentre (position B) to assess measurement-to- calculations agreement at this critical site, where the patient would be located in an MRiPT treatment解説: 空気中での測定された水平方向のビーム偏向と計算されたビーム偏向の比較を示しています。データ点とエラーバー(±3mm)はOCTAVIUS 1500XDRの測定値と位置決め不確かさを表し、実線は100、150、220MeVの陽子線エネルギーに対するRayStationによる計算ビーム軌道を示しています。

さらに、正方形の放射線野の形状変形についても、ガンマ解析(𝛾-analysis)と呼ばれる評価指標を用いて、測定と計算の間で高い一致度を示しました。臨床的に許容される基準である3%/3mm(線量閾値10%)では、すべてのエネルギーと位置で97%を超えるパス率を達成し、線量計算の性能が臨床的に許容できるレベルであることを確認しました。

Table 2
Table 2. Table 2. 𝛾-pass rates (2%/2 mm, 5% dose threshold) for measured and calculated 2D dose distributions at all the detector positions (A–D) and proton beam energies of 100, 150, and 220 MeV, reported with and without auto-alignment. Auto-alignment values (in mm) are listed to indicate the magnitude of corrections applied to centre the measured and calculated dose distributions.解説: 検出器の全位置(A-D)と100、150、220 MeVの陽子線ビームエネルギーにおける、測定および計算された2D線量分布のガンマパス率(2%/2mm、5%線量閾値)を示しています。オートアライメントの有無で報告されており、オートアライメントの値(mm)も記載されています。
Figure 8
Figure 8. Figure 8. Top: 2 %/2 mm 𝛾-distributions illustrating the agreement between in-air measured and calculated 2D dose distributions in the presence of the MF at position B for proton energies of 100, 150, and 220 MeV, shown both with auto-alignment (left). Central horizontal dose profile for 220 MeV, allowing direct comparison of measurements and calculations with auto-alignment (right). Coloured rhomboid markers indicate the measured dose values at individual detector channels, with red and blue denoting the left and right penumbra points, respectively, light blue marking central‐region detector readings, and green representing intermediate off-axis points along the field解説: 空気中での測定と計算された2D線量分布の2%/2mmガンマ分布を示しています。上段は位置Bでの100、150、220MeVの陽子線エネルギーに対する結果で、オートアライメントあり(左)と220MeVの中央水平線量プロファイル(右)が含まれます。下段は100MeVの最低エネルギーでの位置A~Dすべての測定位置に対するガンマ分布を示しています。

臨床的条件での検証と今後の展望

この研究では、単純な空気中での検証だけでなく、水等価ファントム内での陽子線輸送を伴う臨床的に関連性の高い食道がんの治療計画も評価しました。複数の水等価深度(WED)で線量分布を測定し、RayStationの計算結果と比較しました。その結果、浅い深度から深い深度まで、すべてのWEDで測定値と計算値の間に良好な一致が確認されました。ガンマ解析のパス率は、自動アラインメントなしで95%以上を維持し、これはTPSが磁場下での陽子線のエネルギー損失と線量沈着を正確にモデル化できることを示しています。

Table 3
Table 3. Table 3. 𝛾-pass rates (2%/2 mm, 5% dose threshold) for measured and calculated 2D dose distributions at all water equivalent depths (WED) for both fields of the oesophagus treatment plan, reported with and without auto-alignment. Auto-alignment values (in mm) are listed to indicate the magnitude of corrections applied to centre the measured and calculated dose distributions.解説: 食道治療計画の両フィールドにおける、すべての水等価深度(WED)での測定および計算された2D線量分布のガンマパス率(2%/2mm、5%線量閾値)を示しています。オートアライメントの有無で報告されており、オートアライメントの値(mm)も記載されています。
Figure 9
Figure 9. Figure 9. Right: 2 %/2 mm 𝛾-distributions illustrating the agreement between in-phantom measured and calculated 2D dose distributions in the presence of the MF of both the oesophagus radiation treatment fields (Top: Field A, Bottom: Field B) Left: Central horizontal dose profile for the 17.46 water-equivalent depth (WED), allowing direct comparison of measurements and calculations with auto-alignment.解説: ファントム内測定における、食道放射線治療分野(上:Field A、下:Field B)の磁場存在下での測定および計算された2D線量分布の一致を示す2%/2mmガンマ分布と、17.46cm水等価深度(WED)での中央水平線量プロファイルを示しています。

本研究は、実際のインビームMRスキャナーで測定された高解像度3D磁場マップを、商業用TPS(RayStation)の研究バージョンに初めて実装し、その陽子線線量計算機能を検証した点で画期的な成果です。これまでの研究は、シミュレーションツールや自社開発のTPSに限られていたり、理想化された磁場条件を仮定したりすることが多かったため、本研究はMRiPTの臨床実装に向けた重要な基盤を築くものです。このワークフローは、磁場が存在する陽子線治療計画の基礎を確立し、MRiPTの臨床応用への貢献が期待されます。今後の展望としては、デルタ電子の生成と輸送など、現在のRayStation MC陽子線線量計算アルゴリズムでは考慮されていない詳細な物理プロセスを組み込むことで、さらなる精度向上も考えられるでしょう。

用語集

  • 陽子線治療 (PT): 陽子線を用いた放射線治療の一種で、がん組織に高い線量を集中させ、周囲の健康な組織への影響を最小限に抑える治療法。
  • 磁気共鳴イメージング (MRI): 磁場と電波を用いて体内の詳細な画像を生成する医療画像診断法。優れた軟部組織コントラストが特徴。
  • 磁場 (MF): 磁力が発生している空間。MRI装置は強力な磁場を生成し、これが陽子線の軌道に影響を与える。
  • MR統合陽子線治療 (MRiPT): 陽子線治療装置にMRI機能を統合した治療システム。治療中にリアルタイムで画像誘導を行うことを目指す。
  • 治療計画システム (TPS): 放射線治療の計画を作成するためのソフトウェア。どの部位にどれくらいの線量を照射するかなどを計算する。
  • モンテカルロ法 (MC): ランダムな数値シミュレーションを繰り返すことで、陽子線が物質中でどのように振る舞うかを予測する計算手法。
  • ペンシルビームスキャニング (PBS): 細い陽子線ビーム(ペンシルビーム)をスキャンすることで、がん組織の形状に合わせて精密に線量を照射する技術。
  • 線量測定 (Dosimetric validation): 治療計画で計算された放射線量が、実際に装置から出力される量と一致するかを実験的に検証すること。
  • 水等価深度 (WED): 放射線が水中でどれだけの深さに到達するかを示す指標。生体組織の密度差を水と同じ効果に換算する。
  • ガンマ解析 (𝛾-analysis): 二つの線量分布(例:測定値と計算値)の差異を定量的に評価するための手法。位置と線量の両方のずれを考慮する。
  • フリンジフィールド (Fringe field): MRIスキャナーの中心部の均一な磁場から離れた、磁場の強度が徐々に弱くなり、不均一になる周辺領域。
powered by visionary imaging services, inc.
免責事項

発行 Beyond the Pixel 編集部(Visionary Imaging Services, Inc.)

編集責任者 小林 達明(診療放射線技師/博士(情報学)) ORCID

本記事は生成 AI を用いて作成しています。公開後に気づいた誤りは編集部が随時修正します。

論文の引用・解釈の方針は免責事項をご覧ください。

コメントを残す

メールアドレスが公開されることはありません。 ※ が付いている欄は必須項目です