JP7255319B2 - 動態解析装置、動態解析システム及びプログラム - Google Patents
動態解析装置、動態解析システム及びプログラム Download PDFInfo
- Publication number
- JP7255319B2 JP7255319B2 JP2019070992A JP2019070992A JP7255319B2 JP 7255319 B2 JP7255319 B2 JP 7255319B2 JP 2019070992 A JP2019070992 A JP 2019070992A JP 2019070992 A JP2019070992 A JP 2019070992A JP 7255319 B2 JP7255319 B2 JP 7255319B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- measurement points
- movement
- dynamic
- chest
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims description 25
- 238000012933 kinetic analysis Methods 0.000 title 1
- 238000005259 measurement Methods 0.000 claims description 66
- 230000033001 locomotion Effects 0.000 claims description 62
- 238000003384 imaging method Methods 0.000 claims description 61
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000002601 radiography Methods 0.000 claims description 5
- 230000005855 radiation Effects 0.000 description 40
- 238000012545 processing Methods 0.000 description 25
- 230000003434 inspiratory effect Effects 0.000 description 18
- 208000006545 Chronic Obstructive Pulmonary Disease Diseases 0.000 description 15
- 238000004891 communication Methods 0.000 description 15
- 238000010191 image analysis Methods 0.000 description 14
- 238000000034 method Methods 0.000 description 12
- 230000000241 respiratory effect Effects 0.000 description 10
- 238000001514 detection method Methods 0.000 description 9
- 238000003745 diagnosis Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 210000004072 lung Anatomy 0.000 description 7
- 230000004202 respiratory function Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 206010073306 Exposure to radiation Diseases 0.000 description 4
- 239000004065 semiconductor Substances 0.000 description 4
- 210000000988 bone and bone Anatomy 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000002438 flame photometric detection Methods 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 208000028399 Critical Illness Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
胸部の動態を放射線撮影することにより取得された胸部動態画像の任意のフレーム画像の2以上の肋骨領域上に計測点を1点ずつ設定する設定手段と、
前記設定された複数の計測点を前記胸部動態画像の同一区間のフレーム画像においてトラッキングして前記複数の計測点のそれぞれに対応する追跡点を時間方向に複数取得し、取得した追跡点ごとに、時間方向に隣接する追跡点からの移動量及び/又は移動方向を算出して、前記複数の計測点の移動量及び/又は移動方向の時系列データを取得する算出手段と、
を備え、
前記算出手段は、
前記複数の計測点のうちの2つの計測点の組ごとに前記時系列データの類似度を算出する。
コンピューターを、
胸部の動態を放射線撮影することにより取得された胸部動態画像の任意のフレーム画像の2以上の肋骨領域上に計測点を1点ずつ設定する設定手段、
前記設定された複数の計測点を前記胸部動態画像の同一区間のフレーム画像においてトラッキングして前記複数の計測点のそれぞれに対応する追跡点を時間方向に複数取得し、取得した追跡点ごとに、時間方向に隣接する追跡点からの移動量及び/又は移動方向を算出して、前記複数の計測点の移動量及び/又は移動方向の時系列データを取得する算出手段、
として機能させ、
前記算出手段は、
前記複数の計測点のうちの2つの計測点の組ごとに前記時系列データの類似度を算出する。
まず、構成を説明する。
図1に、本実施形態における動態解析システム100の全体構成を示す。
図1に示すように、動態解析システム100は、撮影装置1と、撮影用コンソール2とが通信ケーブル等により接続され、撮影用コンソール2と、診断用コンソール3とがLAN(Local Area Network)等の通信ネットワークNTを介して接続されて構成されている。動態解析システム100を構成する各装置は、DICOM(Digital Image and Communications in Medicine)規格に準じており、各装置間の通信は、DICOMに則って行われる。
撮影装置1は、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性(サイクル)を持つ胸部の動態を撮影する撮影手段である。動態撮影とは、被写体に対し、X線等の放射線をパルス状にして所定時間間隔で繰り返し照射するか(パルス照射)、もしくは、低線量率にして途切れなく継続して照射する(連続照射)ことで、動態を示す複数の画像を取得することをいう。動態撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。なお、以下の実施形態では、パルス照射により胸部の動態撮影を行う場合を例にとり説明する。
放射線照射制御装置12は、撮影用コンソール2に接続されており、撮影用コンソール2から入力された放射線照射条件に基づいて放射線源11を制御して放射線撮影を行う。撮影用コンソール2から入力される放射線照射条件は、例えば、パルスレート、パルス幅、パルス間隔、1撮影あたりの撮影フレーム数、X線管電流の値、X線管電圧の値、付加フィルター種等である。パルスレートは、1秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射1回当たりの放射線照射時間である。パルス間隔は、1回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。
放射線検出部13は、被写体Mを挟んで放射線源11と対向するように設けられている。
撮影用コンソール2は、放射線照射条件や画像読取条件を撮影装置1に出力して撮影装置1による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置1により取得された動態画像を撮影技師等の撮影実施者によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
撮影用コンソール2は、図1に示すように、制御部21、記憶部22、操作部23、表示部24、通信部25を備えて構成され、各部はバス26により接続されている。
)等により構成される。制御部21のCPUは、操作部23の操作に応じて、記憶部22に記憶されているシステムプログラムや各種処理プログラムを読み出してRAM内に展開し、展開されたプログラムに従って後述する撮影制御処理を始めとする各種処理を実行し、撮影用コンソール2各部の動作や、撮影装置1の放射線照射動作及び読み取り動作を集中制御する。
診断用コンソール3は、撮影用コンソール2から動態画像を取得し、取得した動態画像や動態画像の解析結果を表示して医師の診断を支援するための動態解析装置である。
診断用コンソール3は、図1に示すように、制御部31、記憶部32、操作部33、表示部34、通信部35を備えて構成され、各部はバス36により接続されている。
次に、上記動態解析システム100における動作について説明する。
まず、撮影装置1、撮影用コンソール2による撮影動作について説明する。
図2に、撮影用コンソール2の制御部21において実行される撮影制御処理を示す。撮影制御処理は、制御部21と記憶部22に記憶されているプログラムとの協働により実行される。
次に、診断用コンソール3における動作について説明する。
診断用コンソール3においては、通信部35を介して撮影用コンソール2から動態画像の一連のフレーム画像が受信されると、制御部31と記憶部32に記憶されているプログラムとの協働により図3に示す画像解析処理が実行される。
まず、変数nに1が設定される(ステップS11)。
ここでは、n番目のフレーム画像の2以上の肋骨領域上に計測点が1点ずつ設定される。
例えば、図4に示すように、表示部34に1枚目のフレーム画像を表示し、ユーザーによる操作部33の操作により複数の肋骨領域上に1点ずつ指定された複数の点が計測点として設定される。
あるいは、制御部31が自動的に肋骨領域上の複数点を計測点に設定することとしてもよい。例えば、1枚目のフレーム画像から肋骨領域を抽出し、抽出した肋骨領域上の代表点(例えば、第〇肋骨、第△肋骨、第□肋骨の最もカーブしている点等)に自動的に計測点を設定することとしてもよい。肋骨領域の抽出は、例えば、特開平5-176919豪公報に記載されているモデル関数とSobelオペレーターを用いた肋骨抽出方法のように、公知の画像処理技術によって抽出することができる。
なお、本実施形態では、1枚目のフレーム画像に計測点を設定する場合を例にとり説明するが、1枚目に限らず、任意のフレーム画像に計測点を設定することができる。
例えば、n番目のフレーム画像において、n-1番目のフレーム画像で追跡点として抽出された画素(n=2の場合は、n-1番目のフレーム画像で計測点として設定された画素)を中心とするM×N画素(M、Nは正の整数、例えば、3×3)からなるエリアに探索エリアが設定され、探索エリア内で画素値がn-1番目のフレーム画像で追跡点として抽出された画素の画素値に最も近い画素が追跡点として抽出される。
例えば、図5に示すように、追跡点がAtからAt+1に移動した場合、AtからAt+1までの距離が移動量、AtとAt+1を結ぶ線とx軸方向とのなす角度θが移動方向として算出される。これにより、各計測点のフレーム画像間の移動量及び/又は移動方向が算出される。
2つの計測点の動きの同期性を示す指標値とは、2つの計測点の動き(移動量や移動方向)がどれだけ類似しているかを示す指標値であり、例えば、2つの計測点の時系列データの相互相関係数、コサイン類似度が挙げられる。
呼吸位相の特定は、例えば、以下の(1)~(5)の手順により行うことができる。なお、以下の説明では、画像の左上の座標を(0、0)とし、画像の右側、下側へいくほど座標値が増えることとして説明する。
(1)まず、動態画像の各フレーム画像から横隔膜の位置を算出する(図6(a)参照)。
例えば、各フレーム画像の肺野領域の下側のエッジ部分を横隔膜境界部として抽出し、横隔膜境界部の或るx座標の位置に基準点を設定し、設定した基準点のy座標を横隔膜の位置として求めることができる。肺野領域は、公知のいずれの方法を用いて抽出してもよい。例えば、各画素の画素値のヒストグラムから判別分析によって閾値を求め、この閾値より高信号の領域を肺野領域候補として1次抽出する。次いで、1次抽出された肺野領域候補の境界付近でエッジ検出を行い、境界付近の小領域でエッジが最大となる点を境界に沿って抽出すれば肺野領域の境界を抽出することができる。
(2)横隔膜の位置の最高位、最低位を算出する。
(3)(1)の微分値を算出する(図6(b)参照)。
(4)呼気相については、最高位のフレーム画像以降に微分値が“+”→”-”となる最初の変曲点(微分値0の点)を”呼気相開始”とし、呼気相開始後に最初に“-”→”+”となる変曲点を”呼気相終了”とする。ノイズの影響を考慮し、開始/終了条件に「連続してiフレーム(iは正の整数)微分値0が続く」ことを加えてもよい。
(5)吸気相については、最低位のフレーム画像以降に微分値が“-”→”+”となる最初の変曲点を”吸気相開始”とし、吸気相開始後に最初に“+”→”-”となる変曲点を”吸気相終了”とする。ノイズの影響を考慮し、開始/終了条件に「連続してiフレーム微分値0が続く」ことを加えてもよい。
したがって、医師等のユーザーは、胸部動態画像における複数の肋骨の動きの関係を時系列に捉えることができるので、呼吸機能の状態(正常、異常や疾患の程度等)を把握することが可能となる。そのため、呼吸検査による患者の負担を軽減することが可能となる。また、時系列データは、治療効果の効果確認にも使用することができ、治療計画の立案を支援することができる。
図7は、吸気相における呼吸正常者、COPD軽症患者、COPD重症患者の上、中、下の後方肋骨(図9参照)の動き方向の同期関係を示す図であり、図8は、その動き方向を矢印で模式的に示す図である。図7、図8に示すように、呼吸正常者、COPD軽症患者、COPD重症患者では、吸気相の上、中、下の後方肋骨の動き方向の同期性が異なる。呼吸正常者については、上、中の後方肋骨と下の後方肋骨が吸気/呼気相で異なる方向に動く傾向がある。COPDの軽症患者については、吸気相において中の後方肋骨と下の後方肋骨が同方向に同期して動く傾向がある。COPDの重症患者については、上、中、下の後方肋骨が吸気相で同方向に同期して動く傾向がある。
例えば、上記実施形態においては、胸部動態画像の全区間のフレーム画像から複数の計測点をトラッキングして時系列データを取得することとしたが、これに限定されず、胸部動態画像の一部の区間のフレーム画像から複数の計測点をトラッキングして時系列データを取得することしてもよい。この場合、複数の計測点のトラッキングは胸部動態画像の同一区間にて行い、同一区間の時系列データを取得する。なお、本発明では、設定された複数の計測点で同一区間にてトラッキングや時系列データの取得を行えばよく、これに加えて、一部の計測点において他の区間でトラッキングや時系列データを取得することを妨げるものではない。
等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。
1 撮影装置
11 放射線源
12 放射線照射制御装置
13 放射線検出部
14 読取制御装置
2 撮影用コンソール
21 制御部
22 記憶部
23 操作部
24 表示部
25 通信部
26 バス
3 診断用コンソール
31 制御部
32 記憶部
33 操作部
34 表示部
35 通信部
36 バス
Claims (5)
- 胸部の動態を放射線撮影することにより取得された胸部動態画像の任意のフレーム画像の2以上の肋骨領域上に計測点を1点ずつ設定する設定手段と、
前記設定された複数の計測点を前記胸部動態画像の同一区間のフレーム画像においてトラッキングして前記複数の計測点のそれぞれに対応する追跡点を時間方向に複数取得し、取得した追跡点ごとに、時間方向に隣接する追跡点からの移動量及び/又は移動方向を算出して、前記複数の計測点の移動量及び/又は移動方向の時系列データを取得する算出手段と、
を備え、
前記算出手段は、
前記複数の計測点のうちの2つの計測点の組ごとに前記時系列データの類似度を算出する、動態解析装置。 - 前記類似度は、相互相関係数又はコサイン類似度である、請求項1に記載の動態解析装置。
- 前記算出手段による算出結果を表示する表示手段を備える請求項1又は2に記載の動態解析装置。
- 胸部の動態を放射線撮影することにより胸部動態画像を取得する撮影装置と、
請求項1~3のいずれか一項に記載の動態解析装置と、
を備える動態解析システム。 - コンピューターを、
胸部の動態を放射線撮影することにより取得された胸部動態画像の任意のフレーム画像の2以上の肋骨領域上に計測点を1点ずつ設定する設定手段、
前記設定された複数の計測点を前記胸部動態画像の同一区間のフレーム画像においてトラッキングして前記複数の計測点のそれぞれに対応する追跡点を時間方向に複数取得し、取得した追跡点ごとに、時間方向に隣接する追跡点からの移動量及び/又は移動方向を算出して、前記複数の計測点の移動量及び/又は移動方向の時系列データを取得する算出手段、
として機能させ、
前記算出手段は、
前記複数の計測点のうちの2つの計測点の組ごとに前記時系列データの類似度を算出する、プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019070992A JP7255319B2 (ja) | 2019-04-03 | 2019-04-03 | 動態解析装置、動態解析システム及びプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019070992A JP7255319B2 (ja) | 2019-04-03 | 2019-04-03 | 動態解析装置、動態解析システム及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2020168172A JP2020168172A (ja) | 2020-10-15 |
JP7255319B2 true JP7255319B2 (ja) | 2023-04-11 |
Family
ID=72745453
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019070992A Active JP7255319B2 (ja) | 2019-04-03 | 2019-04-03 | 動態解析装置、動態解析システム及びプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7255319B2 (ja) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003529433A (ja) | 2000-04-07 | 2003-10-07 | レモ・イーエイチエフ | 呼吸運動の測定および装置 |
JP2004167109A (ja) | 2002-11-21 | 2004-06-17 | Osaka Industrial Promotion Organization | 3次元計測方法、3次元計測システム、画像処理装置、及びコンピュータプログラム |
JP2010158288A (ja) | 2009-01-06 | 2010-07-22 | Konica Minolta Holdings Inc | 画像表示装置およびプログラム |
JP2015043894A (ja) | 2013-08-28 | 2015-03-12 | 国立大学法人金沢大学 | X線動画像解析装置、x線動画像解析プログラム及びx線動画像撮像装置 |
JP2015229101A (ja) | 2014-06-03 | 2015-12-21 | 株式会社イデアクエスト | 呼吸運動測定装置 |
JP2018183493A (ja) | 2017-04-27 | 2018-11-22 | コニカミノルタ株式会社 | 画像表示システム及び画像処理装置 |
JP2020168173A (ja) | 2019-04-03 | 2020-10-15 | コニカミノルタ株式会社 | 動態解析装置、動態解析システム及びプログラム |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3134009B2 (ja) * | 1990-11-21 | 2001-02-13 | アーチ・デベロップメント・コーポレーション | 画像処理方法及び装置 |
-
2019
- 2019-04-03 JP JP2019070992A patent/JP7255319B2/ja active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003529433A (ja) | 2000-04-07 | 2003-10-07 | レモ・イーエイチエフ | 呼吸運動の測定および装置 |
JP2004167109A (ja) | 2002-11-21 | 2004-06-17 | Osaka Industrial Promotion Organization | 3次元計測方法、3次元計測システム、画像処理装置、及びコンピュータプログラム |
JP2010158288A (ja) | 2009-01-06 | 2010-07-22 | Konica Minolta Holdings Inc | 画像表示装置およびプログラム |
JP2015043894A (ja) | 2013-08-28 | 2015-03-12 | 国立大学法人金沢大学 | X線動画像解析装置、x線動画像解析プログラム及びx線動画像撮像装置 |
JP2015229101A (ja) | 2014-06-03 | 2015-12-21 | 株式会社イデアクエスト | 呼吸運動測定装置 |
JP2018183493A (ja) | 2017-04-27 | 2018-11-22 | コニカミノルタ株式会社 | 画像表示システム及び画像処理装置 |
JP2020168173A (ja) | 2019-04-03 | 2020-10-15 | コニカミノルタ株式会社 | 動態解析装置、動態解析システム及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
JP2020168172A (ja) | 2020-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6436182B2 (ja) | 動態画像解析装置 | |
JP6413927B2 (ja) | 動態解析装置及び動態解析システム | |
JP6217241B2 (ja) | 胸部診断支援システム | |
JP6743662B2 (ja) | 動態画像処理システム | |
JP6958202B2 (ja) | 動態画像処理装置及びプログラム | |
JP6540807B2 (ja) | 撮影用コンソール | |
JP2017200565A (ja) | 動態解析装置及び動態解析システム | |
JP6740910B2 (ja) | 動態画像処理システム | |
JP2019051322A (ja) | 動態解析システム | |
JP2017169830A (ja) | 動態解析装置 | |
JP6690774B2 (ja) | 動態解析システム、プログラム及び動態解析装置 | |
JP2016209267A (ja) | 医用画像処理装置及びプログラム | |
JP2018196693A (ja) | 動態解析システム | |
JP7099086B2 (ja) | 動態画像処理装置及びプログラム | |
JP7047806B2 (ja) | 動態解析装置、動態解析システム、動態解析プログラム及び動態解析方法 | |
JP7255319B2 (ja) | 動態解析装置、動態解析システム及びプログラム | |
JP2014079655A (ja) | 放射線画像撮影システム | |
JP6950483B2 (ja) | 動態撮影システム | |
JP6962030B2 (ja) | 動態解析装置、動態解析システム、動態解析プログラム及び動態解析方法 | |
JP2017217047A (ja) | 画像表示システム | |
JP6888721B2 (ja) | 動態画像処理装置、動態画像処理プログラム及び動態画像処理方法 | |
JP6874484B2 (ja) | 動態画像処理システム | |
JP2022143288A (ja) | 動態画像解析装置及びプログラム | |
JP2020141841A (ja) | 動態解析装置及びプログラム | |
WO2011093221A1 (ja) | 動態画像処理システム及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211224 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20221122 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20221129 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230120 |
|
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: 20230228 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230313 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7255319 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |