JP7313284B2 - 冠血流予備量比シミュレーションパラメータのカスタマイズ、キャリブレーション、及び/又はトレーニング - Google Patents

冠血流予備量比シミュレーションパラメータのカスタマイズ、キャリブレーション、及び/又はトレーニング Download PDF

Info

Publication number
JP7313284B2
JP7313284B2 JP2019554672A JP2019554672A JP7313284B2 JP 7313284 B2 JP7313284 B2 JP 7313284B2 JP 2019554672 A JP2019554672 A JP 2019554672A JP 2019554672 A JP2019554672 A JP 2019554672A JP 7313284 B2 JP7313284 B2 JP 7313284B2
Authority
JP
Japan
Prior art keywords
coronary artery
user
artery tree
tree segmentation
adjusted
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
JP2019554672A
Other languages
English (en)
Other versions
JP2020516347A5 (ja
JP2020516347A (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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2020516347A publication Critical patent/JP2020516347A/ja
Publication of JP2020516347A5 publication Critical patent/JP2020516347A5/ja
Application granted granted Critical
Publication of JP7313284B2 publication Critical patent/JP7313284B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • 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
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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/0044Features 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 heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/14535Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring haematocrit
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/7475User input or interface means, e.g. keyboard, pointing device, joystick
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • 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
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/60ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
    • 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
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Cardiology (AREA)
  • Quality & Reliability (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)

Description

以下は、一般的にはイメージング、特に冠血流予備量比シミュレーション(FFR)のカスタマイズ、キャリブレーション、及び/又はトレーニングに関するものであり、FFR-コンピューター断層撮影(FFR-CT)への特定のアプリケーションで説明され、X線FFRにも対応している。
冠血流予備量比(FFR)は、FFRインデックスを介して、石灰化プラーク又は軟質プラークによる冠動脈病変の血行動態的重要性を定量化するためのカテーテル検査室(Cath Lab)の侵襲性測定である。インデックスは、冠動脈造影中に行われた圧力測定から計算された冠動脈狭窄の機能的重症度を示し、充血状態での近位圧力(口の近く)に対する遠位血圧(狭窄の背後)として定義される。つまり、FFRインデックスは、狭窄が存在しない状態での最大流量と比較して、狭窄が存在する状態での血管を流れる最大流量を表する。 FFR値は0乃至1の絶対値であり、値0.50は、特定の狭窄により血圧が50%低下することを示する。
侵襲性FFRプロシージャでは、カテーテルを大腿動脈又は放射状骨動脈に挿入し、カテーテルを狭窄部まで前進させる必要があり、カテーテルの先端のセンサーは、血管の形状、コンプライアンス及び抵抗、及び/又は他の特性に影響するさまざまな造影剤によって促進される状態の間、狭窄部の間の圧力を検出する。非侵襲性アプローチでは、冠動脈を通る血流と圧力をシミュレートする計算流体力学(CFD)シミュレーションを通じて、心臓のCT画像データ(たとえば、造影冠動脈CT血管造影、CCTA)からFFRインデックスを推定する。これには、CCTA画像データを使用して冠動脈ツリーの幾何モデルを導出し、そこから特徴を抽出し、シミュレーションの特徴から境界条件を決定することが含まれる。
FFR-CTシミュレーションソフトウェアは、CCTA画像から冠動脈をセグメント化するために、工場で調整されたプリセットアルゴリズムパラメータ、又は冠動脈を通じてシミュレートされた血流を制御する末梢抵抗のスケーリングファクターを備えた臨床現場で展開されている。しかしながら、使用パターン、個人の習慣、及び/又は患者集団は、サイトごとに異なる場合がある。たとえば、血管をセグメント化するとき、ある臨床医は血管壁の内面の内側でセグメント化することができ、別の臨床医は血管壁でセグメント化することができる。別の例では、血管壁の実際の周囲は、より低い画質の画像と比較してより高い画質の画像においてより明瞭であり得る。その結果、FFR-CTシミュレーションの結果の精度は、臨床医の間で、施設ごとに異なる場合がある。
本明細書で説明される態様は、上記で参照された問題及びその他に対処する。
一態様では、計算システムは、セグメンタ及び境界条件決定器を備えた生物物理シミュレータを含む、コンピュータ実行可能命令を備えたコンピュータ可読記憶媒体を含む。計算システムは、生物物理シミュレータを実行して、適応冠動脈ツリーセグメンテーション及び適応境界条件の少なくとも一方と心臓撮像データとを用いて冠血流予備量比指数を計算するように構成されたプロセッサをさらに含む。
別の態様では、コンピュータ可読記憶媒体は、コンピュータ可読命令でエンコードされ、計算システムのプロセッサによって実行されるとき、プロセッサに、心臓撮像データを受信させ、生物物理シミュレータを実行させて、適応した冠動脈ツリーセグメンテーション及び適応した境界条件の少なくとも1つ及び心臓撮像データを用いて冠血流予備量比インデックスを計算させる。
別の態様において、方法は、心臓撮像データを受信するステップと、適応した冠動脈ツリーセグメンテーション及び適応した境界条件の少なくとも1つ及び心臓撮像データを用いて冠血流予備量比インデックス計算するステップを含む。
当業者は、添付の説明を読んで理解すると、本出願のさらに他の態様を認識するであろう。
本発明は、様々な構成要素及び構成要素の配置、ならびに様々なステップ及びステップの配置の形をとることができる。図面は、好ましい実施形態を説明するためだけのものであり、本発明を限定するものとして解釈されるべきではない。
計算システム及び撮像システムを含むシステムを概略的に示している。 生物物理シミュレータの例を示している。 冠動脈のセグメント化の例を示している。 冠動脈の解剖学的モデルの例を示している。 本明細書の一実施形態による例示的な方法を示す図である。 本明細書の一実施形態による別の例示的な方法を示す。 本明細書の一実施形態による別の例示的な方法を示す。
図1は、CTスキャナなどの撮像システム102を含むシステム100を概略的に示している。変形形態では、システム100はX線撮像装置を含む。撮像システム102は、ほぼ静止したガントリ104と、静止したガントリ104によって回転可能に支持され、z軸を中心に検査領域108の周りを回転する回転ガントリ106とを含む。カウチなどの被験者支持部110は、検査領域108内の対象又は被験者を支持する。
X線管などの放射線源112は、回転ガントリ106によって回転可能に支持され、回転ガントリ106と共に回転し、検査領域108を横断する放射線を放出する。放射線感受性検出器アレイ114は、検査領域108の間で放射線源112の反対側の角度の弧を定める。放射線感受性検出器アレイ114は、検査領域108を横断する放射線を検出し、それを示す電気信号(投影データ)を生成する。再構成器116は、投影データを再構成し、検査領域108内に位置する被験者又は対象のスキャン部分を示す体積画像データを生成する。
この例では、計算システム118はオペレータコンソールとして機能する。コンソール118は、プロセッサ120(例えば、マイクロプロセッサ、中央処理装置など)と、一時的媒体を除外し、物理メモリデバイスなどの非一時的媒体を含むコンピュータ可読記憶媒体122とを含む。コンソール118は、モニタなどの人間が読み取り可能な出力装置、及びキーボード、マウスなどの入力装置をさらに含む。
コンピュータ可読記憶媒体122は、少なくとも生物物理シミュレータ126のための命令124を含む。プロセッサ120は、オペレータがグラフィカルユーザインターフェース(GUI)又はそれ以外を介してスキャナ102と対話し、及び/又はスキャナ102を操作できるようにする命令124及び/又はソフトウェアを実行するように構成される。プロセッサ120は、追加的又は代替的に、搬送波、信号及び/又は他の一時的な媒体によって運ばれるコンピュータ可読命令を実行してもよい。
変形例では、生物物理シミュレータ126は、コンソール118及びシステム100とは分離される他の計算システムの一部である。この例では、他の計算システムは、プロセッサ、コンピュータ読み取り可能な媒体を含むが、オペレータがスキャナ102と対話し、及び/又はスキャナ102を操作することを可能にするソフトウェアを含まないという点でコンソール118に類似している。
生物物理シミュレータ126は、少なくとも体積画像データを処理して生物物理シミュレーションを実行するように構成される。 FFRに関して、生物物理シミュレータは、FFRインデックスに基づくCCTA画像データ(又はX線血管造影図)を決定する。FFRインデックスは、ディスプレイモニタを介して表示され、保存され、別のデバイスに伝達されることができる。以下により詳細に説明するように、生物物理シミュレータ126は、セグメンテーションフィードバックループ及び/又は境界条件フィードバックループなどの1つ以上のフィードバックループを含む。フィードバックループは、個々の臨床医、施設などのトレーニング、シミュレーションのカスタマイズ及び/又はキャリブレーション、及び/又は予測目的に使用できる。
このように、生物物理シミュレータ126は臨床現場に配備され、例えば、CCTA画像データ(又はX線血管造影図)からの冠動脈ツリーのセグメント化、又は冠動脈を通る血流のシミュレートされた流れを制御する末梢抵抗のためのスケーリングファクタのために、個々のユーザ及び/又は現場に合わせられる。一例では、これにより、使用パターン、個々の習慣、及び/又はサイトごとに異なる患者集団のさまざまな結果が緩和される。結果として、臨床医の間、サイトの間からのFFR-CT(又はX線FFR)シミュレーション結果の精度は、本明細書に記載のアプローチが利用されない構成及び/又は本明細書で説明される生物物理シミュレータ126は省略される構成と比較して改善される可能性がある。
図2は、FFR-CTに関連した生物物理シミュレータ126の例を概略的に示している。この例では、生物物理シミュレータ126は、セグメンテーションフィードバックループ204を備えたセグメンタ202、境界条件フィードバックループ208を備えた境界条件決定器206、及びフローシミュレータ210を含む。変形例では、セグメンテーションフィードバックループ204は省略される。別の変形例では、境界条件フィードバックループ208は省略される。生物物理シミュレータ126は、入力として、撮像システム100、データリポジトリ(例えば、放射線情報システム(RIS)、写真及びアーカイブシステム(PACS)など)、及び/又は他の装置からCCTA撮像データを受け取る。
セグメンタ202は、セグメンテーションアルゴリズムを使用して、CCTA撮像データから冠動脈ツリーをセグメント化する。セグメンテーションは、自動的に(たとえば、機械学習など)又は半自動で(たとえば、ユーザ支援を受けて)実行できる。一例では、セグメンテーションは、そこから冠動脈の中心線及び/又は管腔形状(たとえば、直径、周囲、断面積など)を識別及び/又は抽出することを含む。セグメンテーションは、ボクセルの強度、オブジェクトの形状、及び/又はその他の特性に基づくことができる。図3は、血管管腔の対向する壁302を示す個々の血管の一部のセグメンテーションを示し、図4は、セグメント化された冠動脈ツリー400を示す。
CCTA撮像データから冠動脈ツリーを抽出するための適切なアプローチの例は、Zhengらの「モデル駆動型アプローチとデータ駆動型アプローチとの組み合わせによるCTAでの堅牢かつ正確な冠動脈中心線抽出」(Med Image Comput Assist Interv、 2013; 16(Pt 3):74-81)、Ecabertらの「モデルベースの適応フレームワークを使用したCT画像での心臓と大血管とのセグメンテーション」(Med Image Anal、 2011年 12月; 15(6):863-76)、及びFreimanらの「自動冠動脈内腔セグメンテーションにおける部分ボリュームモデリングを考慮することによるCCTAベースの病変の血行動態の重要性評価の改善」(Med Phys、 2017年3月; 44(3):1040-1049)に開示される。本明細書では他のアプローチも考えられる。
ある例では、CCTA撮像データはトレーニングセットであり、専門家はすでに既知の「グラウンドトゥルース」侵襲性FFR測定値と一致するFFRインデックスのシミュレーションをもたらす正確な冠動脈ツリーセグメンテーションを実行しており、トレーニングセットのための基準セグメンテーションが提供される。セグメンタ202は、トレーニングセットから冠動脈ツリーをセグメント化し、ユーザは、セグメンタ202のツールを使用して冠動脈ツリーセグメンテーションを自由に操作する。セグメンタ202は、ユーザ調整冠動脈ツリーセグメンテーションを基準セグメンテーションと比較する。セグメンタ202は、フィードバックループ204を介して、ユーザ調整された冠動脈ツリーセグメンテーションと基準セグメンテーションとの違いを示すフィードバックを提供する。フィードバックは、後続の調整された冠動脈のセグメンテーションが基準セグメンテーション、セグメンテーションを改善するための推奨、及び/又はセグメンテーションに関する他の情報に近いかを示してもよい。
別の例では、CCTA撮像データはトレーニングセットである。この場合、セグメンタ202は、ユーザ調整された冠動脈のセグメンテーションと基準セグメンテーションとの間の差を学習し、偏差をユーザ固有の較正データとして保存する。セグメンタ202は、例えば、評価又は分析中の患者CCTA撮像データ(非トレーニングデータセット)からの後続のユーザ調整冠動脈ツリーセグメンテーションを自動的に適応させるために、ユーザ固有の較正データを利用する。このために、セグメンタ202は、フィードバックループ204を介して、ユーザ固有の較正データに基づいて、適応された調整冠動脈のセグメンテーションを視覚的に表示する。適応のために、セグメンタ202は、セグメンテーションの1つ以上のポイントをシフトし得る。例えば、セグメンタ202は、セグメント化された冠状血管の直径を増加又は減少させ、セグメント化された冠状血管の壁を移動させることができる。
ユーザは、視覚的に表示され調整された冠動脈ツリーセグメンテーションを受け入れ、拒否及び/又は変更できる。ユーザはセグメンテーションをやり直すこともできる。ユーザは定期的にトレーニングセット(及び/又は別のトレーニングセット)を使用して、ユーザ固有のキャリブレーションデータを更新できる。たとえば、時間の経過とともに、基準に対してより正確な調整された冠動脈のセグメンテーションを作成するユーザの能力は、経験とともに向上する可能性がある。このような場合、現在のユーザ固有のキャリブレーションデータが過剰に修正される可能性がある。このため、ユーザはトレーニングセット(及び/又は他のトレーニングセット)を処理でき、セグメンタ202は、そこからの偏差に基づいてユーザ固有のキャリブレーションデータを更新できる。追加的又は代替的に、セグメンタ202は、起動時又は別の方法で、トレーニングセットをユーザに提示して、ユーザ固有の較正データを更新することができる。
境界条件決定器206は、ユーザ調整冠動脈ツリーセグメンテーション及び/又はセグメンタ適合ユーザ調整冠動脈ツリーセグメンテーション202から血管内の血流の計算流体力学シミュレーションの境界条件を決定する。 1つのアプローチでは、パラメトリック集中モデルが採用される。このモデルは非線形抵抗を使用する中心線表示を含み、要素は流入及び流出境界条件を示し、要素は、血管形状及び/又は油圧効果(たとえば、表1の1乃至6)を反映する一連の線形及び非線形抵抗要素を含むツリーセグメント伝達関数(直径、周囲、断面積など)を表す。
Figure 0007313284000001
集中モデルの例は、Nickisch他、「インタラクティブな冠動脈血流シミュレーションのための患者固有の集中モデルの学習」(医用画像計算及びコンピュータ支援介入-MICCAI 2015:18th International Conference、LNCS巻9350、2015、vol、 9350、pp 433-441)に開示される。境界条件を導出する例は、2014年6月30日に出願された「FFR-CTに対する患者の特定のモデリングの強化」と題されたEP14174891.3に記載されており、本明細書に参照により含まれる。本明細書では他のアプローチも考えられる。
図示の例では、FFR-CTシミュレーション後に収集された患者固有のデータが境界条件決定器206に提供される。そのようなデータの非限定的な例には、ヘマトクリット、糖尿病の有無、急性冠動脈症候群、血圧が含まれるが、これらに限定されない。データは、ユーザが提供したり、電子患者記録から抽出したり、その他の方法で取得したりできる。患者固有のデータは、フィードバックループ208を介して、境界条件を更新するために使用される。たとえば、以前にFFR-CTシミュレーションを受けた患者の侵襲性FFR測定を使用して、将来のシミュレーションを改善できる。さらに、文書化された結果、追加の生物物理学的及び/又は機能的測定を使用して、境界条件を再調整することができる。
以下は、患者固有のデータを利用する方法の非限定的な例を示している。ヘマトクリットは、赤血球の体積割合の尺度である。一般的に、ヘマトクリットの増加は、血液の粘度が高いことを意味する。血液の粘度は、経験曲線、たとえば40%から45%の間の線形スケーリングを使用して調整できる。糖尿病の存在は一般に、硬い壁と高い心筋抵抗を意味する。対応する抵抗境界条件は、たとえばデフォルト値に10%を追加することで調整できる。治療のカットオフFFR閾値が0.80である場合、急性冠症候群(ACS)の患者に対して、この閾値を例えば0.85に上げて、境界患者を確実に治療することができる。
特定の患者に対する追加の侵襲性基準FFR測定を使用して、(現在のセグメンテーションが与えられた場合)基準とモデル予測の間のより良い一致を得るように、何れかのモデルパラメータ、例えば心筋抵抗を適合させることができる。これは、勾配降下、パラメータ検索などによって行うことができる。心臓イベント又は生存データなどの文書化された患者の結果を同様に使用して、以前に取得したモデル予測と治療を考慮してパラメータを調整できる。たとえば、CT-FFR評価で病変が重要でないと見なされ、その病変が主要な心臓イベントを引き起こした場合、評価を再検討し、FFRしきい値などのパラメータを更新して、結果とよりよく一致させることができる。
さらに、又は代わりに、1つ又は複数の入力パラメータについて、境界条件決定器206は、入力パラメータの1つ又は複数の所定の変動の境界条件も決定する。例えば、Xの測定された血圧について、境界条件決定器206は、少なくとも(X)(±0.Z)の血圧値について境界条件を決定することもできる。ここで、Zは1%、2%、5%、10%、25%などの百分率としての許容範囲を表す。これにより、計算システム118は、入力パラメータ及び入力値周辺の値のFFR値を計算することができる。計算されたFFR値は、パラメータの値がFFR-CTシミュレーションの結果にどのように影響するかを示す。
パラメータ値の信頼性(たとえば、測定された血圧)が低く、値の小さな変化がFFR値に大きな影響を与える場合(たとえば、患者を治療するために患者を観察することにより推奨を変更する)、ユーザは、現在の値を使用するのではなく、より正確な値を取得して使用することを選択できる。代わりに、ユーザは、このパラメータ値がFFR値にどのように影響するかを理解することもできる。代わりに、このパラメータ値を使用せずに、及び/又はこのパラメータの影響を受ける境界条件を調整せずに続行することもできる。パラメータ値の信頼度が高い場合、ユーザは値を続行できる。
フローシミュレータ210は、境界条件でフローシミュレーションを実行し、FFR値を生成及び出力する。例えば、計算流体力学(CFD)アプローチ及び/又は他のアプローチを使用して、フローシミュレーションを実行できる。 FFR値の計算例は、2013年5月10日に出願された「血管の狭窄に対する冠血流予備量比(FFR)値の決定」と題されるUS 2015/0092999 A1、及び2013年10月24日に出願された「冠血流予備量比(FFR)インデックス」と題されるUS 2015/0282765 A1に開示されており、本明細書に参照として含まれる。
図5は、本明細書で説明される実施形態による例示的な方法を示す。
502では、CCTA撮像データのトレーニングセットから冠状血管がセグメント化される。
504で、セグメンテーションの調整を示すユーザ入力が受信され、セグメント化された血管に適用され、ユーザ調整された冠状血管セグメンテーションが生成される。
506で、ユーザ調整された冠状血管セグメンテーションは、トレーニングセットの基準血管セグメンテーションと比較される。
508で、セグメンテーション間の違いが視覚的に表示される。違いは、グラフィック及び/又は数値にすることができる。
この方法は、冠動脈ツリーのセグメンテーショントレーニングに適している。
図6は、本明細書で説明される実施形態による別の例示的な方法を示す。
602では、冠動脈ツリーがCCTA撮像データのトレーニングセットからセグメント化される。
604で、セグメンテーションの調整を示すユーザ入力が受信され、セグメント化されたツリーに適用され、ユーザ調整された冠動脈ツリーのセグメンテーションを生成する。
606で、ユーザ調整された冠動脈ツリーのセグメンテーションは、心臓撮像データのトレーニングセットについての基準ツリーセグメンテーションと比較される。
608で、それらの間の偏差がユーザ固有の較正として保存される。
610で、冠動脈ツリーが評価中の患者のCCTA撮像データからセグメント化される。
612で、セグメンテーションに対する調整を示すユーザ入力が受信され、セグメント化されたツリーに適用され、ユーザ調整された冠動脈ツリーセグメンテーションを生成する。
614で、ユーザ固有の較正がユーザ調整の冠状血管セグメンテーションに適用され、患者用に適合されたユーザ調整された冠動脈ツリーセグメンテーションが生成される。
616で、適合されたユーザ調整冠動脈ツリーセグメンテーションを使用してFFRインデックスが生成される。
図7は、本明細書で説明される実施形態による別の例示的な方法を示す。
702で、冠動脈ツリーが患者のCCTA撮像データからセグメント化される。一例では、冠動脈ツリーは、図6に関連して説明したように適合される。
704で、冠動脈ツリーに基づいて境界条件が決定される。
706で、境界条件を用いてFFR値が決定される。
708で、患者固有のデータ及び/又は他のデータに基づいて境界条件が更新される。
上記は、コンピュータ可読記憶媒体にエンコード又は埋め込まれたコンピュータ可読命令により実装され、コンピュータプロセッサによって実行されると、プロセッサに記述された行為を実行させる。それに加えて、又はその代わりに、コンピュータ可読命令の少なくとも1つは、コンピュータ可読記憶媒体ではない信号、搬送波又は他の一時的な媒体によって運ばれる。さらに、行為の順序は限定的ではないことを理解されたい。したがって、本明細書では他の順序付けが考えられる。加えて、1つ以上の行為が省略されてもよく、及び/又は1つ以上の追加の行為が含まれてもよい。
本発明は、図面及び前述の説明において詳細に図示及び説明されたが、そのような図示及び説明は、限定的ではなく例示的又は例示的であると見なされるべきである。本発明は開示された実施形態に限定されない。開示された実施形態に対する他の変形は、図面、開示、及び添付の特許請求の範囲の検討から、請求された発明を実施する際に当業者によって理解及び達成され得る。
請求項において、「含む」という語は他の要素又はステップを除外せず、不定冠詞「a」又は「an」は複数を除外しない。単一のプロセッサ又は他のユニットは、特許請求の範囲に記載されたいくつかのアイテムの機能を果たし得る。特定の測定値が相互に異なる従属請求項に記載されているという単なる事実は、これらの測定値の組み合わせが有利に使用できないことを示すものではありません。
コンピュータプログラムは、他のハードウェアとともに、又はその一部として提供される光学式記憶媒体又はソリッドステート媒体などの適切な媒体に格納/配布できるが、インターネットやその他の有線又は無線の通信システム。請求項中の参照符号は、範囲を限定するものとして解釈されるべきではない。

Claims (23)

  1. 計算システムであって、
    セグメンタ及び境界条件決定器を備える生物物理シミュレータを含む、コンピュータ実行可能命令を備えるコンピュータ可読記憶媒体と、
    前記生物物理シミュレータに、
    第1のユーザ調整冠動脈ツリーセグメンテーションと既知のグラウンドトゥルース侵襲性冠血流予備量比インデックス測定値に一致する冠血流予備量比インデックスをもたらす基準冠動脈ツリーセグメンテーションとの間の偏差を学習させ、
    前記第1のユーザ調整冠動脈ツリーセグメンテーションと前記基準冠動脈ツリーセグメンテーションとの間の前記学習される偏差に基づいて、非トレーニングデータセットからの第2のユーザ調整冠動脈ツリーセグメンテーションを適応させ、
    心臓撮像データと、前記適応冠動脈ツリーセグメンテーション及び適応境界条件の少なくとも一方とを用いて冠血流予備量比インデックスを計算させる
    ように構成される処理回路と
    を有する、
    計算システム。
  2. 前記心臓撮像データはトレーニング撮像を含み、セグメンタはユーザインタラクションで前記心臓撮像データから前記第1のユーザ調整冠動脈ツリーセグメンテーションを生成し、前記第1のユーザ調整冠動脈ツリーセグメンテーションと前記基準冠動脈ツリーセグメンテーションとの間の前記偏差を示す情報を視覚的に表示する、請求項1に記載のシステム。
  3. 前記心臓撮像データはトレーニング撮像を含み、セグメンタはユーザインタラクションで前記心臓撮像データから第1のユーザ調整冠動脈ツリーセグメンテーションを生成し、前記第1のユーザ調整冠動脈ツリーセグメンテーションと前記基準冠動脈ツリーセグメンテーションとの間の前記偏差を較正データとして記憶する、請求項1に記載のシステム。
  4. 前記セグメンタが、後続して患者の心臓画像データを処理し、ユーザインタラクションで前記心臓撮像データから第2のユーザ調整冠動脈ツリーセグメンテーションを生成する、請求項3に記載のシステム。
  5. 前記適応は、前記第2のユーザ調整冠動脈ツリーセグメンテーションの少なくとも1つのポイントをシフトするステップを含み、前記シフトするステップは、セグメント化される冠状血管の直径を変化させるステップを含む、請求項4に記載のシステム。
  6. 前記セグメンタは、確認又は拒否のために前記適応セグメンテーションを視覚的に表示する、請求項4に記載のシステム。
  7. 前記セグメンタは、前記トレーニング撮像又は他のトレーニング撮像の後続のセグメンテーションに基づいて前記較正データを更新する、請求項3に記載のシステム。
  8. 前記セグメンタは、患者の心臓撮像データを処理し、ユーザインタラクションで前記心臓撮像データから前記第1のユーザ調整冠動脈ツリーセグメンテーションを生成し、
    前記境界条件決定器は、前記適応セグメンテーション又は前記第1のユーザ調整冠動脈ツリーセグメンテーションに基づいて境界条件を決定し、前記適応セグメンテーション又は前記第1のユーザ調整冠動脈ツリーセグメンテーションに基づいて前記冠血流予備量比インデックスを計算するように構成されるフローシミュレータをさらに有する、
    請求項1に記載のシステム。
  9. 前記境界条件決定器は、患者固有のデータに基づいて前記境界条件を更新する、請求項8に記載のシステム。
  10. 前記境界条件決定器は、前記患者固有のデータにおける所定の変動のセットに対する境界条件のセットを決定し、前記フローシミュレータは、前記境界条件のセットに基づいて冠血流予備量比値を計算するように構成される、請求項9に記載のシステム。
  11. 前記境界条件決定器は、文書化される結果に基づいて前記境界条件を更新し、前記フローシミュレータは、前記境界条件のセットに基づいて冠血流予備量比値を計算するように構成される、請求項8に記載のシステム。
  12. 前記患者固有のデータはヘマトクリットを含み、前記境界条件決定器は、上昇するヘマトクリットに応じて血液境界条件の粘度の初期値を所定の範囲の間でスケーリングすることにより前記血液境界条件の粘度を増加させる、請求項8に記載のシステム。
  13. 前記患者固有のデータは糖尿病の存在を示し、前記境界条件決定器は、前記糖尿病の存在の表示に応じて抵抗境界条件の初期値に所定の割合を加えることにより、心筋壁の抵抗境界条件を増加させる、請求項8に記載のシステム。
  14. 前記患者固有のデータは急性冠動脈症候群の存在を示し、前記処理回路は、前記急性冠動脈症候群の存在の表示に応じて治療のためのカットオフ冠血流予備量比閾値を増加させる、請求項8に記載のシステム。
  15. コンピュータ可読命令でエンコードされるコンピュータ可読記憶媒体であって、処理回路によって実行されるとき、前記処理回路に、
    心臓撮像データを受信させ、
    第1のユーザ調整冠動脈ツリーセグメンテーションと既知のグラウンドトゥルース侵襲性冠血流予備量比インデックス測定値に一致する冠血流予備量比インデックスをもたらす基準冠動脈ツリーセグメンテーションとの間の偏差を学習させ、
    前記第1のユーザ調整冠動脈ツリーセグメンテーションと前記基準冠動脈ツリーセグメンテーションとの間の前記学習される偏差に基づいて、非トレーニングデータセットからの第2のユーザ調整冠動脈ツリーセグメンテーションを適応させ、
    前記適応冠動脈ツリーセグメンテーション及び適応境界条件の少なくとも1つ並びに前記心臓撮像データを用いて冠血流予備量比インデックスを計算する
    ように生物物理シミュレータを実行させる、
    コンピュータ可読記憶媒体。
  16. 前記心臓撮像データはトレーニング撮像を含み、前記処理回路はユーザインタラクションで前記心臓撮像データから前記第1のユーザ調整冠動脈ツリーセグメンテーションを生成し、前記第1のユーザ調整冠動脈ツリーセグメンテーションと前記基準冠動脈ツリーセグメンテーションとの間の前記偏差を示す情報を視覚的に表示するか、又は前記偏差を較正データとして記憶する、請求項15に記載のコンピュータ可読記憶媒体。
  17. 前記心臓撮像データは、評価中の患者の撮像データを含み、前記処理回路は、ユーザインタラクションで前記撮像データから前記第1のユーザ調整冠動脈ツリーセグメンテーションを生成する、請求項16に記載のコンピュータ可読記憶媒体。
  18. 前記受信心臓撮像データは患者の心臓撮像データを含み、前記処理回路は、ユーザインタラクションで前記心臓撮像データから前記第1のユーザ調整冠動脈ツリーセグメンテーションを生成し、前記適応セグメンテーション又は前記第1のユーザ調整冠動脈ツリーセグメンテーションに基づいて境界条件を決定し、前記境界条件に基づいて前記冠血流予備量比インデックスを計算し、前記冠血流予備量比インデックスの計算後に収集される患者固有のデータに基づいて前記境界条件を更新する、請求項15に記載のコンピュータ可読記憶媒体。
  19. 計算システムの作動方法であって、前記計算システムは処理回路を有し、
    前記処理回路が、心臓撮像データを受信するステップと、
    前記処理回路が、第1のユーザ調整冠動脈ツリーセグメンテーションと既知のグラウンドトゥルース侵襲性冠血流予備量比インデックス測定値に一致する冠血流予備量比インデックスをもたらす基準冠動脈ツリーセグメンテーションとの間の偏差を学習するステップと、
    前記処理回路が、前記第1のユーザ調整冠動脈ツリーセグメンテーションと前記基準冠動脈ツリーセグメンテーションとの間の前記学習される偏差に基づいて、非トレーニングデータセットからの第2のユーザ調整冠動脈ツリーセグメンテーションを適応するステップと、
    前記処理回路が、適応冠動脈ツリーセグメンテーション及び適応境界条件の少なくとも1つ及び前記心臓撮像データを用いて冠血流予備量比インデックス計算するステップと
    を有する、方法。
  20. 前記心臓撮像データはトレーニング撮像を含み、
    前記処理回路が、ユーザインタラクションで前記心臓撮像データから前記第1のユーザ調整冠動脈ツリーセグメンテーションを生成するステップと、
    前記処理回路が、前記第1のユーザ調整冠動脈ツリーセグメンテーションと前記基準冠動脈ツリーセグメンテーションとの間の前記偏差を示す情報を視覚的に表示するステップ、又は前記偏差を較正データとして記憶するステップのうちの少なくとも一つと
    を更に有する、請求項19に記載の方法。
  21. 前記心臓撮像データは、評価中の患者の撮像データを含み、
    前記処理回路が、ユーザインタラクションで前記撮像データから前記第1のユーザ調整冠動脈ツリーセグメンテーションを生成するステップ
    を更に有する、請求項20に記載の方法。
  22. 前記受信心臓撮像データは、患者の心臓撮像データを含み、
    前記処理回路が、ユーザインタラクションで前記心臓撮像データから前記第1のユーザ調整冠動脈ツリーセグメンテーションを生成するステップと、
    前記処理回路が、前記適応セグメンテーション又は前記第1のユーザ調整冠動脈ツリーセグメンテーションに基づいて境界条件を決定するステップと、
    前記処理回路が、前記境界条件に基づいて前記冠血流予備量比インデックスを計算するステップと、
    前記処理回路が、前記冠血流予備量比インデックスの計算後に収集される患者固有のデータに基づいて前記境界条件を更新するステップと
    を更に有する、請求項19に記載の方法。
  23. 前記計算される冠血流予備量比インデックスは、既知のグラウンドトゥルース侵襲性冠血流予備量比インデックス測定値に一致する、請求項1に記載のシステム。
JP2019554672A 2017-04-06 2018-04-06 冠血流予備量比シミュレーションパラメータのカスタマイズ、キャリブレーション、及び/又はトレーニング Active JP7313284B2 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201762482231P 2017-04-06 2017-04-06
US62/482,231 2017-04-06
US201762535264P 2017-07-21 2017-07-21
US62/535,264 2017-07-21
PCT/EP2018/058885 WO2018185298A1 (en) 2017-04-06 2018-04-06 Fractional flow reserve simulation parameter customization, calibration and/or training

Publications (3)

Publication Number Publication Date
JP2020516347A JP2020516347A (ja) 2020-06-11
JP2020516347A5 JP2020516347A5 (ja) 2021-05-13
JP7313284B2 true JP7313284B2 (ja) 2023-07-24

Family

ID=61952704

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019554672A Active JP7313284B2 (ja) 2017-04-06 2018-04-06 冠血流予備量比シミュレーションパラメータのカスタマイズ、キャリブレーション、及び/又はトレーニング

Country Status (5)

Country Link
US (1) US11195278B2 (ja)
EP (1) EP3607558A1 (ja)
JP (1) JP7313284B2 (ja)
CN (1) CN110574121B (ja)
WO (1) WO2018185298A1 (ja)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7036742B2 (ja) 2016-05-16 2022-03-15 キャスワークス リミテッド 血管評価システム
EP4241694A3 (en) 2016-05-16 2023-12-20 Cathworks Ltd. Selection of vascular paths from images
EP3602569A1 (en) 2017-03-31 2020-02-05 Koninklijke Philips N.V. Simulation of transcatheter aortic valve implantation (tavi) induced effects on coronary flow and pressure
US20220015730A1 (en) * 2018-11-28 2022-01-20 Koninklijke Philips N.V. Most relevant x-ray image selection for hemodynamic simulation
US12039685B2 (en) 2019-09-23 2024-07-16 Cathworks Ltd. Methods, apparatus, and system for synchronization between a three-dimensional vascular model and an imaging device
US11386563B2 (en) * 2019-11-22 2022-07-12 The Regents Of The University Of Michigan Anatomical and functional assessment of coronary artery disease using machine learning
DE102020200750A1 (de) * 2020-01-22 2021-07-22 Siemens Healthcare Gmbh Bereitstellen eines Blutflussparametersatzes einer Gefäßmalformation
CN111754522B (zh) * 2020-06-19 2021-08-03 上海杏脉信息科技有限公司 冠状动脉血流动力学数据的获取方法及装置
CN112711831B (zh) * 2020-12-07 2022-10-28 上海联影医疗科技股份有限公司 血管仿真分析方法、设备、装置、计算机设备和存储介质
CN113648059B (zh) * 2021-08-26 2023-09-29 上海联影医疗科技股份有限公司 手术规划评估方法、计算机设备和存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007511267A (ja) 2003-11-13 2007-05-10 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 変形可能表面を使った三次元セグメント化
US20160133015A1 (en) 2010-08-12 2016-05-12 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015509820A (ja) * 2012-03-15 2015-04-02 シーメンス アクティエンゲゼルシャフト 安静および充血時における冠血流のパーソナライゼーション用フレームワーク
CN104321009B (zh) 2012-05-14 2017-06-13 皇家飞利浦有限公司 针对血管狭窄的血流储备分数(ffr)值的确定
US9757073B2 (en) 2012-11-06 2017-09-12 Koninklijke Philips N.V. Fractional flow reserve (FFR) index
US20140236547A1 (en) 2013-02-15 2014-08-21 Siemens Aktiengesellschaft Patient-specific automated tuning of boundary conditions for distal vessel tree
US8824752B1 (en) * 2013-03-15 2014-09-02 Heartflow, Inc. Methods and systems for assessing image quality in modeling of patient anatomic or blood flow characteristics
US10595806B2 (en) 2013-10-22 2020-03-24 Koninklijke Philips N.V. Fractional flow reserve (FFR) index with adaptive boundary condition parameters
US9390232B2 (en) 2014-03-24 2016-07-12 Heartflow, Inc. Systems and methods for modeling changes in patient-specific blood vessel geometry and boundary conditions
US9087147B1 (en) * 2014-03-31 2015-07-21 Heartflow, Inc. Systems and methods for determining blood flow characteristics using flow ratio
US9754082B2 (en) * 2014-05-30 2017-09-05 Heartflow, Inc. Systems and methods for reporting blood flow characteristics
EP3160335B1 (en) 2014-06-30 2023-09-20 Koninklijke Philips N.V. Apparatus for determining a fractional flow reserve value
US9349178B1 (en) * 2014-11-24 2016-05-24 Siemens Aktiengesellschaft Synthetic data-driven hemodynamic determination in medical imaging
WO2016087396A1 (en) 2014-12-02 2016-06-09 Koninklijke Philips N.V. Fractional flow reserve determination
US9785748B2 (en) 2015-07-14 2017-10-10 Heartflow, Inc. Systems and methods for estimating hemodynamic forces acting on plaque and monitoring patient risk
CN112840408A (zh) * 2018-10-17 2021-05-25 哈特弗罗公司 用于根据脂肪组织来评估心血管疾病和治疗效果的***和方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007511267A (ja) 2003-11-13 2007-05-10 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 変形可能表面を使った三次元セグメント化
US20160133015A1 (en) 2010-08-12 2016-05-12 Heartflow, Inc. Method and system for image processing to determine patient-specific blood flow characteristics
JP2016104328A (ja) 2010-08-12 2016-06-09 ハートフロー, インコーポレイテッド 患者固有の血流のモデリングのための方法およびシステム

Also Published As

Publication number Publication date
US11195278B2 (en) 2021-12-07
WO2018185298A1 (en) 2018-10-11
CN110574121B (zh) 2024-03-19
EP3607558A1 (en) 2020-02-12
US20210110543A1 (en) 2021-04-15
CN110574121A (zh) 2019-12-13
JP2020516347A (ja) 2020-06-11

Similar Documents

Publication Publication Date Title
JP7313284B2 (ja) 冠血流予備量比シミュレーションパラメータのカスタマイズ、キャリブレーション、及び/又はトレーニング
US10595806B2 (en) Fractional flow reserve (FFR) index with adaptive boundary condition parameters
JP6297699B2 (ja) 心臓データを処理する処理装置、そのような処理装置を有するイメージングシステム、処理装置の作動方法、イメージングシステムの作動方法及びコンピュータプログラム
JP6484760B2 (ja) 非侵襲的血流予備量比(ffr)に対する側副血流モデル化
JP2017518844A (ja) 血流予備量比値を算出するための機器
CN109996495B (zh) 用于生物物理模拟的血管树标准化和/或用于经修剪的部分的扩展模拟
JP6749917B2 (ja) iFR−CT
JP2016520357A (ja) 狭窄治療計画
JP7300392B2 (ja) 標準化された冠動脈疾患メトリック
US20190076196A1 (en) Vessel geometry and additional boundary conditions for hemodynamic ffr/ifr simulations from intravascular imaging
JP7426824B2 (ja) 非侵襲的イメージングベースのffrのインタラクションモニタリング
CN111033635B (zh) 基于模型和成像数据的冠状动脉健康状态预测
JP2022510879A (ja) 血行力学的シミュレーションのための最も関連のあるx線画像の選択
JPWO2018185040A5 (ja)
EP4197429A1 (en) System and method for determining partial mass of an organ

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210402

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210402

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220322

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220324

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20220607

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220921

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20221020

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221208

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230330

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230608

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230711

R150 Certificate of patent or registration of utility model

Ref document number: 7313284

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150