JP2018158022A - Medical image diagnostic device and medical image processing device - Google Patents

Medical image diagnostic device and medical image processing device Download PDF

Info

Publication number
JP2018158022A
JP2018158022A JP2017057926A JP2017057926A JP2018158022A JP 2018158022 A JP2018158022 A JP 2018158022A JP 2017057926 A JP2017057926 A JP 2017057926A JP 2017057926 A JP2017057926 A JP 2017057926A JP 2018158022 A JP2018158022 A JP 2018158022A
Authority
JP
Japan
Prior art keywords
image
measurement value
display
medical 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.)
Granted
Application number
JP2017057926A
Other languages
Japanese (ja)
Other versions
JP6878074B2 (en
Inventor
吉田 真樹
Maki Yoshida
真樹 吉田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Medical Systems Corp
Original Assignee
Canon Medical Systems Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Canon Medical Systems Corp filed Critical Canon Medical Systems Corp
Priority to JP2017057926A priority Critical patent/JP6878074B2/en
Publication of JP2018158022A publication Critical patent/JP2018158022A/en
Application granted granted Critical
Publication of JP6878074B2 publication Critical patent/JP6878074B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

PROBLEM TO BE SOLVED: To provide a medical image diagnostic device and a medical image processing device capable of generating an image for easy diagnosis.SOLUTION: A medical image diagnostic device includes: an image acquisition part for acquiring a medical image of a plurality of time phases; an extraction part for extracting a predetermined region of interest in the medical image of the plurality of time phases; a calculation part for calculating a predetermined measurement value on the predetermined region of interest; and a generation part for generating an image in which the measurement values are arranged by spatial position and by time series along the predetermined region of interest.SELECTED DRAWING: Figure 1

Description

本発明の実施形態は、医用画像診断装置及び医用画像処理装置に関する。   Embodiments described herein relate generally to a medical image diagnostic apparatus and a medical image processing apparatus.

複数時相の医用画像の観察に関する技術として、TDC(Time Dentsiy Curve)解析がある。TDC解析は、複数時相の医用画像の画素値を解析する手法の一種である。複数時相の医用画像の任意の座標にROI(関心領域)を設定し、ROI上の時間的な画素値変化をグラフとして描画する。このTDC解析により、例えば、時間経過に伴う造影剤の浸透の様子を観察できるようになり、血管の血流解析が可能となる。   As a technique related to observation of medical images of a plurality of time phases, there is TDC (Time Dentsiy Curve) analysis. TDC analysis is a type of technique for analyzing pixel values of a plurality of time-phase medical images. An ROI (region of interest) is set at an arbitrary coordinate of a medical image of a plurality of time phases, and temporal pixel value changes on the ROI are drawn as a graph. By this TDC analysis, for example, it becomes possible to observe the state of infiltration of the contrast agent over time, and blood flow analysis of blood vessels becomes possible.

しかし、TDC解析では、設定するROIの位置や数を任意に設定可能ではあるが、表示されるグラフには各ROI上の画素値変化のみが描画されることになる。このため、ROI間の空間的な位置関係に意味がある場合でも、各ROI間の空間的な情報を確認できない。そのため、空間的な位置関係を確認するためには、元の医用画像を参照する必要がある。   However, in the TDC analysis, although the position and number of ROIs to be set can be arbitrarily set, only the pixel value change on each ROI is drawn on the displayed graph. For this reason, even when the spatial positional relationship between the ROIs is meaningful, the spatial information between the ROIs cannot be confirmed. Therefore, in order to confirm the spatial positional relationship, it is necessary to refer to the original medical image.

また、複数時相の医用画像の比較のために、過去画像で指定した関心領域を、別の時相画像の任意の位置、例えば該当関心領域の隣などに、重ねて表示する技術もある(特許文献1)。これにより少ないデータサイズ(例えば、表示画像1枚分)で時相画像が比較でき、また、比較する関心領域の表示位置が近いので、視線の移動距離が少なく済むというメリットもある。   In addition, there is a technique for displaying a region of interest specified in a past image in an overlapping manner at an arbitrary position of another time phase image, for example, next to the region of interest for comparison of medical images of a plurality of time phases ( Patent Document 1). As a result, the time phase images can be compared with a small data size (for example, one display image), and the display position of the region of interest to be compared is close, so that the moving distance of the line of sight can be reduced.

しかし、2次元画像上に別の2次元画像の関心領域画像を表示するため、2次元画像としてのx、y軸方向の情報と時間軸方向の情報が、2次元画像上で混在して存在してしまう。このため、比較する時相が多い場合、情報の混在は、対象画像の観察が困難になる要因となる。   However, since the region-of-interest image of another two-dimensional image is displayed on the two-dimensional image, information in the x- and y-axis directions and information in the time-axis direction as a two-dimensional image are mixed on the two-dimensional image. Resulting in. For this reason, when there are many time phases to be compared, mixing of information becomes a factor that makes it difficult to observe the target image.

このため、複数時相の医用画像について、読影者に容易に視認可能な態様で表示する医用画像診断装置及び医用画像処理装置が望まれている。   For this reason, a medical image diagnostic apparatus and a medical image processing apparatus that display medical images of a plurality of time phases in a manner that can be easily visually recognized by a radiogram interpreter are desired.

特開2014−23921号公報JP 2014-23921 A 特開2007−151881号公報JP 2007-151881 A 特開2004−517667号公報JP 2004-517667 A 特開2015−217170号公報JP 2015-217170 A

本実施形態の目的は、読影者に容易に視認可能な画像を生成する医用画像診断装置及び医用画像処理装置を提供することにある。   An object of the present embodiment is to provide a medical image diagnostic apparatus and a medical image processing apparatus that generate an image that can be easily visually recognized by a radiogram interpreter.

本実施形態に係る医用画像診断装置は、複数時相の医用画像を取得する、画像取得部と、前記複数時相の医用画像における所定の関心領域の抽出を行う、抽出部と、前記所定の関心領域に関する所定の測定値を算出する、算出部と、前記所定の関心領域に沿って前記測定値を空間的位置及び時系列毎に並べた画像を生成する、生成部と、を備えている。
本実施形態に係る医用画像処理装置は、複数時相の医用画像を取得する、画像取得部と、前記複数時相の医用画像における所定の関心領域の抽出を行う、抽出部と、前記所定の関心領域に関する所定の測定値を算出する、算出部と、前記所定の関心領域に沿って前記測定値を空間的位置及び時系列毎に並べた画像を生成する、生成部と、を備えている。
The medical image diagnostic apparatus according to the present embodiment includes an image acquisition unit that acquires medical images of a plurality of time phases, an extraction unit that extracts a predetermined region of interest in the medical images of the plurality of time phases, and the predetermined A calculation unit that calculates a predetermined measurement value related to the region of interest, and a generation unit that generates an image in which the measurement values are arranged in spatial positions and time series along the predetermined region of interest. .
The medical image processing apparatus according to the present embodiment includes an image acquisition unit that acquires medical images of a plurality of time phases, an extraction unit that extracts a predetermined region of interest in the medical images of the plurality of time phases, and the predetermined A calculation unit that calculates a predetermined measurement value related to the region of interest, and a generation unit that generates an image in which the measurement values are arranged in spatial positions and time series along the predetermined region of interest. .

第1実施形態に係るX線CT装置の全体構成を説明するブロック図。The block diagram explaining the whole structure of the X-ray CT apparatus which concerns on 1st Embodiment. 第1実施形態に係るX線CT装置で実行される画像表示処理の内容を説明するフローチャートを示す図。The figure which shows the flowchart explaining the content of the image display process performed with the X-ray CT apparatus which concerns on 1st Embodiment. 時相毎に血管の芯線を抽出して2次元画像化する処理を概念的に説明する図。The figure which illustrates notionally the process which extracts the core line of a blood vessel for every time phase, and makes it a two-dimensional image. 芯線に沿ったMAP画像を生成する処理を概念的に説明する図。The figure which illustrates notionally the process which produces | generates the MAP image along a core line. MAP画像に基づいて生成された表示画像と、芯線に沿った2次元画像に基づいて生成された参考画像とを、画面に表示した状態の一例を示す図。The figure which shows an example of the state which displayed on the screen the display image produced | generated based on the MAP image, and the reference image produced | generated based on the two-dimensional image along a core line. MAP画像に基づいて生成された表示画像に、曲線ROIを重ねて表示した画面の一例を示す図。The figure which shows an example of the screen which overlapped and displayed the curve ROI on the display image produced | generated based on the MAP image. 図6の画像を生成して表示するX線CT装置の全体構成を説明するブロック図。The block diagram explaining the whole structure of the X-ray CT apparatus which produces | generates and displays the image of FIG. (a)記憶回路に記憶されている、単色のグラデーションでMAP画像に基づく表示画像を生成するためのテーブルにおける、CT値と輝度値との対応関係を説明する図、(b)記憶回路に記憶されている、カラーでMAP画像に基づく表示画像を生成するためのテーブルにおける、CT値と色との対応関係を説明する図。(A) A diagram for explaining the correspondence between CT values and luminance values in a table for generating a display image based on a MAP image with a single color gradation stored in the storage circuit, and (b) storing in the storage circuit. The figure explaining the correspondence of CT value and color in the table for producing | generating the display image based on the MAP image by the color. 第2実施形態に係るX線CT装置で実行される画像表示処理の内容を説明するフローチャートを示す図。The figure which shows the flowchart explaining the content of the image display process performed with the X-ray CT apparatus which concerns on 2nd Embodiment. 肺の画像にROIとして軸を設定した状態と、測定値として用いる肺の断面積の計測値を説明する図。The figure explaining the state which set the axis | shaft as ROI to the image of a lung, and the measured value of the cross-sectional area of the lung used as a measured value. 第3実施形態に係るX線CT装置で実行される画像表示処理の内容を説明するフローチャートを示す図。The figure which shows the flowchart explaining the content of the image display process performed with the X-ray CT apparatus which concerns on 3rd Embodiment. 心筋壁を曲面として捉えた場合における、曲面の法線方向の平均ベクトル方向を2次元平面に射影することで2次元画像に変換する方法を説明する図。The figure explaining the method of converting into the two-dimensional image by projecting the average vector direction of the normal direction of a curved surface to a two-dimensional plane in the case of capturing the myocardial wall as a curved surface. 心筋壁を曲面として捉えた場合における、曲面に格子状の座標系を配置し、設定した2次元平面の座標系へと曲面の座標系を射影することで、2次元画像に変換する方法を説明した図。Explains how to convert a 2D image by arranging a grid coordinate system on the curved surface and projecting the curved coordinate system to the set 2D plane coordinate system when the myocardial wall is captured as a curved surface Figure. 心筋壁を2次元画像化して2次元のMAP画像を生成し、この2次元のMAP画像を時相毎に並べて3次元のMAP画像を生成する処理を概念的に説明する図。The figure which illustrates notionally the process which makes a myocardial wall into a two-dimensional image, produces | generates a two-dimensional MAP image, arranges this two-dimensional MAP image for every time phase, and produces | generates a three-dimensional MAP image. モダリティに医用画像処理装置を接続して、第1実施形態乃至第3実施形態に係る画像表示処理を医用画像処理装置で実行する変形例を説明する図。The figure explaining the modification which connects a medical image processing apparatus to modality, and performs the image display process which concerns on 1st Embodiment thru | or 3rd Embodiment with a medical image processing apparatus.

以下、図面を参照しながら、本実施形態に係る医用画像診断装置及び医用画像処理装置を説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行うこととする。   Hereinafter, a medical image diagnostic apparatus and a medical image processing apparatus according to the present embodiment will be described with reference to the drawings. In the following description, components having substantially the same functions and configurations are denoted by the same reference numerals, and redundant description will be provided only when necessary.

〔第1実施形態〕
図1は、第1実施形態に係るX線CT装置1の全体構成を説明する図である。この図1に示すX線CT装置1は、本実施形態における医用画像診断装置の一例であり、例えば、架台装置10と、寝台装置30と、コンソール装置40とを備えて構成されている。
[First Embodiment]
FIG. 1 is a diagram illustrating the overall configuration of an X-ray CT apparatus 1 according to the first embodiment. The X-ray CT apparatus 1 shown in FIG. 1 is an example of a medical image diagnostic apparatus according to this embodiment, and includes, for example, a gantry device 10, a bed device 30, and a console device 40.

より具体的には、架台装置10は、X線発生装置11と、X線検出器12と、回転体13と、X線高電圧装置14と、架台制御装置15と、データ収集回路18とを、備えて構成されている。   More specifically, the gantry device 10 includes an X-ray generator 11, an X-ray detector 12, a rotating body 13, an X-ray high voltage device 14, a gantry controller 15, and a data collection circuit 18. , And is configured.

X線発生装置11は、例えば、X線高電圧装置14から高電圧の電力供給を受けて、陰極(フィラメントと呼ぶ場合もある)から陽極(ターゲットと呼ぶ場合もある)に向けて熱電子を照射するX線管(真空管)から構成される。   The X-ray generator 11 receives, for example, a high-voltage power supply from the X-ray high-voltage device 14 and directs thermoelectrons from the cathode (sometimes called a filament) toward the anode (sometimes called a target). It consists of an X-ray tube (vacuum tube) for irradiation.

X線発生装置11は、X線管に限定される必要は無い。例えば、X線発生装置11は、X線管に代えて、電子銃から発生した電子ビームを集束させるフォーカスコイルと電磁偏向させる偏向コイルと、被検体Pの半周を囲い偏向した電子ビーム衝突することによってX線を発生させるターゲットリングとによって構成しても構わない。   The X-ray generator 11 need not be limited to an X-ray tube. For example, instead of the X-ray tube, the X-ray generator 11 collides with a focus coil that focuses an electron beam generated from an electron gun, a deflection coil that electromagnetically deflects, and an electron beam that deflects around a half circumference of the subject P. And a target ring that generates X-rays.

X線検出器12は、例えば、X線管の焦点を中心として1つの円弧に沿ってチャネル方向に複数のX線検出素子が配列された複数のX線検出素子列から構成される。X線検出器12は、チャネル方向に複数のX線検出素子が配列されたX線検出素子列がスライス方向に複数配列された構造を有する。X線検出器12は、X線発生装置11から照射され、被検体Pを通過したX線を検出し、当該X線量に対応した電気信号をデータ収集回路18へと出力する。   The X-ray detector 12 includes, for example, a plurality of X-ray detection element arrays in which a plurality of X-ray detection elements are arranged in the channel direction along one circular arc with the focal point of the X-ray tube as the center. The X-ray detector 12 has a structure in which a plurality of X-ray detection element arrays in which a plurality of X-ray detection elements are arranged in the channel direction are arranged in the slice direction. The X-ray detector 12 detects X-rays irradiated from the X-ray generator 11 and passed through the subject P, and outputs an electrical signal corresponding to the X-ray dose to the data collection circuit 18.

また、X線検出器12は、例えば、グリッドと、シンチレータアレイと、光センサアレイとから構成される間接変換型の検出器である。   Moreover, the X-ray detector 12 is an indirect conversion type detector comprised of, for example, a grid, a scintillator array, and an optical sensor array.

シンチレータアレイは、複数のシンチレータから構成され、シンチレータは入射X線量に応じた光子量の光を出力するシンチレータ結晶にて構成される。   The scintillator array is composed of a plurality of scintillators, and the scintillator is composed of a scintillator crystal that outputs a photon amount of light corresponding to the incident X-ray dose.

グリッドは、シンチレータアレイのX線入射側の面に配置され、散乱X線を吸収する機能を有するX線遮蔽板で構成される。   The grid is an X-ray shielding plate that is disposed on the surface on the X-ray incident side of the scintillator array and has a function of absorbing scattered X-rays.

光センサアレイは、シンチレータからの光量に応じた電気信号に変換する機能を有し、例えば、光電子増倍管等の光センサから構成される。   The photosensor array has a function of converting into an electrical signal corresponding to the amount of light from the scintillator, and is composed of a photosensor such as a photomultiplier tube, for example.

なお、X線検出器12は、間接変換型の検出器に限定されず、入射したX線を直接電気信号に変換する半導体素子から構成される直接変換型の検出器であっても構わない。   The X-ray detector 12 is not limited to an indirect conversion type detector, and may be a direct conversion type detector composed of a semiconductor element that converts incident X-rays directly into an electrical signal.

回転体13は、この回転体13の中心を回転軸として回転自在に支持されており、架台制御装置15の制御に基づいて、回転駆動される。   The rotating body 13 is rotatably supported with the center of the rotating body 13 as a rotation axis, and is driven to rotate based on the control of the gantry control device 15.

X線高電圧装置14は、変圧器(トランス)及び整流器等の電気回路から構成され、X線発生装置11に印加する高電圧を発生する機能を有する高電圧発生装置と、X線発生装置11が照射するX線に応じた出力電圧の制御を行うX線制御装置から構成される。高電圧発生装置は、変圧器方式であってもよいし、インバータ方式であっても構わない。   The X-ray high voltage device 14 includes an electric circuit such as a transformer and a rectifier, and has a function of generating a high voltage to be applied to the X-ray generator 11, and the X-ray generator 11. Is constituted by an X-ray control device that controls the output voltage according to the X-rays irradiated. The high voltage generator may be a transformer system or an inverter system.

架台制御装置15は、CPU等によって構成される処理回路とモータ及びアクチュエータ等の駆動機構から構成される。架台制御装置15は、コンソール装置40若しくは架台装置10に取り付けられた入力装置からの入力信号を受けて、架台の動作制御を行う機能を有する。例えば、架台制御装置15は、入力信号を受けて回転体13を回転させる制御や、架台装置10をチルトさせる制御、及び寝台装置30及び天板33を動作させる制御を行う。   The gantry control device 15 includes a processing circuit including a CPU and a driving mechanism such as a motor and an actuator. The gantry control device 15 has a function of receiving an input signal from the console device 40 or an input device attached to the gantry device 10 and controlling the operation of the gantry. For example, the gantry control device 15 receives the input signal and performs control to rotate the rotating body 13, control to tilt the gantry device 10, and control to operate the bed device 30 and the top plate 33.

データ収集回路(DAS:Data Acquisition System)18は、X線検出器12の各X線検出素子から出力される電気信号に対して増幅処理を行う増幅器と、電気信号をデジタル信号に変換するA/D変換器とから少なくとも構成され、検出データ(純生データ)を生成する。データ収集回路18が生成した検出データは、コンソール装置40へと転送される。   A data acquisition circuit (DAS) 18 includes an amplifier that performs amplification processing on an electric signal output from each X-ray detection element of the X-ray detector 12 and an A / A that converts the electric signal into a digital signal. It comprises at least a D converter and generates detection data (pure raw data). The detection data generated by the data collection circuit 18 is transferred to the console device 40.

寝台装置30は、スキャン対象の被検体Pを載置、移動させる装置であり、基台31と、寝台駆動装置32と、天板33と、支持フレーム34とを備えている。   The couch device 30 is a device for placing and moving the subject P to be scanned, and includes a base 31, a couch driving device 32, a top plate 33, and a support frame 34.

基台31は、支持フレーム34を鉛直方向に移動可能に支持する筐体である。寝台駆動装置32は、被検体Pが載置された支持フレームを支持フレーム34の長軸方向に移動するモータあるいはアクチュエータである。支持フレーム34の上面に設けられた天板33は、被検体Pが載置される板である。   The base 31 is a housing that supports the support frame 34 so as to be movable in the vertical direction. The bed driving device 32 is a motor or actuator that moves the support frame on which the subject P is placed in the longitudinal direction of the support frame 34. The top plate 33 provided on the upper surface of the support frame 34 is a plate on which the subject P is placed.

天板33は、天板33だけを移動させてもよいし、寝台装置30の支持フレーム34ごと移動する方式であってもよい。本実施形態を立位CTに応用する場合は、天板33に相当する患者移動機構を移動する方式であってもよい。   The top plate 33 may move only the top plate 33 or may move along with the support frame 34 of the bed apparatus 30. When this embodiment is applied to standing CT, a method of moving a patient moving mechanism corresponding to the top board 33 may be used.

架台装置10の撮像系と天板33の位置関係の相対的な変更を伴うスキャン(ヘリカルスキャンや位置決めスキャン等)を実行する際には、当該位置関係の相対的な変更は天板33の駆動によって行われてもよいし、架台装置10の走行によって行われてもよく、またそれらの複合によって行われてもよい。   When performing a scan (helical scan, positioning scan, etc.) involving a relative change in the positional relationship between the imaging system of the gantry 10 and the top plate 33, the relative change in the positional relationship is driven by the top plate 33. May be performed by running the gantry device 10, or may be performed by combining them.

コンソール装置40は、記憶回路41と、表示装置42と、入力装置43と、処理回路44とを、備えて構成されている。   The console device 40 includes a storage circuit 41, a display device 42, an input device 43, and a processing circuit 44.

記憶回路41は、例えば、RAM(Random Access Memory)、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等により実現される。記憶回路41は、例えば、投影データや再構成画像データを記憶する。本実施形態においては、記憶回路41は、例えば、画像記憶回路411、画素値対応付け記憶回路412、及び、プログラム記憶回路413を有する。詳しくは後述するが、本実施形態においては、画像記憶回路411は、架台装置10で撮像された医用画像の画像データを記憶する回路であり、画素値対応付け記憶回路412は、測定値と輝度値との対応関係をテーブルとして記憶する回路であり、プログラム記憶回路413は、種々のプログラムを記憶する回路である。   The storage circuit 41 is realized by, for example, a RAM (Random Access Memory), a semiconductor memory element such as a flash memory, a hard disk, an optical disk, or the like. The storage circuit 41 stores, for example, projection data and reconstructed image data. In the present embodiment, the storage circuit 41 includes, for example, an image storage circuit 411, a pixel value association storage circuit 412, and a program storage circuit 413. As will be described in detail later, in the present embodiment, the image storage circuit 411 is a circuit that stores image data of a medical image captured by the gantry device 10, and the pixel value association storage circuit 412 is a measurement value and luminance. The correspondence relationship with values is stored as a table, and the program storage circuit 413 is a circuit that stores various programs.

表示装置42は、各種の情報を表示する。例えば、表示装置42は、処理回路44によって生成された医用画像(CT画像)や、操作者からの各種操作を受け付けるためのGUI(Graphical User Interface)等を出力する。例えば、表示装置42は、液晶ディスプレイやCRT(Cathode Ray Tube)ディスプレイ等によって構成される。   The display device 42 displays various information. For example, the display device 42 outputs a medical image (CT image) generated by the processing circuit 44, a GUI (Graphical User Interface) for receiving various operations from the operator, and the like. For example, the display device 42 includes a liquid crystal display, a CRT (Cathode Ray Tube) display, or the like.

入力装置43は、操作者からの各種の入力操作を受け付け、受け付けた入力操作を電気信号に変換して処理回路44に出力する。例えば、入力装置43は、投影データを収集する際の収集条件や、CT画像を再構成する際の再構成条件、CT画像から後処理画像を生成する際の画像処理条件等を操作者から受け付ける。例えば、入力装置43は、マウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック等により実現される。   The input device 43 receives various input operations from the operator, converts the received input operations into electrical signals, and outputs them to the processing circuit 44. For example, the input device 43 receives from the operator collection conditions when collecting projection data, reconstruction conditions when reconstructing a CT image, image processing conditions when generating a post-process image from a CT image, and the like. . For example, the input device 43 is realized by a mouse, a keyboard, a trackball, a switch, a button, a joystick, or the like.

処理回路44は、入力装置43から出力される入力操作の電気信号に応じてX線CT装置1の全体的な動作を制御する。本実施形態においては、処理回路44は、例えば、システム制御機能441、前処理機能442、再構成処理機能443、画像取得機能444、抽出機能445、算出機能446、生成機能447、及び、表示機能448を有する。   The processing circuit 44 controls the overall operation of the X-ray CT apparatus 1 in accordance with an input operation electric signal output from the input device 43. In the present embodiment, the processing circuit 44 includes, for example, a system control function 441, a preprocessing function 442, a reconstruction processing function 443, an image acquisition function 444, an extraction function 445, a calculation function 446, a generation function 447, and a display function. 448.

図1における実施形態では、システム制御機能441、前処理機能442、再構成処理機能443、画像取得機能444、抽出機能445、算出機能446、生成機能447、及び、表示機能448にて行われる各処理機能は、コンピュータによって実行可能なプログラムの形態で記憶回路41のプログラム記憶回路413に格納されている。処理回路44はプログラムを記憶回路41のプログラム記憶回路413から読み出し、実行することで、各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路は、図1の処理回路44内に示された各機能を有することとなる。なお、図1においては単一の処理回路44にてシステム制御機能441、前処理機能442、再構成処理機能443、画像取得機能444、抽出機能445、算出機能446、生成機能447、及び、表示機能448にて行われる処理機能が実現されるものとして説明したが、複数の独立したプロセッサを組み合わせて処理回路44を構成し、各プロセッサがプログラムを実行することにより機能を実現するものとしても構わない。   In the embodiment in FIG. 1, the system control function 441, the preprocessing function 442, the reconstruction processing function 443, the image acquisition function 444, the extraction function 445, the calculation function 446, the generation function 447, and the display function 448 are performed. The processing functions are stored in the program storage circuit 413 of the storage circuit 41 in the form of a program that can be executed by a computer. The processing circuit 44 is a processor that realizes a function corresponding to each program by reading the program from the program storage circuit 413 of the storage circuit 41 and executing the program. In other words, the processing circuit in the state where each program is read has the functions shown in the processing circuit 44 of FIG. In FIG. 1, the system control function 441, the preprocessing function 442, the reconstruction processing function 443, the image acquisition function 444, the extraction function 445, the calculation function 446, the generation function 447, and the display are performed by a single processing circuit 44. Although the processing function performed by the function 448 has been described as being realized, the processing circuit 44 may be configured by combining a plurality of independent processors, and the function may be realized by each processor executing a program. Absent.

これら各機能のうち、システム制御機能441は、入力装置43を介して操作者から受け付けた入力操作に基づいて、処理回路44の各種機能を制御する。   Among these functions, the system control function 441 controls various functions of the processing circuit 44 based on an input operation received from the operator via the input device 43.

前処理機能442は、データ収集回路18から出力されたデータに対して対数変換処理やオフセット補正処理、チャネル間の感度補正処理、ビームハードニング補正等の前処理を施したデータを生成する。なお、前処理前のデータおよび前処理後のデータを総称して投影データと称する場合もある。   The preprocessing function 442 generates data obtained by performing preprocessing such as logarithmic conversion processing, offset correction processing, interchannel sensitivity correction processing, and beam hardening correction on the data output from the data acquisition circuit 18. Note that data before preprocessing and data after preprocessing may be collectively referred to as projection data.

再構成処理機能443は、前処理機能442にて生成された投影データに対して、フィルタ補正逆投影法や逐次近似再構成法等を用いた再構成処理を行ってCT画像データを生成する。このCT画像データは、例えば、記憶回路41の画像記憶回路411に記憶される。CT画像データは、CT値の集合体であり、本実施形態における複数時相の医用画像の一例である。すなわち、CT画像データは、経時的要素を含む画像として、記憶回路41の画像記憶回路411に格納される。   The reconstruction processing function 443 performs CT reconstruction on the projection data generated by the preprocessing function 442 using a filter-corrected back projection method, a successive approximation reconstruction method, or the like. The CT image data is stored in the image storage circuit 411 of the storage circuit 41, for example. The CT image data is an aggregate of CT values, and is an example of a medical image having a plurality of time phases in the present embodiment. That is, the CT image data is stored in the image storage circuit 411 of the storage circuit 41 as an image including temporal elements.

画像取得機能444、抽出機能445、算出機能446、生成機能447、及び、表示機能448については、詳しくは後述するが、概略的には、画像取得機能444は、例えば、記憶回路41の画像記憶回路411に記憶されている、複数時相のCT画像データを取得する。或いは、再構成処理機能443で生成された複数時相のCT画像データを取得する。   The image acquisition function 444, the extraction function 445, the calculation function 446, the generation function 447, and the display function 448 will be described in detail later, but, roughly, the image acquisition function 444 is, for example, the image storage function of the storage circuit 41. A plurality of time-phase CT image data stored in the circuit 411 is acquired. Alternatively, CT image data of a plurality of time phases generated by the reconstruction processing function 443 is acquired.

抽出機能445は、複数時相のCT画像データから、所定の関心領域の抽出を行う。算出機能446は、所定の関心領域に関するCT値を算出する。生成機能447は、所定の関心領域に沿って、例えば、CT値を空間的位置及び時系列毎に並べた画像を生成する。表示機能448は、生成機能447で生成された画像に基づいて表示画像を生成し、表示装置42の画面に表示させる。   The extraction function 445 extracts a predetermined region of interest from a plurality of time phase CT image data. The calculation function 446 calculates a CT value related to a predetermined region of interest. The generation function 447 generates, for example, an image in which CT values are arranged for each spatial position and time series along a predetermined region of interest. The display function 448 generates a display image based on the image generated by the generation function 447 and displays it on the screen of the display device 42.

上述したように、本実施形態においては、処理回路44は、例えば、プロセッサにより構成される。ここで、プロセッサという文言は、例えば、CPU(Central Processing Unit)、GPU(Graphics Processing Unit)、或いは、特定用途向け集積回路(Application Specific Integrated Circuit:ASIC)、プログラマブル論理デバイス(例えば、単純プログラマブル論理デバイス(Simple Programmable Logic Device:SPLD)、複合プログラマブル論理デバイス(Complex Programmable Logic Device:CPLD)、及び、フィールドプログラマブルゲートアレイ(Field Programmable Gate Array:FPGA))等の回路を意味する。プロセッサは、記憶回路41に保存されたプログラムを読み出して実行することにより機能を実現する。なお、記憶回路41にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成して構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、プロセッサは、プロセッサ単一の回路として構成されている場合に限らず、複数の独立した回路を組み合わせて、1つのプロセッサとして構成し、その機能を実現するようにしてもよい。さらに、図1における複数の構成要素を1つのプロセッサへ統合して、その機能を実現するようにしてもよい。   As described above, in the present embodiment, the processing circuit 44 is configured by a processor, for example. Here, the term processor refers to, for example, a CPU (Central Processing Unit), a GPU (Graphics Processing Unit), an application specific integrated circuit (ASIC), a programmable logic device (for example, a simple programmable logic device) It means circuits such as (Simple Programmable Logic Device: SPLD), Complex Programmable Logic Device (CPLD), and Field Programmable Gate Array (FPGA). The processor implements a function by reading and executing a program stored in the storage circuit 41. Instead of storing the program in the storage circuit 41, the program may be directly incorporated into the processor circuit. In this case, the processor realizes the function by reading and executing the program incorporated in the circuit. Note that the processor is not limited to being configured as a single circuit of the processor, but may be configured as a single processor by combining a plurality of independent circuits to realize the function. Further, a plurality of components in FIG. 1 may be integrated into one processor to realize the function.

なお、図1のX線CT装置1においては、コンソール装置40の構成要素は1つの筐体に格納されているが、これらコンソール装置40の構成要素は、必ずしも、1つの筐体に格納されている必要はなく、別筐体に格納されて複数存在してもよい。例えば、前処理機能442、再構成処理機能443等の処理回路を分散して有するように構成しても構わない。換言すれば、コンソール装置40は、単一のコンソールにて複数の機能を実行するのではなく、複数の機能を別々のコンソールが実行することにしても構わない。   In the X-ray CT apparatus 1 of FIG. 1, the components of the console device 40 are stored in one housing, but the components of the console device 40 are not necessarily stored in one housing. It is not necessary to be present, and there may be a plurality of cases stored in separate cases. For example, the processing circuit such as the preprocessing function 442 and the reconstruction processing function 443 may be distributed. In other words, the console device 40 may execute a plurality of functions by separate consoles instead of executing a plurality of functions by a single console.

また、X線CT装置1には、X線発生装置11とX線検出器12とが一体として被検体Pの周囲を回転するRotate/Rotate−Type(第3世代CT) 、リング状にアレイされた多数のX線検出素子が固定され、X線発生装置11のみが被検体Pの周囲を回転するStationary/Rotate−Type(第4世代CT)等様々なタイプがあり、いずれのタイプでも本実施形態へ適用可能である。   Further, in the X-ray CT apparatus 1, an X-ray generator 11 and an X-ray detector 12 are integrally rotated and rotated in the form of a Rotate / Rotate-Type (third generation CT) that rotates around the subject P. There are various types such as Stationary / Rotate-Type (fourth generation CT) in which a large number of X-ray detection elements are fixed and only the X-ray generation device 11 rotates around the subject P. Applicable to form.

次に、本実施形態に係るX線CT装置1で実行される画像表示処理について説明する。図2は、本実施形態に係る画像表示処理の内容を説明するフローチャートを示す図である。この図2に示すように、本実施形態に係る画像表示処理においては、複数時相の医用画像の中で、関心領域に含まれる組織に芯線となる部位が存在する場合において、芯線に関する測定値を算出し、この測定値を空間的位置及び時系列毎に並べた画像を生成し、この画像に基づく表示画像をユーザに提示する。本実施形態においては、例えば、測定値は、CT画像データに基づいて算出されたCT値などである。以下においては、関心領域にある組織が心血管である場合を例に、一実施形態を説明する。   Next, an image display process executed by the X-ray CT apparatus 1 according to the present embodiment will be described. FIG. 2 is a flowchart illustrating the contents of the image display process according to the present embodiment. As shown in FIG. 2, in the image display processing according to the present embodiment, in a plurality of time-phase medical images, a measurement value related to a core line is present when a portion that becomes a core line exists in a tissue included in a region of interest. Is calculated, an image in which the measured values are arranged for each spatial position and time series is generated, and a display image based on this image is presented to the user. In the present embodiment, for example, the measurement value is a CT value calculated based on CT image data. In the following, an embodiment will be described taking the case where the tissue in the region of interest is a cardiovascular as an example.

図2に示すように、本実施形態に係る画像表示処理は、概略的には、注目血管を抽出する処理(ステップS10)、注目血管の芯線を抽出する処理(ステップS12)、芯線形状を2次元画像化する処理(ステップS14)、芯線上の解析対象の測定値を算出する処理(ステップS16)、芯線に沿ったMAP画像を生成する処理(ステップS18)、及び、このMAP画像を表示する処理(ステップS20)を、備えている。注目血管を抽出する処理(ステップS10)、注目血管の芯線を抽出する処理(ステップS12)、及び、芯線形状を2次元画像化する処理(ステップS14)は、上述した画像取得機能444と抽出機能445とにより実現される処理であり、芯線上の解析対象の測定値を算出する処理(ステップS16)は、上述した算出機能446により実現される処理であり、芯線に沿ったMAP画像を生成する処理(ステップS18)は、上述した生成機能447により実現される処理であり、MAP画像を表示する処理(ステップS20)は、上述した表示機能448により実現される処理である。   As shown in FIG. 2, the image display processing according to the present embodiment generally includes processing for extracting a target blood vessel (step S10), processing for extracting a core line of the target blood vessel (step S12), and a core shape of 2 A process for forming a dimensional image (step S14), a process for calculating a measurement value to be analyzed on the core line (step S16), a process for generating a MAP image along the core line (step S18), and displaying this MAP image A process (step S20) is provided. The process of extracting the target blood vessel (step S10), the process of extracting the core line of the target blood vessel (step S12), and the process of converting the core line shape into a two-dimensional image (step S14) are the image acquisition function 444 and the extraction function described above. The process (step S16) for calculating the measurement value of the analysis target on the core line is a process realized by the calculation function 446 described above, and generates a MAP image along the core line. The process (step S18) is a process realized by the generation function 447 described above, and the process of displaying a MAP image (step S20) is a process realized by the display function 448 described above.

各処理をより詳しく説明すると、まず、注目血管を抽出する処理(ステップS10)では、複数時相の医用画像を例えば画像記憶回路411から取得して、この医用画像から、既存の血管抽出アルゴリズムを用いて、心血管など、人体内の血管領域を各時相で抽出する。抽出する血管は、ユーザが観察したい血管を、マウス等のユーザインターフェースを用いて、入力装置43から操作して選択しても良いし、心臓の血管解析で一般的に観察される主要な血管を自動的に選択しても良い。   To describe each process in more detail, first, in the process of extracting the target blood vessel (step S10), a medical image of a plurality of time phases is acquired from, for example, the image storage circuit 411, and an existing blood vessel extraction algorithm is extracted from this medical image. Used to extract blood vessel regions in the human body, such as cardiovascular, at each time phase. The blood vessel to be extracted may be selected by operating the input device 43 using a user interface such as a mouse, or a main blood vessel that is generally observed in the blood vessel analysis of the heart. It may be selected automatically.

次に、注目血管の芯線を抽出する処理(ステップS12)では、注目血管を抽出する処理(ステップS10)で抽出した各時相の血管領域に対して、血管の芯線に該当する曲線を既存の血管芯線抽出アルゴリズムを用いて抽出する。ここで、抽出される血管の長さは各時相によって差がある場合がある。そのため、特定の1つまたは複数の時相で抽出された血管の芯線の座標を他の時相に位置合わせすることで、各時相において血管の芯線上の各座標の該当位置を求めるものとする。これにより、芯線上の各点の座標は各時相間で対応付けられる。   Next, in the process of extracting the core line of the target blood vessel (step S12), a curve corresponding to the core line of the blood vessel is added to the existing blood vessel region in each time phase extracted in the process of extracting the target blood vessel (step S10). Extraction is performed using a blood vessel core line extraction algorithm. Here, the length of the extracted blood vessel may vary depending on each time phase. Therefore, by aligning the coordinates of the blood vessel core line extracted in one or more specific time phases with other time phases, the corresponding position of each coordinate on the blood vessel core line in each time phase is obtained. To do. Thereby, the coordinate of each point on a core line is matched between each time phase.

また、上記の処理で用いる時相は、ユーザが求める血管の長さが最低限抽出できている時相を選択してもよいし、抽出できた芯線の長さが最も長い時相を自動で選択してもよい。あるいは全時相の中で最高画素値となる画素の画素値を各画素の画素値として採用するTimeMIP画像を生成して用いてもよいし、最高画素値の代わりに平均値などを用いるTimeStack画像を生成して用いてもよい。また、採用した血管芯線抽出アルゴリズムによっては、注目血管を抽出する処理(ステップS10)と、注目血管の芯線を抽出する処理(ステップS12)とを、同時に実行してもよい。   In addition, the time phase used in the above processing may be selected as the time phase in which the length of the blood vessel desired by the user can be extracted at the minimum, or the time phase in which the extracted core length is the longest is automatically selected. You may choose. Alternatively, a TimeMIP image that uses the pixel value of the pixel that is the highest pixel value in all time phases as the pixel value of each pixel may be generated and used, or a TimeStack image that uses an average value or the like instead of the highest pixel value May be generated and used. Further, depending on the blood vessel core line extraction algorithm employed, the process of extracting the target blood vessel (step S10) and the process of extracting the core line of the target blood vessel (step S12) may be performed simultaneously.

次に、芯線形状を2次元画像化する処理(ステップS14)では、ステップS10及びステップS12で抽出した3次元の血管の芯線を、2次元で扱うために、抽出した血管の芯線の形状に従って、血管の芯線に対するCPR(Curved Planar Reconstruction)画像またはSPR(Stretched Curved Planar Reconstruction)を生成して、2次元画像に変換する。CPR画像は、任意の曲面で画像を切断した断面を再構成するMPR(Multi Planar Reconstruction)画像の一種であり、本実施形態では、血管の芯線に沿った曲面で画像を切断した断面の再構成画像を指す。   Next, in the process of converting the core line shape into a two-dimensional image (step S14), in order to handle the three-dimensional blood vessel core line extracted in step S10 and step S12 in two dimensions, according to the extracted blood vessel core line shape, A CPR (Curved Planar Reconstruction) image or an SPR (Stretched Curved Planar Reconstruction) for a blood vessel core line is generated and converted into a two-dimensional image. The CPR image is a kind of MPR (Multi Planar Reconstruction) image for reconstructing a cross section obtained by cutting an image with an arbitrary curved surface. In this embodiment, the reconstruction of a cross section obtained by cutting the image with a curved surface along the core line of the blood vessel is used. Refers to the image.

CPR画像とSPR画像のどちらを生成するかはユーザが任意に選択できる。後述する芯線に沿ったMAP画像を生成する処理(ステップS18)において、血管の曲線形状の様子を観察しやすくするためにCPR画像を選択しても良いし、空間的な情報として血管の進行方向の長さのみ観察したい場合はSPR画像を選択してもよい。図3は、抽出した血管芯線を2次元画像としてSPR画像に変換した態様の一例を示している。   The user can arbitrarily select whether to generate a CPR image or an SPR image. In the process of generating a MAP image along the core line (step S18), which will be described later, a CPR image may be selected in order to make it easy to observe the curved shape of the blood vessel, and the direction of blood vessel progression as spatial information When it is desired to observe only the length of the SPR image, the SPR image may be selected. FIG. 3 shows an example of a mode in which the extracted blood vessel core line is converted into an SPR image as a two-dimensional image.

次に、芯線上の解析対象値を算出する処理(ステップS16)では、注目血管の芯線を抽出する処理(ステップS12)、又は、芯線形状を2次元画像化する処理(ステップS14)で得られた芯線上の各点の座標において、解析対象となる測定値を各時相で算出する。解析対象とする測定値は、例えばCT値であってもよいし、血管径など血管の芯線で通常用いられる計測値のどれかでもよい。また解析対象となる測定値の算出では、芯線上の各点の1つの画素から測定値を算出してもよいし、芯線上の各点の周辺にある複数の画素を選択して、選択した複数の画素の平均値や最大値などから測定値を算出してもよい。   Next, the process of calculating the analysis target value on the core line (step S16) is obtained by the process of extracting the core line of the blood vessel of interest (step S12) or the process of converting the core line shape into a two-dimensional image (step S14). At the coordinates of each point on the core line, the measurement value to be analyzed is calculated at each time phase. The measurement value to be analyzed may be, for example, a CT value, or any of the measurement values that are normally used for blood vessel cores, such as the blood vessel diameter. In the calculation of the measurement value to be analyzed, the measurement value may be calculated from one pixel at each point on the core line, or a plurality of pixels around each point on the core line are selected and selected. The measurement value may be calculated from an average value or a maximum value of a plurality of pixels.

次に、芯線に沿ったMAP画像を生成する処理(ステップS18)では、芯線形状を2次元画像化する処理(ステップS14)で生成した複数時相のCPR画像またはSPR画像の中から、代表となる時相を選択する。選択した時相のCPR画像またはSPR画像における血管芯線の進行方向を画像のx軸、時間軸方向をy軸方向に取った2次元画像の各画素の画素値として、芯線上の解析対象の測定値を算出する処理(ステップS16)で算出した解析対象の測定値(CT値など)を割り当てたMAP画像を生成する。つまり、ステップS16で算出した各時相の解析対象の測定値に基づいて、1時相分の芯線に沿った測定値をx軸と平行に配置したものを、y軸と平行な方向に時相順で並べて配置したMAP画像を生成する。結果として上記の処理で生成したMAP画像は、芯線形状を2次元画像化する処理(ステップS14)で生成したSPR画像の芯線上の測定値を抽出して順々に並べたような画像が生成できることになる。   Next, in the process of generating a MAP image along the core line (step S18), the representative and the CPR images or SPR images of a plurality of time phases generated in the process of converting the core line shape into a two-dimensional image (step S14) Select a time phase. Measurement of the analysis target on the core line, using the progression direction of the blood vessel core line in the selected time phase CPR image or SPR image as the pixel value of each pixel of the two-dimensional image taking the time axis direction as the y-axis direction. A MAP image to which the analysis target measurement value (CT value or the like) calculated in the value calculation processing (step S16) is assigned is generated. That is, based on the measurement value of the analysis target of each time phase calculated in step S16, the measurement value along the core for one time phase is arranged in parallel with the x axis in the direction parallel to the y axis. MAP images arranged in phase order are generated. As a result, the MAP image generated by the above process generates an image in which the measurement values on the core line of the SPR image generated by the process of converting the core line shape into a two-dimensional image (step S14) are extracted and arranged in order. It will be possible.

図4は、上記の処理の概要を模式的に説明する図である。この図4に例示するように、ステップS18の処理では、x軸方向に関心領域の測定値が割り当てられた帯状の領域が複数生成され、この複数の帯状の領域がy軸方向に時系列に並んだMAP画像が生成される。本実施形態においては、測定値がCT値であるので、CT値の集合で、MAP画像が構成されることとなる。なお、生成したMAP画像のx軸、y軸はユーザの好みで逆転させてもよい。つまり、x軸を時間軸、y軸を血管の芯線の進行方向としてもよい。   FIG. 4 is a diagram schematically illustrating the outline of the above processing. As illustrated in FIG. 4, in the process of step S <b> 18, a plurality of band-like areas to which the measurement value of the region of interest is assigned in the x-axis direction are generated, and the plurality of band-like areas are time-sequentially in the y-axis direction. A side-by-side MAP image is generated. In the present embodiment, since the measurement value is a CT value, a MAP image is composed of a set of CT values. Note that the x-axis and y-axis of the generated MAP image may be reversed according to the user's preference. That is, the x axis may be the time axis and the y axis may be the blood vessel core line progression direction.

次に、MAP画像を表示する処理(ステップS20)では、芯線に沿ったMAP画像を生成する処理(ステップS18)で生成したMAP画像の各測定値に、所定の輝度値や色を割り当てた上で表示画像を生成し、表示装置42の画面に表示する。図5は、生成したMAP画像を表示装置42の画面に表示した一例を示している。生成したMAP画像だけでも空間的情報を有しているが、図5の例では、より直感的に理解しやすいように代表的な時相におけるSPR画像を参考画像として、画面に並べて表示している。SPR画像の代わりに、処理に用いた元の医用画像である3次元画像を参考画像として、並べて表示してもよい。代表的な時相はユーザが任意に指定してもよいし、抽出できた芯線の長さが最も長い時相を自動で選択してもよい。あるいは、生成したMAP画像上の画素をクリックすることで、その画素が対応する時相をユーザが選択できるようにしてもよい。すなわち、コンソール装置40が取得した医用画像、及び/又は、この取得した医用画像に基づいて生成された画像を、参考画像として生成し、MAP画像に並べて表示して、ユーザの利便性を向上させることも可能である。無論、このSPR画像などの参考画像は、省略することも可能である。   Next, in the process of displaying the MAP image (step S20), a predetermined luminance value or color is assigned to each measurement value of the MAP image generated in the process of generating the MAP image along the core line (step S18). Then, a display image is generated and displayed on the screen of the display device 42. FIG. 5 shows an example in which the generated MAP image is displayed on the screen of the display device 42. The generated MAP image alone has spatial information, but in the example of FIG. 5, the SPR image at a representative time phase is displayed as a reference image side by side on the screen so that it can be understood more intuitively. Yes. Instead of the SPR image, a three-dimensional image that is an original medical image used for processing may be displayed side by side as a reference image. A representative time phase may be arbitrarily designated by the user, or a time phase having the longest extracted core wire may be automatically selected. Alternatively, by clicking a pixel on the generated MAP image, the user may be able to select a time phase corresponding to the pixel. That is, a medical image acquired by the console device 40 and / or an image generated based on the acquired medical image is generated as a reference image and displayed side by side on the MAP image to improve user convenience. It is also possible. Of course, the reference image such as the SPR image can be omitted.

本実施形態に係るX線CT装置1で生成したMAP画像及び表示画像は、通常の2次元画像と同様の画像処理が可能である。そのため、ウィンドウ幅(WW)/ウィンドウレベル(WL)の調整やカラーマップ表示、スムージングなどの画像フィルタ処理、指定した画素の画素値表示、指定した画素を基準点とした領域抽出処理などが可能である。   The MAP image and display image generated by the X-ray CT apparatus 1 according to the present embodiment can be subjected to image processing similar to a normal two-dimensional image. Therefore, adjustment of window width (WW) / window level (WL), color map display, image filter processing such as smoothing, pixel value display of designated pixels, region extraction processing using designated pixels as reference points, and the like are possible. is there.

また、生成したMAP画像上で時相指定用の曲線ROIを設置し、曲線ROIの各座標に該当する時相の画素値を用いた画像を新たに生成することができる。図6は、曲線ROIを画面の表示画像に重畳させて表示させた画像の一例を示している。つまり、曲線という連続性のあるROIで時相を選択することで、ある程度時間軸方向の連続性を保つという制限のもと、TimeMIP画像のような画像が生成できる。通常のTimeMIP画像の場合、体動や心拍等による一時的なノイズやアーティファクトの影響で極端に画素値が高くなる場合があるが、図6に示すように、曲線ROIを指定することで、用いる時相を指定できるようになり、ノイズやアーティファクトの影響を避けることができる。時相指定用の曲線ROIの設置位置は、マウスなどのユーザインターフェースである入力装置43を用いて、ユーザが全ての座標を指定して設定してもよいし、各時相で最高画素となる時相をそれぞれ初期設定位置として設定した後、ユーザが好みで変更可能としてもよい。   Further, a curve ROI for designating a time phase can be set on the generated MAP image, and an image using a pixel value of a time phase corresponding to each coordinate of the curve ROI can be newly generated. FIG. 6 shows an example of an image displayed by superimposing the curve ROI on the display image on the screen. That is, an image such as a TimeMIP image can be generated under the restriction that continuity in the time axis direction is maintained to some extent by selecting a time phase with a continuous ROI called a curve. In the case of a normal TimeMIP image, the pixel value may become extremely high due to the influence of temporary noise or artifacts due to body movement or heartbeat, but it is used by specifying a curve ROI as shown in FIG. You can now specify the time phase and avoid the effects of noise and artifacts. The installation position of the curve ROI for designating the time phase may be set by the user by designating all coordinates using the input device 43 such as a user interface such as a mouse, or becomes the highest pixel in each time phase. After setting each time phase as an initial setting position, the user may be able to change it as desired.

図7は、ユーザの指定に基づいて曲線ROIが表示できるX線CT装置1の全体構成を説明するブロック図であり、上述した図1に対応する図である。この図7に示すように、このX線CT装置1においては、処理回路44に、追加的に、指定機能449と、測定値取得機能450と、グラフ表示機能451とが設けられる。これらの機能も、処理回路44がプログラム記憶回路413に格納されているプログラムを読み込んで実行することにより実現される。   FIG. 7 is a block diagram for explaining the overall configuration of the X-ray CT apparatus 1 capable of displaying a curved line ROI based on the user's designation, and corresponds to FIG. 1 described above. As shown in FIG. 7, in the X-ray CT apparatus 1, the processing circuit 44 is additionally provided with a designation function 449, a measurement value acquisition function 450, and a graph display function 451. These functions are also realized by the processing circuit 44 reading and executing a program stored in the program storage circuit 413.

指定機能449は、ユーザに、入力装置43の操作に基づいて、表示装置42に表示された表示画像における、時相と空間的位置を指定させる機能を実現する。測定値取得機能450は、ユーザが指定した時相と空間的位置に基づいて、表示画像上における、指定された時相と空間的位置に対応するMAP画像の測定値を取得する。グラフ表示機能451は、測定値取得機能450が取得した測定値に基づいて、曲線ROIをグラフとして生成して、再び、表示画面に重畳して表示する。   The designation function 449 realizes a function that allows the user to designate a time phase and a spatial position in the display image displayed on the display device 42 based on the operation of the input device 43. The measurement value acquisition function 450 acquires the measurement value of the MAP image corresponding to the designated time phase and spatial position on the display image based on the time phase and spatial position designated by the user. The graph display function 451 generates a curve ROI as a graph based on the measurement value acquired by the measurement value acquisition function 450, and displays the curve ROI superimposed on the display screen again.

例えば、デフォルトで表示された曲線ROIにノイズが含まれている場合、ユーザは、これに変えて、適正な時相と空間的位置とを指定することにより、ノイズを排除して、指定した時相と空間的位置の測定値を用いてグラフを生成することができるようになる。グラフは、曲線ROIに限られるものでは無く、折れ線グラフや棒グラフ、点による分布図などの種々の態様で生成し、表示画像に重ねて表示させることができる。   For example, when noise is included in the curve ROI displayed by default, the user can specify the appropriate time phase and spatial position instead of the noise to specify the time. Graphs can be generated using phase and spatial position measurements. The graph is not limited to the curved line ROI, but can be generated in various forms such as a line graph, a bar graph, and a distribution diagram by points, and can be displayed on the display image.

なお、このMAP画像に基づく表示画像を、表示装置42の画面に表示させる際は、白黒のモノトーンによる表示でもよいし、各種の色が付いたカラーによる表示でもよい。白黒のモノトーンによる画像表示の場合、MAP画像における各画素の測定値を輝度値に変換して表示画像を生成し、表示装置42の画面に表示させる。   When the display image based on the MAP image is displayed on the screen of the display device 42, it may be displayed in black and white monotone, or may be displayed in colors with various colors. In the case of displaying an image in black and white monotone, a measurement value of each pixel in the MAP image is converted into a luminance value, and a display image is generated and displayed on the screen of the display device 42.

図8(a)は、白黒のモノトーンのように、単色のグラデーションでMAP画像に基づく表示画像を生成する場合における、測定値と輝度値との対応関係を説明する図である。ここでは、一例として、測定値がCT値である場合を示している。この図8(a)の例では、CT値が低いほど輝度値が低く、CT値が高いほど輝度値が高くなるように対応付けられている。例えば、CT値が−1000の場合は最も輝度値が低く、黒が対応付けられており、CT値が+1000の場合は最も輝度値が高く、白が対応付けられている。このため、輝度値が高くて明るい部分はCT値が高いという前提で、ユーザは表示画像を観察できる。本実施形態においては、この測定値と輝度値との対応関係は、テーブルとして、記憶回路41の画素値対応付け記憶回路412に格納されている。このため、処理回路44は、記憶回路41から、このテーブルを読み込んで、このテーブルに基づいて、MAP画像に割り付けられている測定値に対応付けられた輝度値を算出し、単色のグラデーションで表現された表示画像を生成して、表示装置42に表示させる。   FIG. 8A is a diagram for explaining a correspondence relationship between a measurement value and a luminance value when a display image based on a MAP image is generated with a monochrome gradation like a black and white monotone. Here, as an example, the case where the measured value is a CT value is shown. In the example of FIG. 8A, the brightness value is lower as the CT value is lower, and the brightness value is higher as the CT value is higher. For example, when the CT value is −1000, the luminance value is the lowest and associated with black, and when the CT value is +1000, the luminance value is the highest and associated with white. For this reason, the user can observe the display image on the premise that the bright value is high and the bright portion has a high CT value. In the present embodiment, the correspondence relationship between the measurement value and the luminance value is stored in the pixel value association storage circuit 412 of the storage circuit 41 as a table. For this reason, the processing circuit 44 reads this table from the storage circuit 41, calculates the luminance value associated with the measurement value assigned to the MAP image based on this table, and expresses it with a monochrome gradation. The displayed display image is generated and displayed on the display device 42.

各種の色が付いたカラーによる画像表示の場合、MAP画像における測定値を色に変換して表示画像を生成し、表示装置42の画面に表示させる。図8(b)では、カラーのグラデーションで、MAP画像に基づく表示画像を生成する場合における、測定値と色との対応関係を説明する図である。ここでも、一例として、測定値がCT値である場合を示している。この図8(b)の例では、CT値が低い場合は白色が対応付けられており、CT値が高い場合は赤色が対応付けられている。このため、色が赤みを帯びるに従ってCT値が高いという前提で、ユーザは表示画像を観察できる。本実施形態においては、この測定値と色との対応関係は、テーブルとして、記憶回路41の画素値対応付け記憶回路412に格納されている。このため、処理回路44は、記憶回路41から、このテーブルを読み込んで、このテーブルに基づいて、MAP画像に割り付けられている測定値に対応付けられた色を算出し、カラーの表示画像を生成して、表示装置42に表示させる。   In the case of displaying an image in a color with various colors, the measurement value in the MAP image is converted into a color to generate a display image and displayed on the screen of the display device 42. FIG. 8B is a diagram for explaining the correspondence between measured values and colors when a display image based on a MAP image is generated with color gradation. Here, as an example, a case where the measurement value is a CT value is shown. In the example of FIG. 8B, when the CT value is low, white is associated, and when the CT value is high, red is associated. For this reason, the user can observe the display image on the assumption that the CT value increases as the color becomes reddish. In the present embodiment, the correspondence between the measurement value and the color is stored in the pixel value association storage circuit 412 of the storage circuit 41 as a table. Therefore, the processing circuit 44 reads this table from the storage circuit 41, calculates a color associated with the measurement value assigned to the MAP image based on this table, and generates a color display image. Then, it is displayed on the display device 42.

以上のように、本実施形態に係るX線CT装置1によれば、例えば、1枚のMAP画像がもつ空間的位置情報と時系列的時間情報から血流解析が可能となる。このため、血管の狭窄により造影剤の浸透がどのように妨げられるかを観察できるようになる。また、動脈瘤への造影剤の流入がある場合、造影剤の流れが変化するはずであるから、動脈瘤の状態解析や経過観察にも利用できる。また血管内に存在するプラークへの造影剤の浸透の様子から、不安定プラークかどうかなど、プラークの状態解析にも利用可能である。   As described above, according to the X-ray CT apparatus 1 according to the present embodiment, for example, blood flow analysis can be performed from the spatial position information and time-series time information of one MAP image. For this reason, it becomes possible to observe how the penetration of the contrast agent is hindered by the narrowing of the blood vessel. Further, when the contrast agent flows into the aneurysm, the flow of the contrast agent should change, and therefore, it can be used for aneurysm state analysis and follow-up observation. It can also be used to analyze the state of a plaque such as whether it is an unstable plaque from the state of penetration of a contrast medium into a plaque present in a blood vessel.

換言すれば、関心領域の空間的情報を保持したまま、時相の変化に伴う観察対象の測定値の変化を確認できる。これは、従来技術であるTDC解析と異なり、空間的な情報をも含んでいるので、例えば、造影剤の浸透の様子から血流解析や血管プラークの状態解析などが、1つの画像を見るだけで可能となる。つまり、従来のMAP画像の別形態とも言える。また、解析対象である測定値の変化を2次元画像として、TDC解析よりも1次元大きい範囲で捉えるため、TDC解析のようにROIの設定に適した座標を探すための手間が少なくなる。またTDC解析のようにグラフとして表示するよりも、表示画像の輝度や色として確認できるため、直感的に解析しやすくなる。   In other words, it is possible to confirm the change in the measurement value of the observation object accompanying the change in the time phase while maintaining the spatial information of the region of interest. Unlike TDC analysis, which is a conventional technique, this also includes spatial information. For example, blood flow analysis or vascular plaque state analysis from the state of infiltration of a contrast agent only looks at one image. Is possible. That is, it can be said to be another form of the conventional MAP image. In addition, since the change of the measurement value to be analyzed is captured as a two-dimensional image in a range that is one dimension larger than the TDC analysis, the trouble of searching for coordinates suitable for setting the ROI as in the TDC analysis is reduced. Moreover, since it can confirm as the brightness | luminance and color of a display image rather than displaying as a graph like TDC analysis, it becomes easy to analyze intuitively.

なお、本実施形態においては、時系列的な情報として時相を用いたが、それ以外の情報を、時系列的な情報とすることもできる。例えば、心位相や収縮期から拡張期までといった位相情報や、心拍数、撮影時期等を時系列的な情報と置き換えても、同様の処理が可能である。   In the present embodiment, the time phase is used as time-series information, but other information may be time-series information. For example, the same processing can be performed by replacing phase information such as cardiac phase and systole to diastole, heart rate, imaging time, etc. with time-series information.

また、本実施形態では、血管に対して、上述した画像表示処理を行ったが、血管以外にも芯線を抽出可能な組織であれば、上述した画像表示処理を同様に適用可能である。例えば、気管支の解析にも、上述した画像表示処理を適用可能である。この場合、例えば、血管の芯線の代わりに気管支の芯線を抽出し、CT値の代わりに気管支の直径や単位長さ当たりの体積等の計測値を解析対象の測定値として用いれば、上述した画像表示処理が可能となる。   In the present embodiment, the above-described image display processing is performed on a blood vessel. However, the above-described image display processing can be similarly applied to any tissue other than blood vessels that can extract a core line. For example, the image display process described above can also be applied to bronchial analysis. In this case, for example, if the bronchus core line is extracted instead of the blood vessel core line, and the measurement value such as the bronchus diameter or the volume per unit length is used as the analysis target measurement value instead of the CT value, the image described above Display processing is possible.

〔第2実施形態〕
上述した第1実施形態に係るX線CT装置1においては、関心領域に含まれる組織の芯線を抽出し、この芯線に関する解析対象の測定値を算出し、これを空間的位置及び時系列毎に並べてMAP画像を生成することとしたが、第2実施形態では、複数時相の画像の中で、関心領域に含まれる組織に芯線となる部位が存在しない場合でも、観察対象となる測定値を算出し、空間的位置及び時系列毎に並べた画像を生成できるようにしたものである。以下、上述した第1実施形態と異なる部分を説明する。
[Second Embodiment]
In the X-ray CT apparatus 1 according to the first embodiment described above, the core line of the tissue included in the region of interest is extracted, the measurement value of the analysis target related to this core line is calculated, and this is calculated for each spatial position and time series. Although the MAP images are generated side by side, in the second embodiment, the measurement value to be observed is obtained even in the case where there is no core part in the tissue included in the region of interest among the images of a plurality of time phases. It is possible to generate an image calculated and arranged for each spatial position and time series. Hereinafter, a different part from 1st Embodiment mentioned above is demonstrated.

第2実施形態に係る画像診断装置の一例であるX線CT装置1の内部構成は、上述した第1実施形態の図1と同様である。本実施形態においては、芯線が存在しない関心領域である組織として、肺を例に説明する。   The internal configuration of an X-ray CT apparatus 1 that is an example of an image diagnostic apparatus according to the second embodiment is the same as that of FIG. 1 of the first embodiment described above. In the present embodiment, a lung will be described as an example of a tissue that is a region of interest where no core wire exists.

図9は、本実施形態に係るコンソール装置40で実行される画像表示処理の内容を説明するフローチャートを示す図である。この図9に示すように、本実施形態に係る画像表示処理は、概略的には、軸を設定する処理(ステップS30)、軸上の解析対象の測定値を算出する処理(ステップS32)、軸に沿ったMAP画像を生成する処理(ステップS34)、及び、MAP画像を表示する処理(ステップS36)を、備えている。軸を設定する処理(ステップS30)は、上述した画像取得機能444と抽出機能445とにより実現される処理であり、軸上の解析対象の測定値を算出する処理(ステップS32)は、上述した算出機能446により実現される処理であり、軸に沿ったMAP画像を生成する処理(ステップS34)は、上述した生成機能447により実現される処理であり、MAP画像を表示する処理(ステップS36)は、上述した表示機能448により実現される処理である。   FIG. 9 is a flowchart illustrating the contents of image display processing executed by the console device 40 according to the present embodiment. As shown in FIG. 9, the image display processing according to the present embodiment roughly includes processing for setting an axis (step S30), processing for calculating a measurement value of an analysis target on the axis (step S32), A process of generating a MAP image along the axis (step S34) and a process of displaying the MAP image (step S36) are provided. The processing for setting the axis (step S30) is realized by the image acquisition function 444 and the extraction function 445 described above, and the processing for calculating the measurement value of the analysis target on the axis (step S32) is described above. The process realized by the calculation function 446 and the process for generating the MAP image along the axis (step S34) is the process realized by the generation function 447 described above, and the process for displaying the MAP image (step S36). Is a process realized by the display function 448 described above.

各処理をより詳しく説明すると、まず、軸を設定する処理(ステップS30)では、肺が含まれた複数時相の医用画像を例えば画像記憶回路411から取得して、この複数時相の医用画像の中から特定の時相を指定し、2次元画像または3次元画像として、表示装置42の画面に肺領域を表示する。画面に表示された肺領域に対して、ユーザがマウスなどのユーザインターフェースを用いて、入力装置43から線分または曲線のROIを設定し、以降、そのROIを軸とみなす。さらに、特定時相の画像で設定した線分または曲線ROIの各座標に対応する座標を、位置合わせ処理などにより、その他の時相でも算出し、時相間で対応付けて記録しておく。図10は、肺の画像にユーザが設定したROIを表示した画像の一例を示している。この図10の例では、曲線ROI上のドットが、ユーザが指定した座標を表している。   To describe each process in more detail, first, in the process of setting an axis (step S30), a plurality of time-phase medical images including the lungs are acquired from, for example, the image storage circuit 411, and the plurality of time-phase medical images are acquired. A specific time phase is designated from among the images, and the lung region is displayed on the screen of the display device 42 as a two-dimensional image or a three-dimensional image. For a lung region displayed on the screen, a user sets a ROI of a line segment or a curve from the input device 43 using a user interface such as a mouse, and thereafter, the ROI is regarded as an axis. Furthermore, the coordinates corresponding to the coordinates of the line segment or the curve ROI set in the image of the specific time phase are also calculated in other time phases by the alignment processing or the like, and are recorded in association with each other. FIG. 10 shows an example of an image in which the ROI set by the user is displayed on the lung image. In the example of FIG. 10, the dots on the curve ROI represent the coordinates designated by the user.

次に、軸上の解析対象の測定値を算出する処理(ステップS32)では、上述した軸を設定する処理(ステップS34)で設定した軸に沿って、解析対象となる測定値を各時相で算出する。ここで解析対象となる測定値の例として、肺の幅や肺の断面積、単位長さ当たりの体積などの計測値が考えられる。また、軸上の各点における計測値の算出に用いる断面は、軸の接線に対して垂直な面でも良いし、人体の高さ方向に垂直な面でもよい。図10の例では、曲線ROIを軸として肺に設定した場合に、肺の断面積を計測値とし、人体の高さ方向に垂直な面を計測値の算出に用いる断面としている。   Next, in the process of calculating the measurement value of the analysis target on the axis (step S32), the measurement value to be analyzed is set to each time phase along the axis set in the process of setting the axis (step S34) described above. Calculate with Here, as examples of measurement values to be analyzed, measurement values such as lung width, lung cross-sectional area, and volume per unit length are conceivable. Further, the cross section used for calculation of the measured value at each point on the axis may be a plane perpendicular to the axis tangent or a plane perpendicular to the height direction of the human body. In the example of FIG. 10, when the lung is set with the curve ROI as an axis, a cross-sectional area of the lung is used as a measurement value, and a plane perpendicular to the height direction of the human body is used as a cross-section for calculation of the measurement value.

次に、軸に沿ったMAP画像を生成する処理(ステップS34)では、軸を設定する処理(ステップS30)で設定した軸の進行方向を画像のx軸、時間軸方向をy軸方向に取った2次元画像の各画素の画素値として、軸上の解析対象の測定値を算出する処理(ステップS32)で算出した各時相の解析対象である測定値を割り当てた画像を生成する。この軸に沿ったMAP画像を生成する処理(ステップS34)で生成されるMAP画像は、上述した第1実施形態の図2における芯線に沿ったMAP画像を生成する処理(ステップS18)と同様である。すなわち、図4で説明したような画像生成処理が行われる。   Next, in the process of generating the MAP image along the axis (step S34), the axis traveling direction set in the axis setting process (step S30) is taken as the x-axis of the image, and the time axis direction is taken as the y-axis direction. As the pixel value of each pixel of the two-dimensional image, an image to which the measurement value that is the analysis target of each time phase calculated in the process of calculating the measurement value of the analysis target on the axis (step S32) is generated. The MAP image generated by the process of generating the MAP image along the axis (step S34) is the same as the process of generating the MAP image along the core line in FIG. 2 of the first embodiment (step S18). is there. That is, the image generation process described with reference to FIG. 4 is performed.

次に、MAP画像を表示する処理(ステップS36)では、軸に沿ったMAP画像を生成する処理(ステップS34)で生成したMAP画像に基づいて生成した表示画像を、表示装置42の画面に表示する。このMAP画像を表示する処理(ステップS36)は、上述した第1実施形態の図2におけるMAP画像を表示する処理(ステップS20)と同様である。すなわち、図5に示したようなMAP画像が、表示装置42に表示される。なお、上述した第1実施形態と同様に、より直感的に理解しやすいように、MAP画像と並べて、代表的な時相におけるSPR画像や組織の3次元画像を、参考画像として並べて表示するようにしてもよい。   Next, in the process of displaying the MAP image (step S36), the display image generated based on the MAP image generated in the process of generating the MAP image along the axis (step S34) is displayed on the screen of the display device 42. To do. The process of displaying the MAP image (step S36) is the same as the process of displaying the MAP image (step S20) in FIG. 2 of the first embodiment described above. That is, a MAP image as shown in FIG. 5 is displayed on the display device 42. As in the first embodiment described above, an SPR image and a three-dimensional image of a tissue at a representative time phase are displayed side by side as reference images so as to be more intuitively understood. It may be.

以上のように、本実施形態に係るX線CT装置1によれば、肺などの芯線が存在しない組織が関心領域である場合でも、芯線の代わりに軸を設定して、この軸に関する測定値を算出し、この測定値を空間的位置及び時系列毎に並べたMAP画像を生成することとした。このため、芯線が存在しない組織が関心領域である場合でも、MAP画像を生成し、表示することができるようになる。これにより、例えば肺が関心領域の組織である場合には、呼吸下における肺の動きを可視化できるようになる。具体的には、例えば肺に腫瘍が生じた場合、呼吸時の肺の動きに影響を与えることになるが、その変化の様子を1枚のMAP画像として生成して表示することにより、ユーザが容易に確認できるようにすることができる。   As described above, according to the X-ray CT apparatus 1 according to the present embodiment, even when a tissue such as the lung where no core wire is present is a region of interest, an axis is set instead of the core wire, and a measurement value related to this axis is set. And a MAP image in which the measured values are arranged for each spatial position and time series is generated. For this reason, even when a tissue having no core line is a region of interest, a MAP image can be generated and displayed. Thereby, for example, when the lung is the tissue of the region of interest, the movement of the lung under respiration can be visualized. Specifically, for example, when a tumor occurs in the lung, the movement of the lung during breathing will be affected. By generating and displaying the change as a single MAP image, the user can It can be easily confirmed.

〔第3実施形態〕
上述した第1実施形態では、関心領域として組織の芯線を抽出し、第2実施形態では、関心領域として組織の軸を設定することとしたが、医用画像を用いた診断においては、関心領域が心臓や肺などの組織の面である場合もある。そこで、第3実施形態に係るX線CT装置1では、複数時相の医用画像において、関心領域となる面、特に曲面がある場合において、観察対象となる測定値を算出し、この測定値を空間的位置及び時系列毎に並べた画像を生成することとした。以下、上述した第1実施形態及び第2実施形態と異なる部分を説明する。
[Third Embodiment]
In the first embodiment described above, the tissue core line is extracted as the region of interest, and in the second embodiment, the tissue axis is set as the region of interest. However, in the diagnosis using the medical image, the region of interest is It may also be in terms of tissues such as the heart and lungs. Therefore, in the X-ray CT apparatus 1 according to the third embodiment, in a medical image having a plurality of phases, when there is a surface to be a region of interest, particularly a curved surface, a measurement value to be observed is calculated, and this measurement value is calculated. An image arranged for each spatial position and time series is generated. Hereinafter, a different part from 1st Embodiment and 2nd Embodiment which were mentioned above is demonstrated.

第3実施形態に係る画像診断装置の一例であるX線CT装置1の内部構成は、上述した第1実施形態の図1と同様である。以下においては、曲面を有する関心領域として心筋を例に、本実施形態を説明する。   The internal configuration of an X-ray CT apparatus 1 that is an example of an image diagnostic apparatus according to the third embodiment is the same as that of FIG. 1 of the first embodiment described above. In the following, the present embodiment will be described using the myocardium as an example of the region of interest having a curved surface.

図11は、本実施形態に係るコンソール装置40で実行される画像表示処理の内容を説明するフローチャートを示す図である。この図11に示すように、本実施形態に係る画像表示処理は、概略的には、心筋領域を抽出する処理(ステップS40)、心筋形状を2次元画像化する処理(ステップS42)、心筋上の解析対象の測定値を算出する処理(ステップS44)、心筋に沿ったMAP画像を生成する処理(ステップS46)、及び、MAP画像を表示する処理(ステップS48)を、備えている。心筋領域を抽出する処理(ステップS40)、及び、心筋形状を2次元画像化する処理(ステップS42)は、上述した画像取得機能444と抽出機能445とにより実現される処理であり、心筋上の解析対象の測定値を算出する処理(ステップS44)は、上述した算出機能446により実現される処理であり、心筋に沿ったMAP画像を生成する処理(ステップS46)は、上述した生成機能447により実現される処理であり、MAP画像を表示する処理(ステップS48)は、上述した表示機能448により実現される処理である。   FIG. 11 is a flowchart illustrating the contents of image display processing executed by the console device 40 according to the present embodiment. As shown in FIG. 11, the image display processing according to the present embodiment generally includes processing for extracting a myocardial region (step S40), processing for forming a myocardial shape into a two-dimensional image (step S42), (Step S44), a process of generating a MAP image along the myocardium (step S46), and a process of displaying the MAP image (step S48). The process of extracting the myocardial region (step S40) and the process of converting the myocardial shape into a two-dimensional image (step S42) are processes realized by the image acquisition function 444 and the extraction function 445 described above. The process of calculating the measurement value to be analyzed (step S44) is a process realized by the calculation function 446 described above, and the process of generating a MAP image along the myocardium (step S46) is performed by the generation function 447 described above. The process that is realized and the process of displaying the MAP image (step S48) is a process that is realized by the display function 448 described above.

まず、心筋領域を抽出する処理(ステップS40)では、注目対象となる心筋壁領域を、既存の血管抽出アルゴリズムを用いて抽出する。抽出する心筋壁領域は、左心室の心壁など、解剖学的に特定可能な領域であってもよいし、事前に抽出した心臓領域からユーザがマウス操作などにより曲面を指定してもよいし、ユーザが注目したい特定の心血管を選択し、選択した心血管からの距離などから、選択した心血管の影響が強い領域であってもよい。抽出した心筋領域は、既存の位置合わせアルゴリズムにより、領域を示す座標を各時相間で座標に対応付けて記録する。   First, in the process of extracting the myocardial region (step S40), the myocardial wall region as the target of interest is extracted using an existing blood vessel extraction algorithm. The extracted myocardial wall region may be an anatomically identifiable region such as the left ventricular heart wall, or the user may specify a curved surface from the previously extracted heart region by operating the mouse or the like. The region may be a region in which the influence of the selected cardiovascular is strong based on the distance from the selected cardiovascular and the like by selecting a specific cardiovascular desired by the user. The extracted myocardial region is recorded by associating coordinates indicating the region with the coordinates between the time phases by an existing alignment algorithm.

次に、心筋形状を2次元画像化する処理(ステップS42)では、上述した心筋領域を抽出する処理(ステップS40)で求めた心筋領域を、2次元画像に変換する。曲面の2次元画像への変換は、既存のアルゴリズムを用いる。既存のアルゴリズムの例としては、曲面の法線方向の平均ベクトル方向へ射影する方法や、曲面状に格子状の座標系を配置し、設定した座標系へと曲面を射影する方法などがある。図12は、曲面の法線方向の平均ベクトル方向を2次元平面に射影することで、曲面を2次元画像に変換する方法の概要を示した図である。図13は、曲面に格子状の座標系を配置し、2次元平面に設定した座標系へと曲面を射影することで2次元画像に変換する方法の概要を示した図である。   Next, in the process of converting the myocardial shape into a two-dimensional image (step S42), the myocardial area obtained in the above-described process of extracting the myocardial region (step S40) is converted into a two-dimensional image. An existing algorithm is used to convert a curved surface into a two-dimensional image. Examples of existing algorithms include a method of projecting in the average vector direction of the normal direction of the curved surface, and a method of projecting the curved surface into a set coordinate system by arranging a lattice coordinate system on the curved surface. FIG. 12 is a diagram showing an outline of a method for converting a curved surface into a two-dimensional image by projecting the average vector direction of the normal direction of the curved surface onto a two-dimensional plane. FIG. 13 is a diagram showing an outline of a method of converting a two-dimensional image by arranging a grid-like coordinate system on a curved surface and projecting the curved surface onto a coordinate system set on a two-dimensional plane.

次に、心筋上の解析対象の測定値を算出する処理(ステップS44)では、心筋領域を抽出する処理(ステップS40)、及び、心筋形状を2次元画像化する処理(ステップS42)で抽出した心筋領域上の各点の座標について、解析対象となる測定値を各時相で算出する。解析対象となる測定値は、例えばCT値であってもよいし、心壁の移動距離など、心筋解析で通常用いられる計測値でもよい。また解析対象となる測定値の算出では、心筋上の各点の画素1つから算出してもよいし、心筋上の各点の周辺にある複数の画素を選択して平均値や最大値などを算出してもよい。   Next, in the process of calculating the measurement value of the analysis target on the myocardium (step S44), it is extracted by the process of extracting the myocardial region (step S40) and the process of converting the myocardial shape into a two-dimensional image (step S42). For the coordinates of each point on the myocardial region, the measurement value to be analyzed is calculated at each time phase. The measurement value to be analyzed may be, for example, a CT value, or a measurement value that is normally used in myocardial analysis, such as a moving distance of the heart wall. The measurement value to be analyzed may be calculated from one pixel at each point on the myocardium, or a plurality of pixels around each point on the myocardium may be selected to obtain an average value, maximum value, etc. May be calculated.

次に、心筋に沿ったMAP画像を生成する処理(ステップS46)では、心筋形状を2次元画像化する処理(ステップS42)で生成した複数時相の2次元画像(心筋領域を2次元画像として射影したもの)の中から、代表となる時相を選択する。本実施形態では、選択した時相の2次元画像における心筋領域の接線方向(つまり、2次元画像の面に平行な方向)にx軸、y軸を設定し、心筋領域の法線方向(つまり、2次元画像の面の法線方向)にz軸を割り当てた3次元画像を生成するものとする。   Next, in the process of generating a MAP image along the myocardium (step S46), a two-dimensional image having a plurality of phases generated by the process of generating a myocardial shape into a two-dimensional image (step S42) Select a representative time phase from the projected ones). In the present embodiment, the x-axis and y-axis are set in the tangential direction of the myocardial region (that is, the direction parallel to the plane of the two-dimensional image) in the selected two-dimensional image at the time phase, and the normal direction of the myocardial region (that is, It is assumed that a three-dimensional image in which the z axis is assigned to the normal direction of the surface of the two-dimensional image is generated.

生成する3次元画像の各画素の画素値は、心筋上の解析対象の測定値を算出する処理(ステップS44)で求めた解析対象の測定値を、心筋領域の1時相分の測定値分布をx軸、y軸に平行に割り当てて1時相分の画像を生成し、さらに、z軸と平行な方向に時系列に各時相の画像を並べて配置した画像を生成する。結果として、上記の処理で生成した画像は、時相毎に、心筋形状を2次元画像化する処理(ステップS42)で生成した2次元画像に測定値を割り当てて画像を生成し、この時相毎の画像を時系列に並べたような画像が生成できることになる。この処理の概要を示したのが、図14である。   The pixel value of each pixel of the three-dimensional image to be generated is the measurement value distribution for one time phase of the myocardial region obtained from the measurement value of the analysis target obtained in the process of calculating the measurement value of the analysis target on the myocardium (step S44). Are assigned in parallel to the x-axis and the y-axis to generate an image for one time phase, and further, an image in which images of each time phase are arranged in time series in a direction parallel to the z-axis is generated. As a result, the image generated by the above process generates an image by assigning a measurement value to the two-dimensional image generated in the process (step S42) of forming the myocardial shape into a two-dimensional image for each time phase. An image in which each image is arranged in time series can be generated. FIG. 14 shows an outline of this process.

次に、MAP画像を表示する処理(ステップS48)では、心筋に沿ったMAP画像を生成する処理(ステップS46)で生成したMAP画像に基づいて生成した表示画像を、表示装置42の画面に表示する。ここで、本実施形態で生成するMAP画像及び表示画像は、通常の3次元画像と同様の画像処理が適用可能である。このため、ウィンドウ幅(WW)/ウィンドウレベル(WL)の調整やカラーマップ表示、スムージングなどの画像フィルタ処理、指定した画素の画素値表示、指定した画素を基準点とした領域抽出処理などが可能である。通常のMPR画像のようにMAP画像の任意断面を画像に表示して空間的情報と時間的情報を合わせて観察してもよいし、MAP画像上の任意の点を指定し、該当する時相に対応する元画像(処理に用いた3次元画像)を画面に表示してもよい。   Next, in the process of displaying the MAP image (step S48), the display image generated based on the MAP image generated in the process of generating the MAP image along the myocardium (step S46) is displayed on the screen of the display device 42. To do. Here, the MAP image and the display image generated in the present embodiment can be applied with the same image processing as a normal three-dimensional image. Therefore, adjustment of window width (WW) / window level (WL), color map display, image filter processing such as smoothing, pixel value display of specified pixels, region extraction processing using specified pixels as reference points, etc. are possible It is. An arbitrary cross section of the MAP image may be displayed on the image like a normal MPR image, and the spatial information and the temporal information may be observed together, an arbitrary point on the MAP image may be designated, and the corresponding time phase The original image (three-dimensional image used for processing) corresponding to the above may be displayed on the screen.

以上のように、本実施形態に係るX線CT装置1によれば、表示装置42の画面に表示された表示画像に基づいて、時相の変化に伴う、心筋領域に対するCT値や移動量の変化を、空間的情報と時間的情報を合わせた状態で観察することができる。すなわち、観察対象の測定値(CT値や移動量など)の変化を、心臓の拡張期と収縮期といった離散的な解析ではなく、各時相での推移の仕方が表示画像で観察できるため、虚血領域などが存在する場合、観察対象の測定値の異常な動きを観察できることになる。   As described above, according to the X-ray CT apparatus 1 according to the present embodiment, based on the display image displayed on the screen of the display device 42, the CT value and the movement amount for the myocardial region associated with the change in time phase can be obtained. The change can be observed in a state where spatial information and temporal information are combined. In other words, the change in the measurement value (CT value, movement amount, etc.) of the observation target can be observed in the display image, not the discrete analysis of the cardiac diastole and systole, but the transition of each time phase. When an ischemic region or the like exists, an abnormal movement of the measurement value of the observation target can be observed.

なお、上述した例では、心筋に対して図11の画像表示処理を行ったものであるが、心筋以外にも解析対象となる面を設定できる組織であれば、本実施形態を適用することができる。例えば、複数の関節の接続面に対しても、上述した画像表示処理が適用できる。この場合、特定の関節に対して、隣接する関節との接続面を抽出し、隣接関節との距離を観察対象の測定値として用いることで、関節の動きに伴う、関節の接続面における隣接関節との距離を可視化できる。これにより、表示装置42の表示画像に基づいて、例えば、関節痛の原因となる箇所を推定できるようになる。   In the above-described example, the image display processing of FIG. 11 is performed on the myocardium. However, the present embodiment can be applied to any tissue that can set the analysis target surface other than the myocardium. it can. For example, the above-described image display processing can be applied to connection surfaces of a plurality of joints. In this case, for a specific joint, the connection surface with the adjacent joint is extracted, and the distance from the adjacent joint is used as the measurement value of the observation target, so that the adjacent joint on the connection surface of the joint accompanying the movement of the joint Can be visualized. Thereby, based on the display image of the display device 42, for example, a location causing joint pain can be estimated.

本実施形態においても、上述した第1実施形態や第2実施形態と同様に、関心領域の空間的情報を保持したまま、時相の変化に伴う観察対象の測定値の変化を確認できる。また、解析対象の測定値の変化を3次元画像として捉えたことにより、解析対象の測定値の変化を、ユーザは直感的に解析できるようになる。   Also in the present embodiment, as in the first embodiment and the second embodiment described above, it is possible to confirm a change in the measurement value of the observation target accompanying a change in time phase while maintaining the spatial information of the region of interest. In addition, by capturing the change in the measurement value of the analysis target as a three-dimensional image, the user can intuitively analyze the change in the measurement value of the analysis target.

〔第1乃至第3実施形態の変形例〕
上述した第1実施形態乃至第3実施形態においては、本実施形態が適用されるモダリティがX線CT装置1である場合を例として説明したが、本実施形態は、他のモダリティである医用画像診断装置に適用することもできる。他のモダリティに本実施形態を適用した場合、そのモダリティが備えるコンソールの処理回路が、上述した画像表示処理を実行することにより、上記各実施形態が実現される。
[Modifications of First to Third Embodiments]
In the above-described first to third embodiments, the case where the modality to which the present embodiment is applied is the X-ray CT apparatus 1 has been described as an example. However, the present embodiment is a medical image having another modality. It can also be applied to a diagnostic device. When the present embodiment is applied to other modalities, the above-described embodiments are realized by the console processing circuit included in the modality executing the image display process described above.

また、本実施形態は、ワークステーション等の医用画像処理装置に適用することもできる。例えば、図15に示すように、医用画像処理装置WSをモダリティMDに接続し、医用画像処理装置WSが、上述した画像表示処理を実行することにより、本発明を実現できる。すなわち、モダリティがスキャン等により医用画像を生成し、この生成された医用画像を画像処理装置WSが取得した後、上述した画像表示処理を実行することにより、画像処理装置WSは、所定の関心領域に沿って、測定値を空間的位置及び時系列毎に並べた画像を生成することが出来る。   The present embodiment can also be applied to a medical image processing apparatus such as a workstation. For example, as shown in FIG. 15, the medical image processing device WS is connected to the modality MD, and the medical image processing device WS executes the above-described image display processing, thereby realizing the present invention. That is, the modality generates a medical image by scanning or the like, and after the image processing apparatus WS acquires the generated medical image, the image processing apparatus WS executes the above-described image display process, so that the image processing apparatus WS has a predetermined region of interest. Then, it is possible to generate an image in which the measurement values are arranged for each spatial position and time series.

具体的には、図15に示すように、医用画像処理装置WSは、制御部100と、入力インターフェース102と、記憶回路141と、表示装置142と、処理回路144とを備えて構成されている。制御部100は、この医用画像処理装置WSの全体的な制御を行う回路であり、例えば、プロセッサにより構成されており、必要なプログラムを記憶回路141のプログラム記憶回路1413から読み出して実行することで、各種の制御を実現する。   Specifically, as illustrated in FIG. 15, the medical image processing apparatus WS includes a control unit 100, an input interface 102, a storage circuit 141, a display device 142, and a processing circuit 144. . The control unit 100 is a circuit that performs overall control of the medical image processing apparatus WS. The control unit 100 includes, for example, a processor, and reads a necessary program from the program storage circuit 1413 of the storage circuit 141 and executes it. Realize various controls.

入力インターフェース102は、医用画像処理装置WSを、モダリティMDに接続するためのインターフェースである。ここでは、モダリティMDは、スキャン等により複数時相の医用画像を生成可能な医用機器であり、例えば、上述したようなX線CT装置により構成されている。   The input interface 102 is an interface for connecting the medical image processing apparatus WS to the modality MD. Here, the modality MD is a medical device that can generate medical images having a plurality of time phases by scanning or the like, and is configured by, for example, an X-ray CT apparatus as described above.

記憶回路141は、図1を用いて説明した各実施形態における記憶回路41に相当する回路である。また、記憶回路141は、画像記憶回路1411と、画素値対応付け記憶回路1412と、プログラム記憶回路1413とを備えており、画像記憶回路1411は、上述した各実施形態における画像記憶回路411に相当する回路であり、画素値対応付け記憶回路1412は、上述した各実施形態における画素値対応付け記憶回路412に相当する回路であり、プログラム記憶回路1413は、上述した各実施形態におけるプログラム記憶回路413に相当する回路である。   The memory circuit 141 is a circuit corresponding to the memory circuit 41 in each embodiment described with reference to FIG. The storage circuit 141 includes an image storage circuit 1411, a pixel value association storage circuit 1412, and a program storage circuit 1413. The image storage circuit 1411 corresponds to the image storage circuit 411 in each of the above-described embodiments. The pixel value association storage circuit 1412 is a circuit corresponding to the pixel value association storage circuit 412 in each embodiment described above, and the program storage circuit 1413 is a program storage circuit 413 in each embodiment described above. It is a circuit corresponding to.

処理回路144は、例えば、プロセッサにより構成されており、必要なプログラムを記憶回路141のプログラム記憶回路1413から読み出して実行することで、上述した各種の機能を実現する。すなわち、処理回路144は、プログラム記憶回路1413に格納されている画像表示処理プログラムを読み込んで実行することにより、上述した各実施形態における画像表示処理を実現する。   The processing circuit 144 is configured by, for example, a processor, and implements the various functions described above by reading a necessary program from the program storage circuit 1413 of the storage circuit 141 and executing it. That is, the processing circuit 144 reads and executes the image display processing program stored in the program storage circuit 1413, thereby realizing the image display processing in each embodiment described above.

この画像表示処理を処理回路144が実行することにより、処理回路144は、画像取得機能1444と、抽出機能1445と、算出機能1446と、生成機能1447と、表示機能1448とを実現する。画像取得機能1444は、上述した各実施形態における画像取得機能444に相当する機能であり、抽出機能1445は、上述した各実施形態における抽出機能445に相当する機能であり、算出機能1446は、上述した各実施形態における算出機能446に相当する機能であり、生成機能1447は、上述した各実施形態における生成機能447に相当する機能であり、表示機能1448は、上述した各実施形態における表示機能448に相当する機能である。   When the processing circuit 144 executes this image display processing, the processing circuit 144 realizes an image acquisition function 1444, an extraction function 1445, a calculation function 1446, a generation function 1447, and a display function 1448. The image acquisition function 1444 is a function corresponding to the image acquisition function 444 in each embodiment described above, the extraction function 1445 is a function corresponding to the extraction function 445 in each embodiment described above, and the calculation function 1446 is described above. The generation function 1447 is a function corresponding to the generation function 447 in each of the above-described embodiments, and the display function 1448 is a display function 448 in each of the above-described embodiments. It is a function equivalent to.

図15に示す医用画像処理装置WSは、医用画像をモダリティMDから取得する以外は、上述したX線CT装置1と同様の処理を実行する。すなわち、処理回路144が、上述した画像表示処理を実行することにより、上述した第1実施形態乃至第3実施形態と同様の処理が実現され、関心領域に沿って測定値を空間的位置及び時系列毎に並べた画像を生成するとともに、表示装置142に表示させることができる。   The medical image processing apparatus WS shown in FIG. 15 performs the same processing as the X-ray CT apparatus 1 described above, except that a medical image is acquired from the modality MD. That is, the processing circuit 144 executes the above-described image display processing, thereby realizing the same processing as in the first to third embodiments described above, and the measurement value is measured along the region of interest and the spatial position and time. An image arranged for each series can be generated and displayed on the display device 142.

なお、上述した第1実施形態乃至第3実施形態の画像取得機能444は、特許請求の範囲における画像取得部の一例であり、第1実施形態乃至第3実施形態の抽出機能445は、特許請求の範囲における抽出部の一例であり、第1実施形態乃至第3実施形態の算出機能446は、特許請求の範囲における算出部の一例であり、第1実施形態乃至第3実施形態の生成機能447は、特許請求の範囲における生成部の一例であり、第1実施形態乃至第3実施形態の表示機能448は、特許請求の範囲における表示部の一例である。また、第1実施形態乃至第3実施形態の指定機能449は、特許請求の範囲における指定部の一例であり、第1実施形態乃至第3実施形態の測定値取得機能450は、特許請求の範囲における測定値取得部の一例であり、第1実施形態乃至第3実施形態のグラフ表示機能451は、特許請求の範囲におけるグラフ表示部の一例である。   The image acquisition function 444 of the first to third embodiments described above is an example of an image acquisition unit in the scope of claims, and the extraction function 445 of the first to third embodiments is claimed. The calculation function 446 of the first to third embodiments is an example of the calculation unit in the claims, and the generation function 447 of the first to third embodiments. Is an example of a generation unit in the claims, and the display function 448 of the first to third embodiments is an example of a display unit in the claims. The designation function 449 of the first to third embodiments is an example of a designation unit in the scope of claims, and the measurement value acquisition function 450 of the first to third embodiments is claimed. The graph display function 451 of the first to third embodiments is an example of a graph display unit in the claims.

また、上述した第1実施形態乃至第3実施形態の変形例の画像取得機能1444は、特許請求の範囲における画像取得部の一例であり、第1実施形態乃至第3実施形態の変形例の抽出機能1445は、特許請求の範囲における抽出部の一例であり、第1実施形態乃至第3実施形態の変形例の算出機能1446は、特許請求の範囲における算出部の一例であり、第1実施形態乃至第3実施形態の変形例の生成機能1447は、特許請求の範囲における生成部の一例であり、第1実施形態乃至第3実施形態の変形例の表示機能1448は、表示部の一例である。   Moreover, the image acquisition function 1444 of the modification of the first to third embodiments described above is an example of an image acquisition unit in the scope of claims, and the modification of the first to third embodiments is extracted. The function 1445 is an example of an extraction unit in the scope of claims, and the calculation function 1446 of a modification of the first to third embodiments is an example of a calculation section in the scope of claims, and is described in the first embodiment. The generation function 1447 of the modification of the third embodiment is an example of a generation unit in the claims, and the display function 1448 of the modification of the first to third embodiments is an example of a display unit. .

以上、いくつかの実施形態を説明したが、これらの実施形態は、例としてのみ提示したものであり、発明の範囲を限定することを意図したものではない。本明細書で説明した新規な装置および方法は、その他の様々な形態で実施することができる。また、本明細書で説明した装置および方法の形態に対し、発明の要旨を逸脱しない範囲内で、種々の省略、置換、変更を行うことができる。添付の特許請求の範囲およびこれに均等な範囲は、発明の範囲や要旨に含まれるこのような形態や変形例を含むように意図されている。   Although several embodiments have been described above, these embodiments are presented as examples only and are not intended to limit the scope of the invention. The novel apparatus and methods described herein can be implemented in a variety of other forms. In addition, various omissions, substitutions, and changes can be made to the forms of the apparatus and method described in the present specification without departing from the spirit of the invention. The appended claims and their equivalents are intended to include such forms and modifications as fall within the scope and spirit of the invention.

1…X線CT装置、10…架台装置、11…X線発生装置、12…X線検出器、13…回転体、14…X線高電圧装置、15…架台制御装置、データ収集回路、30…寝台装置、40…コンソール装置、41…記憶回路、42…表示装置、43…入力装置、44…処理回路、441…システム制御機能、442…前処理機能、443…再構成処理機能、444…画像取得機能、445…抽出機能、446…算出機能、447…生成機能、448…表示機能 DESCRIPTION OF SYMBOLS 1 ... X-ray CT apparatus, 10 ... Mount apparatus, 11 ... X-ray generator, 12 ... X-ray detector, 13 ... Rotating body, 14 ... X-ray high voltage apparatus, 15 ... Mount control apparatus, Data acquisition circuit, 30 ... Bed device, 40 ... Console device, 41 ... Storage circuit, 42 ... Display device, 43 ... Input device, 44 ... Processing circuit, 441 ... System control function, 442 ... Pre-processing function, 443 ... Reconfiguration processing function, 444 ... Image acquisition function, 445 ... extraction function, 446 ... calculation function, 447 ... generation function, 448 ... display function

Claims (14)

複数時相の医用画像を取得する、画像取得部と、
前記複数時相の医用画像における所定の関心領域の抽出を行う、抽出部と、
前記所定の関心領域に関する所定の測定値を算出する、算出部と、
前記所定の関心領域に沿って前記測定値を空間的位置及び時系列毎に並べた画像を生成する、生成部と、
を備える医用画像診断装置。
An image acquisition unit for acquiring medical images of a plurality of time phases;
An extraction unit that performs extraction of a predetermined region of interest in the medical images of the plurality of time phases;
A calculation unit for calculating a predetermined measurement value related to the predetermined region of interest;
A generation unit that generates an image in which the measurement values are arranged for each spatial position and time series along the predetermined region of interest;
A medical image diagnostic apparatus comprising:
前記所定の測定値は、CT値である、請求項1に記載の医用画像診断装置。   The medical image diagnostic apparatus according to claim 1, wherein the predetermined measurement value is a CT value. 前記所定の測定値は、前記所定の関心領域に含まれる組織の位置に基づく計測値である、請求項1に記載の医用画像診断装置。   The medical image diagnostic apparatus according to claim 1, wherein the predetermined measurement value is a measurement value based on a position of a tissue included in the predetermined region of interest. 前記抽出部は、前記関心領域として組織の芯線を抽出し、前記算出部は、芯線に関する所定の測定値を算出する、請求項1乃至請求項3のいずれかに記載の医用画像診断装置。   The medical image diagnostic apparatus according to claim 1, wherein the extraction unit extracts a core line of a tissue as the region of interest, and the calculation unit calculates a predetermined measurement value related to the core line. 前記抽出部は、芯線抽出アルゴリズムを用いて、前記医用画像を解析することにより、前記芯線を抽出する、請求項4に記載の医用画像診断装置。   The medical image diagnosis apparatus according to claim 4, wherein the extraction unit extracts the core line by analyzing the medical image using a core line extraction algorithm. 前記抽出部は、前記関心領域として組織の軸を抽出し、前記算出部は、前記軸に関する所定の測定値を算出する、請求項1乃至請求項3のいずれかに記載の医用画像診断装置。   The medical image diagnostic apparatus according to claim 1, wherein the extraction unit extracts a tissue axis as the region of interest, and the calculation unit calculates a predetermined measurement value related to the axis. 前記抽出部は、前記医用画像を用いたユーザの設定に基づいて、前記軸を抽出する、請求項6に記載の医用画像診断装置。   The medical image diagnosis apparatus according to claim 6, wherein the extraction unit extracts the axis based on a user setting using the medical image. 前記抽出部は、前記関心領域として組織の面に基づく2次元画像を抽出し、前記算出部は、前記2次元画像に関する所定の測定値を算出し、前記生成部は、前記2次元画像に前記測定値を割り当て、これを時系列に並べた前記画像を生成する、請求項1乃至請求項3のいずれかに記載の医用画像診断装置。   The extraction unit extracts a two-dimensional image based on a tissue surface as the region of interest, the calculation unit calculates a predetermined measurement value regarding the two-dimensional image, and the generation unit adds the two-dimensional image to the two-dimensional image. The medical image diagnostic apparatus according to claim 1, wherein measurement images are assigned and the images in which the measurement values are arranged in time series are generated. 前記抽出部は、前記組織の面を曲面として特定し、前記曲面を2次元化することにより、前記2次元画像を抽出する、請求項8に記載の医用画像診断装置。   The medical image diagnosis apparatus according to claim 8, wherein the extraction unit specifies the surface of the tissue as a curved surface, and extracts the two-dimensional image by making the curved surface two-dimensional. 前記生成部は、前記医用画像及び/又は前記医用画像に基づいて生成された画像を、前記画像に並べて表示する参考画像として生成する、請求項1乃至請求項9のいずれかに記載の医用画像診断装置。   The medical image according to any one of claims 1 to 9, wherein the generation unit generates the medical image and / or an image generated based on the medical image as a reference image to be displayed side by side on the image. Diagnostic device. 前記生成部が生成した前記画像における前記測定値に対応付けられた輝度値に基づく表示画像を生成し、この表示画像を画面に表示させる、表示部をさらに備える、請求項1乃至請求項10のいずれかに記載の医用画像診断装置。   The display unit according to claim 1, further comprising: a display unit configured to generate a display image based on a luminance value associated with the measurement value in the image generated by the generation unit and display the display image on a screen. The medical image diagnostic apparatus according to any one of the above. 前記生成部が生成した前記画像における前記測定値に対応付けられた色に基づく表示画像を生成し、この表示画像を画面に表示させる、表示部をさらに備える、請求項1乃至請求項10のいずれかに記載の医用画像診断装置。   The display unit according to any one of claims 1 to 10, further comprising a display unit that generates a display image based on a color associated with the measurement value in the image generated by the generation unit and displays the display image on a screen. A medical image diagnostic apparatus according to claim 1. 前記表示画像を用いて、ユーザが、前記表示画像上で時相と空間的位置を指定する指定部と、
前記指定部で指定された時相と空間的位置に基づいて、前記表示画像上の当該時相の当該空間的位置に対応する前記測定値を取得する、測定値取得部と、
前記測定値取得部で取得した測定値をグラフ化して前記表示画像に重畳して表示する、グラフ表示部と、
を備える、請求項11又は請求項12に記載の医用画像診断装置。
Using the display image, a user specifies a time phase and a spatial position on the display image;
A measurement value acquisition unit that acquires the measurement value corresponding to the spatial position of the time phase on the display image based on the time phase and spatial position specified by the specification unit;
A graph display unit that displays the measurement value acquired by the measurement value acquisition unit in a graph and superimposed on the display image;
The medical image diagnostic apparatus of Claim 11 or Claim 12 provided with these.
複数時相の医用画像を取得する、画像取得部と、
前記複数時相の医用画像における所定の関心領域の抽出を行う、抽出部と、
前記所定の関心領域に関する所定の測定値を算出する、算出部と、
前記所定の関心領域に沿って前記測定値を空間的位置及び時系列毎に並べた画像を生成する、生成部と、
を備える医用画像処理装置。
An image acquisition unit for acquiring medical images of a plurality of time phases;
An extraction unit that performs extraction of a predetermined region of interest in the medical images of the plurality of time phases;
A calculation unit for calculating a predetermined measurement value related to the predetermined region of interest;
A generation unit that generates an image in which the measurement values are arranged for each spatial position and time series along the predetermined region of interest;
A medical image processing apparatus comprising:
JP2017057926A 2017-03-23 2017-03-23 Medical image diagnostic equipment and medical image processing equipment Active JP6878074B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017057926A JP6878074B2 (en) 2017-03-23 2017-03-23 Medical image diagnostic equipment and medical image processing equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017057926A JP6878074B2 (en) 2017-03-23 2017-03-23 Medical image diagnostic equipment and medical image processing equipment

Publications (2)

Publication Number Publication Date
JP2018158022A true JP2018158022A (en) 2018-10-11
JP6878074B2 JP6878074B2 (en) 2021-05-26

Family

ID=63795891

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017057926A Active JP6878074B2 (en) 2017-03-23 2017-03-23 Medical image diagnostic equipment and medical image processing equipment

Country Status (1)

Country Link
JP (1) JP6878074B2 (en)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03184531A (en) * 1989-10-27 1991-08-12 Arch Dev Corp Method and device for measuring flow amount of blood
JP2004181041A (en) * 2002-12-05 2004-07-02 Ziosoft Inc Image analysis apparatus, function map preparing means for image analysis apparatus, method and program of image analysis
JP2004230086A (en) * 2003-01-31 2004-08-19 Toshiba Corp Image processor, image data processing method, and program
JP2009022455A (en) * 2007-07-18 2009-02-05 Toshiba Corp Medical imaging diagnostic device
CN102028494A (en) * 2011-01-24 2011-04-27 海纳医信(北京)软件科技有限责任公司 Method and system for processing cerebral perfusion image sequence
US20120072190A1 (en) * 2010-09-16 2012-03-22 Siemens Corporation Method and System for Non-Invasive Assessment of Coronary Artery Disease
JP2012187342A (en) * 2011-03-14 2012-10-04 Toshiba Corp Medical image diagnostic apparatus and medical image processing apparatus
US20130261445A1 (en) * 2012-04-02 2013-10-03 Dirk Ertel Representation of blood vessels and tissue in the heart
JP2014128651A (en) * 2012-11-30 2014-07-10 Toshiba Corp Medical-image diagnostic evaluation device
US20140376791A1 (en) * 2013-06-25 2014-12-25 Siemens Medical Solutions Usa, Inc. Reconstruction of time-varying data
JP2015112232A (en) * 2013-12-11 2015-06-22 株式会社東芝 Image analysis apparatus and x-ray diagnostic apparatus
JP2015167790A (en) * 2014-03-10 2015-09-28 株式会社東芝 Medical image processor
US20160317113A1 (en) * 2015-04-29 2016-11-03 Bayer Pharma Aktiengesellschaft Determining the velocity of a fluid using an imaging method

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH03184531A (en) * 1989-10-27 1991-08-12 Arch Dev Corp Method and device for measuring flow amount of blood
JP2004181041A (en) * 2002-12-05 2004-07-02 Ziosoft Inc Image analysis apparatus, function map preparing means for image analysis apparatus, method and program of image analysis
JP2004230086A (en) * 2003-01-31 2004-08-19 Toshiba Corp Image processor, image data processing method, and program
JP2009022455A (en) * 2007-07-18 2009-02-05 Toshiba Corp Medical imaging diagnostic device
US20120072190A1 (en) * 2010-09-16 2012-03-22 Siemens Corporation Method and System for Non-Invasive Assessment of Coronary Artery Disease
CN102028494A (en) * 2011-01-24 2011-04-27 海纳医信(北京)软件科技有限责任公司 Method and system for processing cerebral perfusion image sequence
JP2012187342A (en) * 2011-03-14 2012-10-04 Toshiba Corp Medical image diagnostic apparatus and medical image processing apparatus
US20130261445A1 (en) * 2012-04-02 2013-10-03 Dirk Ertel Representation of blood vessels and tissue in the heart
JP2014128651A (en) * 2012-11-30 2014-07-10 Toshiba Corp Medical-image diagnostic evaluation device
US20140376791A1 (en) * 2013-06-25 2014-12-25 Siemens Medical Solutions Usa, Inc. Reconstruction of time-varying data
JP2015112232A (en) * 2013-12-11 2015-06-22 株式会社東芝 Image analysis apparatus and x-ray diagnostic apparatus
JP2015167790A (en) * 2014-03-10 2015-09-28 株式会社東芝 Medical image processor
US20160317113A1 (en) * 2015-04-29 2016-11-03 Bayer Pharma Aktiengesellschaft Determining the velocity of a fluid using an imaging method

Also Published As

Publication number Publication date
JP6878074B2 (en) 2021-05-26

Similar Documents

Publication Publication Date Title
US11450001B2 (en) Medical image processing apparatus, medical image processing method, and x-ray CT apparatus
US8553955B2 (en) Image processing apparatus, X-ray computed tomography apparatus, and image processing method
JP2016026556A (en) X-ray ct equipment and image forming device
US10111626B2 (en) X-ray CT apparatus
JP6509131B2 (en) X-ray CT apparatus, image processing apparatus, and image reconstruction method
JP7341803B2 (en) X-ray CT equipment, medical image processing equipment, and X-ray CT systems
US10299750B2 (en) Medical image processing apparatus and X-ray CT apparatus
JP6985047B2 (en) X-ray CT device
US10706596B2 (en) X-ray computed tomography apparatus and image generation apparatus
JP6878074B2 (en) Medical image diagnostic equipment and medical image processing equipment
US11321829B2 (en) Medical image processing apparatus, method for processing medical image, and non-transitory computer medium storing computer program
JP6373558B2 (en) X-ray CT system
JP7236239B2 (en) X-ray CT system and CT scan execution program
JP6923414B2 (en) X-ray computed tomography equipment, medical imaging equipment, and programs
US20220383501A1 (en) Medical information display apparatus, medical image processing apparatus, medical information display method, and computer program product
US20230200759A1 (en) Medical image diagnosis apparatus and scanning-range setting method
JP7426310B2 (en) X-ray computerized tomography imaging device
JP2016202886A (en) Medical image processor, x-ray ct apparatus and image processing method
JP2023045878A (en) Medical image processing device, medical image processing method, and program
JP2020058794A (en) Medical image processing device, medical image processing method and medical image processing program
JP2024039641A (en) Medical information quantification device, medical information quantification method, medical information quantification program, and medical image diagnostic device
JP2021037148A (en) Medical diagnostic imaging apparatus, medical diagnostic imaging program, and imaging planning apparatus
CN110881993A (en) X-ray CT apparatus, medical image processing apparatus, and X-ray CT system
JP2020005760A (en) X-ray ct apparatus
JP2019050860A (en) Medical image diagnostic device and image processing device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20191216

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20201127

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210105

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210226

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210428

R150 Certificate of patent or registration of utility model

Ref document number: 6878074

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150