JP6671482B2 - Ｃｐｒ画像生成装置、方法およびプログラム - Google Patents
Ｃｐｒ画像生成装置、方法およびプログラム Download PDFInfo
- Publication number
- JP6671482B2 JP6671482B2 JP2018536988A JP2018536988A JP6671482B2 JP 6671482 B2 JP6671482 B2 JP 6671482B2 JP 2018536988 A JP2018536988 A JP 2018536988A JP 2018536988 A JP2018536988 A JP 2018536988A JP 6671482 B2 JP6671482 B2 JP 6671482B2
- Authority
- JP
- Japan
- Prior art keywords
- cross
- section
- region
- image
- interest
- 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
Links
- 238000000034 method Methods 0.000 title claims description 29
- 230000003902 lesion Effects 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 5
- 210000004351 coronary vessel Anatomy 0.000 description 46
- 210000004165 myocardium Anatomy 0.000 description 21
- 238000010586 diagram Methods 0.000 description 14
- 210000000621 bronchi Anatomy 0.000 description 12
- 238000003745 diagnosis Methods 0.000 description 10
- 210000004204 blood vessel Anatomy 0.000 description 9
- 206010028980 Neoplasm Diseases 0.000 description 8
- 238000000605 extraction Methods 0.000 description 7
- 239000000284 extract Substances 0.000 description 6
- 210000002429 large intestine Anatomy 0.000 description 6
- 238000002591 computed tomography Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000002600 positron emission tomography Methods 0.000 description 5
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 208000010125 myocardial infarction Diseases 0.000 description 3
- 230000036262 stenosis Effects 0.000 description 3
- 208000037804 stenosis Diseases 0.000 description 3
- 238000001574 biopsy Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000004195 computer-aided diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000001839 endoscopy Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000004217 heart function Effects 0.000 description 1
- 210000000936 intestine Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 201000005202 lung cancer Diseases 0.000 description 1
- 208000020816 lung neoplasm Diseases 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 238000009751 slip forming Methods 0.000 description 1
- 230000002966 stenotic effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/466—Displaying means of special interest adapted to display 3D data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2215/00—Indexing scheme for image rendering
- G06T2215/06—Curved planar reformation of 3D line structures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
- G06T2219/008—Cut plane or projection plane definition
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Vascular Medicine (AREA)
- Pulmonology (AREA)
- Human Computer Interaction (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Architecture (AREA)
- Physiology (AREA)
- Dentistry (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
本発明は、血管等の構造物の長手方向に沿った断面を表すＣＰＲ（Curved Planer Reconstruction）画像を生成するＣＰＲ画像生成装置、方法およびプログラムに関するものである。
従来、被写体の血管、腸、気管支、動脈等の管腔を有する管状構造物の病変を観察するために、ＣＴ（Computed Tomography）装置およびＭＲＩ（Magnetic Resonance Imaging）装置等により得られた物体の３次元画像から３次元構造のイメージをボリュームレンダリングにより取得し、３次元画像から内視鏡を模してあたかもその管腔内に入って見ているように管腔内面を観察できるモードを有する医用画像表示装置が知られている。さらに、３次元画像処理としてＭＩＰ（Maximum Intensity Projection）処理、ＭｉｎＩＰ（Minimum Intensity Projection）処理、ＭＰＲ（Multi Planar Reconstruction）処理、ＣＰＲ（Curved Planer Reconstruction）処理等が知られている。特に、ＣＰＲ処理は３次元画像における任意の曲面を指定し、指定した曲面に沿った３次元画像を２次元画像に再構成してＣＰＲ画像を生成するものである。これにより、管状構造物の内壁面の形態について、長手方向の断面を一画面で表示することができる（例えば特許文献１参照）。
また、上記のＣＴ装置等から得られた３次元画像の他に、被写体の機能を示す機能画像も用いられている。機能画像としては、例えば単光子放出断層撮影（Single Photon Emission Tomography）により取得されたＳＰＥＣＴ画像、および陽電子放出断層撮影（Positron Emission Tomography）により取得されたＰＥＴ画像等が知られている。このような機能画像をＣＰＲ画像と併せて表示することにより、例えば冠動脈を診断対象とする場合、冠動脈の状態と心臓の機能の状態とを同時に確認することができる。
一方、血管等のように蛇行しながら延びる管状構造物においては、ＣＰＲ画像により表された管状構造物の縦断面は、管状構造物の蛇行に沿った湾曲形状またはねじれ形状を有する切断曲面となる。このため、ＣＰＲ画像に慣れないユーザには、ＣＰＲ画像が管状構造物の各位置をどの向きに切断した曲面を表しているのかを判断することが難しかった。このため、縦断面を表す指標をＣＰＲ画像とともに表示することにより、切断曲面と管状構造物との位置関係を示す手法が提案されている（特許文献２参照）。また、特許文献２には、血管および気管支のように分岐を含む管状構造物において、全ての分岐を含むマルチパスＣＰＲ画像を生成する手法が提案されている。特許文献２に記載されたマルチパスＣＰＲ画像を生成する手法は、分岐間の区画毎に切断曲面を設定し、各区画のＣＰＲ画像を１つの画像上に生成する手法である。なお、切断曲面の設定は、切断曲面の各位置の法線方向等の情報に基づいて行われる。これにより、管状構造物の各分岐の構造を含むＣＰＲ画像を生成することができるため、管状構造物の分岐の構造を一画面で表示することができる。
しかしながら、特許文献２に記載された手法では、管状構造物の分岐位置において切断曲面が連続しない。このため、生成されるマルチパスＣＰＲ画像において、分岐位置の画像に境界が生じてしまう。
一方、心筋梗塞等の疾患は、冠動脈の狭窄部分に血栓が詰まり、詰まった冠動脈の下流側に血液が流れなくなる結果、心筋が活動しなくなることにより発生する。心筋梗塞等が発生する可能性は、血管の狭窄率が高いほど高くなる。血管が狭窄する原因は血管壁に形成されるプラークである。このため、血管においてプラークがどの程度発生しているかをＣＰＲ画像において観察することは、極めて重要である。また、肺癌等の肺における腫瘍等の病変の確認のために、気管支に内視鏡を挿入して病変の生検を行う場合があるが、この場合、病変に到達するまでの経路および気管支と病変の位置との関係をＣＰＲ画像において観察することが重要である。
しかしながら、これらの心筋、プラークおよび病変等の関心領域は、管状構造物の１つの縦断面上にのみ存在するものではなく、管状構造物の各断面の中心または重心を結んだ管状構造物の中心線（以下、芯線とする）に垂直な断面において、芯線を中心とする放射状の任意の位置に存在する。このため、従来のように１つの縦切断面のみを設定してＣＰＲ画像を生成したのでは、ＣＰＲ画像において関心領域を良好に観察することができない場合がある。
本発明は上記事情に鑑みなされたものであり、分岐、心筋、プラークおよび病変等の関心領域を含むようにＣＰＲ画像を生成することを目的とする。
本発明によるＣＰＲ画像生成装置は、撮影により取得された３次元画像から、対象構造物を抽出する構造物抽出手段と、
対象構造物の基準線上の各点において、基準線に垂直な断面を設定する断面設定手段と、
各断面上に関心領域が存在する場合、基準線および関心領域を含む、各断面における１つの切断面を決定する切断面決定手段と、
各断面における切断面を含むＣＰＲ画像を生成する画像生成手段とを備えたことを特徴とするものである。
対象構造物の基準線上の各点において、基準線に垂直な断面を設定する断面設定手段と、
各断面上に関心領域が存在する場合、基準線および関心領域を含む、各断面における１つの切断面を決定する切断面決定手段と、
各断面における切断面を含むＣＰＲ画像を生成する画像生成手段とを備えたことを特徴とするものである。
ＣＰＲ画像は、対象構造物の芯線に沿った断面を表す画像であるが、対象構造物が例えば血管または気管支のように複数の分岐を含む場合、分岐毎に複数の芯線が存在する。このような場合、ＣＰＲ画像を生成するためには、１つの芯線を決定する必要がある。「基準線」とは、対象構造物が複数の芯線を含む場合に、ＣＰＲ画像を生成するための１つの芯線を意味する。なお、対象構造物が１つの芯線のみを含む場合、その１つの芯線が基準線となる。
なお、本発明によるＣＰＲ画像生成装置においては、対象構造物は管状構造物であってもよい。
また、本発明によるＣＰＲ画像生成装置においては、関心領域は、対象構造物の分岐、対象構造物に隣接する他の構造物、対象構造物内の他の構造物、および病変の少なくとも１つであってもよい。
「対象構造物に隣接する他の構造物」としては、例えば冠動脈を対象構造物とした場合における心筋等を挙げることができる。「対象構造物内の他の構造物」とは、例えば冠動脈を対象構造物とした場合における狭窄部分に形成されたプラーク等を挙げることができる。
また、本発明によるＣＰＲ画像生成装置においては、画像生成手段は、各断面における切断面が連続する切断曲面を生成し、切断曲面を含むＣＰＲ画像を生成するものであってもよい。
また、本発明によるＣＰＲ画像生成装置においては、各断面に複数の関心領域が存在する場合、切断面決定手段は、１つの関心領域を選択するものであってもよい。
また、本発明によるＣＰＲ画像生成装置においては、ＣＰＲ画像を表示する表示手段をさらに備えるものとしてもよい。
この場合、表示手段は、各断面における切断面の基準となる断面の切断面からの基準線を中心とした角度を表す情報を表示するものであってもよい。
本発明によるＣＰＲ画像生成方法は、撮影により取得された３次元画像から、対象構造物を抽出し、
対象構造物の基準線上の各点において、基準線に垂直な断面を設定し、
各断面上に関心領域が存在する場合、基準線および関心領域を含む、各断面における１つの切断面を決定し、
各断面における切断面を含むＣＰＲ画像を生成することを特徴とするものである。
対象構造物の基準線上の各点において、基準線に垂直な断面を設定し、
各断面上に関心領域が存在する場合、基準線および関心領域を含む、各断面における１つの切断面を決定し、
各断面における切断面を含むＣＰＲ画像を生成することを特徴とするものである。
なお、本発明によるＣＰＲ画像生成方法をコンピュータに実行させるためのプログラムとして提供してもよい。
本発明によれば、３次元画像から対象構造物および対象構造物の基準線が抽出され、基準線上の各点において、基準線に垂直な断面が設定される。そして、各断面において関心領域が存在する場合には、基準線および対象構造物に関連する関心領域を含む、各断面における１つの切断面が決定され、各断面における１つの切断面を含むＣＰＲ画像が生成される。このため、各断面において基準線を中心とする放射状の任意の方向に存在する関心領域を含むようにＣＰＲ画像を生成することができる。したがって、ＣＰＲ画像において、対象構造物に含まれる複数の関心領域を良好に観察することができる。また、各断面において１つの切断面を決定しているため、例えば関心領域が分岐である場合において、ＣＰＲ画像に複数の分岐を含ませることができ、かつ分岐間においてＣＰＲ画像が連続するようになる。このため、境界がない自然な印象のＣＰＲ画像を生成することができる。
以下、図面を参照して本発明の実施形態について説明する。図１は、本発明の第１の実施形態によるＣＰＲ画像生成装置を適用した、診断支援システムの概要を示すハードウェア構成図である。図１に示すように、診断支援システムでは、本実施形態によるＣＰＲ画像生成装置１、３次元画像撮影装置２、および画像保管サーバ３が、ネットワーク４を経由して通信可能な状態で接続されている。そして、診断支援システムにおいては、被検体の診断対象となる部位についてのＣＰＲ画像を、ＣＰＲ画像生成装置１において生成するものである。
３次元画像撮影装置２は、被検体の診断対象となる部位を撮影することにより、その部位を表す３次元画像を生成する装置であり、具体的には、ＣＴ装置、ＭＲＩ装置、およびＰＥＴ（Positron Emission Tomography)装置等である。この３次元画像撮影装置２により生成された３次元画像は画像保管サーバ３に送信され、保存される。なお、本実施形態においては、被検体の診断対象部位は冠動脈であり、３次元画像撮影装置２はＣＴ装置であり、被検体の胸部の３次元画像が生成されるものとする。
画像保管サーバ３は、各種データを保存して管理するコンピュータであり、大容量外部記憶装置およびデータベース管理用ソフトウェアを備えている。画像保管サーバ３は、有線あるいは無線のネットワーク４を介して他の装置と通信を行い、画像データ等を送受信する。具体的には３次元画像撮影装置２において生成された３次元画像等の画像データをネットワーク経由で取得し、大容量外部記憶装置等の記録媒体に保存して管理する。なお、画像データの格納形式およびネットワーク４経由での各装置間の通信は、ＤＩＣＯＭ（Digital Imaging and COmmunication in Medicine）等のプロトコルに基づいている。
ＣＰＲ画像生成装置１は、１台のコンピュータに、本発明のＣＰＲ画像生成プログラムをインストールしたものである。コンピュータは、診断を行う医師が直接操作するワークステーションまたはパーソナルコンピュータでもよいし、それらとネットワークを介して接続されたサーバコンピュータでもよい。ＣＰＲ画像生成プログラムは、ＤＶＤ（Digital Versatile Disc）あるいはＣＤ−ＲＯＭ（Compact Disk Read Only Memory）等の記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。または、ネットワークに接続されたサーバコンピュータの記憶装置、もしくはネットワークストレージに、外部からアクセス可能な状態で記憶され、要求に応じて医師が使用するコンピュータにダウンロードされ、インストールされる。
図２は、コンピュータにＣＰＲ画像生成プログラムをインストールすることにより実現されるＣＰＲ画像生成装置の概略構成を示す図である。図２に示すように、ＣＰＲ画像生成装置１は、標準的なワークステーションの構成として、ＣＰＵ（Central Processing Unit）１１、メモリ１２およびストレージ１３を備えている。また、ＣＰＲ画像生成装置１には、ディスプレイ１４および、マウス等の入力部１５が接続されている。
ストレージ１３には、ネットワーク４を経由して画像保管サーバ３から取得した、同一被検体についての最新および前回の診断時における過去の３次元画像、並びに処理に必要な情報を含む各種情報が記憶されている。なお、本実施形態においては、同一被検体についての胸部を対象部位とする３次元画像Ｇ０が記憶されているものとする。
また、メモリ１２には、ＣＰＲ画像生成プログラムが記憶されている。ＣＰＲ画像生成プログラムは、ＣＰＵ１１に実行させる処理として、３次元画像撮影装置２が取得した、３次元画像Ｇ０を取得する画像取得処理、３次元画像Ｇ０から対象構造物を抽出する構造物抽出処理、対象構造物の芯線の１つである基準線上の各点において、基準線に垂直な断面を設定する断面設定処理、各断面上に関心領域が存在する場合、基準線および関心領域を含む、各断面における１つの切断面を決定する切断面決定処理、並びに各断面における切断面を含むＣＰＲ画像を生成する画像生成処理を規定する。
そして、ＣＰＵ１１がプログラムに従いこれらの処理を実行することにより、コンピュータは、画像取得部２１、構造物抽出部２２、断面設定部２３、切断面決定部２４、および画像生成部２５として機能する。なお、ＣＰＲ画像生成装置１は、画像取得処理、構造物抽出処理、断面設定処理、切断面決定処理、および画像生成処理をそれぞれ行う複数のプロセッサまたは処理回路を備えるものであってもよい。
画像取得部２１は、診断の対象部位である冠動脈すなわち心臓を含む胸部の３次元画像Ｇ０を画像保管サーバ３から取得する。なお、３次元画像Ｇ０が既にストレージ１３に記憶されている場合には、画像取得部２１は、ストレージ１３から３次元画像Ｇ０を取得するようにしてもよい。
構造物抽出部２２は、例えば、特開２０１０−２００９２５号公報および特開２０１０−２２０７４２号公報等に記載された手法により、３次元画像Ｇ０から冠動脈領域を抽出する。この方法では、まず、ボリュームデータを構成するボクセルデータの値に基づいて、冠動脈の芯線を構成する複数の候補点の位置と主軸方向を算出する。もしくは、ボリュームデータについてヘッセ行列を算出し、算出されたヘッセ行列の固有値を解析することにより、冠動脈の芯線を構成する複数の候補点の位置情報と主軸方向を算出する。そして、候補点周辺のボクセルデータについて冠動脈らしさを表す特徴量を算出し、算出された特徴量に基づいてそのボクセルデータが冠動脈領域を表すものであるか否かを判別する。特徴量に基づく判別は、マシンラーニングにより予め取得された評価関数に基づいて行う。これにより、ボリュームデータから冠動脈領域３０が抽出される。図３に、抽出された冠動脈領域３０の一部を示す。
断面設定部２３は、冠動脈領域３０の基準線上の各点において、基準線に垂直な断面を設定する。上記特開２０１０−２００９２５号公報および特開２０１０−２２０７４２号公報に記載された手法では、冠動脈領域３０を抽出する過程において、冠動脈の芯線が設定される。また、芯線を構成する候補点のそれぞれについて、位置と主軸方向が算出される。このため、各候補点において、主軸方向と垂直な断面（直交断面）を、芯線に垂直な断面として設定することができる。
ここで、本実施形態において抽出される冠動脈領域３０は、図３に示すように分岐を有するため、芯線は分岐において２つに分かれる。ＣＰＲ画像を生成するためには、１つの芯線を決定する必要がある。このため、断面設定部２３は、２つの芯線のうち、ＣＰＲ画像を生成するための芯線を基準線に決定し、基準線において基準線に垂直な断面を設定する。なお、基準線の決定は、抽出された冠動脈領域をディスプレイ１４に表示して、操作者の入力部１５からの指示を受け付けることにより決定する。本実施形態においては、図３において上側から左側に延びる冠動脈領域３０Ａの芯線が基準線Ｌ０に決定されるものとする。図３において基準線Ｌ０を実線で示す。また、図３において分岐から下方に延びる冠動脈領域３０Ｂの芯線Ｃ０を破線で示す。また、図３においては、基準線Ｌ０上の候補点を黒丸で示す。
断面設定部２３は、基準線Ｌ０上の各点、すなわち各候補点において、基準線Ｌ０に垂直な断面を設定する。なお、図３においては、説明のために候補点を間引いて５つの断面Ｐｉ、Ｐｉ＋ａ、Ｐｉ＋ｂ、Ｐｉ＋ｃ、Ｐｉ＋ｄを設定した状態を示している。
切断面決定部２４は、断面設定部２３が基準線Ｌ０に設定した各断面Ｐｋ（ｋ＝１〜ｎ：ｎは候補点の数）上に関心領域が存在する場合、基準線Ｌ０および関心領域を含む、各断面Ｐｋにおける１つの切断面Ｃｋを決定する。このため、切断面決定部２４は、各断面Ｐｋにおいて関心領域を設定する。なお、関心領域としては、各断面Ｐｋにおいてプラーク３２が形成されている領域および心筋３３の領域とすることができる。また、冠動脈領域３０における分岐３４の位置においては、分岐３４を関心領域とすることができる。なお、関心領域の設定のために、切断面決定部２４は、３次元画像Ｇ０から各断面Ｐｋの断面画像ＰＧｋを生成する。
ここで、切断面決定部２４は、図４に示すように、各断面画像ＰＧｋにおいて、基準線Ｌ０を中心とした複数（ここでは４方向）の冠動脈領域３０の内腔３１の径ｄ１〜ｄ４を算出する。ここで、冠動脈の内腔３１にプラーク３２が形成されていると、内腔３１の径は狭くなる。このため、切断面決定部２４は、径ｄ１〜ｄ４のうちの最小となる径となる内腔３１の位置をプラーク３２が存在する関心領域に設定する。図４においては、径ｄ１となる内腔３１の位置を関心領域に設定する。
また、冠動脈は心臓の表面に存在するため、断面画像ＰＧｋには、図５に示すように心筋３３が含まれる。心筋３３は、冠動脈に血流が十分に供給されるか否かにより、活動状況が変化するため、心筋梗塞等の診断上、重要な領域である。切断面決定部２４は、各断面Ｐｋの断面画像において心筋３３が存在する位置を関心領域に決定する。
また、切断面決定部２４は、断面Ｐｋが分岐に存在する場合、図６に示すように、分岐３４を関心領域に設定する。
そして、切断面決定部２４は、関心領域を含むように各断面Ｐｋにおける切断面を決定する。例えば、図４に示すようにプラーク３２が関心領域である場合、径ｄ１を含むように切断面Ｃｋを決定する。また、図５に示すように心筋３３が関心領域である場合、基準線Ｌ０を通り、かつ心筋３３を含むように切断面Ｃｋを決定する。また、図６に示すように、分岐３４が関心領域の場合、基準線Ｌ０としなかった芯線Ｃ０を通るように切断面Ｃｋを決定する。
なお、断面Ｐｋの位置によっては、プラークと心筋、またはプラークと分岐というように、複数の関心領域が存在する場合がある。特に本実施形態においては、冠動脈を診断の対象の構造物としているため、各断面Ｐｋにおいて必ず心筋が含まれる。本実施形態においては、切断面決定部２４は、関心領域について予め優先度を設定しておき、１つの断面に複数の関心領域が存在する場合は、優先度が高い方の関心領域を選択し、選択された関心領域を含むように切断面を決定する。本実施形態においては、プラーク、分岐および心筋の順で優先度を設定するものとする。このため、冠動脈領域３０におけるプラーク３２が存在する領域においては、心筋３３が存在していてもプラーク３２が関心領域として選択される。また、冠動脈領域３０における分岐３４が存在する領域においては、心筋３３が存在していても、分岐３４が関心領域として選択される。なお、各断面ｐｋの断面画像ＰＧｋをディスプレイ１４に表示し、操作者が切断面に含める関心領域を選択してもよい。
画像生成部２５は、切断面決定部２４が決定した各断面Ｐｋの切断面Ｃｋを、補間する等して滑らかに接続して切断曲面を生成する。図７〜図８は切断曲面の生成を説明するための図である。本実施形態においては、冠動脈領域３０の背後に心筋３３が存在し、図７に示すように、冠動脈領域３０Ａの破線で囲んだ部分Ｂ１にプラーク３２が存在しているものとする。この場合、図８に示すように、各断面Ｐｋにおける複数の関心領域が切断されるように、切断面決定部２４により切断曲面ＣＭ０が決定される。なお、図８においては、矢印の方向が各断面Ｐｋにおける切断面Ｃｋの方向である。これにより、部分Ｂ１においてはプラーク３２を切断し、分岐３４においては冠動脈領域３０Ｂを切断し、その他の領域においては心筋３３を切断するように切断曲面ＣＭ０が決定される。
そして、画像生成部２５は、冠動脈領域３０を切断曲面ＣＭ０により切断してＣＰＲ画像を生成する。なお、本実施形態において生成されるＣＰＲ画像は、ストレートＣＰＲ(Straightened CPR)画像とするが、ストレッチＣＰＲ(Stretched CPR)画像や、プロジェクトＣＰＲ(Projected CPR)画像でもよい。図９はＣＰＲ画像を示す図である。図９に示すように、本実施形態において生成されるＣＰＲ画像Ｇ１０においては、冠動脈領域３０において、プラーク３２および分岐３４が含まれるものとなる。また、分岐３４を含むため、基準線Ｌ０としなかった芯線Ｃ０についての冠動脈領域３０Ｂが含まれるものとなる。なお、ＣＰＲ画像Ｇ１０において、プラーク３２および冠動脈領域３０Ｂ以外の部分は、心筋３３である。
次いで、本実施形態において行われる処理について説明する。図１０は本実施形態にお手行われる処理を示すフローチャートである。まず、画像取得部２１が、３次元画像Ｇ０を取得し、構造物抽出部２２が、３次元画像Ｇ０から冠動脈領域３０を抽出する（ステップＳＴ１）。次いで、断面設定部２３が、冠動脈領域３０における基準線Ｌ０上の各点について、基準線Ｌ０に垂直な断面Ｐｋを設定する（ステップＳＴ２）。そして、切断面決定部２４が、各断面Ｐｋ上に関心領域が存在する場合、基準線Ｌ０および関心領域を含む、各断面における１つの切断面Ｃｋを決定する（ステップＳＴ３）。さらに、画像生成部２５が、切断面Ｃｋを含むＣＰＲ画像Ｇ１０を生成し（ステップＳＴ４）、生成されたＣＰＲ画像Ｇ１０をディスプレイ１４に表示し（ステップＳＴ５）、処理を終了する。
このように、本実施形態によれば、各断面において基準線を中心とする放射状の任意の方向に存在する関心領域を含むようにＣＰＲ画像Ｇ１０を生成することができる。このため、ＣＰＲ画像Ｇ１０において冠動脈領域３０に含まれる複数の関心領域を良好に観察することができる。また、各断面Ｐｋにおいて１つの切断面Ｃｋを決定しているため、例えば関心領域が分岐である場合において、ＣＰＲ画像に複数の分岐を含ませることができ、かつ分岐位置においてＣＰＲ画像Ｇ１０が連続するようになる。このため、境界がない自然な印象のＣＰＲ画像Ｇ１０を生成することができる。
また、各断面Ｐｋにおける切断面Ｃｋが連続する切断曲面ＣＭ０を生成し、切断曲面ＣＭ０を含むＣＰＲ画像Ｇ１０を生成することにより、切断面Ｃｋが滑らかに変化する。このため、断面間において、ＣＰＲ画像Ｇ１０に境界が発生してしまうことを防止できる。
また、各断面に複数の関心領域が存在する場合、１つの関心領域を選択することにより、選択された関心領域をＣＰＲ画像Ｇ１０に含めることができる。
なお、上記実施形態において、ＣＰＲ画像Ｇ１０を表示する際に、各断面Ｐｋにおける切断面Ｃｋの基準となる断面の切断面からの基準線Ｌ０を中心とした角度を表す情報を表示してもよい。例えば、図１１の上側に示す基準となる断面ＰＢにおける切断面ＣＢｋを基準となる切断面とする。そして、図１１の下側に示すある断面Ｐｋの切断面Ｃｋについて、基準となる切断面ＣＢｋからの基準線Ｌ０を中心とした角度θｋを求める。そして、求めた角度θｋを表す角度情報をＣＰＲ画像Ｇ１０におけるその切断面Ｃｋを得た断面Ｐｋと対応づけて表示する。図１２は角度情報４０を表示したＣＰＲ画像を示す図である。なお、角度情報４０は適宜サンプリングした断面Ｐｋと対応づけて表示される。これにより、ＣＰＲ画像Ｇ１０において注目している領域が、基準となる切断面からどの程度角度がずれているかを認識することができる。また、角度θｋの大きさに応じて、ＣＰＲ画像Ｇ１０における冠動脈領域３０を色分けして表示してもよい。
なお、上記実施形態においては、冠動脈を対象構造物としているが、他の管状構造物、例えば気管支および大腸等を対象構造物としてもよい。対象構造物が気管支の場合、構造物抽出部２２は、３次元画像Ｇ０から気管支の構造を気管支領域として抽出する。具体的には、例えば特開２０１０−２２０７４２号公報等に記載された手法を用いて、入力された３次元画像Ｇ０に含まれる気管支領域のグラフ構造を、３次元の気管支領域として抽出する。また、対象構造物が大腸の場合、構造物抽出部２２は、３次元画像Ｇ０における大腸の画素値となる領域を大腸領域として抽出する。
また、気管支を対象構造物とする場合、診断に重要なのは気管支の近傍に存在する腫瘍である。また、大腸を対象構造物とする場合、診断に重要なのは大腸内腔壁に存在する腫瘍である。このため、構造物抽出部２２において、ＣＡＤ（Computer-Aided Diagnosis）により腫瘍を抽出することが好ましい。この場合、断面設定部２３は腫瘍を関心領域として切断面Ｃｋを設定すればよい。
とくに、気管支を対象構造物とした場合、本実施形態によれば、気管支の入り口から複数の分岐および腫瘍を含むようにＣＰＲ画像Ｇ１０を生成できる。このため、腫瘍の生検を行うための内視鏡検査を行う際に、内視鏡を気管支のどの分岐を通って腫瘍に到達させればよいかを、ＣＰＲ画像Ｇ１０により、容易に確認することができる。
また、上記実施形態においては、冠動脈等の管状構造物を対象としてＣＰＲ画像を生成しているが、これに限定されるものではなく、脊柱、四肢等の長手方向に延びる構造物を対象構造物としてＣＰＲ画像を生成してもよい。
また、上記実施形態においては、ＣＰＲ画像Ｇ１０と同時に、心筋の機能画像を表示してもよい。これにより、冠動脈の状態と心臓の機能の状態とを同時に確認することができる。
以下、本実施形態の作用効果について説明する。
各断面における切断面が連続する切断曲面を生成し、切断曲面を含むＣＰＲ画像を生成することにより、切断面が滑らかに変化するため、各断面の間において、ＣＰＲ画像に境界が発生してしまうことをより確実に防止できる。
各断面に複数の関心領域が存在する場合、１つの関心領域を選択することにより、選択された関心領域をＣＰＲ画像に含めることができる。
各断面における切断面の、基準となる断面の切断面からの基準線を中心とした角度を表す情報を表示することにより、ＣＰＲ画像において注目している領域が、基準となる断面の切断面からどの程度角度がずれているかを認識することができる。
１ ＣＰＲ画像生成装置
２ ３次元画像撮影装置
３ 画像保管サーバ
４ ネットワーク
１１ ＣＰＵ
１２ メモリ
１３ ストレージ
１４ ディスプレイ
１５ 入力部
２１ 画像取得部
２２ 構造物抽出部
２３ 断面設定部
２４ 切断面決定部
２５ 画像生成部
３０，３０Ａ，３０Ｂ 冠動脈領域
３１ 内腔
３２ プラーク
３３ 心筋
３４ 分岐
４０ 角度情報
Ｃｋ 切断面
ＣＭ０ 切断曲面
Ｇ１０ ＣＰＲ画像
Ｐｋ 断面
ＰＧｋ 断面画像
２ ３次元画像撮影装置
３ 画像保管サーバ
４ ネットワーク
１１ ＣＰＵ
１２ メモリ
１３ ストレージ
１４ ディスプレイ
１５ 入力部
２１ 画像取得部
２２ 構造物抽出部
２３ 断面設定部
２４ 切断面決定部
２５ 画像生成部
３０，３０Ａ，３０Ｂ 冠動脈領域
３１ 内腔
３２ プラーク
３３ 心筋
３４ 分岐
４０ 角度情報
Ｃｋ 切断面
ＣＭ０ 切断曲面
Ｇ１０ ＣＰＲ画像
Ｐｋ 断面
ＰＧｋ 断面画像
Claims (7)
- 撮影により取得された３次元画像から、対象構造物を抽出する構造物抽出手段と、
前記対象構造物の基準線上の各点において、該基準線に垂直な断面を設定する断面設定手段と、
前記各断面上に関心領域が存在する場合、前記基準線および前記関心領域を含む、前記各断面における１つの切断面を決定する切断面決定手段と、
前記各断面における切断面を含むＣＰＲ画像を生成する画像生成手段と、
前記ＣＰＲ画像を表示する表示手段と、
を備え
前記表示手段は、前記各断面における切断面の、基準となる断面の切断面からの前記基準線を中心とした角度を表す情報を含む前記ＣＰＲ画像を表示することを特徴とするＣＰＲ画像生成装置。 - 前記対象構造物は、管状構造物である請求項１記載のＣＰＲ画像生成装置。
- 前記関心領域は、前記対象構造物の分岐、該対象構造物に隣接する他の構造物、該対象構造物内の他の構造物、および病変の少なくとも１つである請求項１または２記載のＣＰＲ画像生成装置。
- 前記画像生成手段は、前記各断面における切断面が連続する切断曲面を生成し、該切断曲面を含む前記ＣＰＲ画像を生成する請求項１から３のいずれか１項記載のＣＰＲ画像生成装置。
- 前記各断面に複数の関心領域が存在する場合、前記切断面決定手段は、１つの関心領域を選択する請求項１から４のいずれか１項記載のＣＰＲ画像生成装置。
- 撮影により取得された３次元画像から、対象構造物を抽出し、
前記対象構造物の基準線上の各点において、該基準線に垂直な断面を設定し、
前記各断面上に関心領域が存在する場合、前記基準線および前記関心領域を含む、前記各断面における１つの切断面を決定し、
前記各断面における切断面を含むＣＰＲ画像を生成し、
前記ＣＰＲ画像、および前記各断面における切断面の、基準となる断面の切断面からの前記基準線を中心とした角度を表す情報を表示することを特徴とするＣＰＲ画像生成方法。 - 撮影により取得された３次元画像から、対象構造物を抽出する手順と、
前記対象構造物の基準線上の各点において、該基準線に垂直な断面を設定する手順と、
前記各断面上に関心領域が存在する場合、前記基準線および前記関心領域を含む、前記各断面における１つの切断面を決定する手順と、
前記各断面における切断面を含むＣＰＲ画像を生成する手順と、
前記ＣＰＲ画像、および前記各断面における切断面の、基準となる断面の切断面からの前記基準線を中心とした角度を表す情報を表示する手順をコンピュータに実行させることを特徴とするＣＰＲ画像生成プログラム。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016169176 | 2016-08-31 | ||
| JP2016169176 | 2016-08-31 | ||
| PCT/JP2017/024637 WO2018042870A1 (ja) | 2016-08-31 | 2017-07-05 | Ｃｐｒ画像生成装置、方法およびプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2018042870A1 JPWO2018042870A1 (ja) | 2019-01-24 |
| JP6671482B2 true JP6671482B2 (ja) | 2020-03-25 |
Family
ID=61301469
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018536988A Active JP6671482B2 (ja) | 2016-08-31 | 2017-07-05 | Ｃｐｒ画像生成装置、方法およびプログラム |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20190197762A1 (ja) |
| JP (1) | JP6671482B2 (ja) |
| WO (1) | WO2018042870A1 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109920027B (zh) * | 2019-03-01 | 2021-02-09 | 数坤（北京）网络科技有限公司 | 一种曲面重建方法、设备及计算机可读介质 |
| WO2020217807A1 (ja) * | 2019-04-25 | 2020-10-29 | 富士フイルム株式会社 | 画像の向き設定装置、方法およびプログラム |
| CN111754506B (zh) * | 2020-07-01 | 2024-02-06 | 杭州脉流科技有限公司 | 基于腔内影像的冠脉狭窄率计算方法、装置、系统和计算机存储介质 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4515615B2 (ja) * | 2000-09-14 | 2010-08-04 | 株式会社日立メディコ | 画像表示装置 |
| GB2395880B (en) * | 2002-11-27 | 2005-02-02 | Voxar Ltd | Curved multi-planar reformatting of three-dimensional volume data sets |
| JP4717935B2 (ja) * | 2009-03-23 | 2011-07-06 | 富士フイルム株式会社 | 画像処理装置および方法並びにプログラム |
| JP5572470B2 (ja) * | 2010-07-28 | 2014-08-13 | 富士フイルム株式会社 | 診断支援装置および方法並びにプログラム |
| US9189847B2 (en) * | 2013-11-12 | 2015-11-17 | Siemens Aktiengesellschaft | 2D visualization for rib analysis |
-
2017
- 2017-07-05 JP JP2018536988A patent/JP6671482B2/ja active Active
- 2017-07-05 WO PCT/JP2017/024637 patent/WO2018042870A1/ja active Application Filing
-
2019
- 2019-02-26 US US16/286,424 patent/US20190197762A1/en not_active Abandoned
Also Published As
| Publication number | Publication date |
|---|---|
| WO2018042870A1 (ja) | 2018-03-08 |
| US20190197762A1 (en) | 2019-06-27 |
| JPWO2018042870A1 (ja) | 2019-01-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4662766B2 (ja) | 最適化視野マップを生成する方法および画像化システム並びにコンピューターワークステーション及びコンピューター読取可能媒体 | |
| JP5312801B2 (ja) | 医用画像ビューイングプロトコル | |
| JP6080248B2 (ja) | ３次元画像表示装置および方法並びにプログラム | |
| US9466117B2 (en) | Segmentation highlighter | |
| JP2006075599A (ja) | 管状の解剖学的構造の可視化方法 | |
| CN104050711B (zh) | 医用图像处理装置以及医用图像处理方法 | |
| JP2023139022A (ja) | 医用画像処理方法、医用画像処理装置、医用画像処理システム及び医用画像処理プログラム | |
| JP5934071B2 (ja) | 管状構造の最短経路探索装置および方法並びにプログラム | |
| JP2007135858A (ja) | 画像処理装置 | |
| JP6671482B2 (ja) | Ｃｐｒ画像生成装置、方法およびプログラム | |
| JP3989896B2 (ja) | 医用画像処理装置、関心領域抽出方法、ならびに、プログラム | |
| CN102548480B (zh) | 医用图像显示设备和方法 | |
| JP6066197B2 (ja) | 手術支援装置、方法およびプログラム | |
| US20200188027A1 (en) | Virtual stent placement apparatus, virtual stent placement method, and virtual stent placement program | |
| JP6734111B2 (ja) | 所見情報作成装置及びシステム | |
| US20160157726A1 (en) | Projection image generating device, method and program | |
| JP5707457B2 (ja) | 心筋運動の可視化 | |
| EP3828836B1 (en) | Method and data processing system for providing a two-dimensional unfolded image of at least one tubular structure | |
| JP5992853B2 (ja) | 手術支援装置、方法およびプログラム | |
| JP7277652B2 (ja) | 修正指示領域表示装置、方法およびプログラム | |
| JPWO2019064791A1 (ja) | 仮想ステント設置装置、方法およびプログラム | |
| JP6494942B2 (ja) | 画像処理装置、及び、画像処理方法 | |
| WO2022270150A1 (ja) | 画像処理装置、方法およびプログラム | |
| JP6662580B2 (ja) | 医用画像処理装置 | |
| JPWO2020090445A1 (ja) | 領域修正装置、方法およびプログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180919 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190723 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190919 |
|
| 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: 20200212 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200303 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6671482 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |