JP2019202154A - Medical image processing apparatus, medical image processing system, and medical image processing method - Google Patents
Medical image processing apparatus, medical image processing system, and medical image processing method Download PDFInfo
- Publication number
- JP2019202154A JP2019202154A JP2019119706A JP2019119706A JP2019202154A JP 2019202154 A JP2019202154 A JP 2019202154A JP 2019119706 A JP2019119706 A JP 2019119706A JP 2019119706 A JP2019119706 A JP 2019119706A JP 2019202154 A JP2019202154 A JP 2019202154A
- Authority
- JP
- Japan
- Prior art keywords
- unit
- responsible
- stenosis
- blood vessel
- region
- 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.)
- Granted
Links
- 238000012545 processing Methods 0.000 title claims abstract description 55
- 238000003672 processing method Methods 0.000 title description 2
- 208000031481 Pathologic Constriction Diseases 0.000 claims abstract description 93
- 230000036262 stenosis Effects 0.000 claims abstract description 93
- 208000037804 stenosis Diseases 0.000 claims abstract description 93
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 92
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 90
- 230000000302 ischemic effect Effects 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 208000010125 myocardial infarction Diseases 0.000 claims abstract description 24
- 230000002966 stenotic effect Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 description 63
- 239000003550 marker Substances 0.000 description 63
- 230000002107 myocardial effect Effects 0.000 description 49
- 238000002591 computed tomography Methods 0.000 description 35
- 230000017531 blood circulation Effects 0.000 description 30
- 238000011282 treatment Methods 0.000 description 28
- 238000010586 diagram Methods 0.000 description 26
- 206010061216 Infarction Diseases 0.000 description 23
- 230000007574 infarction Effects 0.000 description 23
- 238000000034 method Methods 0.000 description 21
- 239000000284 extract Substances 0.000 description 18
- 210000004165 myocardium Anatomy 0.000 description 18
- 238000000605 extraction Methods 0.000 description 15
- 238000004891 communication Methods 0.000 description 13
- 238000004088 simulation Methods 0.000 description 12
- 239000002872 contrast media Substances 0.000 description 11
- 230000010412 perfusion Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 238000011144 upstream manufacturing Methods 0.000 description 8
- 238000002595 magnetic resonance imaging Methods 0.000 description 7
- 210000003484 anatomy Anatomy 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000009206 nuclear medicine Methods 0.000 description 5
- 238000001356 surgical procedure Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 3
- 239000008280 blood Substances 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000002651 drug therapy Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 208000031225 myocardial ischemia Diseases 0.000 description 3
- 235000016709 nutrition Nutrition 0.000 description 3
- 230000035764 nutrition Effects 0.000 description 3
- 208000024891 symptom Diseases 0.000 description 3
- 230000002308 calcification Effects 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 230000001364 causal effect Effects 0.000 description 2
- HVYWMOMLDIMFJA-DPAQBDIFSA-N cholesterol Chemical compound C1C=C2C[C@@H](O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2 HVYWMOMLDIMFJA-DPAQBDIFSA-N 0.000 description 2
- 239000011248 coating agent Substances 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 150000002632 lipids Chemical class 0.000 description 2
- 230000004060 metabolic process Effects 0.000 description 2
- 230000001338 necrotic effect Effects 0.000 description 2
- 235000015097 nutrients Nutrition 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 238000012502 risk assessment Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 229910052688 Gadolinium Inorganic materials 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 206010057469 Vascular stenosis Diseases 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- UIWYJDYFSGRHKR-UHFFFAOYSA-N gadolinium atom Chemical compound [Gd] UIWYJDYFSGRHKR-UHFFFAOYSA-N 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012966 insertion method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Abstract
Description
本発明の実施形態は、医用画像処理装置に関する。 Embodiments described herein relate generally to a medical image processing apparatus.
一般に、虚血性心疾患では、冠動脈の閉塞や狭窄等により、心筋への血流が阻害され、血液の供給が不足もしくは途絶えることにより、心臓に障害が生じる。症状としては、主に前胸部、時に左腕や背中に痛みや圧迫感を感じる。虚血性心疾患の患者に対しての治療法は、大きく分けると、バイパス手術、PCI(カテーテル手術)、薬物療法の3つに分類される。 In general, in ischemic heart disease, blood flow to the myocardium is inhibited due to occlusion or stenosis of the coronary arteries, and the supply of blood is insufficient or disrupted, resulting in damage to the heart. Symptoms include pain and pressure mainly in the front chest and sometimes in the left arm and back. The treatment methods for patients with ischemic heart disease are roughly classified into three types: bypass surgery, PCI (catheter surgery), and drug therapy.
バイパス手術は、図20に示すように、狭くなっていたり閉塞していたりする血管に対して他の血管を繋げることにより、その血管を介して虚血となっている部位により多くの血液が流れるようにしてやる治療法である。 In the bypass operation, as shown in FIG. 20, by connecting another blood vessel to a blood vessel that is narrowed or occluded, more blood flows through the blood vessel through the blood vessel. This is a treatment method.
PCIは、図21,図22に示すように、細い管状構造の治療器具を閉塞や狭窄を起こしている血管に直接挿入して強制的に血管を広げる治療法である。 As shown in FIGS. 21 and 22, PCI is a treatment method for forcibly expanding a blood vessel by directly inserting a treatment device having a thin tubular structure into a blood vessel that is obstructed or narrowed.
薬物療法は、心臓の虚血を改善したり、血栓ができるのを予防したりする治療法である。 Drug therapy is a treatment that improves cardiac ischemia and prevents the formation of blood clots.
医師がこれら3つの治療法のうちのどの治療法を選択するかの指標として、FFR(Fractional Flow Reserve)がある。 FFR (Fractional Flow Reserve) is used as an index of which of the three treatment methods a doctor selects.
狭窄の進行度の評価は、一般的に、プレッシャーワイヤーを直接血管に挿入して行う。
プレッシャーワイヤーは、図23に示すように挿入され、狭窄部の前後の圧力Pin,Poutを計測する。
Evaluation of the degree of progression of stenosis is generally performed by inserting a pressure wire directly into a blood vessel.
The pressure wire is inserted as shown in FIG. 23 and measures pressures P in and P out before and after the constriction.
ここで、FFRはPout/Pinで定義され、この値が0.8より低いと医師はPCIを治療法として選択し、この値が0.8より高いと医師は薬物療法を治療法として選択する。しかしながら、プレッシャーワイヤーを用いた圧力Pin,Poutの計測は侵襲的であるため、非侵襲的な計測法及びFFRの推定法が望まれている。 Here, FFR is defined by P out / P in , and if this value is lower than 0.8, the doctor selects PCI as the treatment method, and if this value is higher than 0.8, the doctor sets the drug therapy as the treatment method. select. However, since the measurement of the pressures P in and P out using a pressure wire is invasive, a non-invasive measurement method and an FFR estimation method are desired.
そこで、近年、流体解析を用いたシミュレーションベースのFFRの推定法が考案されている。既存のシミュレーションは3D画像を用いてのシミュレーションである。このようなシミュレーションベースのFFRの推定法の基本概念としては、モダリティから取得される血管の形状と、血液等のもつ粘性値等の物理パラメータとをインプットとして、例えばCFD(Computational Fluid Dynamics)で用いられるナビエストークスの式を用いてFFRを推定(計算)する。 Therefore, in recent years, a simulation-based FFR estimation method using fluid analysis has been devised. The existing simulation is a simulation using a 3D image. As a basic concept of such a simulation-based FFR estimation method, for example, CFD (Computational Fluid Dynamics) is used as an input of a blood vessel shape acquired from a modality and physical parameters such as viscosity values of blood and the like. FFR is estimated (calculated) using the Naviestokes equation.
このような3Dシミュレーションでの問題点は多大な計算時間を要することである。そのため、FFRを用いて治療法を選択するまでにも時間を要し、一刻を争うときには不向きであるという不都合が生じる。そこで、改善策として、3D画像を用いてのシミュレーションを2D近似することで、シミュレーションに要する時間を大幅に削減する手法がある。 The problem with such 3D simulation is that it requires a lot of calculation time. Therefore, it takes time to select a treatment method using FFR, and there arises a disadvantage that it is unsuitable when competing for a moment. Therefore, there is a method for greatly reducing the time required for the simulation by improving the simulation using the 3D image by 2D approximation as an improvement measure.
これによりFFRをシミュレーションベースで素早く計算することが可能となり、医師はFFRをより有効的な指標として用いることができる。 This makes it possible to quickly calculate FFR on a simulation basis, and doctors can use FFR as a more effective index.
しかしながら、現在、虚血心筋と治療すべき狭窄部との因果関係やリスク評価を行えていない、つまり、虚血心筋と治療すべき狭窄部との因果関係やリスク評価にFFRを適切に反映させることができておらず、どの狭窄部を最優先に治療すべきか等の判断は医師の経験則に依る所が大きい。このため、不要な治療や見落とし等、ヒューマンエラーが発生する恐れがあるという不都合がある。 However, at present, the causal relationship between the ischemic myocardium and the stenosis to be treated and the risk assessment have not been performed, that is, the FFR is appropriately reflected in the causal relationship and risk assessment between the ischemic myocardium and the stenosis to be treated. The decision as to which stenosis should be treated with the highest priority largely depends on the doctor's empirical rules. For this reason, there is an inconvenience that a human error may occur such as unnecessary treatment or oversight.
また、例えば心筋梗塞に寄与する冠動脈では血流量及び圧力が低下するため、上記したFFRは見かけ上高くなる。このため、心筋梗塞といった重大な症状を引き起こしているにも関わらず、FFRが見かけ上高い。この様な場合においても、症状の見落としや治療法の選択ミス等、ヒューマンエラーが発生する恐れがあるという不都合がある。 In addition, for example, in the coronary artery that contributes to myocardial infarction, the blood flow rate and pressure are reduced, and thus the above-described FFR is apparently high. For this reason, FFR is apparently high in spite of causing serious symptoms such as myocardial infarction. Even in such a case, there is an inconvenience that a human error may occur, such as an oversight of a symptom or an incorrect selection of a treatment method.
目的は、ヒューマンエラーの可能性を低減し得る医用画像処理装置を提供することである。 An object is to provide a medical image processing apparatus that can reduce the possibility of human error.
本実施形態に係る医用画像処理装置は、第1取得部と、計算部と、第2取得部と、第1特定部と、第2特定部と、表示部とを含む。第1取得部は、心臓に関する画像のデータに描出された複数の冠動脈を取得し、当該取得された各冠動脈に描出された少なくとも1つの狭窄部位を取得する。計算部は、前記取得された複数の冠動脈の形状に基づいて、前記取得された各冠動脈の圧較差を計算する。第2取得部は、前記画像に描出された心臓に関する心筋梗塞領域または虚血領域を取得する。第1特定部は、前記取得された心筋梗塞領域または虚血領域を前記取得された各冠動脈と支配域とを関係付ける支配マップに照会することにより当該心筋梗塞領域または虚血領域の責任血管を特定する。第2特定部は、前記特定された責任血管内の狭窄部位に対応した前記圧較差に基づいて責任狭窄を特定する。表示部は、前記特定された責任狭窄が描出された画像を、当該責任狭窄を示す情報と共に表示する。 The medical image processing apparatus according to the present embodiment includes a first acquisition unit, a calculation unit, a second acquisition unit, a first specifying unit, a second specifying unit, and a display unit. The first acquisition unit acquires a plurality of coronary arteries depicted in image data relating to the heart, and obtains at least one stenosis site depicted in each of the acquired coronary arteries. The calculation unit calculates a pressure difference of each acquired coronary artery based on the acquired shapes of the plurality of coronary arteries. The second acquisition unit acquires a myocardial infarction region or an ischemic region related to the heart depicted in the image. The first specifying unit refers to the control map that associates each acquired coronary artery and control region with the acquired myocardial infarction region or ischemic region, thereby determining a responsible blood vessel of the myocardial infarction region or ischemic region. Identify. The second identifying unit identifies the responsible stenosis based on the pressure difference corresponding to the identified stenotic site in the responsible blood vessel. The display unit displays an image depicting the identified responsible stenosis together with information indicating the responsible stenosis.
[第1の実施形態]
図1は、第1の実施形態に係る医用画像処理装置を含む医用画像処理システムの構成例を示す模式図であり、図2は同実施形態に係るテリトリーマップ記憶部に記憶されるテリトリーマップの一例を示す模式図である。図1に示す医用画像処理システム1は、医用画像処理装置10と、CT(Computed Tomography)装置20及びPACS(Picture Archiving and Communication System)30とが、例えばLAN(Local Area Network)や公衆電子通信回線等のネットワーク40を介して通信可能に接続されたシステムである。このため、医用画像処理装置10には、CT装置20及びPACS50との通信を可能にする通信インターフェース103が設けられている。
[First Embodiment]
FIG. 1 is a schematic diagram illustrating a configuration example of a medical image processing system including a medical image processing apparatus according to the first embodiment, and FIG. 2 illustrates a territory map stored in a territory map storage unit according to the first embodiment. It is a schematic diagram which shows an example. A medical image processing system 1 shown in FIG. 1 includes a medical image processing apparatus 10, a CT (Computed Tomography) apparatus 20, and a PACS (Picture Archiving and Communication System) 30, for example, a LAN (Local Area Network) or a public electronic communication line. It is the system connected so that communication was possible via the networks 40, such as. For this reason, the medical image processing apparatus 10 is provided with a communication interface 103 that enables communication with the CT apparatus 20 and the PACS 50.
医用画像処理装置10は、図1に示すように、画像記憶部101、テリトリーマップ記憶部102、通信インターフェース103、制御部104、心臓領域抽出部105、心筋解析部106、冠動脈解析部107、責任血管特定部108、FFR計算部109、責任狭窄特定部110、マーカ発生部111及び表示部112を備えている。以下に、医用画像処理装置10を構成する各部101乃至112の機能について詳細に説明する。 As shown in FIG. 1, the medical image processing apparatus 10 includes an image storage unit 101, a territory map storage unit 102, a communication interface 103, a control unit 104, a cardiac region extraction unit 105, a myocardial analysis unit 106, a coronary artery analysis unit 107, and a responsibility. A blood vessel specifying unit 108, an FFR calculating unit 109, a responsible stenosis specifying unit 110, a marker generating unit 111, and a display unit 112 are provided. Hereinafter, functions of the respective units 101 to 112 constituting the medical image processing apparatus 10 will be described in detail.
画像記憶部101は、制御部104の下でCT装置20又はPACS50から送信された処理画像としての当該被検体の心臓を含む胸部領域に関する複数時相にわたる時系列の3次元造影CT画像データ(以下、単にボリュームデータと表記)を記憶する記憶装置である。 The image storage unit 101 includes time-series three-dimensional contrast CT image data (hereinafter referred to as a plurality of time phases) related to a chest region including the heart of the subject as a processed image transmitted from the CT apparatus 20 or the PACS 50 under the control unit 104. , Simply a volume data).
テリトリーマップ記憶部102は、図2に示すように、冠動脈と各冠動脈により栄養供給がなされる支配域との関係を定義するテリトリーマップ(以下、支配マップと表記)を記憶する記憶装置である。 As shown in FIG. 2, the territory map storage unit 102 is a storage device that stores a territory map (hereinafter referred to as a “dominance map”) that defines the relationship between the coronary arteries and the control areas in which nutrition is supplied by each coronary artery.
心臓領域抽出部105は、ボリュームデータから心輪郭抽出処理等により心臓領域を抽出する。 The heart region extraction unit 105 extracts a heart region from the volume data by a heart contour extraction process or the like.
心筋解析部106は、造影剤濃度に対応するCT値による閾値処理等により、心臓領域抽出部105により抽出された心臓領域から心筋領域を抽出する。また、心筋解析部106は、心筋パフュージョン解析、つまり、抽出した心筋領域内の画素毎又は局所毎に造影剤に関する時間濃度曲線を生成し、その時間濃度曲線に基づいて画素毎又は局所毎に造影剤が流入してから流出するまでの期間に移動する血流の量を算出する。 The myocardial analysis unit 106 extracts a myocardial region from the heart region extracted by the heart region extracting unit 105 by threshold processing using a CT value corresponding to the contrast agent concentration. The myocardial analysis unit 106 generates myocardial perfusion analysis, that is, generates a time concentration curve for the contrast agent for each pixel or each region in the extracted myocardial region, and for each pixel or each region based on the time concentration curve. The amount of blood flow that moves during the period from the inflow of the contrast agent to the outflow is calculated.
例えば、CT装置を用いた撮影では非イオン性造影剤を患者へ注入し、CT値の変化から臓器の灌流情報を描出することができる。このため、CTパフュージョン解析では、例えば、512×512ピクセルで構成されたCT画像(ボリュームデータ)の経時変化を各画素におけるCT値の変化から測定し、血流量等を数値化することができる。このようにして、複数時相のCT画像から臓器の灌流情報(例えば、血流量)を表す1枚のカラーマップが生成される。 For example, in imaging using a CT apparatus, a nonionic contrast agent can be injected into a patient, and organ perfusion information can be depicted from changes in CT values. Therefore, in CT perfusion analysis, for example, a temporal change of a CT image (volume data) composed of 512 × 512 pixels can be measured from a change in CT value in each pixel, and a blood flow rate or the like can be quantified. . In this way, one color map representing organ perfusion information (for example, blood flow) is generated from a plurality of time-phase CT images.
更に、心筋解析部106は、算出した血流量の空間分布から閾値処理により虚血領域を特定する。 Furthermore, the myocardial analyzer 106 identifies an ischemic region by threshold processing from the calculated spatial distribution of blood flow.
冠動脈解析部107は、心臓領域抽出部105により抽出された心臓領域から複数の冠動脈を抽出し、更に、当該抽出した各冠動脈から少なくとも1つの狭窄部位を抽出する。
具体的には、冠動脈解析部107は、冠動脈の血管芯線や血管内壁等に沿って、冠動脈の解剖学的構造やプラーク性状の分析を行い、冠動脈に関するボリュームデータを抽出する、つまり、冠動脈ならびにその冠動脈の内壁に位置する狭窄部位を抽出する。なお、プラーク性状の具体例としては、脂質量、血清コレステロール濃度、硬さ、石灰化度、繊維性被膜(Thin-cap)の厚さ、及びFFR値(但し、本実施形態では、FFR値はFFR計算部109により算出されるものとする)等が挙げられる。
The coronary artery analysis unit 107 extracts a plurality of coronary arteries from the heart region extracted by the heart region extraction unit 105, and further extracts at least one stenosis site from each of the extracted coronary arteries.
Specifically, the coronary artery analysis unit 107 analyzes the anatomical structure and plaque characteristics of the coronary artery along the coronary artery blood vessel core line and the inner wall of the blood vessel, and extracts volume data related to the coronary artery. A stenosis site located on the inner wall of the coronary artery is extracted. Specific examples of the plaque properties include lipid amount, serum cholesterol concentration, hardness, calcification degree, fibrous coating (Thin-cap) thickness, and FFR value (in this embodiment, the FFR value is And the like are calculated by the FFR calculation unit 109).
責任血管特定部108は、心筋解析部106により特定された虚血領域に対して、テリトリーマップ記憶部102に記憶された支配マップを照会することにより、虚血領域へ本来的に栄養供給責任を有している血管(以下、責任血管と表記)を特定する。 The responsible blood vessel specifying unit 108 inherently takes responsibility for supplying nutrients to the ischemic region by inquiring the control map stored in the territory map storage unit 102 for the ischemic region specified by the myocardial analysis unit 106. Identify the blood vessels you have (hereinafter referred to as responsible blood vessels).
FFR計算部109は、冠動脈解析部107により抽出された各狭窄部位に対応するFFRの値をシミュレーションベースで計算する。具体的には、まずFFR計算部109は、冠動脈解析部107により抽出された狭窄部位毎に、心筋解析部106により生成されたカラーマップに基づいて、各狭窄部位の下流の少なくとも一つの位置において組織血流量と、各狭窄部位の上流の少なくとも一つの位置において組織血流量とを算出する。そして、FFR計算部109は、算出した狭窄部位の下流の組織血流量を、算出した狭窄部位の上流の組織血流量で除算することで、少なくとも狭窄部位を含む各位置におけるFFR値を算出する。なお、本実施形態では、FFR計算部109は、上記のような計算手法によりFFR値を算出するものとしたが、FFR値の計算手法はこれに限定されるものでなく、各狭窄部位に対応したFFR値を算出可能であれば、FFR計算部109で用いるFFR値の計算手法として適宜適用可能である。 The FFR calculation unit 109 calculates the FFR value corresponding to each stenosis site extracted by the coronary artery analysis unit 107 on a simulation basis. Specifically, first, the FFR calculation unit 109 performs, for each stenosis site extracted by the coronary artery analysis unit 107, at least one position downstream of each stenosis site based on the color map generated by the myocardial analysis unit 106. A tissue blood flow and a tissue blood flow are calculated at at least one position upstream of each stenosis site. Then, the FFR calculation unit 109 calculates the FFR value at each position including at least the stenosis site by dividing the calculated tissue blood flow downstream of the stenosis site by the calculated tissue blood flow upstream of the stenosis site. In this embodiment, the FFR calculation unit 109 calculates the FFR value by the calculation method as described above, but the FFR value calculation method is not limited to this, and corresponds to each stenosis site. If the calculated FFR value can be calculated, it can be appropriately applied as a method for calculating the FFR value used in the FFR calculation unit 109.
責任狭窄特定部110は、冠動脈解析部107により抽出された狭窄部位のうち、責任血管特定部108により特定された責任血管の内壁に位置する狭窄部位(以下、責任狭窄と表記)、つまり、責任狭窄候補の中からFFR値が閾値未満である狭窄部位を責任狭窄として特定する。 The responsible stenosis specifying unit 110 is a stenosis site (hereinafter referred to as responsible stenosis) located on the inner wall of the responsible blood vessel specified by the responsible blood vessel specifying unit 108 among the stenosis sites extracted by the coronary artery analyzing unit 107, that is, the responsible Among the stenosis candidates, a stenosis site whose FFR value is less than the threshold is specified as a responsible stenosis.
マーカ発生部111は、責任血管特定部108により特定された責任血管や、責任狭窄特定部110により特定された責任狭窄、FFR計算部109により算出されたFFR値、ならびに冠動脈解析部107により抽出された責任狭窄候補等を表すマーカのデータを発生する。これらマーカは、ボリュームデータからレンダリング等により生成された3次元画像や、断面変換(Multi-Planar Reconstruction)により生成された2次元画像に重畳させて、表示部112に表示される。なお、マーカ発生部111から発生したマーカを重畳させる画像としては、CT装置20によるボリュームデータに由来する画像には限定されず、例えばX線診断装置等、他のモダリティから取得される画像であってもよい。 The marker generating unit 111 is extracted by the responsible blood vessel identified by the responsible blood vessel identifying unit 108, the responsible stenosis identified by the responsible stenosis identifying unit 110, the FFR value calculated by the FFR calculating unit 109, and the coronary artery analyzing unit 107. Marker data representing the candidate responsible stenosis is generated. These markers are displayed on the display unit 112 by being superimposed on a three-dimensional image generated by rendering or the like from volume data or a two-dimensional image generated by cross-sectional transformation (Multi-Planar Reconstruction). Note that the image on which the marker generated from the marker generation unit 111 is superimposed is not limited to an image derived from the volume data obtained by the CT apparatus 20, and may be an image acquired from another modality such as an X-ray diagnostic apparatus. May be.
ここで、本実施形態に係る医用画像処理装置10の動作の一例について、図2,図4乃至図7の模式図と、図3に示すフローチャートとを参照しながら説明する。 Here, an example of the operation of the medical image processing apparatus 10 according to the present embodiment will be described with reference to the schematic diagrams of FIGS. 2, 4 to 7, and the flowchart illustrated in FIG. 3.
始めに、制御部104は、通信インターフェース103を介して、CT装置20又はPACS50から胸部領域に関する複数時相にわたる時系列のボリュームデータの入力を受け付けると、当該入力を受け付けたボリュームデータを画像記憶部101に書込む(ステップS1)。 First, when receiving the input of time-series volume data over a plurality of time phases related to the chest region from the CT apparatus 20 or the PACS 50 via the communication interface 103, the control unit 104 receives the input volume data from the image storage unit. Write to 101 (step S1).
続いて、心臓領域抽出部105は、制御部104の下で拍動が比較的少ない特定時相のボリュームデータを処理画像として画像記憶部101から読出し、当該ボリュームデータから心臓領域を抽出する(ステップS2)。 Subsequently, the heart region extraction unit 105 reads volume data of a specific time phase with relatively few pulsations under the control unit 104 from the image storage unit 101 as a processed image, and extracts a heart region from the volume data (step). S2).
次に、冠動脈解析部107は、心臓領域抽出部105により抽出された心臓領域を対象にして冠動脈解析処理を実行する(ステップS3,S4)。具体的には、冠動脈解析部107は、冠動脈の血管芯線や血管内壁等に沿って、冠動脈の解剖学的構造やプラーク性状の分析を行い、冠動脈ならびにその冠動脈の内壁に位置する狭窄部位を抽出する。その後、冠動脈解析部107は、例えば、図4(a),図4(b)に示すように、冠動脈の解剖学的構造を心臓形態画像に重畳させて3次元画像g1又は2次元画像g2として表示部112に表示させる。なお、図4(a),図4(b)に示す画像g1,g2を表示部112に表示させるタイミングは、操作者が任意に設定可能である、つまり、処理途中に表示させてもよいし、処理結果と共に表示させてもよい。 Next, the coronary artery analysis unit 107 performs a coronary artery analysis process on the heart region extracted by the heart region extraction unit 105 (steps S3 and S4). Specifically, the coronary artery analysis unit 107 analyzes the anatomical structure and plaque characteristics of the coronary artery along the coronary vascular core line and the inner wall of the blood vessel, and extracts the stenosis portion located on the coronary artery and the inner wall of the coronary artery To do. Thereafter, as shown in FIGS. 4A and 4B, for example, the coronary artery analysis unit 107 superimposes the anatomical structure of the coronary artery on the heart morphological image to form a three-dimensional image g1 or a two-dimensional image g2. It is displayed on the display unit 112. The timing at which the images g1 and g2 shown in FIGS. 4A and 4B are displayed on the display unit 112 can be arbitrarily set by the operator, that is, they may be displayed during the processing. , It may be displayed together with the processing result.
続いて、心筋解析部106は、造影剤濃度に対応するCT値による閾値処理により、心臓領域抽出部105により抽出された心臓領域から心筋領域を抽出する(ステップS5)。 Subsequently, the myocardial analysis unit 106 extracts a myocardial region from the heart region extracted by the heart region extracting unit 105 by threshold processing using a CT value corresponding to the contrast agent concentration (step S5).
次に、心筋解析部106は、抽出した心筋に限定してCTパフュージョン解析処理を実行する(ステップS6,S7,S8)。具体的には、心筋解析部106は、時系列のボリュームデータに基づいて、抽出した心筋領域内の画素毎又は局所毎に造影剤に関する時間濃度曲線を生成する。その後、心筋解析部106は、それら時間濃度曲線に基づいて画素毎又は局所毎に造影剤が流入してから流出するまでの期間に移動する血流量を算出する。
これにより、例えば図5に示すように、血流量の空間分布を示すカラーマップg3が生成される。そして、心筋解析部106は、生成したカラーマップg3、つまり、算出した血流量の空間分布に基づいて、所定の血流量未満の領域を虚血領域として特定する。
Next, the myocardial analysis unit 106 executes CT perfusion analysis processing only for the extracted myocardium (steps S6, S7, and S8). Specifically, the myocardial analyzer 106 generates a time concentration curve related to a contrast agent for each pixel or each local area in the extracted myocardial region based on time-series volume data. Thereafter, the myocardial analyzer 106 calculates the blood flow volume that moves during the period from the inflow of the contrast medium to the outflow for each pixel or each region based on these time density curves.
Thereby, for example, as shown in FIG. 5, a color map g3 indicating the spatial distribution of the blood flow is generated. Then, the myocardial analysis unit 106 identifies a region less than the predetermined blood flow as an ischemic region based on the generated color map g3, that is, the calculated spatial distribution of blood flow.
続いて、責任血管特定部108は、心筋解析部106により特定された虚血領域に対して、図2に示すように、テリトリーマップ記憶部102に記憶されている支配マップを照会することにより、責任血管を特定する(ステップS9)。 Subsequently, the responsible blood vessel specifying unit 108 refers to the control map stored in the territory map storage unit 102 for the ischemic region specified by the myocardial analysis unit 106, as shown in FIG. A responsible blood vessel is specified (step S9).
次に、FFR計算部109は、責任血管特定部108により特定された責任血管の内壁に位置する狭窄部位毎に、心筋解析部106により生成されたカラーマップg3に基づいて、各狭窄部位の下流の組織血流量と、各狭窄部位の上流の組織血流量とを算出する。そして、FFR計算部109は、算出した狭窄部位の下流の組織血流量を、算出した狭窄部位の上流の組織血流量で除算することで、少なくとも狭窄部位を含む各位置におけるFFR値を算出する(ステップS10)。 Next, for each stenosis site located on the inner wall of the responsible blood vessel identified by the responsible blood vessel identification unit 108, the FFR calculation unit 109 downstream of each stenosis site based on the color map g3 generated by the myocardial analysis unit 106. And the tissue blood flow upstream of each stenosis site are calculated. Then, the FFR calculation unit 109 calculates the FFR value at each position including at least the stenosis site by dividing the calculated tissue blood flow downstream of the stenosis site by the calculated tissue blood flow upstream of the stenosis site ( Step S10).
続いて、責任狭窄特定部110は、FFR計算部109により算出されたFFR値が閾値未満である狭窄部位を責任狭窄として特定する(ステップS11)。 Subsequently, the responsible stenosis specifying unit 110 specifies a stenosis site in which the FFR value calculated by the FFR calculating unit 109 is less than the threshold value as the responsible stenosis (step S11).
しかる後、表示部112は、例えば図6,7に示すように、マーカ発生部111により発生した責任血管、責任狭窄、ならびにFFR値を表すマーカをボリュームデータに由来した3次元画像g4又は2次元画像g5に重畳させて表示する(ステップS12)。 Thereafter, as shown in FIGS. 6 and 7, for example, the display unit 112 displays the responsible blood vessel, the responsible stenosis, and the marker representing the FFR value generated by the marker generating unit 111 as a three-dimensional image g4 or two-dimensional image derived from the volume data. The image is superimposed and displayed on the image g5 (step S12).
以上説明した一実施形態によれば、CT装置20によるボリュームデータから心臓領域、心筋領域、冠動脈及び狭窄部位を抽出可能な心臓領域抽出部105、心筋解析部106及び冠動脈解析部107と、心筋解析部106による処理結果に基づいて責任血管を特定する責任血管特定部108と、冠動脈解析部107、責任血管特定部108及びFFR計算部109による処理結果に基づいて責任狭窄を特定する責任狭窄特定部110と、責任血管や責任狭窄に関するマーカをボリュームデータに由来した3次元画像又は2次元画像に重畳表示させる表示部112とを備えた構成により、例えば、図6に示すように、責任狭窄と支配域との対応関係を医師に対して視覚的に示すことができる、ひいては、ヒューマンエラーの可能性を低減させることができる。 According to the embodiment described above, the heart region extraction unit 105, the myocardial analysis unit 106, and the coronary artery analysis unit 107 that can extract the heart region, myocardial region, coronary artery, and stenosis from the volume data obtained by the CT apparatus 20, and the myocardial analysis. Responsible blood vessel identifying unit 108 that identifies the responsible blood vessel based on the processing result by the unit 106, Responsible stenosis identifying unit that identifies the responsible stenosis based on the processing results by the coronary artery analyzing unit 107, the responsible blood vessel identifying unit 108, and the FFR calculating unit 109 110 and a display unit 112 that superimposes and displays a responsible blood vessel or responsible stenosis marker on a volume data-derived 3D image or 2D image, for example, as shown in FIG. The correspondence with the community can be shown visually to the doctor, which reduces the possibility of human error. It is possible.
また、本実施形態では、FFR計算部109がシミュレーションベースでFFR値を算出する、つまり、プレッシャーワイヤー等、侵襲性をもつものを用いないため、検査時に患者にかかる負担を低減させることができる。 Further, in this embodiment, the FFR calculation unit 109 calculates the FFR value on a simulation basis, that is, since an invasive device such as a pressure wire is not used, the burden on the patient during the examination can be reduced.
なお、本実施形態では、表示部112が表示する画像の一例として図6,7を示したが、表示部112が表示する画像はこれらに限定されるものでなく、例えば、図8A,図8Bに示すような画像g6,g7が表示されるとしてもよい。図8A,図8Bは、FFR計算部109により算出されたFFR値に基づいて、狭窄部位の治療の優先度を順位付し、この処理結果をマーカとして重畳した画像g6,g7の一例を示す。ここでは、治療の優先度が高い狭窄部位ほど狭窄部位を示すマーカ(丸印)が大きくなっている。また、図8Aでは、優先度の他に、冠動脈解析部107によるプラーク性状の分析結果を示す棒グラフをマーカとして重畳した画像g6の一例を示す。なお、図8Aに示す画像g6には、プラーク性状の分析結果の他に、血液の性質等をマーカとして更に重畳させてもよい。また、責任血管が複数存在する場合、複数の責任血管にわたって狭窄部位の治療の優先度を順位付した後に、この処理結果をマーカとして重畳した画像を表示部112に表示させることもできる。 In this embodiment, FIGS. 6 and 7 are shown as examples of images displayed on the display unit 112. However, the images displayed on the display unit 112 are not limited to these, and for example, FIGS. 8A and 8B. Images g6 and g7 as shown in FIG. 8A and 8B show an example of images g6 and g7 in which priorities of treatments for stenotic sites are ranked based on the FFR values calculated by the FFR calculation unit 109, and the processing results are superimposed as markers. Here, the stenosis part (circle) which shows a stenosis part is large, so that a treatment priority is high. FIG. 8A shows an example of an image g6 in which a bar graph indicating the analysis result of the plaque property by the coronary artery analysis unit 107 is superimposed as a marker in addition to the priority. In addition to the plaque property analysis result, blood properties or the like may be further superimposed on the image g6 shown in FIG. 8A as a marker. Further, when there are a plurality of responsible blood vessels, the priority of the treatment of the stenosis site is ranked over the plurality of responsible blood vessels, and then an image in which the processing results are superimposed as a marker can be displayed on the display unit 112.
また、本実施形態に係る医用画像処理装置10によれば、責任血管特定部108による支配マップ照会時に、ある動脈が供給すべき領域の一部に対して別の動脈の領域が張り出しているか否かを検出し、別の動脈の領域が張り出していることが検出された場合に、側副血管による血液供給の可能性があることを示唆することもできる。一般に、側副血管が存在するとFFR値の信頼性は低下するため、例えば、図9に示すように、側副血管の存在を示唆する画像g8を表示部112に表示させることで、よりヒューマンエラーの可能性を低減させることができる。 In addition, according to the medical image processing apparatus 10 according to the present embodiment, whether or not another arterial region protrudes from a part of a region to be supplied by a certain artery when the responsible blood vessel reference unit 108 makes a reference to the control map. If it is detected that another arterial region is overhanging, it may be suggested that there is a possibility of blood supply by the collateral blood vessel. In general, since the reliability of the FFR value decreases when a collateral blood vessel exists, for example, as shown in FIG. 9, by displaying an image g8 suggesting the presence of a collateral blood vessel on the display unit 112, more human error The possibility of this can be reduced.
更に、本実施形態に係る医用画像処理装置10によれば、責任狭窄特定部110により責任狭窄を特定した後に、ボリュームデータに由来する2次元画像から責任狭窄の断面の径や長さを測定し、この測定結果に基づいて最適なカテーテル・ステントのサイズを示唆することもできる。具体的には、例えば、図10に示すように、最適なカテーテル・ステントのサイズを示唆する画像g9を表示部112に表示させることで、よりヒューマンエラーの可能性を低減させることができる。 Furthermore, according to the medical image processing apparatus 10 according to the present embodiment, after the responsible stenosis identifying unit 110 identifies the responsible stenosis, the diameter and length of the cross section of the responsible stenosis are measured from the two-dimensional image derived from the volume data. Based on this measurement result, an optimal catheter / stent size can be suggested. Specifically, for example, as shown in FIG. 10, by displaying an image g9 suggesting an optimal catheter / stent size on the display unit 112, the possibility of human error can be further reduced.
また、本実施形態に係る医用画像処理装置10によれば、心筋解析部106により虚血領域が複数個所存在する旨が示された場合に、ボリュームデータから心筋の厚み等を測定し、この測定結果に基づいて壊死している心筋の支配領域に対応する血管は責任血管として特定しないといった設定を付加することも可能である。 Further, according to the medical image processing apparatus 10 according to the present embodiment, when the myocardial analysis unit 106 indicates that there are a plurality of ischemic regions, the thickness of the myocardium is measured from the volume data, and this measurement is performed. Based on the result, it is possible to add a setting such that the blood vessel corresponding to the dominant region of the myocardium that is necrotic is not identified as the responsible blood vessel.
更に、本実施形態に係る医用画像処理装置10によれば、表示部112により所望の3次元画像を表示したときに、マウスやキーボード、タッチパネル等の図示しない入力インターフェースから責任血管や責任狭窄を選択する旨の入力を受け付けると、当該選択された責任血管や責任狭窄が観察しやすい角度に3次元画像を自動で回転させることも可能である。 Furthermore, according to the medical image processing apparatus 10 according to the present embodiment, when a desired three-dimensional image is displayed on the display unit 112, a responsible blood vessel or responsible stenosis is selected from an input interface (not shown) such as a mouse, a keyboard, or a touch panel. When the input to the effect is accepted, the three-dimensional image can be automatically rotated to an angle at which the selected responsible blood vessel or responsible stenosis is easily observed.
[第2の実施形態]
図11は第2の実施形態に係る医用画像処理装置を含む医用画像処理システムの構成例を示す模式図である。図11に示す医用画像処理システム1は、医用画像処理装置10と、CT(Computed Tomography)装置20、MRI(Magnetic Resonance Imaging)装置30、核医学診断装置40及びPACS(Picture Archiving and Communication System)50とが、例えばLAN(Local Area Network)や公衆電子通信回線等のネットワーク60を介して通信可能に接続されたシステムである。このため、医用画像処理装置10には、CT装置20、MRI装置30、核医学診断装置40及びPACS50との通信を可能にする通信インターフェース103が設けられている。
[Second Embodiment]
FIG. 11 is a schematic diagram illustrating a configuration example of a medical image processing system including a medical image processing apparatus according to the second embodiment. A medical image processing system 1 shown in FIG. 11 includes a medical image processing apparatus 10, a CT (Computed Tomography) apparatus 20, an MRI (Magnetic Resonance Imaging) apparatus 30, a nuclear medicine diagnostic apparatus 40, and a PACS (Picture Archiving and Communication System) 50. Are systems that are communicably connected via a network 60 such as a LAN (Local Area Network) or a public electronic communication line. Therefore, the medical image processing apparatus 10 is provided with a communication interface 103 that enables communication with the CT apparatus 20, the MRI apparatus 30, the nuclear medicine diagnosis apparatus 40, and the PACS 50.
医用画像処理装置10は、図11に示すように、画像記憶部101、テリトリーマップ記憶部102、通信インターフェース103、制御部104、心臓領域抽出部105、心筋解析部106、冠動脈解析部107、責任血管特定部108、FFR計算部109、マーカ発生部111及び表示部111を備えている。第1の実施形態で示した医用画像処理装置1と異なる構成についてのみ、以下説明する。 As shown in FIG. 11, the medical image processing apparatus 10 includes an image storage unit 101, a territory map storage unit 102, a communication interface 103, a control unit 104, a cardiac region extraction unit 105, a myocardial analysis unit 106, a coronary artery analysis unit 107, and a responsibility. A blood vessel specifying unit 108, an FFR calculating unit 109, a marker generating unit 111, and a display unit 111 are provided. Only the configuration different from the medical image processing apparatus 1 shown in the first embodiment will be described below.
画像記憶部101は、制御部104の下でCT装置20又はPACS50から送信された処理画像としての当該被検体の心臓を含む胸部領域に関する複数時相にわたる時系列の3次元造影CT画像データを記憶する記憶装置である。 The image storage unit 101 stores time-series three-dimensional contrast CT image data over a plurality of time phases related to the chest region including the heart of the subject as a processed image transmitted from the CT apparatus 20 or the PACS 50 under the control unit 104. Storage device.
テリトリーマップ記憶部102は、図2に示すように、冠動脈と各冠動脈により栄養供給がなされる支配域との関係を定義するテリトリーマップを記憶する記憶装置である。 As shown in FIG. 2, the territory map storage unit 102 is a storage device that stores a territory map that defines the relationship between the coronary arteries and the dominating area where nutrition is supplied by each coronary artery.
心臓領域抽出部105は、ボリュームデータから心輪郭抽出処理等により心臓領域を抽出する。 The heart region extraction unit 105 extracts a heart region from the volume data by a heart contour extraction process or the like.
心筋解析部106は、造影剤濃度に対応するCT値による閾値処理等により、心臓領域抽出部105により抽出された心臓領域から心筋領域を抽出する。また、心筋解析部106は、MRI装置30による遅延造影(late gadolinium enhancement)や、核医学診断装置40による糖代謝測定等により、抽出した心筋領域のうち壊死している心筋領域、つまり、心筋梗塞に寄与する心筋梗塞領域を特定する。 The myocardial analysis unit 106 extracts a myocardial region from the heart region extracted by the heart region extracting unit 105 by threshold processing using a CT value corresponding to the contrast agent concentration. The myocardial analysis unit 106 also performs a necrotic myocardial region, that is, myocardial infarction among the extracted myocardial regions by delayed contrast enhancement (late gadolinium enhancement) by the MRI apparatus 30 or sugar metabolism measurement by the nuclear medicine diagnostic apparatus 40. The myocardial infarction region that contributes to is identified.
更に、心筋解析部106は、心筋パフュージョン解析、つまり、抽出した心筋領域内の画素毎又は局所毎に造影剤が流入してから流出するまでの期間に移動する血流の量を算出する。例えば、CT装置20を用いた撮影では非イオン性造影剤を患者へ注入し、CT値の変化から臓器の灌流情報を描出することができる。このため、CTパフュージョン解析では、例えば、512×512ピクセルで構成されたCT画像(ボリュームデータ)の経時変化を各画素におけるCT値の変化から測定し、血流量等を数値化することができる。
このようにして、複数時相のCT画像から臓器の灌流情報(例えば、血流量)を表す1枚のカラーマップが生成される。
Furthermore, the myocardial analyzer 106 calculates myocardial perfusion analysis, that is, the amount of blood flow that moves during the period from the inflow of the contrast medium to the outflow for each pixel or local area in the extracted myocardial region. For example, in imaging using the CT apparatus 20, a nonionic contrast agent can be injected into a patient, and organ perfusion information can be depicted from changes in CT values. Therefore, in CT perfusion analysis, for example, a temporal change of a CT image (volume data) composed of 512 × 512 pixels can be measured from a change in CT value in each pixel, and a blood flow rate or the like can be quantified. .
In this way, one color map representing organ perfusion information (for example, blood flow) is generated from a plurality of time-phase CT images.
即ち、心筋解析部106は心筋梗塞領域を特定するだけでなく、例えば、算出した血流量の空間分布から閾値処理により血流低下部位、つまり、虚血領域を特定することもできる。 That is, the myocardial analysis unit 106 not only specifies the myocardial infarction region, but also can specify a blood flow lowering portion, that is, an ischemic region by threshold processing from the calculated spatial distribution of blood flow.
冠動脈解析部107は、心臓領域抽出部105により抽出された心臓領域から複数の冠動脈を抽出し、更に、当該抽出した各冠動脈から少なくとも1つの狭窄部位を抽出する。
具体的には、冠動脈解析部107は、冠動脈の血管芯線や血管内壁等に沿って、冠動脈の解剖学的構造やプラーク性状の分析を行い、冠動脈に関するボリュームデータを抽出する、つまり、冠動脈ならびにその冠動脈の内壁に位置する狭窄部位を抽出する。なお、プラーク性状には、脂質量、血清コレステロール濃度、硬さ、石灰度及び繊維性被膜(Thin-cap)の厚さ等が含まれる。
The coronary artery analysis unit 107 extracts a plurality of coronary arteries from the heart region extracted by the heart region extraction unit 105, and further extracts at least one stenosis site from each of the extracted coronary arteries.
Specifically, the coronary artery analysis unit 107 analyzes the anatomical structure and plaque characteristics of the coronary artery along the coronary artery blood vessel core line and the inner wall of the blood vessel, and extracts volume data related to the coronary artery. A stenosis site located on the inner wall of the coronary artery is extracted. The plaque properties include lipid amount, serum cholesterol concentration, hardness, calcification, and thickness of fibrous coating (Thin-cap).
責任血管特定部108は、心筋解析部106により特定された心筋梗塞領域に対して、テリトリーマップ記憶部102に記憶された支配マップを照会することにより、心筋梗塞領域へ本来的に栄養供給責任を有している血管(以下、梗塞責任血管と表記)を特定する。 The responsible blood vessel specifying unit 108 inherently assumes the nutrition supply responsibility to the myocardial infarction region by inquiring the control map stored in the territory map storage unit 102 for the myocardial infarction region specified by the myocardial analysis unit 106. A blood vessel (hereinafter referred to as an infarct responsible blood vessel) is specified.
なお、心筋解析部106により心筋梗塞領域に代えて虚血領域が特定されている場合、責任血管特定部108は、心筋解析部106により特定された虚血領域に対して、テリトリーマップ記憶部102に記憶された支配マップを照会することにより、虚血領域へ本来的に栄養供給責任を有している血管(以下、虚血責任血管と表記)を特定する。 When the myocardial analysis unit 106 specifies the ischemic region instead of the myocardial infarction region, the responsible blood vessel specifying unit 108 performs the territory map storage unit 102 for the ischemic region specified by the myocardial analysis unit 106. The blood vessel that is inherently responsible for supplying nutrients to the ischemic region (hereinafter, referred to as ischemic responsible blood vessel) is identified by querying the control map stored in FIG.
FFR計算部109は、冠動脈解析部107により抽出された各冠動脈上の狭窄部位を少なくとも含む複数の位置に関するFFR値をシミュレーションベースで計算する。具体的には、まずFFR計算部109は、心筋解析部106により生成されたカラーマップに基づいて、冠動脈内の狭窄部位の下流側の少なくとも一つの位置において組織血流量と、冠動脈内の狭窄部位の上流の少なくとも一つの位置において組織血流量とを算出する。そして、FFR計算部109は、算出した狭窄部位の下流の組織血流量を、算出した狭窄部位の上流の組織血流量で除算することで、少なくとも狭窄部位を含む各位置におけるFFR値を算出する。なお、ここでは、冠動脈内の1つの狭窄部位に対応したFFR値を算出する場合について説明したが、例えば、冠動脈全体のFFR値を算出する場合、冠動脈内の最も上流に位置する狭窄部位の上流の組織血流量を基準値として固定し、他の複数箇所の組織血流量(但し、FFR計算部109が冠動脈内の複数箇所の組織血流量を算出しているものとする)を変数とすることで、FFR計算部109は、当該複数箇所に対応したFFR値、つまり、冠動脈全体のFFR値を算出することができる。 The FFR calculation unit 109 calculates FFR values regarding a plurality of positions including at least a stenosis site on each coronary artery extracted by the coronary artery analysis unit 107 on a simulation basis. Specifically, the FFR calculation unit 109 first determines the tissue blood flow rate and the stenosis site in the coronary artery at at least one position downstream of the stenosis site in the coronary artery based on the color map generated by the myocardial analysis unit 106. Tissue blood flow rate is calculated at at least one position upstream of. Then, the FFR calculation unit 109 calculates the FFR value at each position including at least the stenosis site by dividing the calculated tissue blood flow downstream of the stenosis site by the calculated tissue blood flow upstream of the stenosis site. Here, although the case where the FFR value corresponding to one stenosis site in the coronary artery is calculated has been described, for example, when the FFR value of the entire coronary artery is calculated, the upstream of the stenosis site located most upstream in the coronary artery The tissue blood flow volume is fixed as a reference value, and the tissue blood flow volume at other multiple locations (provided that the FFR calculation unit 109 calculates the tissue blood flow volume at multiple locations within the coronary artery) is a variable. Thus, the FFR calculation unit 109 can calculate the FFR values corresponding to the plurality of locations, that is, the FFR values of the entire coronary artery.
マーカ発生部111は、責任血管特定部108により特定された梗塞責任血管(又は、虚血責任血管)や、FFR計算部109により算出されたFFR値を表すマーカ(マーク)のデータを発生する。これらマーカは、ボリュームデータからレンダリング等により生成された3次元画像や、断面変換(Multi-Planar Reconstruction)により生成された2次元画像に重畳させて、表示部111に表示される。なお、マーカ発生部111から発生したマーカを重畳させる画像としては、CT装置20によるボリュームデータに由来する画像には限定されず、例えばX線診断装置等、他のモダリティから取得される画像であってもよい。 The marker generation unit 111 generates data of an infarct responsible blood vessel (or ischemic responsible blood vessel) specified by the responsible blood vessel specifying unit 108 or a marker (mark) representing the FFR value calculated by the FFR calculation unit 109. These markers are displayed on the display unit 111 by being superimposed on a three-dimensional image generated by rendering or the like from volume data or a two-dimensional image generated by cross-sectional transformation (Multi-Planar Reconstruction). Note that the image on which the marker generated from the marker generation unit 111 is superimposed is not limited to an image derived from the volume data obtained by the CT apparatus 20, and may be an image acquired from another modality such as an X-ray diagnostic apparatus. May be.
ここで、本実施形態に係る医用画像処理装置10の動作の一例について、図2,図13乃至図16の模式図と、図12のフローチャートとを参照しながら説明する。 Here, an example of the operation of the medical image processing apparatus 10 according to the present embodiment will be described with reference to the schematic diagrams of FIGS. 2, 13 to 16, and the flowchart of FIG. 12.
但し、ここでは、画像記憶部101には、CT装置20又はPACS50からの胸部領域に関する複数時相にわたる時系列のボリュームデータが予め記憶されているものとする。また、ここでは、心臓領域抽出部105が、制御部104の下で拍動が比較的少ない特定時相のボリュームデータを処理画像として画像記憶部101から読出し、当該ボリュームデータから心臓領域を既に抽出しているものとする。 However, here, it is assumed that time series volume data over a plurality of time phases related to the chest region from the CT apparatus 20 or the PACS 50 is stored in the image storage unit 101 in advance. Further, here, the heart region extraction unit 105 reads volume data of a specific time phase with relatively few pulsations under the control unit 104 from the image storage unit 101 as a processed image, and already extracts the heart region from the volume data. Suppose you are.
始めに、冠動脈解析部107は、心臓領域抽出部105により抽出された心臓領域を対象にして冠動脈解析処理を実行する(ステップS21)。具体的には、冠動脈解析部107は、冠動脈の血管芯線や血管内壁等に沿って、冠動脈の解剖学的構造やプラーク性状の分析を行い、冠動脈ならびにその冠動脈の内壁に位置する狭窄部位を抽出する。 First, the coronary artery analysis unit 107 executes a coronary artery analysis process for the heart region extracted by the heart region extraction unit 105 (step S21). Specifically, the coronary artery analysis unit 107 analyzes the anatomical structure and plaque characteristics of the coronary artery along the coronary vascular core line and the inner wall of the blood vessel, and extracts the stenosis portion located on the coronary artery and the inner wall of the coronary artery To do.
続いて、心筋解析部106は、造影剤濃度に対応するCT値による閾値処理により、心臓領域抽出部105により抽出された心臓領域から心筋領域を抽出する。その後、心筋解析部106は、MRI装置30による遅延造影や、核医学診断装置40による糖代謝測定等により、抽出した心筋領域のうち壊死している心筋領域、つまり、心筋梗塞に寄与する心筋梗塞領域を特定する(ステップS22,S23,S24)。 Subsequently, the myocardial analyzer 106 extracts a myocardial region from the heart region extracted by the heart region extracting unit 105 by threshold processing using a CT value corresponding to the contrast agent concentration. After that, the myocardial analysis unit 106 performs the delayed contrast imaging by the MRI apparatus 30, the sugar metabolism measurement by the nuclear medicine diagnosis apparatus 40, and the like, and the myocardial infarction that contributes to the myocardial infarction. An area is specified (steps S22, S23, S24).
次に、責任血管特定部108は、心筋解析部106により特定された心筋梗塞領域に対して、テリトリーマップ記憶部102に記憶された支配マップを照会することにより、冠動脈解析部107により抽出された各冠動脈から梗塞責任血管を特定する(ステップS25)。 Next, the responsible blood vessel specifying unit 108 is extracted by the coronary artery analyzing unit 107 by referring to the dominant map stored in the territory map storage unit 102 for the myocardial infarction region specified by the myocardial analyzing unit 106. An infarct responsible blood vessel is identified from each coronary artery (step S25).
続いて、マーカ発生部111は、責任血管特定部108により特定された梗塞責任血管を表すマーカのデータを発生させる(ステップS26)。具体的には、マーカ発生部111は、例えば図13に示すように、梗塞責任血管の輪郭を表すマーカm1を発生させる。 Subsequently, the marker generating unit 111 generates marker data representing the infarct responsible blood vessel specified by the responsible blood vessel specifying unit 108 (step S26). Specifically, the marker generator 111 generates a marker m1 representing the contour of the infarct responsible blood vessel, for example, as shown in FIG.
次に、FFR計算部109は、冠動脈解析部107により抽出された冠動脈毎に、冠動脈全体のFFR値を算出する(ステップS27)。 Next, the FFR calculation unit 109 calculates the FFR value of the entire coronary artery for each coronary artery extracted by the coronary artery analysis unit 107 (step S27).
続いて、マーカ発生部111は、FFR計算部109により算出された各冠動脈に対応したFFR値を表すマーカを発生させる(ステップS28)。具体的には、マーカ発生部111は、例えば図14に示すように、各冠動脈のFFR値の遷移を表すマーカm2を発生させる。 Subsequently, the marker generation unit 111 generates a marker representing an FFR value corresponding to each coronary artery calculated by the FFR calculation unit 109 (step S28). Specifically, for example, as shown in FIG. 14, the marker generating unit 111 generates a marker m2 representing a transition of the FFR value of each coronary artery.
しかる後、表示部111は、例えば図15に示すように、マーカ発生部111により発生した梗塞責任血管を表すマーカm1と各冠動脈のFFR値の遷移を表すマーカm2とをボリュームデータに由来した2次元画像g1に重畳させて表示する(ステップS29)。
なお、表示部111が表示する画像は、上記2次元画像g1だけでなく、例えば図16に示すように、ボリュームデータに由来した3次元画像g2や、他のモダリティから取得される画像にマーカを重畳させた画像であってもよい。
Thereafter, as shown in FIG. 15, for example, the display unit 111 generates a marker m1 representing the infarct responsible blood vessel generated by the marker generating unit 111 and a marker m2 representing the transition of the FFR value of each coronary artery from the volume data 2 The two-dimensional image g1 is superimposed and displayed (step S29).
Note that the image displayed by the display unit 111 is not limited to the two-dimensional image g1, but for example, as shown in FIG. 16, a marker is added to a three-dimensional image g2 derived from volume data or an image acquired from another modality. A superimposed image may be used.
なお、上記動作例の説明では、責任血管特定部108が梗塞責任血管を特定する場合について説明したが、例えば図17のフローチャートに示すように、責任血管特定部108が虚血責任血管を特定する場合においても、医用画像処理装置10は、CT装置20又はMRI装置30による心筋パフュージョン解析や、核医学診断装置40によるSPECT検査等により虚血領域を特定して(ステップS22’,S23’,S24’)、当該特定した虚血領域から虚血責任血管を特定し(ステップS25’)、当該特定した虚血責任血管を表すマーカを発生させる(ステップS26’)こと以外は、上記動作例と同様に動作する。 In the above description of the operation example, the case where the responsible blood vessel specifying unit 108 specifies the infarct responsible blood vessel has been described. For example, as shown in the flowchart of FIG. 17, the responsible blood vessel specifying unit 108 specifies the ischemic responsible blood vessel. Even in such a case, the medical image processing apparatus 10 specifies the ischemic region by myocardial perfusion analysis by the CT apparatus 20 or the MRI apparatus 30, SPECT examination by the nuclear medicine diagnosis apparatus 40, and the like (steps S22 ′, S23 ′, S24 ′), specifying the ischemic responsible blood vessel from the specified ischemic region (step S25 ′), and generating a marker representing the specified ischemic responsible blood vessel (step S26 ′). It operates in the same way.
以上説明した第2の実施形態によれば、CT装置20によるボリュームデータから心臓領域、心筋領域、冠動脈、狭窄部位及び心筋梗塞領域を抽出可能な心臓領域抽出部105、心筋解析部106及び冠動脈解析部107と、心筋解析部106による処理結果に基づいて梗塞責任血管を特定する責任血管特定部108と、梗塞責任血管や、FFR計算部109により算出されたFFR値を表すマーカをボリュームデータに由来した2次元画像又は3次元画像に重畳表示させる表示部111とを備えた構成により、梗塞責任血管を表すマーカが重畳表示された部分のFFR値の信頼度が低い旨を医師に対して提示することができる。 According to the second embodiment described above, the heart region extraction unit 105, the myocardial analysis unit 106, and the coronary artery analysis that can extract the heart region, the myocardial region, the coronary artery, the stenosis site, and the myocardial infarction region from the volume data obtained by the CT apparatus 20. Unit 107, responsible blood vessel specifying unit 108 for identifying the infarct responsible blood vessel based on the processing result by the myocardial analysis unit 106, the infarct responsible blood vessel, and a marker representing the FFR value calculated by the FFR calculating unit 109 are derived from the volume data. The configuration including the display unit 111 that superimposes and displays the two-dimensional image or the three-dimensional image is provided to the doctor that the reliability of the FFR value of the portion where the marker representing the infarct responsible blood vessel is superimposed and displayed is low. be able to.
また、本実施形態では、FFR計算部109がシミュレーションベースでFFR値を算出する、つまり、プレッシャーワイヤー等、侵襲性をもつものを用いないため、検査時に患者にかかる負担を低減させることができる。 Further, in this embodiment, the FFR calculation unit 109 calculates the FFR value on a simulation basis, that is, since an invasive device such as a pressure wire is not used, the burden on the patient during the examination can be reduced.
[第3の実施形態]
次に、第3の実施形態に係る医用画像処理装置について、上記した図11を用いて説明する。本実施形態は、第2の実施形態に示す医用画像処理装置10に、冠動脈内の狭窄部位が治療対象狭窄であるか又は治療非対象狭窄であるかを判定する機能を付加したものである。なお、以下では、図18のフローチャートと図19の模式図とを参照しながら、第2の実施形態とは異なる機能についてのみ説明する。つまり、ステップS21乃至S28の処理は、上記した第2の実施形態と同様であるため、ここでは詳細な説明は省略し、以下では、ステップS30乃至S36の処理について主に説明する。
[Third Embodiment]
Next, a medical image processing apparatus according to the third embodiment will be described with reference to FIG. In this embodiment, the medical image processing apparatus 10 shown in the second embodiment is added with a function of determining whether a stenosis site in a coronary artery is a treatment target stenosis or a non-treatment stenosis. Hereinafter, only functions different from those in the second embodiment will be described with reference to the flowchart in FIG. 18 and the schematic diagram in FIG. 19. That is, the processes in steps S21 to S28 are the same as those in the above-described second embodiment, and thus detailed description thereof is omitted here. The processes in steps S30 to S36 are mainly described below.
冠動脈解析部107は、ステップS21の処理において抽出した狭窄部位毎に、当該狭窄部位が梗塞責任血管内に位置する狭窄であるか否かを判定する(ステップS30)。 The coronary artery analysis unit 107 determines, for each stenosis site extracted in step S21, whether or not the stenosis site is a stenosis located in the infarct responsible blood vessel (step S30).
ステップS30の処理による判定の結果が梗塞責任血管内に位置する狭窄である旨を示す場合(ステップS30のYes)には、心筋解析部106は、ステップS24の処理において特定した心筋梗塞領域内に生存心筋が有るか否かを判定する(ステップS31)。
具体的には、心筋解析部106は、MRI装置30による遅延造影等から当該心筋梗塞領域が心筋の厚さの半分にまで到達しているか否かを判定し、当該判定の結果が到達している旨を示す場合には生存心筋が無いとみなし、当該判定の結果が否を示す場合には生存心筋が有るとみなす。なお、ステップS31の処理による判定の結果が否を示す場合(ステップS31のNo)には、後述するステップS35の処理に進む。
If the result of the determination in step S30 indicates that the stenosis is located in the infarct responsible blood vessel (Yes in step S30), the myocardial analysis unit 106 is within the myocardial infarction region identified in the process in step S24. It is determined whether or not there is a surviving myocardium (step S31).
Specifically, the myocardial analyzer 106 determines whether or not the myocardial infarction region has reached half the thickness of the myocardium from delayed contrast imaging by the MRI apparatus 30 and the result of the determination has reached. If the result indicates that there is no living myocardium, if the result of the determination indicates no, it is considered that there is a living myocardium. If the result of the determination in step S31 indicates NO (No in step S31), the process proceeds to step S35 described later.
ステップS31の処理による判定の結果が生存心筋が有る旨を示す場合(ステップS31のYes)には、マーカ発生部111は、生存心筋が有る梗塞責任血管を表すマーカと治療対象狭窄を表すマーカとを発生させる(ステップS32)。具体的には、マーカ発生部111は、例えば図19に示すように、生存心筋が有る梗塞責任血管の輪郭を表すマーカm3と、治療対象狭窄を表すマーカm4とを発生させる。 When the result of the determination in step S31 indicates that there is a surviving myocardium (Yes in step S31), the marker generating unit 111 includes a marker indicating the infarct responsible blood vessel having the surviving myocardium and a marker indicating the treatment target stenosis. Is generated (step S32). Specifically, for example, as shown in FIG. 19, the marker generating unit 111 generates a marker m3 that represents the contour of the infarct responsible blood vessel having a living myocardium and a marker m4 that represents the treatment target stenosis.
ここで、ステップS30の処理による判定の結果が否を示す場合(ステップS30のNo)には、冠動脈解析部107は、ステップS27の処理において算出された各冠動脈のFFR値のうち当該狭窄部位に対応したFFR値が閾値以下であるか否かを判定する(ステップS33)。 Here, when the result of the determination in the process of step S30 indicates NO (No in step S30), the coronary artery analysis unit 107 determines that the stenosis site is included in the FFR value of each coronary artery calculated in the process of step S27. It is determined whether or not the corresponding FFR value is equal to or less than a threshold value (step S33).
ステップS33の処理による判定の結果が閾値以下である旨を示す場合(ステップS33のYes)には、マーカ発生部111は、治療対象狭窄を表すマーカを発生させる(ステップS34)。即ち、マーカ発生部111は、図19に示したマーカm4に相当するマーカを発生させる。 When the result of the determination in step S33 indicates that the result is equal to or less than the threshold (Yes in step S33), the marker generating unit 111 generates a marker representing the treatment target stenosis (step S34). That is, the marker generator 111 generates a marker corresponding to the marker m4 shown in FIG.
ステップS33の処理による判定の結果が否を示す場合(ステップS33のNo)には、マーカ発生部111は、治療非対象狭窄を表すマーカを発生させる(ステップS35)。具体的には、マーカ発生部111は、例えば図19に示すように、治療非対象狭窄を表すマーカm5を発生させる。 If the result of the determination in step S33 indicates NO (No in step S33), the marker generating unit 111 generates a marker indicating treatment non-target stenosis (step S35). Specifically, for example, as shown in FIG. 19, the marker generating unit 111 generates a marker m5 representing non-target stenosis.
しかる後、表示部111は、マーカ発生部111により発生した梗塞責任血管を表すマーカm1、各冠動脈のFFR値の遷移を表すマーカm2、生存心筋が有る梗塞責任血管を表すマーカm3、治療対象狭窄を表すマーカm4、及び治療非対象狭窄を表すマーカm5をボリュームデータに由来した3次元画像g3に重畳させて表示する(ステップS36)。なお、表示部111が表示する画像は、上記3次元画像g3だけでなく、ボリュームデータに由来した2次元画像や、他のモダリティから取得される画像にマーカを重畳させた画像であってもよい。また、梗塞責任血管を表すマーカm1と生存心筋が有る梗塞責任血管を表すマーカm3との両方が画像上の同じ位置に重畳される場合、マーカm3が優先的に表示されるものとする。 Thereafter, the display unit 111 displays a marker m1 representing the infarct responsible blood vessel generated by the marker generating unit 111, a marker m2 representing the transition of the FFR value of each coronary artery, a marker m3 representing the infarct responsible blood vessel having a living myocardium, and a treatment target stenosis. And a marker m5 representing non-target stenosis are superimposed and displayed on the three-dimensional image g3 derived from the volume data (step S36). Note that the image displayed by the display unit 111 is not limited to the three-dimensional image g3, but may be a two-dimensional image derived from volume data, or an image in which a marker is superimposed on an image acquired from another modality. . In addition, when both the marker m1 representing the infarct responsible blood vessel and the marker m3 representing the infarct responsible blood vessel having the surviving myocardium are superimposed at the same position on the image, the marker m3 is preferentially displayed.
以上説明した第3の実施形態によれば、心筋梗塞領域に生存心筋が有るか否かや、冠動脈内の狭窄部位が治療対象狭窄であるか又は治療非対象狭窄であるかを判定する心筋解析部106及び冠動脈解析部107と、生存心筋が有る梗塞責任血管を表すマーカや、治療対象狭窄を表すマーカ、ならびに治療非対象狭窄を表すマーカをボリュームデータに由来した3次元画像や2次元画像に重畳表示させる表示部111とを備えた構成により、第1の実施形態に比べて、より多くの情報を医師に対して提示することができる。 According to the third embodiment described above, myocardial analysis for determining whether or not there is a living myocardium in the myocardial infarction region, and whether the stenotic site in the coronary artery is a treatment target stenosis or a non-treatment stenosis 106 and coronary artery analysis unit 107, a marker representing an infarct responsible blood vessel having a surviving myocardium, a marker representing treatment target stenosis, and a marker representing treatment non-target stenosis in a three-dimensional image or two-dimensional image derived from volume data Compared to the first embodiment, more information can be presented to the doctor by the configuration including the display unit 111 that displays the images in a superimposed manner.
以上説明した第2及び第3の実施形態の少なくとも一方によれば、FFR値の信頼度が低い旨や、生存心筋の有無、ならびに狭窄部位が治療するに適した狭窄であるか否かを医師に対して提示することができるため、ヒューマンエラーの可能性を低減させることができる。 According to at least one of the second and third embodiments described above, doctors can determine whether the reliability of the FFR value is low, whether there is a surviving myocardium, and whether the stenotic site is a stenosis suitable for treatment. Therefore, the possibility of human error can be reduced.
なお、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。 In addition, although several embodiment of this invention was described, these embodiment is shown as an example and is not intending limiting the range of invention. These novel embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the scope of the invention. These embodiments and modifications thereof are included in the scope and gist of the invention, and are included in the invention described in the claims and the equivalents thereof.
1…医用画像処理システム、10…医用画像処理装置、20…CT装置、30…PACS、40…ネットワーク、101…画像記憶部、102…テリトリーマップ記憶部、103…通信インターフェース、104…制御部、105…心臓領域抽出部、106…心筋解析部、107…冠動脈解析部、108…責任血管特定部、109…FFR計算部、110…責任狭窄特定部、111…マーカ発生部、112…表示部。 DESCRIPTION OF SYMBOLS 1 ... Medical image processing system, 10 ... Medical image processing apparatus, 20 ... CT apparatus, 30 ... PACS, 40 ... Network, 101 ... Image storage part, 102 ... Territory map storage part, 103 ... Communication interface, 104 ... Control part, DESCRIPTION OF SYMBOLS 105 ... Heart region extraction part, 106 ... Myocardial analysis part, 107 ... Coronary artery analysis part, 108 ... Responsible blood vessel specification part, 109 ... FFR calculation part, 110 ... Responsible stenosis specification part, 111 ... Marker generation part, 112 ... Display part.
本発明の実施形態は、医用画像処理装置、医用画像処理システムおよび医用画像処理方法に関する。 Embodiments described herein relate generally to a medical image processing apparatus , a medical image processing system, and a medical image processing method .
Claims (1)
前記取得された複数の冠動脈の形状に基づいて、前記取得された各冠動脈の圧較差を計算する計算部と、
前記画像に描出された心臓に関する心筋梗塞領域または虚血領域を取得する第2取得部と、
前記取得された心筋梗塞領域または虚血領域を前記取得された各冠動脈と支配域とを関係付ける支配マップに照会することにより当該心筋梗塞領域または虚血領域の責任血管を特定する第1特定部と、
前記特定された責任血管内の狭窄部位に対応した前記圧較差に基づいて責任狭窄を特定する第2特定部と、
前記特定された責任狭窄が描出された画像を、当該責任狭窄を示す情報と共に表示する表示部と、
を具備する医用画像処理装置。 A first acquisition unit that acquires a plurality of coronary arteries depicted in image data relating to the heart, and obtains at least one stenosis site depicted in each of the obtained coronary arteries;
Based on the acquired shapes of the plurality of coronary arteries, a calculation unit that calculates a pressure difference of each acquired coronary artery,
A second acquisition unit for acquiring a myocardial infarction region or an ischemic region related to the heart depicted in the image;
A first specifying unit that specifies a responsible blood vessel of the myocardial infarction region or ischemic region by referring to the control map that associates the acquired coronary artery and the control region with the acquired myocardial infarction region or ischemic region When,
A second identifying unit that identifies the responsible stenosis based on the pressure difference corresponding to the identified stenotic site in the responsible blood vessel;
A display unit that displays an image depicting the identified responsible stenosis together with information indicating the responsible stenosis;
A medical image processing apparatus comprising:
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012263565 | 2012-11-30 | ||
JP2012263565 | 2012-11-30 | ||
JP2012263566 | 2012-11-30 | ||
JP2012263566 | 2012-11-30 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017195297A Division JP6553147B2 (en) | 2012-11-30 | 2017-10-05 | Medical image processing apparatus and medical image processing method |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020136235A Division JP7024023B2 (en) | 2012-11-30 | 2020-08-12 | Medical image processing equipment, medical image processing system and medical image processing method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2019202154A true JP2019202154A (en) | 2019-11-28 |
JP6751184B2 JP6751184B2 (en) | 2020-09-02 |
Family
ID=61248060
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017195297A Active JP6553147B2 (en) | 2012-11-30 | 2017-10-05 | Medical image processing apparatus and medical image processing method |
JP2019119706A Active JP6751184B2 (en) | 2012-11-30 | 2019-06-27 | Medical image processing apparatus, medical image processing system, and medical image processing method |
JP2020136235A Active JP7024023B2 (en) | 2012-11-30 | 2020-08-12 | Medical image processing equipment, medical image processing system and medical image processing method |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017195297A Active JP6553147B2 (en) | 2012-11-30 | 2017-10-05 | Medical image processing apparatus and medical image processing method |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020136235A Active JP7024023B2 (en) | 2012-11-30 | 2020-08-12 | Medical image processing equipment, medical image processing system and medical image processing method |
Country Status (1)
Country | Link |
---|---|
JP (3) | JP6553147B2 (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6553147B2 (en) * | 2012-11-30 | 2019-07-31 | キヤノンメディカルシステムズ株式会社 | Medical image processing apparatus and medical image processing method |
JP7314183B2 (en) * | 2018-06-15 | 2023-07-25 | パイ メディカル イメージング ビー ヴイ | Method and apparatus for quantitative hemodynamic flow analysis |
US20220015730A1 (en) * | 2018-11-28 | 2022-01-20 | Koninklijke Philips N.V. | Most relevant x-ray image selection for hemodynamic simulation |
WO2021141135A1 (en) * | 2020-01-09 | 2021-07-15 | 学校法人東京女子医科大学 | Functional ischemia detection technique by coronary artery ct4d flow imaging |
CN112801999B (en) * | 2021-02-05 | 2023-06-23 | 深圳睿心智能医疗科技有限公司 | Method and device for determining heart coronary artery dominance |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009028515A (en) * | 2007-07-05 | 2009-02-12 | Toshiba Corp | Medical image processor and medical image diagnostic apparatus |
US20100298719A1 (en) * | 2007-10-31 | 2010-11-25 | Samuel Alberg Kock | Method for calculating pressures in a fluid stream through a tube section, especially a blood vessel with atherosclerotic plaque |
JP2011115481A (en) * | 2009-12-07 | 2011-06-16 | Hitachi Medical Corp | Image analyzing device and method, and image analyzing program |
JP2011182899A (en) * | 2010-03-08 | 2011-09-22 | Fujifilm Corp | Diagnosis assisting apparatus, coronary artery analyzing program, and coronary artery analyzing method |
JP2012005636A (en) * | 2010-06-24 | 2012-01-12 | Toshiba Corp | Medical image diagnosis apparatus and medical diagnosis support method |
JP2012081254A (en) * | 2010-09-15 | 2012-04-26 | Toshiba Corp | Medical image processing apparatus and medical image processing method |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8315812B2 (en) * | 2010-08-12 | 2012-11-20 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
JP6553147B2 (en) | 2012-11-30 | 2019-07-31 | キヤノンメディカルシステムズ株式会社 | Medical image processing apparatus and medical image processing method |
WO2014084382A1 (en) | 2012-11-30 | 2014-06-05 | 株式会社 東芝 | Medical-image processing device |
-
2017
- 2017-10-05 JP JP2017195297A patent/JP6553147B2/en active Active
-
2019
- 2019-06-27 JP JP2019119706A patent/JP6751184B2/en active Active
-
2020
- 2020-08-12 JP JP2020136235A patent/JP7024023B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009028515A (en) * | 2007-07-05 | 2009-02-12 | Toshiba Corp | Medical image processor and medical image diagnostic apparatus |
US20100298719A1 (en) * | 2007-10-31 | 2010-11-25 | Samuel Alberg Kock | Method for calculating pressures in a fluid stream through a tube section, especially a blood vessel with atherosclerotic plaque |
JP2011115481A (en) * | 2009-12-07 | 2011-06-16 | Hitachi Medical Corp | Image analyzing device and method, and image analyzing program |
JP2011182899A (en) * | 2010-03-08 | 2011-09-22 | Fujifilm Corp | Diagnosis assisting apparatus, coronary artery analyzing program, and coronary artery analyzing method |
JP2012005636A (en) * | 2010-06-24 | 2012-01-12 | Toshiba Corp | Medical image diagnosis apparatus and medical diagnosis support method |
JP2012081254A (en) * | 2010-09-15 | 2012-04-26 | Toshiba Corp | Medical image processing apparatus and medical image processing method |
Also Published As
Publication number | Publication date |
---|---|
JP2021000455A (en) | 2021-01-07 |
JP6553147B2 (en) | 2019-07-31 |
JP2018027323A (en) | 2018-02-22 |
JP6751184B2 (en) | 2020-09-02 |
JP7024023B2 (en) | 2022-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6334902B2 (en) | Medical image processing device | |
US20220296106A1 (en) | Medical image processing apparatus | |
JP6751184B2 (en) | Medical image processing apparatus, medical image processing system, and medical image processing method | |
EP3169237B1 (en) | Stenosis assessment | |
EP3281135B1 (en) | System and method for vascular tree generation using patient-specific structural and functional data, and joint prior information | |
US10052032B2 (en) | Stenosis therapy planning | |
JP2018509971A (en) | Method and apparatus for quantitative flow analysis | |
JP2015231524A (en) | Image processor, image processing method, and program | |
JP6139116B2 (en) | Medical image processing device | |
JP6362851B2 (en) | Blood vessel analysis device, blood vessel analysis program, and blood vessel analysis device operating method | |
JP2021104337A (en) | Estimating endoluminal path of endoluminal device along lumen | |
JP5972768B2 (en) | Medical image processing device | |
US20120046725A1 (en) | Method for Providing an Aid for Use in the Therapeutic Treatment of a Somatic Subject | |
JP6751178B2 (en) | Medical image processing apparatus, medical image processing method and recording medium | |
JP6530043B2 (en) | Medical image processing apparatus, medical image processing method and recording medium | |
JP7041446B2 (en) | Medical image processing methods, medical image processing equipment and medical image processing systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190723 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190723 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200630 |
|
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: 20200714 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200813 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6751184 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |