JP2009153677A - Kinetic image processing system - Google Patents
Kinetic image processing system Download PDFInfo
- Publication number
- JP2009153677A JP2009153677A JP2007334515A JP2007334515A JP2009153677A JP 2009153677 A JP2009153677 A JP 2009153677A JP 2007334515 A JP2007334515 A JP 2007334515A JP 2007334515 A JP2007334515 A JP 2007334515A JP 2009153677 A JP2009153677 A JP 2009153677A
- Authority
- JP
- Japan
- Prior art keywords
- image
- dynamic
- region
- determined
- anatomical structure
- 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.)
- Pending
Links
- 238000012545 processing Methods 0.000 title claims abstract description 64
- 230000008859 change Effects 0.000 claims abstract description 76
- 210000004072 lung Anatomy 0.000 claims abstract description 59
- 210000003484 anatomy Anatomy 0.000 claims abstract description 58
- 230000002159 abnormal effect Effects 0.000 claims abstract description 53
- 238000009423 ventilation Methods 0.000 claims abstract description 34
- 238000003384 imaging method Methods 0.000 claims description 80
- 230000002123 temporal effect Effects 0.000 claims description 5
- 230000005856 abnormality Effects 0.000 abstract description 17
- 238000003745 diagnosis Methods 0.000 abstract description 9
- 210000000481 breast Anatomy 0.000 abstract 1
- 238000000034 method Methods 0.000 description 29
- 238000010191 image analysis Methods 0.000 description 28
- 230000008569 process Effects 0.000 description 21
- 230000000241 respiratory effect Effects 0.000 description 21
- 238000004891 communication Methods 0.000 description 11
- 238000001514 detection method Methods 0.000 description 11
- 238000006243 chemical reaction Methods 0.000 description 8
- 230000029058 respiratory gaseous exchange Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 230000003902 lesion Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000012790 confirmation Methods 0.000 description 3
- 210000000621 bronchi Anatomy 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 230000004217 heart function Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000003434 inspiratory effect Effects 0.000 description 1
- 210000005240 left ventricle Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000005245 right atrium Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
本発明は、動態画像処理システムに関する。 The present invention relates to a dynamic image processing system.
従来、X線画像等を画像解析し、病変部と推定される画像領域を検出する装置が開発されている(例えば、特許文献1参照)。このような検出装置では、X線画像等において正常組織とは異なる濃度特性を有する病変部を検出し、検出された病変部の候補領域を指し示す表示を行って、医師がX線画像等を観察して異常か正常かの判断を行う際の参考情報を提供している。
しかしながら、肺の換気能については、呼吸のサイクルを通しての観察が必要であり、1枚の画像だけで正常又は異常の判断をすることはできない。また、換気能低下の原因となる局所的な病変部分だけでなく、その病変部分によって影響を受ける肺区域が診断に重要となる場合がある。よって、従来のように病変部のみを検出して表示する方法では、換気能について十分な診断支援ができるとはいえない。 However, lung ventilation requires observation through the respiratory cycle, and normal or abnormal judgment cannot be made with only one image. Moreover, not only the local lesion part that causes a decrease in ventilation capacity, but also the lung area affected by the lesion part may be important for diagnosis. Therefore, it cannot be said that the conventional method of detecting and displaying only the lesioned part can provide sufficient diagnosis support for the ventilation ability.
本発明の課題は、肺の換気能の診断に有効な情報を提供することである。 An object of the present invention is to provide information effective for diagnosis of lung ventilation ability.
請求項1に記載の発明によれば、
被写体の胸部を動態撮影し、複数の時間位相における動態画像を生成する撮影手段と、
前記複数の時間位相における動態画像を用いて、各動態画像に含まれる肺野領域を複数の領域に分割した領域毎に時間の経過に伴う面積変化率及び/又は信号変化率を算出し、当該算出した面積変化率及び/又は信号変化率に基づいて換気能の正常又は異常を領域毎に判定する第1判定手段と、
肺における解剖学的構造の位置及び名称が予め定められた参照画像を用いて、前記複数の時間位相における動態画像につき、解剖学的構造の位置及び名称を判定する第2判定手段と、
前記複数の動態画像のうち少なくとも1つを表示する表示手段と、
前記表示された動態画像において、前記異常と判定された領域につき、前記判定された解剖学的構造の位置及び名称の情報を表示させる制御手段と、
を備える動態画像処理システムが提供される。
According to the invention of
An imaging means for capturing a dynamic image of the subject's chest and generating dynamic images at a plurality of time phases;
Using the dynamic images in the plurality of temporal phases, calculating the area change rate and / or signal change rate over time for each region obtained by dividing the lung field region included in each dynamic image into a plurality of regions, First determination means for determining normal or abnormal ventilation capacity for each region based on the calculated area change rate and / or signal change rate;
Second determination means for determining the position and name of the anatomical structure for the dynamic images in the plurality of temporal phases using a reference image in which the position and name of the anatomical structure in the lung are predetermined;
Display means for displaying at least one of the plurality of dynamic images;
In the displayed dynamic image, control means for displaying information on the position and name of the determined anatomical structure for the area determined to be abnormal,
A dynamic image processing system is provided.
請求項2に記載の発明によれば、
複数の分類レベルで定められる解剖学的構造の位置及び名称に応じて、前記参照画像は複数の分類レベル毎に作成され、
前記第2判定手段は、前記複数の分類レベル毎に作成された参照画像を用いて、各分類レベルにおける解剖学的構造の位置及び名称を判定する請求項1に記載の動態画像処理システムが提供される。
According to invention of
Depending on the position and name of the anatomical structure defined at multiple classification levels, the reference image is created for each of multiple classification levels,
The dynamic image processing system according to
請求項3に記載の発明によれば、
前記参照画像は被写体の属性毎に作成され、
前記動態撮影の被写体の属性に関する情報を入力する入力手段を備え、
前記第2判定手段は、前記被写体の属性毎に作成された参照画像のうち、前記入力された情報に対応する属性の参照画像を用いる請求項1又は2に記載の動態画像処理システムが提供される。
According to invention of
The reference image is created for each attribute of the subject,
Input means for inputting information on the attributes of the subject of the dynamic shooting,
3. The dynamic image processing system according to
請求項4に記載の発明によれば、
前記制御手段は、異常と判定された領域を示す指標情報を表示させる請求項1〜3の何れか一項に記載の動態画像処理システムが提供される。
According to invention of Claim 4,
The dynamic image processing system according to any one of
請求項5に記載の発明によれば、
前記制御手段は、異常と判定された領域を正常と判定された領域とは異なる表示形態で表示させる請求項1〜4の何れか一項に記載の動態画像処理システムが提供される。
According to the invention of
The dynamic image processing system according to any one of
請求項1に記載の発明によれば、肺の換気能に異常がある可能性が高い領域の情報だけでなく、その異常によって影響を受ける可能性のある解剖学的構造の位置及び名称の情報を医師に提供することができる。従って、肺の換気能の診断に有効な情報を医師に提供することができる。 According to the first aspect of the present invention, not only information on a region where there is a high possibility of abnormality in lung ventilation, but also information on the position and name of an anatomical structure that may be affected by the abnormality. Can be provided to doctors. Therefore, it is possible to provide information useful for diagnosis of lung ventilation ability to a doctor.
請求項2に記載の発明によれば、各分類レベルでの解剖学的構造の位置又は名称の情報を医師に提供することができる。医師は各分類レベルの視点から診察が可能となる。 According to the second aspect of the present invention, information on the position or name of the anatomical structure at each classification level can be provided to the doctor. Doctors can see from the viewpoint of each classification level.
請求項3に記載の発明によれば、被写体の特徴に近い参照画像を用いて解剖学的構造の位置及び名称の判定を行うことができ、精度の高い判定が可能となる。 According to the third aspect of the present invention, the position and name of the anatomical structure can be determined using the reference image close to the feature of the subject, and the determination can be made with high accuracy.
請求項4に記載の発明によれば、医師は異常と判定された領域の位置を容易に把握することができる。 According to the fourth aspect of the present invention, the doctor can easily grasp the position of the area determined to be abnormal.
請求項5に記載の発明によれば、医師は異常と判定された領域を容易に識別することが可能となる。
According to the invention described in
まず、構成を説明する。
図1に、本実施形態における動態画像処理システム1を示す。
図1に示すように、動態画像処理システム1は、撮影装置10、撮影用コンソール20、診断用コンソール30と、画像処理装置40と、サーバ50とを備えて構成されている。各構成装置10〜50はネットワークNを介して接続されている。
First, the configuration will be described.
FIG. 1 shows a dynamic
As illustrated in FIG. 1, the dynamic
図2を参照して、撮影装置10、撮影用コンソール20及び診断用コンソール30についてさらに説明する。撮影装置10、撮影用コンソール20、診断用コンソール30は、被写体のX線画像を撮影する際に用いられるものである。
図2に示すように、撮影装置10はX線源11、検出器12、読取部13、サイクル検出部14を備えて構成されている。一方、撮影用コンソール20は、制御部21、記憶部22、操作部23、表示部24、通信部25を備えて構成されている。診断用コンソール30も同様に、制御部31、記憶部32、操作部33、表示部34、通信部35を備えて構成されている。
With reference to FIG. 2, the
As shown in FIG. 2, the
まず、撮影装置10について説明する。
撮影装置10は、被写体Wに対しX線を照射し、検出器12からX線画像を読み取る。撮影装置10では動態撮影が可能である。動態撮影とは撮影を連続的に行い、複数の時間位相における動態画像を得る撮影方法である。動態画像とは動態撮影により得られた撮影画像をいい、本実施形態では動態画像はX線画像である。
First, the
The
X線源11は、撮影用コンソール20の制御部21の制御に従って、X線を照射する。制御されるX線照射条件としては、例えば動態撮影における連続撮影時のパルスレート、パルス幅、パルス間隔、照射開始/終了タイミング、X線管電流、X線管電圧、フィルタ値等が挙げられる。パルスレートとは単位時間あたりの撮影回数をいい、パルス幅は撮影1回あたりのX線照射時間である。パルス間隔は、連続撮影においてあるX線照射開始から次の撮影でX線照射を開始するまでの時間である。
The
検出器12は、被写体Wを挟んでX線源11と対向する位置に配置される。検出器12はX線の検出センサがマトリクス状に配置されたFPD(Flat Panel Detector)等である。すなわち、X線をその強度に応じた電気信号に変換して画素(検出センサ)毎に蓄積するので、検出器12にはX線画像が記録されることとなる。
The
読取部13は、検出器12からX線画像を読み取る処理を行い、読み取られたX線画像を撮影用コンソール20に送信する。なお、読取動作は制御部21により制御される。制御される画像読取条件としては、フレームレート、フレーム間隔、画素サイズ、画像サイズ等がある。フレームレート、フレーム間隔は上記パルスレート、パルス間隔と同意義である。
The
サイクル検出部14は、被写体Wの撮影部位について生体反応のサイクルを検出する。例えば、撮影部位が肺野を含む胸部である場合には呼吸モニタベルト、CCDカメラ、光学カメラ、スパイロメータ等を適用して呼吸サイクルを検出する。また、撮影部位が心臓部である場合、心拍計や心電計等を用いて心拍サイクルを検出する。
サイクル検出部14は、検出したサイクルの情報を、撮影用コンソール20の制御部21に出力する。
The
The
次に、撮影用コンソール20及び診断用コンソール30について説明する。
撮影用コンソール20は技師の撮影操作用として用いられ、撮影条件等の入力を受け付けたり、撮影装置10のX線画像を技師の確認用に表示したりする。診断用コンソール30は、医師の操作用として用いられ、撮影用コンソール20から送信されるX線画像を医師の確認用に表示したりする。
Next, the
The
診断用コンソール30の各部(制御部31、記憶部32、操作部33、表示部34、通信部35)の機能は、撮影用コンソール20の各部(制御部21、記憶部22、操作部23、表示部24、通信部25)と基本的に同一である。よって、ここでは撮影用コンソール20の各部を代表として説明し、診断用コンソール30の各部の説明を省略する。
The functions of each part of the diagnostic console 30 (the
制御部21は、CPU(Central Processing Unit)、RAM(Random Access Memory)等により構成されている。制御部21は記憶部22に記憶されている各種プログラムをCPUにより読み出してRAM内に展開し、展開されたプログラムとの協働により各種演算を行ったり、各部の動作を集中制御したりして処理を実行する。
なお、制御部21はCPUクロックを利用して時間を計時するタイマ機能を有する。
The control unit 21 includes a CPU (Central Processing Unit), a RAM (Random Access Memory), and the like. The control unit 21 reads various programs stored in the
The control unit 21 has a timer function for measuring time using the CPU clock.
記憶部22は、ハードディスク等のメモリであり、制御部21により用いられる各種プログラムやプログラムの実行に必要なパラメータ等を記憶している。例えば、撮影部位毎に最適化された撮影条件(X線照射条件やX線画像の画像読取条件等)を記憶している。
The
操作部23は、キーボードやマウス等を備えて構成され、これらの操作に応じて操作信号を生成して制御部21に出力する。操作部23は、被写体Wの属性に関する情報を含む患者情報を入力するための入力手段の1つである。
The
表示部24はディスプレイを備え、制御部21の表示制御に従って各種操作画面や撮影により得られたX線画像等を表示する。
通信部25は、通信用のインターフェイスを備え、ネットワークNに接続された外部装置と通信を行う。
The
The
次に、画像処理装置40及びサーバ50について説明する。
画像処理装置40及びサーバ50は、撮影により得られたX線画像を提供するために用いられる。
図3を参照して、画像処理装置40について説明する。
画像処理装置40は、X線画像に対し、医師が観察しやすい画質となるように画像処理を施すものである。図3に示すように、画像処理装置40は、制御部41、操作部42、表示部43、記憶部44、通信部45、画像処理部46、画像解析部47を備えて構成されている。
Next, the
The
The
The
制御部41〜通信部45については、上記説明した撮影用コンソール20の制御部21〜通信部25と基本的な機能は同じであるので、ここでは詳細な説明を省略する。
Since the basic functions of the
画像処理部46は、X線画像に対し、階調変換処理、周波数調整処理等の各種画像処理を施す。画像処理は撮影部位に応じた種類のものを、撮影部位に応じた画像処理条件により施す。
The
画像解析部47は、CPUと判定処理に係るプログラムとの協働により、胸部を動態撮影して得られた複数の時間位相におけるX線画像を解析し、換気能の正常又は異常を判定する。また、画像解析部47は、X線画像について解剖学的構造の位置及び名称を判定する。具体的な判定方法については後述する。
The
サーバ50は、大容量メモリを備えてこのメモリに画像処理装置40によって画像処理されたX線画像を保存し、管理する。サーバ50に保存されたX線画像は診断用コンソール30からの要求に応じて配信され、診察に供される。
The
次に、動作について説明する。
本実施形態に係る動態画像処理システム1は、胸部の動態撮影を行い、得られた複数の時間位相における動態画像を用いて画像解析し、換気能の正常又は異常を判定するとともに、解剖学的構造の位置及び名称を判定し、その判定結果を表示する。
図4は、その際に主に機能する撮影装置10、診断用コンソール30、画像処理装置40における処理の流れを示すフローチャートである。
Next, the operation will be described.
The dynamic
FIG. 4 is a flowchart showing the flow of processing in the
図4に示すように、まず撮影装置10において動態撮影を行い、複数の時間位相における動態画像を生成する(ステップS1)。
撮影にあたっては、撮影技師が撮影用コンソール20の操作部23を介して、被写体Wに関する患者情報の入力、撮影部位の指定操作等を行う。患者情報には被写体Wつまり患者の氏名の他、年齢、性別、体重、身長等の患者の属性を示す情報が含まれる。なお、患者情報の入力手段として操作部23の例を挙げたが、診断用コンソール30やサーバ50等で管理されている患者情報を取得することとしてもよく、この場合は通信部25が入力手段として機能する。
As shown in FIG. 4, first, dynamic imaging is performed by the
At the time of imaging, the imaging engineer inputs patient information regarding the subject W, specifies an imaging region, and the like via the
撮影用コンソール20では、制御部21が指定された撮影部位に応じた撮影条件を記憶部22から読み出し、撮影装置10のX線源11におけるX線照射条件、読取部13における画像読取条件として設定する。以下、撮影技師により「肺(換気)」の撮影部位が指定されたとして説明する。換気能を見るために肺を撮影する場合、呼吸サイクルは平均0.3回/秒程度であることから、これを考慮して少なくとも1つの呼吸位相について動態撮影を行うため、例えば下記の撮影条件が設定される。
フレームレート(パルスレート):3フレーム/秒(つまり1秒あたり3回の撮影)
画素サイズ:400μm
画像サイズ:40cm×30cm
管電圧:120kV
管電流:50mA
撮影タイミング:吸気から呼気への変換点のタイミング(撮影開始タイミング)からフレーム間隔時間毎
In the
Frame rate (pulse rate): 3 frames / second (that is, 3 shots per second)
Pixel size: 400μm
Image size: 40cm x 30cm
Tube voltage: 120 kV
Tube current: 50 mA
Shooting timing: Every frame interval from the timing of the conversion point from inspiration to expiration (shooting start timing)
なお、制御部21はサイクル検出部14により検出された呼吸サイクルの情報に基づいて、フレームレート等の条件を修正する。例えば、検出された呼吸サイクルに基づいて、1サイクルが所定のフレーム数(例えば、10フレーム)で撮影されるように制御部21がフレームレートを算出し、設定し直す。上記のフレームレートの条件例でいえば、サイクル検出部14により検出された呼吸サイクル数が0.25回/秒であった場合、フレームレートは2.5フレーム/秒に修正される。
The control unit 21 corrects conditions such as a frame rate based on the information on the respiratory cycle detected by the
撮影条件を設定後、制御部21はサイクル検出部14により検出される呼吸サイクルの情報に基づいて、撮影開始タイミングかどうか、つまり1つの呼吸サイクルの動態が始まるタイミング(例えば、吸気→呼気の変換点)かどうかを判断する。撮影開始タイミングであれば、制御部21はX線源11及び読取部13を制御して動態撮影を開始させる。また、制御部21は動態撮影開始に合わせて撮影の開始から終了までに要した撮影時間を経時する。
After setting the imaging conditions, the control unit 21 determines whether or not the imaging start timing, that is, the timing at which the dynamics of one respiratory cycle starts (for example, conversion from inspiration to expiration) based on the respiratory cycle information detected by the
撮影装置10では、設定されたX線照射条件に従ってX線源11から所定のパルスレートでX線を照射する。同様に読取部13は設定された画像読取条件に従って、所定フレームレートで検出器12からX線画像の読取処理を行う。このX線照射動作と画像読取動作は制御部21が同期させる。これにより、複数の時間位相における動態画像が生成され、撮影用コンソール20に出力される。
In the
撮影用コンソール20では、制御部21の表示制御により動態撮影により得られた各時間位相の動態画像を表示部24に表示する。撮影技師が画質等を確認するためである。撮影技師により操作部23を介して承認操作がなされると、制御部21は各時間位相の動態画像に一連の撮影を識別するためのIDや、患者情報、撮影時間の情報を付帯して診断用コンソール30に送信する。診断用コンソール30でも同様に確認用の表示を行う(ステップS2)。そして、承認操作がなされると、各時間位相の動態画像を画像処理装置40に送信する。
In the
画像処理装置40では、各時間位相の動態画像に対し、画像処理部46により肺(換気)の撮影部位に応じた画像処理を施した後、画像解析部47により換気能の正常又は異常の判定を行う(ステップS3)。換気能の正常又は異常の判定は、肺野領域の面積変化率及び/又は信号変化率を算出することにより行う。
In the
正常又は異常の判定に係る処理については、図5を参照して説明する。
図5に示すように、画像解析部47は各時間位相の動態画像のうち任意の動態画像を基準画像として設定する(ステップS31)。次いで、画像解析部47はこの基準画像から肺野領域を検出し、検出した肺野領域に外接する略矩形の領域において図6の点線で示すように複数の小領域に分割する(ステップS32)。分割する小領域は、例えば0.4〜4cm角とすればよい。
Processing related to normality or abnormality determination will be described with reference to FIG.
As shown in FIG. 5, the
なお、肺野領域の検出方法は何れの方法であってもよい。例えば、基準画像の信号値のヒストグラムから判別分析によって閾値を求め、この閾値より高信号の領域を肺野領域として1次検出する。次いで、1次検出された領域の境界付近でエッジ検出を行い、境界付近の小領域でエッジが最大となる点を境界に沿って抽出すれば肺野領域の境界を検出することができる。 Note that the lung field region detection method may be any method. For example, a threshold value is obtained by discriminant analysis from a histogram of signal values of the reference image, and a region having a signal higher than the threshold value is primarily detected as a lung field region. Next, edge detection is performed in the vicinity of the boundary of the first detected region, and a point in the edge in the small region near the boundary is extracted along the boundary, so that the boundary of the lung field region can be detected.
次いで、ローカルマッチングにより、他の時間位相の動態画像において基準画像の各分割領域と対応する領域を求める(ステップS33)。具体的には、基準画像と時間位相が隣接する(つまり時間位相が最も近い)動態画像を、例えば基準画像の2倍の大きさの領域に分割する。そして、この2倍の大きさをもつ分割領域の中で基準画像の対応する分割領域を移動させて、移動毎にマッチング度を算出し、マッチング度が最大となる位置を求める。マッチング度とは画像の一致性を示す度合いをいい、最小二乗法や相互相関によって求めることができる。この処理を時間位相が隣接する各動態画像間で繰り返し行うことにより、基準画像の各分割領域が他の時間位相の動態画像のどの領域に対応するのかを求めることができる。 Next, regions corresponding to the divided regions of the reference image in the dynamic images of other time phases are obtained by local matching (step S33). Specifically, a dynamic image in which the reference image and the time phase are adjacent (that is, the closest time phase) is divided into, for example, an area twice as large as the reference image. Then, the corresponding divided area of the reference image is moved in the divided area having the double size, the matching degree is calculated for each movement, and the position where the matching degree is maximized is obtained. The degree of matching refers to the degree of image matching, and can be obtained by the least square method or cross-correlation. By repeating this process between the dynamic images adjacent in time phase, it is possible to determine which region of the dynamic image in another time phase corresponds to each divided region of the reference image.
次いで、画像解析部47は対応する領域毎に面積変化率及び/又は信号変化率を算出する(ステップS34)。まず、画像解析部47は基準画像の分割領域の中心点を頂点とする領域を生成する。図6の実線で示す領域が新たに生成した領域である。そして、新たに生成した領域の面積について時間の経過に伴う面積変化率を求める。面積は領域の画素数で求め、面積変化率は時間位相が隣接する動態画像間で求める。つまり、時間位相が隣接する2つの動態画像において基準画像を元にそれぞれ対応するとされた4つの領域の中心点により定まる領域の面積の割合を面積変化率として、隣接する動態画像間で順次算出する。
Next, the
一方、新たに生成した領域毎に、その領域内の各画素の信号値をその画素数で除算し、平均信号値を求める。そして、この平均信号値の時間の経過に伴う変化率を信号変化率として求める。つまり、時間位相が隣接する2つの動態画像において基準画像を元にそれぞれ対応するとされた4つの領域の中心点により定まる領域の平均信号値の割合を信号変化率として、隣接する動態画像間で順次算出する。なお、平均信号値ではなく、領域内の各画素の信号値の総和を算出し、この信号値の総和の時間の経過に伴う変化率を信号変化率として採用してもよい。 On the other hand, for each newly generated area, the signal value of each pixel in the area is divided by the number of pixels to obtain an average signal value. Then, the rate of change of the average signal value with the passage of time is obtained as the signal rate of change. That is, the ratio of the average signal value in the region determined by the center point of the four regions that correspond to each other based on the reference image in the two dynamic images whose temporal phases are adjacent to each other is sequentially used as the signal change rate between the adjacent dynamic images. calculate. Instead of the average signal value, the sum of the signal values of each pixel in the region may be calculated, and the rate of change of the sum of the signal values over time may be adopted as the signal rate of change.
面積変化率及び/又は信号変化率を算出すると、画像解析部47は対応する領域毎に求められた面積変化率及び/又は信号変化率を元に、領域毎に換気能が正常か異常かを判定する(ステップS35)。ここでは、多面的に換気能の異常性を判定するため、時間による変化と、各分割領域の周辺における変化という観点から判定を行う。なお、以下の判定は面積変化率又は信号変化率の一方のみを用いて行ってもよいし、双方を用いて行うことも可能である。
After calculating the area change rate and / or signal change rate, the
時間による変化の観点から判定する場合、肺野を左右と上下に分け、各時間位相の動態画像において算出された面積変化率及び/又は信号変化率から時間の経過によって肺野の左右と上下でどのように面積や信号値が変化しているのかを見る。肺野の上下方向とは図6の動態画像でいう副走査方向であり、左右方向とは図6の動態画像でいう主走査方向である。 When judging from the viewpoint of changes due to time, the lung field is divided into left and right and up and down, and the area change rate and / or signal change rate calculated in the dynamic image of each time phase is divided into the left and right and top and bottom of the lung field over time. See how the area and signal values change. The vertical direction of the lung field is the sub-scanning direction referred to in the dynamic image in FIG. 6, and the horizontal direction is the main scanning direction referred to in the dynamic image in FIG.
左右での変化を見る場合は、肺野領域において副走査方向の位置が略同一である領域の面積変化率及び/又は信号変化率の平均値を左右それぞれの肺野領域で求める。肺野は呼吸によって上下方向に伸張と収縮を繰り返し、その伸張と収縮の程度は左右の肺野で略同一であることから、換気能が正常であれば面積や信号値の変化も左右の肺野で略同一であるのが一般的である。よって、左右の肺野領域で求められた平均値の差が閾値以上であれば、左右で面積変化率及び/又は信号変化率が大きく異なり換気能の異常が考えられるため、その平均値を求めた領域を異常と判定する。平均値の差が閾値を下回る場合には正常と判定する。この判定を肺野領域の副走査方向において順次行う。なお、閾値は経験的、実験的に異常と判定すべき値を予め求めて設定すればよい。 When viewing the change in the left and right, the average value of the area change rate and / or the signal change rate of the region where the position in the sub-scanning direction is substantially the same in the lung field region is obtained in each of the left and right lung field regions. The lung field repeatedly expands and contracts in the vertical direction by breathing, and the extent of the expansion and contraction is almost the same in the left and right lung fields. Generally, it is almost the same in the field. Therefore, if the difference between the average values obtained in the left and right lung field regions is greater than or equal to the threshold value, the area change rate and / or signal change rate is significantly different between the left and right, and abnormal ventilation may be considered. Is determined to be abnormal. When the difference between the average values is below the threshold value, it is determined as normal. This determination is sequentially performed in the sub-scanning direction of the lung field region. The threshold value may be set by obtaining in advance a value that should be determined empirically or experimentally as abnormal.
上下での変化を見る場合は、肺野領域において副走査方向の位置が略同一である領域の面積変化率及び/又は信号変化率の平均値を求める。これを肺野領域の副走査方向において順次求める。そして、各副走査方向において算出した平均値の関係が、下記条件1、2を満たさない場合、異常と判定する。逆に条件1、2を満たす場合には正常と判定する。
条件1:上の肺野(肺先)であるほど小さい
条件2:上の肺野(肺先)であるほど、平均値の変化の開始タイミングが遅い
通常、肺野は吸気によって肺野下部が伸張し、呼気によって肺野下部が収縮して元の大きさに戻る。つまり、換気能が正常であれば肺野上部の伸張や収縮の程度は肺野下部に比べて小さいとともに、伸張や収縮を始めるタイミングも肺野下部に比べて遅い。この上下方向における特徴を示す上記条件1、2を満たさない場合には、換気能の異常が考えられる。
In the case of looking at the vertical change, the average value of the area change rate and / or the signal change rate of the region in the lung field region where the position in the sub-scanning direction is substantially the same is obtained. This is sequentially obtained in the sub-scanning direction of the lung field region. If the relationship between the average values calculated in each sub-scanning direction does not satisfy the following
Condition 1: The smaller the upper lung field (lung tip), the smaller. Condition 2: The higher the lung field (lung tip), the slower the start of the mean value change. It expands, and the exhalation causes the lower lung field to contract and return to its original size. In other words, if the ventilation ability is normal, the upper lung field is less stretched and contracted than the lower lung field, and the timing to start stretching and contracting is slower than the lower lung field. If the
一方、各分割領域の周辺での変化という観点から判定する場合、肺野領域において副走査方向が略同一である領域間で面積変化率及び/又は信号変化率の差を求め、この差を所定の閾値と比較する。通常、肺野は吸気によって肺野下部が伸張し、呼気によって肺野下部が収縮して元の大きさに戻る。つまり肺野の上下方向に伸張と収縮を繰り返すのであり、肺野の左右方向では上下方向ほど大きな変化が見られない。よって、求めた面積変化率及び/又は信号変化率の差が閾値より大きい場合には、主走査方向での面積変化率及び/又は信号変化率が異常に大きいと推定し、その領域を異常と判定する。閾値を下回る場合には当該領域を正常と判定する。 On the other hand, when judging from the viewpoint of the change in the vicinity of each divided region, the difference in area change rate and / or signal change rate is obtained between regions in which the sub-scanning direction is substantially the same in the lung field region, and this difference is determined in advance. Compare with the threshold value. Usually, the lower lung field expands by inspiration, and the lower lung field contracts by exhalation and returns to the original size. That is, expansion and contraction are repeated in the vertical direction of the lung field, and the change in the horizontal direction of the lung field is not as great as in the vertical direction. Therefore, when the difference between the obtained area change rate and / or signal change rate is larger than the threshold value, it is estimated that the area change rate and / or signal change rate in the main scanning direction is abnormally large, and the region is regarded as abnormal. judge. If it falls below the threshold, the area is determined to be normal.
なお、基準画像に対する面積変化率及び/又は信号変化率を求めて、上記と同様に換気能の正常又は異常をさらに判定することとしてもよい。基準画像に対する面積変化率、信号変化率の算出方法は基準画像と他の時間位相の動態画像との間で求めるという点で対象とする画像が異なるが、上述した方法と同様である。また、判定についても求めた基準画像に対する面積変化率及び/又は信号変化率を用いて、上述した方法と同様に行えばよい。これにより、隣接する時間位相での面積や信号値の変化という面に加えて、基準となる時間位相からの変化という面からも正常又は異常の判定を行うことができ、多面的な判定が可能となる。 In addition, it is good also as calculating | requiring the area change rate and / or signal change rate with respect to a reference | standard image, and further determining the normality or abnormality of ventilation capacity similarly to the above. The calculation method of the area change rate and the signal change rate with respect to the reference image is the same as the above-described method, although the target image is different in that it is obtained between the reference image and the dynamic image of another time phase. The determination may be performed in the same manner as described above using the area change rate and / or signal change rate with respect to the obtained reference image. As a result, in addition to the area and signal value change at the adjacent time phase, normal or abnormal judgment can be performed from the aspect of change from the reference time phase, and multi-faceted judgment is possible. It becomes.
次いで、画像解析部47は換気能の正常又は異常の判定結果の情報と、異常と判定された場合にはその画像部分の面積変化率及び/又は信号変化率の情報とを各時間位相の動態画像に付帯させる(ステップS36)。
以上のように、換気能について正常又は異常の判定を終えると、図4に示すように各時間位相の動態画像について解剖学的構造の位置及び名称を判定する処理を画像解析部47において行う(ステップS4)。
Next, the
As described above, when the determination of normality or abnormality of the ventilation ability is completed, the
図7を参照して、解剖学的構造の位置及び名称を判定する処理を説明する。
図7に示すように、まず画像解析部47は各時間位相における動態画像を解析し、呼吸位相を決定する(ステップS41)。
With reference to FIG. 7, the process for determining the position and name of the anatomical structure will be described.
As shown in FIG. 7, first, the
図8は、1つの呼吸サイクルにおいて動態撮影された複数の時間位相T(T=t0〜t6)の動態画像を示す図である。図8に示すように呼気時には肺の空気が排気されるため、肺野が収縮している。吸気時には肺に空気が流入するため、肺野が伸張している。よって、画像解析部47は各時間位相Tに係る動態画像の肺野領域の面積(画素数)を算出し、この面積が最大となる時間位相t0から最小となる時間位相t3までの動態画像を呼気位相、面積が最小の時間位相t3から最大の時間位相t6までの動態画像を吸気位相のものとして決定する。なお、横隔膜の高さから呼吸位相を特定することとしてもよい。横隔膜の高さが最大のところが最大呼気位であり、高さが最小のところが最大吸気位である。
FIG. 8 is a diagram showing dynamic images of a plurality of time phases T (T = t 0 to t 6 ) that are dynamic images in one respiratory cycle. As shown in FIG. 8, since the lung air is exhausted during expiration, the lung field is contracted. When inhaling, air flows into the lungs, so the lung field is stretched. Therefore, the
次いで、画像解析部47は、動態画像に付帯されている患者情報に基づいて、複数の参照画像のうち、後述のマッチング処理に用いる参照画像を決定する(ステップS42)。参照画像とは予め肺野について解剖学的構造の名称及び位置が定められた画像をいう。
Next, the
解剖学的構造の位置及び名称は複数の分類レベルで定められており、参照画像はその分類レベル毎に作成されて保存されている。一例として、3つの分類レベル1〜3に応じて作成された参照画像を図9A、図9B、図9Cに示す。図9A、図9B、図9Cに示すように、参照画像においては肺野領域を複数の解剖学的構造に分類しており、各解剖学構造の位置及び名称が定められている。図9A、図9B、図9Cにおいて、各解剖学的構造の画像領域に示されているのはその解剖学的構造の名称である。
The position and name of the anatomical structure are determined at a plurality of classification levels, and a reference image is created and stored for each classification level. As an example, reference images created according to three
分類レベル1は3つの分類レベルの中で最も大きい分類で解剖学的構造の名称及び位置が定められたものであり、分類レベル2は分類レベル1の解剖学的構造をさらに細かく分類したものである。分類レベル3は肺の気管支単位で解剖学的構造を定めたものであり、分類レベル1、2よりもさらに細かい分類となっている。すなわち、図9A、図9B、図9Cに示すように、同じ位置にある領域(斜線で示す領域)であっても、分類レベルが異なればその解剖学的構造の名称は分類レベル1の参照画像では右中葉領域であり、分類レベル2の参照画像では右第4領域、分類レベル3の参照画像では右第○気管支領域と異なることになる。
上記参照画像は1つの呼吸サイクルの時間位相に応じて複数作成されている。つまり、各時間位相に係る複数の参照画像を1セットとして、記憶部44に保存している。
さらに、参照画像は患者の属性に応じて作成されている。患者の属性とは患者の性別、年齢、体重、身長等の身体的な属性をいう。例えば、20代〜60代までの年代別、男女別に参照画像を作成し、記憶部44に保存している。
A plurality of the reference images are created according to the time phase of one respiratory cycle. That is, a plurality of reference images relating to each time phase are stored in the
Furthermore, the reference image is created according to the patient's attribute. Patient attributes refer to physical attributes such as patient gender, age, weight, and height. For example, reference images are created for each age group from the 20s to 60s and for each gender, and stored in the
よって、ステップS42では動態画像に付帯されている患者情報に基づいて被写体Wとなった患者の属性を判別し、記憶部44に記憶されている参照画像のうち、判別した属性に対応する参照画像を決定し、読み出す。
Therefore, in step S42, the attribute of the patient who has become the subject W is determined based on the patient information attached to the dynamic image, and the reference image corresponding to the determined attribute among the reference images stored in the
次に、参照画像は上述のように複数の時間位相に応じて作成されているので、画像解析部47はステップS41で決定した呼吸位相に基づいて、動態画像の呼吸位相と参照画像の呼吸位相とが対応するように、各時間位相の動態画像と各時間位相の参照画像との対応付けを行う(ステップS43)。対応付けが終了すると、画像解析部47は対応付けた動態画像と参照画像とを略一致させるマッチング処理を行う(ステップS44)。
Next, since the reference image is created according to a plurality of time phases as described above, the
図10を参照して、マッチング処理を説明する。
図10に示すように、画像解析部47は、各時間位相の動態画像について信号値の補正を行う(ステップS441)。X線画像はX線照射条件等の違いによって信号値にばらつきが生じる場合があり、整合をとる必要があるためである。具体的には、動態画像の信号値のヒストグラムから最大信号値及び最小信号値を求め、これを予め定めた基準信号値に一致させるように信号値の変換を行う。
同様に、画像解析部47は参照画像についても信号値の補正を行い(ステップS442)、X線画像及び参照画像を同じ基準信号値に合わせる。
The matching process will be described with reference to FIG.
As shown in FIG. 10, the
Similarly, the
補正を終えると、画像解析部47は呼吸位相が対応する動態画像と参照画像に対し、回転、平行移動を繰り返す処理(アフィン変換等)を施して初期位置合わせを行う(ステップS445)。これは、X線画像と参照画像に含まれる被写体Wの画像の位置がずれていることが多いので、後段の処理を行う前になるべく位置が一致するように合わせておくためである。初期位置合わせは、対応する動態画像と参照画像の組合せ全てについて行う。
When the correction is completed, the
次いで、画像解析部47は初期位置合わせをした動態画像を、複数の小領域に分割する(ステップS443)。そして、画像解析部47は、動態画像の分割領域毎に参照画像の画像領域との対応付けを行う(ステップS445)。対応付けの方法は、図5のステップS33のローカルマッチングの手法と同様である。つまり、参照画像をX線画像の分割領域の約2倍の大きさの領域毎に分割する。この参照画像の分割領域の中で、動態画像の分割領域を移動させ、最もマッチング度が高い領域をその動態画像の分割領域と対応付ける。
Next, the
次いで、対応付けた動態画像の分割領域と参照画像の領域とが略一致するように、動態画像に非線形ワーピング処理を施す(ステップS446)。非線形ワーピング処理では、動態画像の各分割領域の中心位置から、対応付けられた参照画像の領域の中心位置へのシフト値を求め、各分割領域についてのシフト値を用いて2次元10次多項式による近似処理することにより、動態画像の全画素について参照画像の全画素に対するシフト値を求める。そして、求めたシフト値分だけ動態画像の全画素をシフトさせる。 Next, nonlinear warping processing is performed on the dynamic image so that the divided region of the associated dynamic image and the region of the reference image substantially match (step S446). In the nonlinear warping process, a shift value from the center position of each divided area of the dynamic image to the center position of the associated reference image area is obtained, and a two-dimensional 10th order polynomial is used using the shift value for each divided area. By performing approximation processing, shift values for all pixels of the reference image are obtained for all pixels of the dynamic image. Then, all the pixels of the dynamic image are shifted by the calculated shift value.
図11を参照して、非線形ワーピング処理について説明する。図11に示すのは呼吸位相が対応する動態画像g1と参照画像r1である。
動態画像g1を複数領域に分割したものが動態画像g2である。この動態画像g2の分割領域と参照画像r2とを対応付けたのが参照画像r2である。動態画像g2は各分割領域と参照画像r2において対応付けられた領域とが略一致するように非線形ワーピング処理される。非線形ワーピング処理後の画像が動態画像g3である。図11から分かるように、動態画像g3と参照画像r2における被写体Wの画像部分が略一致している。
The non-linear warping process will be described with reference to FIG. FIG. 11 shows a dynamic image g1 and a reference image r1 corresponding to the respiratory phase.
A dynamic image g2 is obtained by dividing the dynamic image g1 into a plurality of regions. The reference image r2 associates the divided area of the dynamic image g2 with the reference image r2. The dynamic image g2 is subjected to nonlinear warping processing so that each divided region and the region associated with the reference image r2 substantially coincide with each other. The image after the nonlinear warping process is the dynamic image g3. As can be seen from FIG. 11, the image portion of the subject W in the dynamic image g3 and the reference image r2 substantially coincide.
以上のマッチング処理を、複数の分類レベルについて行う。すなわち、各分類レベルの参照画像を用いてそれぞれマッチング処理を行う。そして、マッチング処理を終えると、図7のステップS45の処理に移行する。
ステップS45では、分類レベル毎に、マッチング処理の結果、それぞれ対応する参照画像に略一致させた動態画像において解剖学的構造の位置及び名称を判別する(ステップS45)。つまり、参照画像とのマッチング処理した場合、参照画像において名称及び位置が定められている解剖学的構造と一致する画像領域についてはその名称の解剖学的構造であると判別する。例えば、分類レベル1の参照画像とのマッチングの結果、図12に示すように動態画像g1の肺野領域において分類レベル1に係る解剖学的構造の名称及び位置を判別できる。
The above matching process is performed for a plurality of classification levels. That is, matching processing is performed using the reference images of the respective classification levels. When the matching process is completed, the process proceeds to step S45 in FIG.
In step S45, the position and name of the anatomical structure are determined in the dynamic image substantially matched with the corresponding reference image as a result of the matching process for each classification level (step S45). That is, when matching processing with a reference image is performed, an image region that matches an anatomical structure whose name and position are determined in the reference image is determined to be an anatomical structure of that name. For example, as a result of matching with the reference image of the
次いで、画像解析部47は解剖学的構造の位置及び名称の判別結果と、正常又は異常の判定情報に基づいて、各時間位相の動態画像において異常と判定された画像領域が属する解剖学的構造の位置及び名称を判別する。この場合も各分類レベルで判別する(ステップS46)。例えば、図9A、図9B、図9Cに示す斜線部分が異常と判定されている場合、この異常と判定された部分が属する解剖学的構造の名称は、分類レベル1では右中葉領域、分類レベル2では右第4領域、分類レベル3では右第○気管支領域と判別する。
Next, the
画像解析部47は、判定結果として、各時間位相の動態画像において判別された解剖学的構造の位置及び名称、さらに異常と判定された画像領域がある場合にはこの画像領域が属する解剖学的構造の位置及び名称の情報を、各時間位相の動態画像に付帯させる(ステップS47)。
その後、判定結果の情報が付帯された各時間位相の動態画像は通信部44を介してサーバ50に送信される。
As a determination result, the
Thereafter, the dynamic image of each time phase accompanied with the information of the determination result is transmitted to the
サーバ50では、付帯情報とともに各時間位相の動態画像をデータベース化してメモリに保存する。診断用コンソール30から要求があれば、サーバ50は要求に係る患者の動態画像群を送信する。
The
図4に示すように、診断用コンソール30では、制御部31の表示制御により、サーバ50から取得した各時間位相の動態画像を表示部34に表示する(ステップS5)。このとき、制御部31は時間位相に応じて各動態画像を連続的に切り替えて動画として表示させる。医師は動的な変化を把握することが可能である。
As shown in FIG. 4, the
次いで、制御部31は表示した動態画像の付帯情報に基づいて、換気能の判定結果及び解剖学的構造の判定結果を表示させる(ステップS6)。具体的には、換気能の判定結果として、異常と判定された画像領域を指し示す指標情報を表示させるとともに、正常組織と識別できるように表示形態を異ならせて表示させる。また、異常と判定された画像領域については面積変化率及び/又は信号変化率の情報を表示させる。一方、解剖学的構造の判定結果としては、異常と判定された画像領域が属する解剖学的構造の名称を分類レベル毎に表示させる。
Next, the
図13A、図13B図13Cに、その表示画面例を示す。
図13Aに示すように、制御部31は動態画像g1において異常と判定された画像領域g11の位置を指し示す指標情報a1を表示させるとともに、当該画像領域g11を正常と判定された画像領域と異なる色で表示させる等、表示形態を異なるものとして識別可能に表示させる。また、異常と判定された画像領域が属する解剖学的構造(分類レベル1)の境界線g12を表示させてその位置を示すとともに、解剖学的構造の名称を示す指標情報a2を表示させる。
FIG. 13A, FIG. 13B, and FIG. 13C show examples of the display screen.
As shown in FIG. 13A, the
図13Bに示すのは分類レベル2、図13Cは分類レベル3の場合の表示例である。この場合、異常と判定された画像領域g11が属する解剖学的構造の境界線g12は分類レベル2又は3に応じたものとする。また、指標情報a2の表示内容を、分類レベル2又は3に応じた名称に変更する。
このように、各分類レベルの表示画面を自動的に切り替えて表示してもよいし、切替操作があったときのみ切り替えて表示することとしてもよい。
FIG. 13B shows a display example in the case of
As described above, the display screens of the respective classification levels may be automatically switched and displayed, or may be switched and displayed only when a switching operation is performed.
図14は、異常と判定された画像領域について算出された面積変化率の表示画面例である。
例えば異常と判定された画像領域が2つある場合、制御部31は図14に示すようにこれらを異常領域A、Bと表してそれぞれ算出された面積変化率を時間位相(又は呼吸位相)に応じてグラフ化して表示させる。1つの呼吸サイクルの中で面積変化率をグラフ化することにより、異常領域Aでは呼気時に異常があり、異常領域Bでは一定間隔毎に異常が生じている等、異常の傾向を把握することが可能となる。
FIG. 14 is a display screen example of the area change rate calculated for the image area determined to be abnormal.
For example, when there are two image regions determined to be abnormal, the
この際、呼気と吸気の位相を把握するための参考情報として呼吸位相のグラフを表示させることとしてもよい。呼吸位相のグラフとしては、例えば図14に示すように横隔膜の高さを各時間位相の動態画像から求めて時間位相によってプロットしたものでもよいし、スパイロメータ等によって得られた測定値をプロットしたスパイログラムであってもよい。
なお、面積変化率の表示例のみ示したが、同様に信号変化率の表示を行ってもよい。
At this time, a breathing phase graph may be displayed as reference information for grasping the phase of expiration and inspiration. As a graph of the respiratory phase, for example, as shown in FIG. 14, the diaphragm height may be obtained from the dynamic image of each time phase and plotted by the time phase, or the measurement value obtained by a spirometer or the like is plotted. It may be a spirogram.
Although only the display example of the area change rate is shown, the signal change rate may be displayed similarly.
以上のように、本実施形態によれば、撮影装置10において被写体Wの胸部を動態撮影し、複数の時間位相における動態画像を生成する。画像処理装置40の画像解析部47により、各動態画像を用いて各動態画像に含まれる肺野領域を複数の領域に分割したその領域毎に面積変化率及び/又は信号変化率を算出し、当該面積変化率及び/又は信号変化率を元に換気能の正常又は異常の判定を行う。さらに、参照画像を用いて各動態画像につき、解剖学的構造の位置及び名称を判定し、その判定結果の情報を動態画像に付帯させる。診断用コンソール30の制御部31により、表示部34に動態画像を表示させるとともに、動態画像の付帯情報に基づいて、異常と判定された画像領域についてはその解剖学的構造の位置及び名称を表示させる。
As described above, according to the present embodiment, the
これにより、肺の換気能の診断に有効な情報を医師に提供することができる。すなわち、換気能の異常箇所の判定情報だけでなく、その異常箇所が影響する可能性のある解剖学的構造の位置及び名称の情報を医師に提供することができる。 This makes it possible to provide doctors with information that is effective in diagnosing the ventilation ability of the lungs. That is, it is possible to provide not only the determination information of the abnormal part of the ventilation ability but also the information of the position and name of the anatomical structure that the abnormal part may influence.
参照画像は解剖学的構造の分類レベルに応じて作成され、解剖学的構造の位置及び名称は複数の分類レベル毎に判定して表示するので、医師は各分類レベルで解剖学的構造を把握することができる。 Reference images are created according to the classification level of the anatomical structure, and the position and name of the anatomical structure are determined and displayed for each classification level, so the doctor grasps the anatomical structure at each classification level. can do.
また、参照画像は患者の属性毎に作成され、患者情報に対応する属性の参照画像を用いて解剖学的構造の位置又は名称の判定を行うので、精度の高い判定を行うことができる。 Moreover, since the reference image is created for each patient attribute and the position or name of the anatomical structure is determined using the reference image having the attribute corresponding to the patient information, a highly accurate determination can be performed.
また、異常と判定された画像領域を表示する際には、その画像領域を指し示す指標情報を表示する。さらに、正常と判定された領域とは異なる表示形態で表示するので、医師は異常と判定された画像領域を容易に把握することができる。 Further, when displaying an image area determined to be abnormal, index information indicating the image area is displayed. Furthermore, since the display is performed in a display form different from the area determined to be normal, the doctor can easily grasp the image area determined to be abnormal.
なお、上述した説明は本発明の好適な一例であり、これに限定されるものではない。
例えば、参照画像は通常呼吸と深呼吸の呼吸方法別に作成して用意することとしてもよい。
The above description is a preferred example of the present invention, and the present invention is not limited to this.
For example, the reference image may be prepared and prepared for each breathing method of normal breathing and deep breathing.
また、換気能の正常又は異常の判定は、全ての時間位相の動態画像を用いるのではなく、指定された時間位相の動態画像のみを用いて判定することとしてもよい。例えば、医師が呼気から吸気への変換期に異常があると判断しており、その変換期における判定のみを希望する場合には、診断用コンソール30において判定を行う時間位相として呼気から吸気への変換期の指定操作を行う。この指定操作を受けて画像処理装置40の画像解析部47では、呼気から吸気への変換期周辺の時間位相に係る動態画像のみを用いて正常又は異常の判定の処理を実行すればよい。これにより、処理する画像枚数が減り処理効率が向上する。
Further, the determination of normal or abnormal ventilation ability may be made by using only a dynamic image of a specified time phase instead of using a dynamic image of all time phases. For example, if the doctor determines that there is an abnormality in the conversion period from exhalation to inspiration, and only wants to make a determination in the conversion period, the time phase for determination in the
また、上記実施形態では胸部を動態撮影し、肺の換気能の正常又は異常を判定するとともに、肺野の解剖学的構造の位置及び名称を判定する例を説明したが、例えば心臓部等の他の部位に適用することとしてもよい。つまり、心臓部を動態撮影し、得られた動態画像を用いて上述した方法と同様に心臓部の時間の経過に伴う面積変化率及び/又は信号変化率を算出して、この面積変化率及び/又は信号変化率から心臓部の心機能の正常又は異常を判定する。正常又は異常の判定基準は心臓部に応じた基準とする。また、心臓部について予め右心房や左心室等の解剖学的構造の位置及び名称を定めた参照画像を用いて撮影した動態画像について解剖学的構造の位置及び名称を判定する。これによれば、心機能についても医師に情報提供することが可能となる。 In the above embodiment, the chest is dynamically photographed to determine normality or abnormality of the lung ventilation and the position and name of the anatomical structure of the lung field are described. It may be applied to other parts. That is, by taking a dynamic image of the heart, and using the obtained dynamic image, the area change rate and / or signal change rate with the passage of time of the heart is calculated in the same manner as described above. The normal or abnormal cardiac function of the heart is determined from the signal change rate. Normal or abnormal criteria are determined according to the heart. Further, the position and name of the anatomical structure are determined with respect to the dynamic image taken using the reference image in which the position and name of the anatomical structure such as the right atrium and the left ventricle are determined in advance for the heart. According to this, it becomes possible to provide information to the doctor regarding the cardiac function.
また、診断用コンソール30において判定結果を表示する構成を説明したが、撮影用コンソール20や他の装置(診断用に用いられるPC等)において表示することとしてもよい。また、画像解析を行う画像処理装置40を設けて画像処理装置40において判定を行う構成を説明したが、診断用コンソール30や他の装置において上記換気能の正常又は異常の判定、解剖学的構造の判定を行うプログラムをインストールし、算出を行うこととしてもよい。
Moreover, although the structure which displays a determination result in the
また、上記説明した処理に係るプログラムを記憶するコンピュータ読み取り可能な媒体としては、ROM等のメモリの他、DVD等の可搬型のものも適用可能である。また、婦ログラムのデータをネットワークを介して提供する媒体として、キャリアウェーブ(搬送波)も適用可能である。 Further, as a computer-readable medium for storing the program related to the above-described processing, a portable type such as a DVD can be applied in addition to a memory such as a ROM. Also, a carrier wave can be used as a medium for providing female program data via a network.
1 動態画像処理システム
10 撮影装置
11 X線源
12 検出器
13 読取部
14 サイクル検出部
20 撮影用コンソール
21 制御部
30 診断用コンソール
31 制御部
34 表示部
40 画像処理装置
46 画像処理部
47 画像解析部
DESCRIPTION OF
Claims (5)
前記複数の時間位相における動態画像を用いて、各動態画像に含まれる肺野領域を複数の領域に分割した領域毎に時間の経過に伴う面積変化率及び/又は信号変化率を算出し、当該算出した面積変化率及び/又は信号変化率に基づいて換気能の正常又は異常を領域毎に判定する第1判定手段と、
肺における解剖学的構造の位置及び名称が予め定められた参照画像を用いて、前記複数の時間位相における動態画像につき、解剖学的構造の位置及び名称を判定する第2判定手段と、
前記複数の動態画像のうち少なくとも1つを表示する表示手段と、
前記表示された動態画像において、前記異常と判定された領域につき、前記判定された解剖学的構造の位置及び名称の情報を表示させる制御手段と、
を備える動態画像処理システム。 An imaging means for capturing a dynamic image of the subject's chest and generating dynamic images at a plurality of time phases;
Using the dynamic images in the plurality of temporal phases, calculating the area change rate and / or signal change rate over time for each region obtained by dividing the lung field region included in each dynamic image into a plurality of regions, First determination means for determining normal or abnormal ventilation capacity for each region based on the calculated area change rate and / or signal change rate;
Second determination means for determining the position and name of the anatomical structure for the dynamic images in the plurality of temporal phases using a reference image in which the position and name of the anatomical structure in the lung are predetermined;
Display means for displaying at least one of the plurality of dynamic images;
In the displayed dynamic image, control means for displaying information on the position and name of the determined anatomical structure for the area determined to be abnormal,
A dynamic image processing system comprising:
前記第2判定手段は、前記複数の分類レベル毎に作成された参照画像を用いて、各分類レベルにおける解剖学的構造の位置及び名称を判定する請求項1に記載の動態画像処理システム。 Depending on the position and name of the anatomical structure defined at multiple classification levels, the reference image is created for each of multiple classification levels,
The dynamic image processing system according to claim 1, wherein the second determination unit determines a position and a name of an anatomical structure at each classification level using a reference image created for each of the plurality of classification levels.
前記動態撮影の被写体の属性に関する情報を入力する入力手段を備え、
前記第2判定手段は、前記被写体の属性毎に作成された参照画像のうち、前記入力された情報に対応する属性の参照画像を用いる請求項1又は2に記載の動態画像処理システム。 The reference image is created for each attribute of the subject,
Input means for inputting information on the attributes of the subject of the dynamic shooting,
The dynamic image processing system according to claim 1, wherein the second determination unit uses a reference image having an attribute corresponding to the input information among reference images created for each attribute of the subject.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007334515A JP2009153677A (en) | 2007-12-26 | 2007-12-26 | Kinetic image processing system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007334515A JP2009153677A (en) | 2007-12-26 | 2007-12-26 | Kinetic image processing system |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2009153677A true JP2009153677A (en) | 2009-07-16 |
Family
ID=40958344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007334515A Pending JP2009153677A (en) | 2007-12-26 | 2007-12-26 | Kinetic image processing system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2009153677A (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012026145A1 (en) * | 2010-08-27 | 2012-03-01 | コニカミノルタエムジー株式会社 | Diagnosis assistance system and program |
JP2012115581A (en) * | 2010-12-03 | 2012-06-21 | Konica Minolta Medical & Graphic Inc | Thoracic diagnosis support information generation method and thoracic diagnosis support system |
CN103079466A (en) * | 2010-08-27 | 2013-05-01 | 柯尼卡美能达医疗印刷器材株式会社 | Thoracic diagnosis assistance system and program |
JP2013169400A (en) * | 2012-02-22 | 2013-09-02 | Konica Minolta Inc | Dynamic state diagnosis support information generating system, dynamic state diagnosis support information generating method, and dynamic state analyzer |
WO2013140749A1 (en) * | 2012-03-23 | 2013-09-26 | 富士フイルム株式会社 | Image analysis device, method, and program |
WO2014162911A1 (en) * | 2013-04-05 | 2014-10-09 | 東芝メディカルシステムズ株式会社 | Medical image processing apparatus and medical image processing method |
JP2014528348A (en) * | 2011-10-12 | 2014-10-27 | ザ・ジョンズ・ホプキンス・ユニバーシティ | Evaluation of regional cardiac function and dyssynchrony from dynamic imaging modalities using intracardiac motion |
JP2015226710A (en) * | 2014-06-02 | 2015-12-17 | 株式会社東芝 | Medical image processor and medical image processing method |
JP2016116867A (en) * | 2014-12-22 | 2016-06-30 | 株式会社東芝 | Medical image processing apparatus, medical image diagnostic apparatus and medical image processing program |
JP2017018681A (en) * | 2016-10-17 | 2017-01-26 | コニカミノルタ株式会社 | Kinetic analysis system |
US9801555B2 (en) | 2013-08-28 | 2017-10-31 | Konica Minolta, Inc. | Thoracic diagnosis assistance system |
JP2018532515A (en) * | 2015-11-09 | 2018-11-08 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | X-ray image intake quality monitoring |
CN115299956A (en) * | 2022-08-19 | 2022-11-08 | 山东大学 | Myocardial ischemia detection method and system based on deterministic learning and electrocardiogram |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6214590A (en) * | 1985-07-12 | 1987-01-23 | Toshiba Corp | Image diagnosing device |
JP2004000410A (en) * | 2002-04-03 | 2004-01-08 | Canon Inc | Processing radiographic image, program, computer-readable medium, and method and system for supporting image diagnosis |
WO2005011501A1 (en) * | 2003-08-01 | 2005-02-10 | Hitachi Medical Corporation | Medical image diagnosis support device and method |
JP2005218796A (en) * | 2004-02-09 | 2005-08-18 | Matsushita Electric Ind Co Ltd | Medical image processor and medical image processing method |
WO2006011545A1 (en) * | 2004-07-30 | 2006-02-02 | Hitachi Medical Corporation | Medical image diagnosis assisting system, device and image processing program |
JP2006239195A (en) * | 2005-03-04 | 2006-09-14 | Fuji Photo Film Co Ltd | Computer-aided image diagnosis support method, image diagnosis support device and program |
JP2007068554A (en) * | 2005-09-02 | 2007-03-22 | Fujifilm Corp | Positioning device, positioning method and its program |
JP2007216007A (en) * | 2006-01-25 | 2007-08-30 | Siemens Medical Solutions Usa Inc | Method and system for labeling and identifying lymph node in medical image |
WO2007129493A1 (en) * | 2006-05-02 | 2007-11-15 | National University Corporation Nagoya University | Medical image observation support device |
-
2007
- 2007-12-26 JP JP2007334515A patent/JP2009153677A/en active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6214590A (en) * | 1985-07-12 | 1987-01-23 | Toshiba Corp | Image diagnosing device |
JP2004000410A (en) * | 2002-04-03 | 2004-01-08 | Canon Inc | Processing radiographic image, program, computer-readable medium, and method and system for supporting image diagnosis |
WO2005011501A1 (en) * | 2003-08-01 | 2005-02-10 | Hitachi Medical Corporation | Medical image diagnosis support device and method |
JP2005218796A (en) * | 2004-02-09 | 2005-08-18 | Matsushita Electric Ind Co Ltd | Medical image processor and medical image processing method |
WO2006011545A1 (en) * | 2004-07-30 | 2006-02-02 | Hitachi Medical Corporation | Medical image diagnosis assisting system, device and image processing program |
JP2006239195A (en) * | 2005-03-04 | 2006-09-14 | Fuji Photo Film Co Ltd | Computer-aided image diagnosis support method, image diagnosis support device and program |
JP2007068554A (en) * | 2005-09-02 | 2007-03-22 | Fujifilm Corp | Positioning device, positioning method and its program |
JP2007216007A (en) * | 2006-01-25 | 2007-08-30 | Siemens Medical Solutions Usa Inc | Method and system for labeling and identifying lymph node in medical image |
WO2007129493A1 (en) * | 2006-05-02 | 2007-11-15 | National University Corporation Nagoya University | Medical image observation support device |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103068312B (en) * | 2010-08-27 | 2015-07-15 | 柯尼卡美能达医疗印刷器材株式会社 | Diagnosis assistance system and program |
CN103068312A (en) * | 2010-08-27 | 2013-04-24 | 柯尼卡美能达医疗印刷器材株式会社 | Diagnosis assistance system and program |
CN103079466A (en) * | 2010-08-27 | 2013-05-01 | 柯尼卡美能达医疗印刷器材株式会社 | Thoracic diagnosis assistance system and program |
JP5874636B2 (en) * | 2010-08-27 | 2016-03-02 | コニカミノルタ株式会社 | Diagnosis support system and program |
CN104887258A (en) * | 2010-08-27 | 2015-09-09 | 柯尼卡美能达医疗印刷器材株式会社 | Diagnosis assistance system |
WO2012026145A1 (en) * | 2010-08-27 | 2012-03-01 | コニカミノルタエムジー株式会社 | Diagnosis assistance system and program |
US9064302B2 (en) | 2010-08-27 | 2015-06-23 | Konica Minolta, Inc. | Diagnosis assistance system and computer readable storage medium |
JP2012115581A (en) * | 2010-12-03 | 2012-06-21 | Konica Minolta Medical & Graphic Inc | Thoracic diagnosis support information generation method and thoracic diagnosis support system |
US10448901B2 (en) | 2011-10-12 | 2019-10-22 | The Johns Hopkins University | Methods for evaluating regional cardiac function and dyssynchrony from a dynamic imaging modality using endocardial motion |
JP2014528348A (en) * | 2011-10-12 | 2014-10-27 | ザ・ジョンズ・ホプキンス・ユニバーシティ | Evaluation of regional cardiac function and dyssynchrony from dynamic imaging modalities using intracardiac motion |
JP2013169400A (en) * | 2012-02-22 | 2013-09-02 | Konica Minolta Inc | Dynamic state diagnosis support information generating system, dynamic state diagnosis support information generating method, and dynamic state analyzer |
US9117287B2 (en) | 2012-03-23 | 2015-08-25 | Fujifilm Corporation | Image analysis apparatus, method, and program |
JP2013192912A (en) * | 2012-03-23 | 2013-09-30 | Fujifilm Corp | Image analysis device and method, and program |
WO2013140749A1 (en) * | 2012-03-23 | 2013-09-26 | 富士フイルム株式会社 | Image analysis device, method, and program |
JP2014210171A (en) * | 2013-04-05 | 2014-11-13 | 東芝メディカルシステムズ株式会社 | Medical image processor and medical image processing method |
WO2014162911A1 (en) * | 2013-04-05 | 2014-10-09 | 東芝メディカルシステムズ株式会社 | Medical image processing apparatus and medical image processing method |
CN105101878A (en) * | 2013-04-05 | 2015-11-25 | 东芝医疗系统株式会社 | Medical image processing apparatus and medical image processing method |
US9730666B2 (en) | 2013-04-05 | 2017-08-15 | Toshiba Medical Systems Corporation | Medical image processing apparatus and medical image processing method |
US9801555B2 (en) | 2013-08-28 | 2017-10-31 | Konica Minolta, Inc. | Thoracic diagnosis assistance system |
JP2015226710A (en) * | 2014-06-02 | 2015-12-17 | 株式会社東芝 | Medical image processor and medical image processing method |
JP2016116867A (en) * | 2014-12-22 | 2016-06-30 | 株式会社東芝 | Medical image processing apparatus, medical image diagnostic apparatus and medical image processing program |
JP2018532515A (en) * | 2015-11-09 | 2018-11-08 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | X-ray image intake quality monitoring |
JP7004648B2 (en) | 2015-11-09 | 2022-01-21 | コーニンクレッカ フィリップス エヌ ヴェ | X-ray image Intake quality monitoring |
JP2017018681A (en) * | 2016-10-17 | 2017-01-26 | コニカミノルタ株式会社 | Kinetic analysis system |
CN115299956A (en) * | 2022-08-19 | 2022-11-08 | 山东大学 | Myocardial ischemia detection method and system based on deterministic learning and electrocardiogram |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2009153677A (en) | Kinetic image processing system | |
JP6436182B2 (en) | Dynamic image analyzer | |
JP5874636B2 (en) | Diagnosis support system and program | |
JP5962237B2 (en) | Chest diagnosis support information generation method | |
JP5919717B2 (en) | Dynamic medical image generation system | |
JP6597548B2 (en) | Dynamic analysis system | |
JP5200656B2 (en) | Dynamic shooting system | |
JP6805918B2 (en) | Radiation image processing equipment and radiation imaging system | |
JP2012110400A (en) | Dynamic diagnosis support information generation system | |
JP6743662B2 (en) | Dynamic image processing system | |
JPWO2012026146A1 (en) | Chest diagnosis support system and program | |
JP2009153678A (en) | Kinetic image processing system | |
JP2018000281A (en) | Dynamic state analysis system | |
JP5617559B2 (en) | Radiation imaging system | |
JP5617577B2 (en) | Chest diagnosis support information generation method | |
JP5125750B2 (en) | Dynamic image diagnosis support system | |
JP2017169830A (en) | Dynamic analysis apparatus | |
JP2012110399A (en) | Dynamic diagnosis support information generation system | |
JP5625799B2 (en) | Dynamic diagnosis support information generation system | |
JP2018148965A (en) | Dynamic analysis system | |
JP2020044445A (en) | Dynamic analysis system, program, and dynamic analysis apparatus | |
JP2014147844A (en) | Dynamic diagnosis support information generation system | |
JP2012245395A (en) | Dynamic image-capturing system and diagnosis-assisting information-generating method | |
JP2012110398A (en) | Diagnosis support information generation system | |
JP2012115582A (en) | Thoracic part diagnosis aiding system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100308 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120201 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20121106 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20130305 |