JP6772908B2 - 動態解析システム及びプログラム - Google Patents
動態解析システム及びプログラム Download PDFInfo
- Publication number
- JP6772908B2 JP6772908B2 JP2017045532A JP2017045532A JP6772908B2 JP 6772908 B2 JP6772908 B2 JP 6772908B2 JP 2017045532 A JP2017045532 A JP 2017045532A JP 2017045532 A JP2017045532 A JP 2017045532A JP 6772908 B2 JP6772908 B2 JP 6772908B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- analysis result
- dynamic
- result image
- feature amount
- 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 149
- 230000006870 function Effects 0.000 claims description 77
- 230000008859 change Effects 0.000 claims description 75
- 210000004072 lung Anatomy 0.000 claims description 33
- 230000017531 blood circulation Effects 0.000 claims description 26
- 230000000737 periodic effect Effects 0.000 claims description 20
- 238000009423 ventilation Methods 0.000 claims description 20
- 238000003745 diagnosis Methods 0.000 claims description 14
- 238000012933 kinetic analysis Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 description 49
- 238000003384 imaging method Methods 0.000 description 39
- 230000008569 process Effects 0.000 description 39
- 230000005855 radiation Effects 0.000 description 35
- 238000012545 processing Methods 0.000 description 31
- 238000004891 communication Methods 0.000 description 17
- 210000000038 chest Anatomy 0.000 description 12
- 238000001514 detection method Methods 0.000 description 10
- 210000003109 clavicle Anatomy 0.000 description 8
- 210000001562 sternum Anatomy 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000003247 decreasing effect Effects 0.000 description 4
- 230000029058 respiratory gaseous exchange Effects 0.000 description 4
- 239000004065 semiconductor Substances 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000010349 pulsation Effects 0.000 description 3
- 238000010835 comparative analysis Methods 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 238000005206 flow analysis Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000003434 inspiratory effect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000000241 respiratory effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 125000002066 L-histidyl group Chemical group [H]N1C([H])=NC(C([H])([H])[C@](C(=O)[*])([H])N([H])[H])=C1[H] 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000037433 frameshift Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/38—Registration of image sequences
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
- A61B6/468—Arrangements for interfacing with the operator or the patient characterised by special input means allowing annotation or message recording
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Physiology (AREA)
- Quality & Reliability (AREA)
- Human Computer Interaction (AREA)
- Cardiology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
生体の動態を放射線撮影することにより得られた複数の動態画像のそれぞれから診断対象の機能に係る特徴量の時間変化を示すグラフを生成し、生成したグラフに基づいて前記特徴量の時間変化の周期を取得する周期取得手段と、
前記複数の動態画像から生成した複数の前記特徴量の時間変化を示すグラフの少なくとも一つの周期を調整することにより、前記複数の周期調整データとして、前記特徴量の時間変化の周期が等しい複数のグラフを生成する周期調整データ生成手段と、
前記複数の周期調整データの同一位相ごとの差分情報を生成する差分情報生成手段と、
前記差分情報を出力する出力手段と、
を備える。
前記周期調整データ生成手段は、開始タイミングの位相が一致した前記複数の周期調整データを生成し、
前記差分情報生成手段は、前記複数の周期調整データにおける同一タイミングの前記特徴量の差分値を算出する。
生体の動態を放射線撮影することにより得られた第一の動態画像を解析して第一の解析結果画像を生成し、前記生体の動態を放射線撮影することにより得られた第二の動態画像を解析して第二の解析結果画像を生成する解析手段と、
前記第一の解析結果画像と前記第二の解析結果画像との同一位相ごとの差分情報を生成する差分情報生成手段と、
前記差分情報を出力する出力手段と、
を備える。
前記第一の解析結果画像から診断対象の機能に係る第一の特徴量の時間変化の周期と、
前記第二の解析結果画像から前記診断対象の機能に係る第二の特徴量の時間変化の周期とを取得する周期取得手段と、
前記周期取得手段により取得された周期を調整して、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい、前記第一の解析結果画像と前記第二の解析結果画像とを生成する周期調整データ生成手段と、
を備える。
前記周期調整データ生成手段は、前記第一の解析結果画像と前記第二の解析結果画像との少なくとも一方においてフレーム画像を追加又は削除することにより、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい、前記第一の解析結果画像と前記第二の解析結果画像とを生成する。
前記周期調整データ生成手段は、前記第一の解析結果画像と前記第二の解析結果画像との少なくとも一方に追加するフレーム画像における各画素の信号値を、追加される解析結果画像の複数のフレーム画像の同じ位置の画素の信号値を用いて補間することにより算出する。
前記周期調整データ生成手段は、前記第一の解析結果画像と前記第二の解析結果画像との少なくとも一方から選択された複数のフレーム画像の信号値に基づいて補間画像を生成することにより、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい、前記第一の解析結果画像と前記第二の解析結果画像とを生成する。
前記周期調整データ生成手段は、開始タイミングの位相が一致した、前記第一の解析結果画像と前記第二の解析結果画像とを生成し、
前記差分情報生成手段は、前記第一の解析結果画像と前記第二の解析結果画像とにおける同一タイミングのフレーム画像の同一座標の画素同士の信号値の差分値を算出する。
前記周期取得手段は、前記第一の解析結果画像から前記診断対象の機能に係る特徴量の時間変化を示すグラフと、前記第二の解析結果画像から前記診断対象の機能に係る特徴量の時間変化を示すグラフとを生成し、生成したグラフに基づいて前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とを取得し、
前記周期調整データ生成手段は、前記第一の解析結果画像から生成した前記第一の特徴量の時間変化を示すグラフと、前記第二の解析結果画像から生成した前記第二の特徴量の時間変化を示すグラフとの少なくとも一方の周期を調整することにより、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい複数のグラフを生成する。
前記周期調整データ生成手段は、開始タイミングの位相が一致した、前記第一の解析結果画像と前記第二の解析結果画像とを生成し、
前記差分情報生成手段は、前記第一の解析結果画像と前記第二の解析結果画像とにおける同一タイミングの信号値の差分値を算出する。
前記第一の動態画像と前記第二の動態画像は、胸部の動態画像である。
前記第一の動態画像と前記第二の動態画像は、胸部の動態画像であり、
前記診断対象の機能は、肺野の血流機能又は換気機能であり、
前記第一の特徴量と前記第二の特徴量は、関心領域の画素の信号値である。
前記第一の動態画像と前記第二の動態画像は、胸部の動態画像であり、
前記診断対象の機能は、横隔膜の機能であり、
前記第一の特徴量と前記第二の特徴量は、横隔膜位置である。
前記出力手段は、前記差分情報を、前記第一の解析結果画像と前記第二の解析結果画像の少なくとも一方と合成して出力する。
コンピューターを、
生体の動態を放射線撮影することにより得られた第一の動態画像を解析して第一の解析結果画像を生成し、前記生体の動態を放射線撮影することにより得られた第二の動態画像を解析して第二の解析結果画像を生成する解析手段、
前記第一の解析結果画像と前記第二の解析結果画像との同一位相ごとの差分情報を生成する差分情報生成手段、
前記差分情報を出力する出力手段、
として機能させる。
前記コンピューターを、さらに、
前記第一の解析結果画像から診断対象の機能に係る第一の特徴量の時間変化の周期と、前記第二の解析結果画像から前記診断対象の機能に係る第二の特徴量の時間変化の周期とを取得する周期取得手段、
前記周期取得手段により取得された周期を調整して、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい、前記第一の解析結果画像と前記第二の解析結果画像とを生成する周期調整データ生成手段、
として機能させる。
前記第一の動態画像と前記第二の動態画像は、胸部の動態画像である。
前記第一の動態画像と前記第二の動態画像は、胸部の動態画像であり、
前記診断対象の機能は、肺野の血流機能又は換気機能であり、
前記第一の特徴量と前記第二の特徴量は、関心領域の画素の信号値である。
前記第一の動態画像と前記第二の動態画像は、胸部の動態画像であり、
前記診断対象の機能は、横隔膜の機能であり、
前記第一の特徴量と前記第二の特徴量は、横隔膜位置である。
前記出力手段は、前記差分情報を、前記第一の解析結果画像と前記第二の解析結果画像の少なくとも一方と合成して出力する。
〔動態解析システム100の構成〕
まず、第1の実施形態の構成を説明する。
図1に、第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に示す差分表示処理Aが実行される。
まず、受信した動態画像と比較する過去の動態画像が選択される(ステップS10)。
ステップS10においては、例えば、記憶部32に記憶されている被写体Mの過去の動態画像の一覧が表示部34に表示され、表示された動態画像の中からユーザーが所望する動態画像を操作部33により選択することとしてもよいし、記憶部32に記憶されている被写体Mの過去の動態画像のうち、検査日が直近の動態画像を制御部31が自動的に選択することとしてもよい。なお、受信した動態画像を現在の動態画像、現在の動態画像と比較する動態画像を過去の動態画像と呼ぶ。
例えば、ステップS11においては、まず、現在の動態画像を構成するフレーム画像と過去の動態画像を構成するフレーム画像のそれぞれから肺野領域の輪郭を検出する。例えば、各フレーム画像において、各画素の信号値(濃度値)のヒストグラムから判別分析によって閾値を求め、この閾値より高信号の領域を肺野領域候補として1次抽出する。1次抽出された肺野領域候補の境界付近でエッジ検出を行い、境界付近の小領域でエッジが最大となる点を境界に沿って抽出すれば肺野領域の境界を抽出することができる。次いで、現在の動態画像と過去の動態画像の全フレーム画像の中から一枚のフレーム画像を選択して基準画像とし、他のフレーム画像の肺野形状が基準画像の肺野形状に一致するようにワーピング処理を施す。基準画像は、現在の動態画像又は過去の動態画像の一フレーム目のフレーム画像としてもよいし、最大呼気位や最大吸気位のフレーム画像としてもよい。あるいは、基準となる肺野形状の画像を予め用意しておき、双方の動態画像の肺野形状を基準となる肺野形状に合わせることとしてもよい。これにより、動態画像を比較するときの肺野形状が常に統一されるので、ユーザーは、同じ環境下で動態画像の比較評価を行うことが可能となる。
次いで、ローパスフィルター処理後の2つの動態画像の各フレーム画像にROI(関心領域)を設定する。ここで、ROIの位置は、診断対象の機能の特徴が最もよく表れる位置に設定することが好ましい。診断対象の機能が換気である場合、ROIは心臓と背骨を除く肺野領域内に設定することが好ましい。ROIの設定は、画像処理により自動的に行ってもよいし、ユーザーによる操作部33の操作に応じて(手動で)行ってもよい。
次いで、現在の動態画像と過去の動態画像のそれぞれにおいて、各フレーム画像のROI内の各画素の信号値の代表値(例えば、平均値、中央値等)を算出し、算出した代表値を時系列に(フレーム画像順に)プロットすることにより信号値の時間変化を示す波形のグラフ(波形グラフと呼ぶ)を生成する。
そして生成した波形グラフから、現在の動態画像と過去の動態画像のそれぞれにおける信号値の時間変化の周期を取得する。周期は、波形グラフの極大点(又は極小点)から次の極大点(又は極小点)までの時間を周期として算出することができる。
次いで、フィルター処理後の2つの動態画像の各フレーム画像にROIを設定する。ここで、ROIの位置は、診断対象の機能の特徴が最もよく表れる位置に設定することが好ましい。診断対象の機能が血流である場合、心臓領域にROIを設定することが好ましい。ROIの設定は、画像処理により自動的に行ってもよいし、ユーザーによる操作部33の操作に応じて(手動で)行ってもよい。
次いで、現在の動態画像と過去の動態画像のそれぞれにおいて、各フレーム画像のROI内の各画素の信号値の代表値(例えば、平均値、中央値等)を算出し、算出した代表値を時系列に(フレーム番号順に)プロットすることにより信号値の時間変化を示す波形のグラフ(波形グラフ)を生成する。
そして生成した波形グラフから、現在の動態画像と過去の動態画像のそれぞれにおける信号値の時間変化の周期を取得する。
ステップS13では、図4に示すように、現在の動態画像から取得した周期と過去の動態画像から取得した周期が異なる場合に、周期の調整を行い、周期が等しい2つの周期調整データを生成する。本実施形態において、周期調整データとは、周期を合わせた複数の動態画像(ここでは、周期を合わせた現在の動態画像と過去の動態画像)のそれぞれである。
次いで、調整後の周期にそれぞれの動態画像の周期を合わせるために、それぞれの動態画像において周期毎に追加又は削除するフレーム画像数を決定する。追加又は削除するフレーム画像が0である場合は、その動態画像に対して周期を調整する処理は行わない。そして、調整が必要な動態画像の周期毎に、決定されたフレーム画像数のフレーム画像の追加又は削除を行って2つの動態画像の周期を合わせることにより、周期が等しい2つの動態画像、すなわち、2つの周期調整データが生成される。
ステップS15においては、例えば、いずれかの周期調整データ(例えば、過去の動態画像の周期調整データ)の各画素に、算出された差分値(符号及び絶対値)に応じた色を重畳することにより差分画像を生成する。
上記第1の実施形態においては、被写体の動態を撮影することにより得られた2つの動態画像を比較する場合について説明したが、上記差分画像表示処理Aは、動態画像を画素毎又は複数画素のブロック毎に動態解析を行うことにより取得された解析結果画像同士を比較する場合(例えば、制御部31が動態画像を解析して解析結果画像を生成する機能を備え、現在の動態画像を解析した解析結果画像と記憶部32に記憶されている同一患者の過去の解析結果画像を比較する場合)にも適用できる。すなわち、比較する2つの解析結果画像に対して上記差分表示処理AのステップS11〜S16と同様の処理を行うことで、それらの差異を可視化してユーザーが容易に把握しやすくすることができる。
血流信号波形は、一連のフレーム画像の各小領域毎に、時間方向のハイパスフィルター処理(例えば、低域カットオフ周波数0.8Hz)を施してから小領域内の各画素の信号値の代表値(平均値、最大値等)を算出し、算出した代表値の時間変化を示す波形を取得することにより求めることができる。
拍動信号波形としては、以下のいずれかを用いることができる。
(a)心臓領域(又は大動脈領域)にROI(関心領域)を定め、そのROIにおける信号値の時間変化を示す波形
(b)(a)の波形を反転させた信号波形
(c)心電検知センサーより得られた心電信号波形
(d)心壁の動き(位置の変化)を示す信号波形
また、相互相関係数は、以下の[数1]により求めることができる。
次に、本発明の第2の実施形態について説明する。
第2の実施形態における構成は、診断用コンソール3の記憶部32に差分表示処理Bを実行するためのプログラムが記憶されている他は、第1の実施形態で説明したものと同様であるので説明を省略し、以下第2の実施形態の動作について説明する。
図6は、第2の実施形態において、診断用コンソール3により実行される差分表示処理Bを示すフローチャートである。差分表示処理Bは、制御部31と記憶部32に記憶されているプログラムとの協働により実行される。
まず、受信した動態画像と比較する過去の動態画像が選択される(ステップS20)。
ステップS20の処理は、図3のステップS10の処理と同様であるので説明を援用する。
ステップS21の処理は、図3のステップS12の処理と同様であるので説明を援用する。なお、波形グラフにプロットされる各データを信号値データと呼ぶ。
ステップS22では、図4に示すように、現在の動態画像に基づいて生成した波形グラフの周期と過去の動態画像に基づいて生成した波形グラフの周期が異なる場合に、周期の調整を行い、周期を統一する。本実施形態において、周期調整データとは、周期を合わせた複数の波形グラフのそれぞれである。
次いで、調整後の周期にそれぞれの波形グラフの周期を合わせるために、それぞれの波形グラフにおいて周期毎に追加又は削除する信号値データの数を決定する。追加又は削除する信号値データの数が0である場合は、その波形グラフに対して周期を調整する処理は行わない。そして、調整が必要な波形グラフの周期毎に、決定された数の信号値データの追加又は削除を行って2つの波形グラフの周期を合わせることにより、周期が等しい2つの波形グラフ、すなわち、2つの周期調整データが生成される。
ステップS24においては、例えば、2つの周期調整データ(グラフ)を一つにまとめ、そのグラフ上に差分値の波形を追加し、差分値の波形に差分値に応じた色を付して差分グラフを生成する。
上記第2の実施形態においては、被写体の動態を撮影することにより得られた2つの動態画像を比較する場合について説明したが、上記差分画像表示処理Bは、動態画像を画素毎又は複数画素のブロック毎に動態解析を行うことにより取得された解析結果画像同士を比較する場合(例えば、制御部31が動態画像を解析して解析結果画像を生成する機能を備え、現在の動態画像を解析した解析結果画像と記憶部32に記憶されている同一患者の過去の解析結果画像を比較する場合)にも適用できる。すなわち、比較する2つの解析結果画像に対して上記差分表示処理BのステップS21〜S25と同様の処理を行うことで、それらの差異を可視化してユーザーが容易に把握しやすくすることができる。
解析結果画像の例は、第1の実施形態で説明したものと同様であるので説明を援用する。
次に、本発明の第3の実施形態について説明する。
第3の実施形態における構成は、診断用コンソール3の記憶部32に差分表示処理Cを実行するためのプログラムが記憶されている他は、第1の実施形態で説明したものと同様であるので説明を省略し、以下第3の実施形態の動作について説明する。
図8は、第3の実施形態において、診断用コンソール3により実行される差分表示処理Cを示すフローチャートである。差分表示処理Cは、制御部31と記憶部32に記憶されているプログラムとの協働により実行される。
まず、受信した動態画像と比較する過去の動態画像が選択される(ステップS30)。
ステップS30の処理は、図3のステップS10の処理と同様であるので説明を援用する。
例えば、ステップS31においては、まず、現在の動態画像を構成するフレーム画像と過去の動態画像を構成するフレーム画像のそれぞれから鎖骨や胸郭を抽出する。鎖骨や胸郭の抽出は、例えば、各フレーム画像において、予め用意した鎖骨テンプレート、肋骨テンプレート、胸骨のテンプレートを用いてテンプレートマッチングを行う手法や、エッジ検出後にカーブフィッティング関数を当てはめる等の手法により行うことができる。また、鎖骨や胸郭の骨の構造の前知識に基づく位置、形状、サイズ、濃度勾配、方向等の特徴をもとに、抽出した領域が鎖骨や胸骨か否かの精査を行い、過剰抽出されている部分を判別して取り除くこととしてもよい。次いで、現在の動態画像と過去の動態画像の全フレーム画像の中から一枚のフレーム画像を選択して基準画像とし、他のフレーム画像の鎖骨や胸郭が基準画像の鎖骨や胸骨に一致するようにワーピング処理を施す。基準画像は、現在の動態画像又は過去の動態画像の一フレーム目のフレーム画像としてもよいし、最大呼気位や最大吸気位のフレーム画像としてもよい。
ステップS32においては、まず、現在の動態画像と過去の動態画像のそれぞれについで、各フレーム画像から横隔膜の位置を特定する。胸部正面の単純X線画像において、横隔膜は肺野下部に接して見える。そこで、例えば、各フレーム画像において、肺野領域を抽出し、抽出した肺野領域の下部の輪郭を横隔膜の位置として特定することができる。
次いで、現在の動態画像と過去の動態画像のそれぞれにおいて、各フレーム画像のから特定した横隔膜位置のy座標の代表値(例えば、平均値、中央値等)を算出し、算出した代表値を時系列に(フレーム画像順に)プロットすることにより横隔膜位置の時間変化を示す波形のグラフ(波形グラフと呼ぶ)を生成する。なお、画像の水平方向をx方向、垂直方向をy方向とする。また、画像の左上を原点とし、y方向が下に行くにつれて座標値が大きくなるものとする。
そして、生成した波形グラフから、現在の動態画像と過去の動態画像のそれぞれにおける横隔膜位置の時間変化の周期を取得する。周期は、横隔膜位置の時間変化を示す波形の極大点(又は極小点)から次の極大点(又は極小点)までの時間を周期として算出することができる。
ステップS33では、図4に示すように、現在の動態画像から取得した周期と過去の動態画像から取得した周期が異なる場合に、周期の調整を行い、周期が等しい2つの周期調整データを生成する。本実施形態において、周期調整データとは、周期を合わせた複数の動態画像(ここでは、周期を合わせた現在の動態画像と過去の動態画像)のそれぞれである。
次いで、調整後の周期にそれぞれの動態画像の周期を合わせるために、それぞれの動態画像において周期毎に追加又は削除するフレーム画像数を決定する。追加又は削除するフレーム画像が0である場合は、その動態画像に対して周期を調整する処理は行わない。そして、調整が必要な動態画像の周期毎に、決定されたフレーム画像数のフレーム画像の追加又は削除を行って2つの動態画像の周期を合わせることにより、周期が等しい2つの動態画像、すなわち、2つの周期調整データが生成される。
ステップS35においては、例えば、いずれかの周期調整データ(例えば、過去の動態画像の周期調整データ)の各フレーム画像において、現在と過去の横隔膜位置及びその間の領域(差分領域)にそれぞれ異なる色を重畳することにより差分画像を生成する。差分領域には、例えば、差分値の符号に応じた色を重畳する。
1 撮影装置
11 放射線源
12 放射線照射制御装置
13 放射線検出部
14 読取制御装置
2 撮影用コンソール
21 制御部
22 記憶部
23 操作部
24 表示部
25 通信部
26 バス
3 診断用コンソール
31 制御部
32 記憶部
33 操作部
34 表示部
35 通信部
36 バス
Claims (20)
- 生体の動態を放射線撮影することにより得られた複数の動態画像のそれぞれから診断対象の機能に係る特徴量の時間変化を示すグラフを生成し、生成したグラフに基づいて前記特徴量の時間変化の周期を取得する周期取得手段と、
前記複数の動態画像から生成した複数の前記特徴量の時間変化を示すグラフの少なくとも一つの周期を調整することにより、前記複数の周期調整データとして、前記特徴量の時間変化の周期が等しい複数のグラフを生成する周期調整データ生成手段と、
前記複数の周期調整データの同一位相ごとの差分情報を生成する差分情報生成手段と、
前記差分情報を出力する出力手段と、
を備える動態解析システム。 - 前記周期調整データ生成手段は、開始タイミングの位相が一致した前記複数の周期調整データを生成し、
前記差分情報生成手段は、前記複数の周期調整データにおける同一タイミングの前記特徴量の差分値を算出する請求項1に記載の動態解析システム。 - 生体の動態を放射線撮影することにより得られた第一の動態画像を解析して第一の解析結果画像を生成し、前記生体の動態を放射線撮影することにより得られた第二の動態画像を解析して第二の解析結果画像を生成する解析手段と、
前記第一の解析結果画像と前記第二の解析結果画像との同一位相ごとの差分情報を生成する差分情報生成手段と、
前記差分情報を出力する出力手段と、
を備える動態解析システム。 - 前記第一の解析結果画像から診断対象の機能に係る第一の特徴量の時間変化の周期と、前記第二の解析結果画像から前記診断対象の機能に係る第二の特徴量の時間変化の周期とを取得する周期取得手段と、
前記周期取得手段により取得された周期を調整して、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい、前記第一の解析結果画像と前記第二の解析結果画像とを生成する周期調整データ生成手段と、
を備える請求項3に記載の動態解析システム。 - 前記周期調整データ生成手段は、前記第一の解析結果画像と前記第二の解析結果画像との少なくとも一方においてフレーム画像を追加又は削除することにより、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい、前記第一の解析結果画像と前記第二の解析結果画像とを生成する請求項4に記載の動態解析システム。
- 前記周期調整データ生成手段は、前記第一の解析結果画像と前記第二の解析結果画像との少なくとも一方に追加するフレーム画像における各画素の信号値を、追加される解析結果画像の複数のフレーム画像の同じ位置の画素の信号値を用いて補間することにより算出する請求項5に記載の動態解析システム。
- 前記周期調整データ生成手段は、前記第一の解析結果画像と前記第二の解析結果画像との少なくとも一方から選択された複数のフレーム画像の信号値に基づいて補間画像を生成することにより、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい、前記第一の解析結果画像と前記第二の解析結果画像とを生成する請求項4に記載の動態解析システム。
- 前記周期調整データ生成手段は、開始タイミングの位相が一致した、前記第一の解析結果画像と前記第二の解析結果画像とを生成し、
前記差分情報生成手段は、前記第一の解析結果画像と前記第二の解析結果画像とにおける同一タイミングのフレーム画像の同一座標の画素同士の信号値の差分値を算出する請求項5〜7のいずれか一項に記載の動態解析システム。 - 前記周期取得手段は、前記第一の解析結果画像から前記診断対象の機能に係る特徴量の時間変化を示すグラフと、前記第二の解析結果画像から前記診断対象の機能に係る特徴量の時間変化を示すグラフとを生成し、生成したグラフに基づいて前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とを取得し、
前記周期調整データ生成手段は、前記第一の解析結果画像から生成した前記第一の特徴量の時間変化を示すグラフと、前記第二の解析結果画像から生成した前記第二の特徴量の時間変化を示すグラフとの少なくとも一方の周期を調整することにより、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい複数のグラフを生成する請求項4に記載の動態解析システム。 - 前記周期調整データ生成手段は、開始タイミングの位相が一致した、前記第一の解析結果画像と前記第二の解析結果画像とを生成し、
前記差分情報生成手段は、前記第一の解析結果画像と前記第二の解析結果画像とにおける同一タイミングの信号値の差分値を算出する請求項9に記載の動態解析システム。 - 前記第一の動態画像と前記第二の動態画像は、胸部の動態画像である請求項3〜10のいずれか一項に記載の動態解析システム。
- 前記第一の動態画像と前記第二の動態画像は、胸部の動態画像であり、
前記診断対象の機能は、肺野の血流機能又は換気機能であり、
前記第一の特徴量と前記第二の特徴量は、関心領域の画素の信号値である請求項4〜10のいずれか一項に記載の動態解析システム。 - 前記第一の動態画像と前記第二の動態画像は、胸部の動態画像であり、
前記診断対象の機能は、横隔膜の機能であり、
前記第一の特徴量と前記第二の特徴量は、横隔膜位置である請求項4〜10のいずれか一項に記載の動態解析システム。 - 前記出力手段は、前記差分情報を、前記第一の解析結果画像と前記第二の解析結果画像の少なくとも一方と合成して出力する請求項3〜13のいずれか一項に記載の動態解析システム。
- コンピューターを、
生体の動態を放射線撮影することにより得られた第一の動態画像を解析して第一の解析結果画像を生成し、前記生体の動態を放射線撮影することにより得られた第二の動態画像を解析して第二の解析結果画像を生成する解析手段、
前記第一の解析結果画像と前記第二の解析結果画像との同一位相ごとの差分情報を生成する差分情報生成手段、
前記差分情報を出力する出力手段、
として機能させるためのプログラム。 - 前記コンピューターを、さらに、
前記第一の解析結果画像から診断対象の機能に係る第一の特徴量の時間変化の周期と、前記第二の解析結果画像から前記診断対象の機能に係る第二の特徴量の時間変化の周期とを取得する周期取得手段、
前記周期取得手段により取得された周期を調整して、前記第一の特徴量の時間変化の周期と前記第二の特徴量の時間変化の周期とが等しい、前記第一の解析結果画像と前記第二の解析結果画像とを生成する周期調整データ生成手段、
として機能させる請求項15に記載のプログラム。 - 前記第一の動態画像と前記第二の動態画像は、胸部の動態画像である請求項15又は16に記載のプログラム。
- 前記第一の動態画像と前記第二の動態画像は、胸部の動態画像であり、
前記診断対象の機能は、肺野の血流機能又は換気機能であり、
前記第一の特徴量と前記第二の特徴量は、関心領域の画素の信号値である請求項16に記載のプログラム。 - 前記第一の動態画像と前記第二の動態画像は、胸部の動態画像であり、
前記診断対象の機能は、横隔膜の機能であり、
前記第一の特徴量と前記第二の特徴量は、横隔膜位置である請求項16に記載のプログラム。 - 前記出力手段は、前記差分情報を、前記第一の解析結果画像と前記第二の解析結果画像の少なくとも一方と合成して出力する請求項15〜19のいずれか一項に記載のプログラム。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017045532A JP6772908B2 (ja) | 2017-03-10 | 2017-03-10 | 動態解析システム及びプログラム |
US15/915,801 US10540767B2 (en) | 2017-03-10 | 2018-03-08 | Dynamic analysis system |
US16/706,144 US10810741B2 (en) | 2017-03-10 | 2019-12-06 | Dynamic analysis system |
US17/015,951 US11410312B2 (en) | 2017-03-10 | 2020-09-09 | Dynamic analysis system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017045532A JP6772908B2 (ja) | 2017-03-10 | 2017-03-10 | 動態解析システム及びプログラム |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020164856A Division JP2020203191A (ja) | 2020-09-30 | 2020-09-30 | 動態解析システム及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018148964A JP2018148964A (ja) | 2018-09-27 |
JP6772908B2 true JP6772908B2 (ja) | 2020-10-21 |
Family
ID=63444838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017045532A Active JP6772908B2 (ja) | 2017-03-10 | 2017-03-10 | 動態解析システム及びプログラム |
Country Status (2)
Country | Link |
---|---|
US (3) | US10540767B2 (ja) |
JP (1) | JP6772908B2 (ja) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6772908B2 (ja) * | 2017-03-10 | 2020-10-21 | コニカミノルタ株式会社 | 動態解析システム及びプログラム |
JP6950483B2 (ja) * | 2017-11-20 | 2021-10-13 | コニカミノルタ株式会社 | 動態撮影システム |
JP7147507B2 (ja) * | 2018-11-22 | 2022-10-05 | コニカミノルタ株式会社 | 画像処理装置及びプログラム |
JP7188051B2 (ja) * | 2018-12-18 | 2022-12-13 | コニカミノルタ株式会社 | 医用画像管理システム |
JP2020141841A (ja) * | 2019-03-06 | 2020-09-10 | コニカミノルタ株式会社 | 動態解析装置及びプログラム |
JP7211172B2 (ja) | 2019-03-08 | 2023-01-24 | コニカミノルタ株式会社 | 動態画像解析システム及び動態画像処理装置 |
CA3149577A1 (en) * | 2019-08-27 | 2021-03-04 | 4DMedical Limited | Methods, systems, and apparatus for assessing an effect of a medical treatment on organ function |
NL2023710B1 (en) * | 2019-08-27 | 2021-05-11 | Thirona B V | Computer implemented method for estimating lung perfusion from lung images |
JP7435242B2 (ja) * | 2020-05-15 | 2024-02-21 | コニカミノルタ株式会社 | 動態画像解析装置、動態画像解析方法及びプログラム |
JP2020203191A (ja) * | 2020-09-30 | 2020-12-24 | コニカミノルタ株式会社 | 動態解析システム及びプログラム |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US540767A (en) * | 1895-06-11 | Combination can-opener | ||
JP4938480B2 (ja) * | 2007-02-02 | 2012-05-23 | 富士フイルム株式会社 | 放射線撮像装置、放射線撮像方法、及びプログラム |
WO2009093693A1 (ja) * | 2008-01-25 | 2009-07-30 | Konica Minolta Holdings, Inc. | 画像生成装置、画像生成方法およびプログラム |
CN104703539B (zh) | 2012-10-04 | 2018-04-10 | 柯尼卡美能达株式会社 | 图像处理装置和程序 |
JP6772908B2 (ja) * | 2017-03-10 | 2020-10-21 | コニカミノルタ株式会社 | 動態解析システム及びプログラム |
-
2017
- 2017-03-10 JP JP2017045532A patent/JP6772908B2/ja active Active
-
2018
- 2018-03-08 US US15/915,801 patent/US10540767B2/en active Active
-
2019
- 2019-12-06 US US16/706,144 patent/US10810741B2/en active Active
-
2020
- 2020-09-09 US US17/015,951 patent/US11410312B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US20200410684A1 (en) | 2020-12-31 |
US11410312B2 (en) | 2022-08-09 |
US20180260955A1 (en) | 2018-09-13 |
US20200118270A1 (en) | 2020-04-16 |
US10540767B2 (en) | 2020-01-21 |
JP2018148964A (ja) | 2018-09-27 |
US10810741B2 (en) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6772908B2 (ja) | 動態解析システム及びプログラム | |
JP6772873B2 (ja) | 動態解析装置及び動態解析システム | |
JP6743662B2 (ja) | 動態画像処理システム | |
JP6958202B2 (ja) | 動態画像処理装置及びプログラム | |
JP6217241B2 (ja) | 胸部診断支援システム | |
JP6812685B2 (ja) | 動態解析装置 | |
JP6690774B2 (ja) | 動態解析システム、プログラム及び動態解析装置 | |
JP6740910B2 (ja) | 動態画像処理システム | |
JP7020125B2 (ja) | 動態画像解析装置及びプログラム | |
JP6743730B2 (ja) | 動態解析システム | |
JP6848393B2 (ja) | 動態画像処理装置 | |
JP2017169830A (ja) | 動態解析装置 | |
JP2021194140A (ja) | 画像処理装置及び画像処理方法 | |
JP6950507B2 (ja) | 動態画像処理装置 | |
JP2018187310A (ja) | 動態画像処理システム | |
JP2018196693A (ja) | 動態解析システム | |
JP7099086B2 (ja) | 動態画像処理装置及びプログラム | |
JP2019005417A (ja) | 動態画像処理装置及び動態画像処理システム | |
JP6962030B2 (ja) | 動態解析装置、動態解析システム、動態解析プログラム及び動態解析方法 | |
JP2020203191A (ja) | 動態解析システム及びプログラム | |
JP6950483B2 (ja) | 動態撮影システム | |
JP6888721B2 (ja) | 動態画像処理装置、動態画像処理プログラム及び動態画像処理方法 | |
JP2020168173A (ja) | 動態解析装置、動態解析システム及びプログラム | |
JP6874484B2 (ja) | 動態画像処理システム | |
JP2017217047A (ja) | 画像表示システム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190620 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200417 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200616 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200812 |
|
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: 20200901 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200914 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6772908 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |