JP4797173B2 - X-ray diagnosis support apparatus, program, and recording medium - Google Patents
X-ray diagnosis support apparatus, program, and recording medium Download PDFInfo
- Publication number
- JP4797173B2 JP4797173B2 JP2007522243A JP2007522243A JP4797173B2 JP 4797173 B2 JP4797173 B2 JP 4797173B2 JP 2007522243 A JP2007522243 A JP 2007522243A JP 2007522243 A JP2007522243 A JP 2007522243A JP 4797173 B2 JP4797173 B2 JP 4797173B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- moving image
- ray
- chest
- processing unit
- 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
- 238000003745 diagnosis Methods 0.000 title claims description 47
- 238000012545 processing Methods 0.000 claims description 67
- 238000009423 ventilation Methods 0.000 claims description 36
- 238000000034 method Methods 0.000 claims description 18
- 230000010354 integration Effects 0.000 claims description 17
- 238000011976 chest X-ray Methods 0.000 claims description 14
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 11
- 230000033001 locomotion Effects 0.000 claims description 9
- 238000002591 computed tomography Methods 0.000 description 75
- 210000000038 chest Anatomy 0.000 description 34
- 210000004072 lung Anatomy 0.000 description 29
- 230000006870 function Effects 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 7
- 230000000241 respiratory effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000000877 morphologic effect Effects 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 3
- 210000002216 heart Anatomy 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000009613 pulmonary function test Methods 0.000 description 3
- 238000011158 quantitative evaluation Methods 0.000 description 3
- 208000023504 respiratory system disease Diseases 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 208000006545 Chronic Obstructive Pulmonary Disease Diseases 0.000 description 2
- 210000000481 breast Anatomy 0.000 description 2
- 210000000621 bronchi Anatomy 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 230000002685 pulmonary effect Effects 0.000 description 2
- 230000009325 pulmonary function Effects 0.000 description 2
- 206010006458 Bronchitis chronic Diseases 0.000 description 1
- 206010014561 Emphysema Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000003915 air pollution Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000008827 biological function Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 206010006451 bronchitis Diseases 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 208000007451 chronic bronchitis Diseases 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 210000004884 grey matter Anatomy 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000004199 lung function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 230000036412 respiratory physiology Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000000391 smoking effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- 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
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/24—Indexing scheme for image data processing or generation, in general involving graphical user interfaces [GUIs]
-
- 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/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biophysics (AREA)
- Pulmonology (AREA)
- Physiology (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
本発明は、CT(Computed Tomography)画像及びX線動画像を用いてコンピュータ解析を行い、動態画像診断のための情報を得るX線診断支援技術に関する。 The present invention relates to an X-ray diagnosis support technique for performing computer analysis using a CT (Computed Tomography) image and an X-ray moving image to obtain information for dynamic image diagnosis.
従来、呼吸器系疾患の罹患率が世界的に上昇する傾向にあり、特に、慢性閉塞性肺疾患(COPD)の罹患率が高い状況にある。この慢性閉塞性肺疾患は、大気汚染や喫煙が主な原因の肺気腫や慢性気管支炎等の疾患であり、高齢者に多く発症している。このため、呼吸器系疾患の罹患率の上昇に伴い、胸部の疾患に関する定量的な機能評価を簡便に繰り返して行う方法が必要とされている。 Conventionally, the prevalence of respiratory diseases has been increasing worldwide, and in particular, the prevalence of chronic obstructive pulmonary disease (COPD) is high. This chronic obstructive pulmonary disease is a disease such as pulmonary emphysema and chronic bronchitis mainly caused by air pollution and smoking, and occurs frequently in the elderly. For this reason, with an increase in the prevalence of respiratory diseases, a method for simply and repeatedly performing quantitative functional evaluation on chest diseases is required.
一般に、呼吸器系疾患の機能評価は、肺機能検査により行われる。この肺機能検査は、体内に酸素を取り入れ、体外に二酸化炭素を排出する基本的な機能を評価するために、換気、ガスと血流の分布、拡散等が正常に行われているか否かを、定性的及び定量的に検査するものである。しかしながら、肺機能検査では、局所的にまたは左右の肺を別々に検査することができず、最大努力呼吸でないと正確な検査ができないという問題があった。また、受検者には、検査のために多大な体力が必要になるという負担もあった。 In general, functional evaluation of respiratory diseases is performed by a pulmonary function test. This pulmonary function test determines whether ventilation, gas and blood flow distribution, diffusion, etc. are performed normally in order to evaluate the basic function of taking oxygen into the body and discharging carbon dioxide outside the body. Qualitative and quantitative inspection. However, in the pulmonary function test, there is a problem that the left and right lungs cannot be examined locally or separately, and accurate examination cannot be performed unless maximum effort breathing is performed. The examinee also had a burden that a great deal of physical strength was required for the examination.
このような肺機能検査の問題を解決し、受検者に与える負担を軽減するために、肺機能を局所的に定量評価するX線診断装置が開発されている(例えば、特許文献1を参照)。このX線診断装置は、胸部をX線撮影して得た静止画像を記憶手段に記憶し、この静止画像における胸部の任意の部位、濃度値の計算方法、及び結果の表示方法が設定されると、その設定された部位のピクセル値(濃度値)に基づいて濃度に関する解析処理を行い、その解析処理結果を画面に表示する。また、X線診断装置は、胸部をX線撮影して得た時系列のフレームの静止画像を記憶手段に記憶し、この静止画像における胸部の任意の部位、濃度値の計算方法、解析対象のフレーム、及び結果の表示方法が設定されると、その設定された部位のピクセル値に基づいて、フレーム間の濃度に関する解析処理を時系列に行い、その解析処理結果を画面に表示する。このような静止画像の解析処理、及び動画像の解析処理により、読影や診断のために必要な情報を表示するようにしたから、人間の目に頼ることのない診断を実現することができ、読影医によって診断結果がばらつく可能性が低くなる。また、X線診断により得た画像に関する様々な処理技術が開示されている(例えば、特許文献2〜4を参照。) In order to solve such a problem of pulmonary function testing and reduce the burden on the examinee, an X-ray diagnostic apparatus for locally quantitatively evaluating pulmonary function has been developed (for example, see Patent Document 1). . In this X-ray diagnostic apparatus, a still image obtained by X-ray imaging of the chest is stored in a storage unit, and an arbitrary part of the chest in this still image, a density value calculation method, and a result display method are set. Then, an analysis process related to the density is performed based on the pixel value (density value) of the set part, and the analysis process result is displayed on the screen. In addition, the X-ray diagnostic apparatus stores a still image of a time-series frame obtained by X-ray imaging of the chest in a storage unit, an arbitrary part of the chest in this still image, a method for calculating a density value, and an analysis target When the frame and the result display method are set, the analysis processing regarding the density between frames is performed in time series based on the pixel value of the set portion, and the analysis processing result is displayed on the screen. Since the information necessary for interpretation and diagnosis is displayed by such still image analysis processing and moving image analysis processing, it is possible to realize diagnosis that does not rely on human eyes, The possibility that the diagnostic results vary by the interpretation doctor is reduced. Various processing techniques relating to images obtained by X-ray diagnosis are disclosed (see, for example, Patent Documents 2 to 4).
ところで、従来、X線照射によって被検体の断層像(CT画像)を得るX線CT装置が利用されている。このX線CT装置は、被検体にX線を照射し、臓器、血液、灰白質等の人体組織のX線吸収率の差を検出器により検出し、これをコンピュータ処理して再構成することにより、検査部位の断層面(スライス面)の画像を得るものである。このようにして得られたCT画像に基づいて、医師は患者の病状等を精度高く診断することができる。 Conventionally, an X-ray CT apparatus that obtains a tomographic image (CT image) of a subject by X-ray irradiation has been used. This X-ray CT apparatus irradiates a subject with X-rays, detects differences in X-ray absorption rates of human tissues such as organs, blood, and gray matter with a detector, and reconfigures them by computer processing. Thus, an image of the tomographic plane (slice plane) of the examination site is obtained. Based on the CT image thus obtained, a doctor can diagnose a patient's medical condition or the like with high accuracy.
また、従来、被検体の胸部をスクリーニング検査する装置の開発が行われている。例えば、FPD(Flat Panel Detector)は、X線を電気信号に変換することにより画像を得る機器であり、画像を直接デジタル化する平面検出器を使用した撮影機器である。このFPDを使用することにより、画像の確認、濃度調整及び幾何学的調整を行い、安定した動画像を得ることができる。例えば、吸気から呼気までの一連の呼吸動作の動画像を得て、肺尖部から横隔膜までの距離を計測して動画像を解析処理し、その結果を読影や診断のために必要な情報として表示することができる。しかしながら、FPDを使用して得た動画像を読影するには、呼吸生理学の知識が必要である。 Conventionally, an apparatus for screening examination of a subject's chest has been developed. For example, an FPD (Flat Panel Detector) is a device that obtains an image by converting X-rays into an electrical signal, and is an imaging device that uses a flat panel detector that directly digitizes the image. By using this FPD, image confirmation, density adjustment and geometric adjustment can be performed to obtain a stable moving image. For example, a moving image of a series of breathing movements from inspiration to expiration is obtained, the distance from the lung apex to the diaphragm is measured, the moving image is analyzed, and the results are used as information necessary for interpretation and diagnosis. Can be displayed. However, in order to interpret moving images obtained using FPD, knowledge of respiratory physiology is required.
〔特許文献1〕特開平7−194583号公報
〔特許文献2〕特開平5−7579号公報
〔特許文献3〕特開平4−134568号公報
〔特許文献4〕特開昭64−17631号公報[Patent Document 1] Japanese Patent Laid-Open No. 7-194583 [Patent Document 2] Japanese Patent Laid-Open No. 5-7579 [Patent Document 3] Japanese Patent Laid-Open No. 4-134568 [Patent Document 4] Japanese Patent Laid-Open No. 64-17631
このような状況の下で、X線CT装置により得られたCT画像や動画像の形態情報と、動画像を解析処理して得られた動態情報とを関連付け、例えば胸部に関する形態情報と呼吸性動態情報とを重ね合わせた新たな情報とし、この重ね合わせた情報を読影や診断のために利用することが望まれている。 Under such circumstances, the CT image or moving image morphological information obtained by the X-ray CT apparatus is associated with the dynamic information obtained by analyzing the moving image. It is desired that new information is superimposed on the dynamic information, and the superimposed information is used for interpretation and diagnosis.
そこで、本発明は、このような状況を鑑みてなされたものであり、その目的は、被検体である肺等の機能について局所的な定量評価を行い、読影や診断のために有効に利用できる評価結果を得ることが可能なX線診断支援装置、プログラム及び記憶媒体を提供することにある。 Therefore, the present invention has been made in view of such a situation, and the purpose thereof is to perform local quantitative evaluation on the function of the subject lung or the like, and can be used effectively for interpretation and diagnosis. An object is to provide an X-ray diagnosis support apparatus, a program, and a storage medium that can obtain an evaluation result.
本発明によるX線診断支援装置は、診断部位の一連の動作のX線動画像が格納された動画像格納部と、前記診断部位のX線静止画像が格納された静止画像格納部と、前記動画像格納部から動画像を読み出し、該動画像に基づいて前記診断部位に関する機能を示す動態情報を生成する動画像処理部と、前記静止画像格納部から静止画像を読み出し、該静止画像から前記診断部位に関する新たな画像を生成する静止画像処理部と、前記動画像処理部により生成された動態情報と、前記動画像及び前記静止画像処理部により生成された新たな画像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示する統合部とを備えたことを特徴とする。 An X-ray diagnosis support apparatus according to the present invention includes a moving image storage unit that stores an X-ray moving image of a series of operations of a diagnostic region, a still image storage unit that stores an X-ray still image of the diagnostic region, A moving image is read from the moving image storage unit, a moving image processing unit that generates dynamic information indicating a function related to the diagnosis part based on the moving image, a still image is read from the still image storage unit, and the moving image is read from the still image At least one of a still image processing unit that generates a new image relating to a diagnostic region, dynamic information generated by the moving image processing unit, and a new image generated by the moving image and the still image processing unit And an integration unit that integrates the form information and displays the integrated information on a screen.
ここで、形態情報とは、例えば、血管や気管支の太さや分布、心臓の大きさや位置、横隔膜の位置や形、骨の浸潤や骨折、肺の性質、異常物(癌、水、異物、カビ、空気の塊等)の存在を認識できる静止画像または動画像に関する情報をいう。また、動態情報とは、動画像を解析処理して得られる生態の機能情報であり、例えば、血管や気管支の拡散または収縮、心臓の拍動、横隔膜の動きまたは可動領域、肺の換気、異常物の可動性に関する情報をいう。この動態情報には、呼吸性動態情報及び心拍動静動態情報等がある。 Here, the morphological information includes, for example, the thickness and distribution of blood vessels and bronchi, the size and position of the heart, the position and shape of the diaphragm, bone infiltration and fracture, the nature of the lungs, abnormal substances (cancer, water, foreign matter, mold, etc. , Information on still images or moving images capable of recognizing the presence of air masses. The dynamic information is biological function information obtained by analyzing and processing moving images, for example, diffusion or contraction of blood vessels or bronchi, heart pulsation, diaphragm movement or movable region, lung ventilation, abnormalities, etc. Information about the mobility of things. This dynamic information includes respiratory dynamic information and heartbeat dynamic information.
また、本発明によるX線診断支援装置は、前記動画像格納部に格納されたX線動画像を、胸部の後面から撮影し、吸気から呼気までの一連の呼吸動作を示す胸部X線動画像とし、前記静止画像格納部に格納されたX線静止画像を、胸部のスライス面を示すCT画像とし、前記動画像処理部が、胸部X線動画像の各フレームについて横隔膜の位置を計測し、該横隔膜の位置から最大吸気のフレーム及び最大呼気のフレームを特定し、予め設定された胸部のエリア毎に、前記最大吸気のフレームと最大呼気のフレームとの間のピクセル差分値を算出して換気情報とし、該換気情報を動態情報として生成し、前記静止画像処理部が、CT画像に基づいて胸部の後面から見たレイサム画像、コロナル像、及びサジタル像のうちの少なくとも一つを新たな画像として生成することを特徴とする。 Further, the X-ray diagnosis support apparatus according to the present invention captures an X-ray moving image stored in the moving image storage unit from the rear surface of the chest, and shows a series of breathing operations from inspiration to expiration. The X-ray still image stored in the still image storage unit is a CT image showing a slice plane of the chest, and the moving image processing unit measures the position of the diaphragm for each frame of the chest X-ray moving image, A frame of maximum inspiration and a frame of maximum exhalation are specified from the position of the diaphragm, and a pixel difference value between the frame of maximum inspiration and the frame of maximum exhalation is calculated for each predetermined chest area, and ventilation is performed. Information, and the ventilation information is generated as dynamic information, and the still image processing unit newly generates at least one of a laysum image, a coronal image, and a sagittal image viewed from the back of the chest based on the CT image. And generating an image.
また、本発明によるX線診断支援装置は、前記動画像格納部に格納されたX線動画像を、胸部の後面から撮影し、吸気から呼気までの一連の呼吸動作を示す胸部X線動画像とし、前記静止画像格納部に格納されたX線静止画像を、胸部のスライス面を示すCT画像とし、前記動画像処理部が、胸部X線動画像の各フレームについて横隔膜の位置を計測し、予め設定された胸部のエリア毎に、胸部X線動画像の各フレームについて、時系列の前後におけるフレーム間のピクセル差分値を算出して換気情報とし、該換気情報を動態情報として生成し、前記静止画像処理部が、CT画像に基づいて胸部の後面から見たレイサム画像、コロナル像、及びサジタル像のうちの少なくとも前記レイサム画像を含む新たな画像を生成し、該生成したレイサム画像について横隔膜の位置を計測し、前記統合部が、静止画像処理部により計測された横隔膜の位置と、動画像処理部により各フレームについて計測された横隔膜の位置とを比較し、両位置が一致する動画像のフレームを特定し、動画像処理部により生成された動態情報と、前記動画像、前記特定した動画像のフレーム、静止画像処理部により生成されたレイサム画像、コロナル像、及びサジタル像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示することを特徴とする。 Further, the X-ray diagnosis support apparatus according to the present invention captures an X-ray moving image stored in the moving image storage unit from the rear surface of the chest, and shows a series of breathing operations from inspiration to expiration. The X-ray still image stored in the still image storage unit is a CT image showing a slice plane of the chest, and the moving image processing unit measures the position of the diaphragm for each frame of the chest X-ray moving image, For each frame of the chest X-ray moving image that is set in advance, for each frame of the chest X-ray moving image, the pixel difference value between the frames before and after the time series is calculated as ventilation information, and the ventilation information is generated as dynamic information, A still image processing unit generates a new image including at least the above-described ray-thum image among a ray-thum image, a coronal image, and a sagittal image viewed from the back of the chest based on the CT image, and the generated ray-sum The position of the diaphragm is measured for the image, and the integration unit compares the position of the diaphragm measured by the still image processing unit with the position of the diaphragm measured for each frame by the moving image processing unit. Moving image information generated by the moving image processing unit, the moving image, the specified moving image frame, the laysum image generated by the still image processing unit, the coronal image, and the sagittal image Is integrated with at least one of the form information, and the integrated information is displayed on a screen.
また、本発明によるX線診断支援装置は、前記統合部が、特定した動画像のフレームと、静止画像処理部により生成されたレイサム画像とを、シフト及び回転させながら、各位置におけるピクセル差分値を算出し、該ピクセル差分値の絶対値の総和が最も小さいシフト及び回転位置を求め、該位置をマッチングした位置に決定し、該マッチング位置により、前記動態情報と形態情報とを統合して画面に表示することを特徴とする。 In the X-ray diagnosis support apparatus according to the present invention, the integration unit shifts and rotates the identified moving image frame and the laysum image generated by the still image processing unit, and the pixel difference value at each position. Calculating the shift and rotation position with the smallest sum of absolute values of the pixel difference values, determining the position as a matching position, and integrating the dynamic information and the form information by the matching position. It is characterized by displaying.
また、本発明によるX線診断支援プログラムは、診断部位の一連の動作のX線動画像が格納された動画像格納部と、前記診断部位のX線静止画像が格納された静止画像格納部とを備え、前記X線動画像及び静止画像を用いて診断を支援するX線診断支援装置が実行するプログラムであって、該X線診断支援装置を構成するコンピュータに、前記動画像格納部から動画像を読み出し、該動画像に基づいて前記診断部位に関する機能を示す動態情報を生成する処理と、前記静止画像格納部から静止画像を読み出し、該静止画像から前記診断部位に関する新たな画像を生成する処理と、前記生成された動態情報と、動画像及び新たな画像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示する処理とを実行させることを特徴とする。 An X-ray diagnosis support program according to the present invention includes a moving image storage unit that stores an X-ray moving image of a series of operations of a diagnostic region, and a still image storage unit that stores an X-ray still image of the diagnostic region. A program executed by an X-ray diagnosis support apparatus that supports diagnosis using the X-ray moving image and the still image, the computer configuring the X-ray diagnosis support apparatus from the moving image storage unit A process of reading an image, generating dynamic information indicating a function related to the diagnostic region based on the moving image, reading a still image from the still image storage unit, and generating a new image related to the diagnostic region from the still image And a process of integrating the generated dynamic information and at least one form information of the moving image and the new image and displaying the integrated information on a screen. To.
本発明によれば、被検体である肺等の機能の局所的な定量評価を行い、例えば、X線CT装置により得られたCT画像である形態情報と、動画像を解析処理して得られた換気に関する動態情報とを関連付けることができるから、胸部に関する形態情報と動態情報とを統合した新たな情報を、読影や診断のために有効に利用することが可能となる。 According to the present invention, a local quantitative evaluation of the function of a subject lung or the like is performed, and for example, obtained by analyzing morphological information, which is a CT image obtained by an X-ray CT apparatus, and a moving image. Therefore, new information obtained by integrating morphological information and dynamic information regarding the chest can be effectively used for interpretation and diagnosis.
1 X線診断支援装置
2 X線動画像装置
3 X線CT装置
11 CPU
12 RAM
13 ROM
14 HD
15,16 I/F
17 表示器
18 マウス
19 キーボード
21 動画像格納部
22 CT画像格納部
23 動画像処理部
24 CT画像処理部
25 統合部1 X-ray diagnosis support device 2 X-ray moving image device 3 X-ray CT device 11 CPU
12 RAM
13 ROM
14 HD
15,16 I / F
17 Display 18 Mouse 19 Keyboard 21 Moving Image Storage Unit 22 CT Image Storage Unit 23 Moving Image Processing Unit 24 CT Image Processing Unit 25 Integration Unit
以下、本発明の実施の形態について、図面を用いて詳細に説明する。
〔構成〕
図1は、本発明の実施の形態によるX線診断支援装置を含む全体システムの構成を示す概略図である。このシステムは、X線診断支援装置1、X線動画像装置2、及びX線CT装置3を備えており、X線診断支援装置1の内部構成はハードウェア資源を表している。このX線診断支援装置1は、プログラムに従って各処理を実行するCPU11、プログラムやデータを一時的に格納するRAM12、OS等のシステムプログラムやシステムデータが格納されているROM13、各処理を実行するためのプログラムやデータ、動画像及びCT画像が格納されるHD14、X線動画像装置2及びX線CT装置3との間の情報の入出力を中継するI/F15、動画像及びCT画像等を画面に表示する表示器17、オペレータの操作を入力するマウス18、キーボード19、及び、表示器17等を中継するI/F16を備えている。CPU11は、I/F15を介してX線動画像装置2及びX線CT装置3からそれぞれ動画像及びCT画像を入力し、HD14に格納する。また、CPU11は、各処理を実行するためのプログラム及びデータをHD14から読み出してRAM12に格納し、プログラムに従って、I/F16を介してオペレータによるマウス18及びキーボード19の操作により、HD14から動画像及びCT画像を読み出し、各処理を実行し、実行結果をI/F16を介して表示器17に表示する。Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
〔Constitution〕
FIG. 1 is a schematic diagram showing a configuration of an entire system including an X-ray diagnosis support apparatus according to an embodiment of the present invention. This system includes an X-ray diagnosis support apparatus 1, an X-ray moving image apparatus 2, and an X-ray CT apparatus 3, and the internal configuration of the X-ray diagnosis support apparatus 1 represents hardware resources. The X-ray diagnosis support apparatus 1 includes a CPU 11 that executes each process according to a program, a RAM 12 that temporarily stores programs and data, a ROM 13 that stores system programs and system data such as an OS, and each process. Programs and data, HD 14 for storing moving images and CT images, I / F 15 for relaying information input / output between the X-ray moving image apparatus 2 and the X-ray CT apparatus 3, moving images and CT images, etc. A display 17 displayed on the screen, a mouse 18 for inputting an operator's operation, a keyboard 19, and an I / F 16 for relaying the display 17 and the like are provided. The CPU 11 inputs a moving image and a CT image from the X-ray moving image apparatus 2 and the X-ray CT apparatus 3 via the I / F 15 and stores them in the HD 14. In addition, the CPU 11 reads a program and data for executing each process from the HD 14 and stores them in the RAM 12. According to the program, the CPU 11 operates the mouse 18 and the keyboard 19 through the I / F 16 to move the moving image and the data from the HD 14. The CT image is read out, each process is executed, and the execution result is displayed on the display unit 17 via the I / F 16.
図2は、図1に示したX線診断支援装置1の機能構成を示すブロック図である。このX線診断支援装置1は、動画像格納部21、CT画像格納部22、動画像処理部23、CT画像処理部24、及び、統合部25を備えている。動画像格納部21には、前述したFPD等のX線動画像装置2により撮影された胸部X線動画像が格納されている。胸部X線動画像は、吸気から呼気までの一連の呼吸動作における動画像であり、受検者毎に格納されている。CT画像格納部22には、X線CT装置3により撮影された胸部のCT画像が格納されている。CT画像は、胸部の各部を水平にスライスした静止画像であり、受検者毎に格納されている。ここで、動画像格納部21及びCT画像格納部22は、図1におけるHD14に相当する。 FIG. 2 is a block diagram showing a functional configuration of the X-ray diagnosis support apparatus 1 shown in FIG. The X-ray diagnosis support apparatus 1 includes a moving image storage unit 21, a CT image storage unit 22, a moving image processing unit 23, a CT image processing unit 24, and an integration unit 25. The moving image storage unit 21 stores a chest X-ray moving image captured by the X-ray moving image apparatus 2 such as the FPD described above. The chest X-ray moving image is a moving image in a series of breathing motions from inspiration to expiration, and is stored for each examinee. The CT image storage unit 22 stores CT images of the chest imaged by the X-ray CT apparatus 3. The CT image is a still image obtained by slicing each part of the chest horizontally, and is stored for each examinee. Here, the moving image storage unit 21 and the CT image storage unit 22 correspond to the HD 14 in FIG.
動画像処理部23は、オペレータの操作により、動画像格納部21から動画像を読み出し、横隔膜の位置を計測して移動量を求め、最大吸気時及び最大呼気時の動画像のフレームを特定してこれらのフレームのピクセル値の差であるピクセル差分値や、フレーム間のピクセル差分値を用いて、分割した胸部エリア毎に相対的な換気情報(換気量)を求める。また、オペレータの操作により、動画像格納部21から読み出した動画像、及び、計測した横隔膜の位置情報を画面に表示する。 The moving image processing unit 23 reads out a moving image from the moving image storage unit 21 by an operator's operation, measures the position of the diaphragm to obtain a moving amount, and specifies a frame of the moving image at the time of maximum inspiration and maximum expiration. Then, relative ventilation information (ventilation volume) is obtained for each divided chest area using a pixel difference value that is a difference between pixel values of these frames and a pixel difference value between frames. Further, the moving image read from the moving image storage unit 21 and the measured position information of the diaphragm are displayed on the screen by the operation of the operator.
CT画像処理部24は、オペレータの操作により、CT画像格納部22からCT画像を読み出し、読み出したCT画像間で線形補間を行い、等方的なCTデータを作成し、そこから、表示のためのMPR(Multi Planer Reconstruction)画像(コロナル像(胸部を前面から見た断層像)やサジタル像(胸部を側面から見た断層像)、及びX線写真のようなDRR(Digital Reconstruction Radiographs:レイサム(Ray Sum)画像)を作成し、このレイサム画像から横隔膜の位置を計測する。また、オペレータの操作により、CT画像格納部22から読み出したCT画像、及び、作成したコロナル像及びサジタル像を画面に表示する。 The CT image processing unit 24 reads out the CT image from the CT image storage unit 22 by the operator's operation, performs linear interpolation between the read CT images, creates isotropic CT data, and from there for display MPR (Multi Planer Reconstruction) images (coronal images (tomographic images of the chest viewed from the front), sagittal images (tomographic images of the chest viewed from the side), and DRR (Digital Reconstruction Radiographs: Latham ( Ray Sum) image) is measured, and the position of the diaphragm is measured from the ray sum image, and the CT image read from the CT image storage unit 22 and the created coronal image and sagittal image are displayed on the screen by the operation of the operator. indicate.
統合部25は、CT画像と呼吸レベルが一致する動画像のフレームを特定し、その後、CT画像と呼吸レベルが一致する動画像のフレームと、CT画像から作成されたレイサム画像との間で位置合わせを行う。また、統合部25は、動画像処理部23から、胸部を分割したエリア情報を示す胸部分割情報及び胸部エリア毎の相対的な換気情報を入力し、当該エリア毎の換気情報をコロナル像及び動画像に重ね合わせる。また、統合部25は、オペレータの操作により、換気情報を重ね合わせたコロナル像及び動画像を表示する。 The integration unit 25 identifies a moving image frame having a respiratory level that matches the CT image, and then positions the moving image frame having a respiratory level that matches the CT image and a laysum image created from the CT image. Align. Further, the integration unit 25 inputs from the moving image processing unit 23 chest division information indicating area information obtained by dividing the chest and relative ventilation information for each chest area, and converts the ventilation information for each area into a coronal image and a moving image. Superimpose on the image. Further, the integration unit 25 displays a coronal image and a moving image in which ventilation information is superimposed by an operator's operation.
〔動作〕
次に、図1及び図2に示したX線診断支援装置1の動作について説明する。図3は、X線診断支援装置1の動作を示すフローチャート図である。図3に示すように、X線診断支援装置1は、胸部X線動画像及びCT画像を読み出し、前処理を行い、動画像とCT画像との間の位置合わせ、動画像の動態解析、重ね合わせの一連の処理を行う。以下、具体的に説明する。動画像格納部21には、胸部X線動画像が、1344×1344ピクセルの画像サイズ(解像度)、0.32mm/ピクセル、12ビット/4096の階調にて、時系列に30フレーム分格納されているものとする。CT画像格納部22には、512×512ピクセルの画像サイズ、0.64mm/ピクセル、16ビット/65536の階調にて、スライス厚5mmのCT画像が66枚分格納されているものとする。[Operation]
Next, the operation of the X-ray diagnosis support apparatus 1 shown in FIGS. 1 and 2 will be described. FIG. 3 is a flowchart showing the operation of the X-ray diagnosis support apparatus 1. As shown in FIG. 3, the X-ray diagnosis support apparatus 1 reads a chest X-ray moving image and a CT image, performs preprocessing, alignment between the moving image and the CT image, dynamic analysis of the moving image, and superposition. A series of processing is performed. This will be specifically described below. The moving image storage unit 21 stores chest X-ray moving images for 30 frames in time series with an image size (resolution) of 1344 × 1344 pixels, 0.32 mm / pixel, and gradation of 12 bits / 4096. It shall be. It is assumed that the CT image storage unit 22 stores 66 CT images having a slice thickness of 5 mm at an image size of 512 × 512 pixels, a gradation of 0.64 mm / pixel, and 16 bits / 65536.
動画像処理部23は、動画像格納部21から前記30フレーム分の動画像を読み出す(ステップ301)。そして、CT画像の画像サイズに合わせるために、1344×1344ピクセルの画像サイズを512×512ピクセルの画像サイズに調整し、30フレーム分の当該サイズの動画像を作成する(ステップ302)。 The moving image processing unit 23 reads the 30 frames worth of moving images from the moving image storage unit 21 (step 301). Then, in order to match the image size of the CT image, the image size of 1344 × 1344 pixels is adjusted to the image size of 512 × 512 pixels, and a moving image of the size corresponding to 30 frames is created (step 302).
CT画像処理部24は、CT画像格納部22から前記66枚分のCT画像を読み出す(ステップ303)。そして、5mmのスライス厚の2枚のCT画像に対して、線形補間により画像データのボクセル化を行う(ステップ304)。具体的には、前記2枚のCT画像を用いて、1スライスのCT画像を7.8125スライスのCT画像に補間する。つまり、読み出した512ピクセル×512ピクセル×66スライス(枚)のCT画像をそれぞれ補間することにより、512ピクセル×512ピクセル×512スライス(枚)のCT画像を生成する。この場合、全体で、327.68mm×327.68mm×327.68mmのサイズである。そして、CT画像処理部24は、内挿補間したCT画像を用いて、MPR表示のためのコロナル像及びサジタル像を作成すると共に、DRR(レイサム画像)を作成する(ステップ305)。この場合、コロナル像、サジタル像及びレイサム画像は、それぞれ512×512ピクセルの画像サイズ及び0.64mm/ピクセルの画像となる。 The CT image processing unit 24 reads the 66 CT images from the CT image storage unit 22 (step 303). Then, voxelization of image data is performed by linear interpolation on two CT images having a slice thickness of 5 mm (step 304). Specifically, using the two CT images, a 1-slice CT image is interpolated into a 7.8125-slice CT image. That is, a CT image of 512 pixels × 512 pixels × 512 slices (sheets) is generated by interpolating the read CT images of 512 pixels × 512 pixels × 66 slices (sheets). In this case, the total size is 327.68 mm × 327.68 mm × 327.68 mm. Then, the CT image processing unit 24 creates a coronal image and a sagittal image for MPR display by using the interpolated CT image, and also creates a DRR (latham image) (step 305). In this case, the coronal image, the sagittal image, and the laysum image have an image size of 512 × 512 pixels and an image of 0.64 mm / pixel, respectively.
動画像処理部23は、ステップ302により作成された全フレームの動画像について、肺野領域を認識する(ステップ306)。また、CT画像処理部24は、閾値処理により、肺野領域を認識する(ステップ306)。具体的には、CT画像処理部24は、空気のCT値を閾値として設定し、各ピクセルにおけるCT値と閾値とを比較し、CT値が閾値以下のピクセルを抽出し、その抽出したピクセルの領域を肺野領域として認識する。 The moving image processing unit 23 recognizes a lung field region for all frames of moving images created in step 302 (step 306). Further, the CT image processing unit 24 recognizes the lung field region by threshold processing (step 306). Specifically, the CT image processing unit 24 sets the CT value of air as a threshold value, compares the CT value of each pixel with the threshold value, extracts pixels whose CT value is equal to or less than the threshold value, and extracts the extracted pixels. Recognize the region as a lung field region.
動画像処理部23は、ステップ302により作成された全フレームの動画像について、肺尖部及び横隔膜の位置を検出し、肺尖部から横隔膜までの距離を算出する(ステップ307)。また、CT画像処理部24は、ステップ305により作成されたレイサム画像について、肺尖部及び横隔膜の位置を検出し、肺尖部から横隔膜までの距離を算出する(ステップ307)。具体的には、動画像処理部23は1フレーム目の動画像について、CT画像処理部24はレイサム画像について、ピクセル値に基づき、その値が大きく変化する境界を検出し、肺野領域の上部において肺尖部の座標(位置)を決定し、肺野領域の下部において横隔膜の座標を決定する。そして、動画像処理部23は、1フレーム目の動画像から決定した肺尖部の座標の周辺を関心領域(ROI)として設定すると共に、同様に横隔膜の座標の周辺を関心領域として設定し、2フレーム目以降の動画像について、その関心領域に対して追跡し、ピクセル値が大きく変化する境界を検出し、肺尖部及び横隔膜の座標をそれぞれ決定する。 The moving image processing unit 23 detects the positions of the lung apex and the diaphragm in the moving image of all frames created in step 302, and calculates the distance from the lung apex to the diaphragm (step 307). In addition, the CT image processing unit 24 detects the positions of the lung apex and the diaphragm for the latham image created in step 305, and calculates the distance from the apex to the diaphragm (step 307). Specifically, the moving image processing unit 23 detects a boundary where the value of the moving image of the first frame and the CT image processing unit 24 greatly change the value of the latsum image based on the pixel value, and the upper part of the lung field region. The coordinates (position) of the pulmonary apex are determined, and the coordinates of the diaphragm are determined below the lung field region. The moving image processing unit 23 sets the periphery of the lung apex coordinates determined from the moving image of the first frame as a region of interest (ROI), and similarly sets the periphery of the coordinates of the diaphragm as a region of interest. The moving image after the second frame is tracked with respect to the region of interest, the boundary where the pixel value changes greatly is detected, and the coordinates of the lung apex and the diaphragm are determined.
統合部25は、ステップ307において算出した肺尖部から横隔膜までの距離を動画像処理部23及びCT画像処理部24から入力し、レイサム画像における当該距離と動画像の各フレームにおける距離とを比較し、距離が等しいフレームを特定する(ステップ308)。すなわち、CT画像の呼吸レベルと動画像の呼吸レベルとが一致するフレームを、複数のフレームから特定する。 The integration unit 25 inputs the distance from the lung apex to the diaphragm calculated in step 307 from the moving image processing unit 23 and the CT image processing unit 24, and compares the distance in the laysum image with the distance in each frame of the moving image. Then, frames having the same distance are specified (step 308). That is, a frame in which the respiratory level of the CT image matches the respiratory level of the moving image is specified from a plurality of frames.
統合部25は、レイサム画像をCT画像処理部24から入力し、特定した動画像のフレームとレイサム画像との間で位置合わせを行う(ステップ309)。具体的には、入力した特定フレームとレイサム画像とを少しずつシフト及び回転させながら、各ピクセルにおける両画像のピクセル値の差分を算出し、この差分値の絶対値の総和が最も小さいシフト及び回転位置を求める。この位置を、両画像がマッチングした位置に決定する。つまり、統合部25は、以下の数式のRを、両画像を少しずつシフト及び回転させてそれぞれ求め、Rが最小となるシフト及び回転位置をマッチング位置に決定する。
ここで、Rはピクセル差分値の絶対値の総和、F(x,y)は動画像の特定フレーム、G(x,y)はレイサム画像、Rotationは回転関数(入力は画像及び回転角度、出力は回転された画像)、xは横方向の座標(0<x<M)、yは縦方向の座標(0<y<N)、Mは横方向の画像サイズ(単位はピクセル)、Nは縦方向の画像サイズ(単位はピクセル)、dxは横方向のシフト量(0<dx<10,単位はピクセル)、dyは縦方向のシフト量(0<dy<10,単位はピクセル)、drは回転量(0<dr<5,単位はdeg)とする。The integration unit 25 inputs the laysum image from the CT image processing unit 24, and performs alignment between the frame of the identified moving image and the laysum image (step 309). Specifically, while shifting and rotating the input specific frame and the laysum image little by little, the difference between the pixel values of both images in each pixel is calculated, and the shift and rotation with the smallest sum of absolute values of the difference values are calculated. Find the position. This position is determined as a position where both images are matched. That is, the integration unit 25 obtains R in the following mathematical formula by shifting and rotating both images little by little, and determines the shift and rotation position where R is the minimum as the matching position.
Here, R is the sum of absolute values of pixel difference values, F (x, y) is a specific frame of a moving image, G (x, y) is a ray-sum image, Rotation is a rotation function (input is image and rotation angle, output) Is a horizontal coordinate (0 <x <M), y is a vertical coordinate (0 <y <N), M is a horizontal image size (unit is pixel), and N is Vertical image size (unit is pixel), dx is horizontal shift amount (0 <dx <10, unit is pixel), dy is vertical shift amount (0 <dy <10, unit is pixel), dr Is a rotation amount (0 <dr <5, unit is deg).
尚、統合部25は、入力したフレーム及びレイサム画像の512×512ピクセルの画像サイズを、それぞれ128×128ピクセルの画像サイズに縮小し、この縮小したフレームとレイサム画像とを用いて、マッチング位置を決定するようにしてもよい。これにより、マッチング位置を決定するための時間を短縮することができる。 Note that the integration unit 25 reduces the image size of 512 × 512 pixels of the input frame and the laysum image to an image size of 128 × 128 pixels, respectively, and uses this reduced frame and the laysum image to determine the matching position. It may be determined. Thereby, the time for determining a matching position can be shortened.
動画像処理部23は、ステップ308において特定した動画像のフレームについて、左肺野領域及び右肺野領域に対してx軸方向に水平なエリア毎に所定の数だけ分割する(ステップ310)。ここで、肺機能における生態上の追従精度の観点からすると、左右の肺野領域に対してそれぞれ8分割するのが好適である。 The moving image processing unit 23 divides the frame of the moving image identified in step 308 by a predetermined number for each area horizontal in the x-axis direction with respect to the left lung field and the right lung field (step 310). Here, from the viewpoint of ecological follow-up accuracy in lung function, it is preferable to divide the left and right lung field regions into 8 parts.
動画像処理部23は、ステップ310において分割したエリア毎に、ピクセル値の平均を算出し、時系列に隣り合う動画像のフレーム間のピクセル値の差分をそれぞれ算出する。また、最大吸気時のピクセル値と最大呼気時のピクセル値との差分も算出する(ステップ311)。具体的には、動画像処理部23は、時系列の動画像フレームにおいて、第1のフレームと第2のフレームとの間のピクセル差分値(ピクセル値の平均の差分値)、第2フレームと第3フレームとの間のピクセル差分値、・・・をそれぞれエリア毎に算出する。また、動画像処理部23は、ステップ307において算出した全ての動画像フレームの肺尖部から横隔膜までの距離から、その距離が最大となるフレームを特定すると共に、その距離が最小となるフレームを特定する。ここで、距離が最大のフレームは、最大吸気時の形態を示し、距離が最小のフレームは、最大呼気時の形態を示している。そして、動画像処理部23は、最大吸気時のフレームのピクセル値と、最大呼気時のフレームのピクセル値との差分であるピクセル差分値をエリア毎に算出する。ここで、動画像処理部23が算出したピクセル差分値は、換気情報となる。 The moving image processing unit 23 calculates an average of pixel values for each area divided in step 310, and calculates a difference in pixel values between frames of moving images adjacent in time series. Also, the difference between the pixel value at the time of maximum inspiration and the pixel value at the time of maximum expiration is calculated (step 311). Specifically, the moving image processing unit 23, in a time-series moving image frame, the pixel difference value (average difference value of pixel values) between the first frame and the second frame, the second frame, A pixel difference value from the third frame,... Is calculated for each area. Further, the moving image processing unit 23 specifies the frame having the maximum distance from the distances from the lung apex to the diaphragm of all the moving image frames calculated in step 307, and selects the frame having the minimum distance. Identify. Here, the frame with the maximum distance indicates the form at the time of maximum inspiration, and the frame with the minimum distance indicates the form at the time of maximum expiration. Then, the moving image processing unit 23 calculates, for each area, a pixel difference value that is a difference between the pixel value of the frame at the time of the maximum inspiration and the pixel value of the frame at the time of the maximum expiration. Here, the pixel difference value calculated by the moving image processing unit 23 is ventilation information.
統合部25は、胸部分割情報、エリア毎のフレーム間ピクセル差分値、及び最大吸気/呼気時ピクセル差分値を動画像処理部23から入力し、エリア毎のフレーム間ピクセル差分値に応じた濃淡色を動画像に重ね合わせて表示し、また、エリア毎の最大吸気/呼気時ピクセル差分値に応じた濃淡色をコロナル像及び動画像に重ね合わせて表示する(ステップ312)。画面表示の詳細については後述する。 The integration unit 25 inputs the chest division information, the inter-frame pixel difference value for each area, and the maximum inspiration / expiration pixel difference value from the moving image processing unit 23, and the gray color corresponding to the inter-frame pixel difference value for each area. Are superimposed and displayed on the moving image, and the shade color corresponding to the maximum inspiration / expiration pixel difference value for each area is superimposed and displayed on the coronal image and the moving image (step 312). Details of the screen display will be described later.
〔画面表示〕
図4は、図1に示した表示器17に表示される画面の一例を示す図である。この画面は、オペレータがマウス18やキーボード19を操作することにより表示される。図4において、画面左上の領域401にはCT画像、画面中央上の領域402には動画像、画面左下の領域403にはコロナル像、その右隣の領域404にはサジタル像、その右隣の領域405には領域402と同じ動画像、その右隣の領域406には肺尖部から横隔膜までの距離を示すグラフ、画面右側にはオペレータによる操作領域407〜409が表示される。以下、オペレータの操作から画面表示までの動作について説明する。〔Screen display〕
FIG. 4 is a diagram showing an example of a screen displayed on the display unit 17 shown in FIG. This screen is displayed when the operator operates the mouse 18 or the keyboard 19. In FIG. 4, a CT image is in a region 401 at the upper left of the screen, a moving image is in a region 402 at the center of the screen, a coronal image is in a region 403 at the lower left of the screen, a sagittal image is in a region 404 to the right of the screen. The area 405 displays the same moving image as the area 402, the area 406 on the right side of the moving image displays a graph indicating the distance from the lung apex to the diaphragm, and the operation areas 407 to 409 by the operator are displayed on the right side of the screen. Hereinafter, the operation from the operation of the operator to the screen display will be described.
まず、オペレータが領域407のホルダーまたはファイルを選択し、領域409の「Open images」ボタンを押下すると、動画像がHD14(動画像格納部21)から読み出され、領域402,405に表示される。この動画像は、吸気から呼気までの一連の呼吸動作の画像である。同様に、オペレータが領域407において同じ受検者のホルダーまたはファイルを選択し、領域409の「Open CT」ボタンを押下すると、CT画像がHD14(CT画像格納部22)から読み出され、領域401に表示される。領域408には、x軸方向の画像サイズ、y軸方向の画像サイズ、ピクセルサイズ、スライス厚、階調ビット数が設定及び表示される。 First, when the operator selects a folder or file in the area 407 and presses the “Open images” button in the area 409, a moving image is read from the HD 14 (moving image storage unit 21) and displayed in the areas 402 and 405. . This moving image is a series of breathing motion images from inspiration to expiration. Similarly, when the operator selects the same examinee's holder or file in the area 407 and presses the “Open CT” button in the area 409, the CT image is read out from the HD 14 (CT image storage unit 22) and stored in the area 401. Is displayed. In the area 408, the image size in the x-axis direction, the image size in the y-axis direction, the pixel size, the slice thickness, and the number of gradation bits are set and displayed.
オペレータが領域409の「Lung area」ボタンを押下すると、動画像の肺野認識が行われるとともに、動画像における横隔膜の位置を計測して移動量が求められ、動画像の各フレームにおける肺尖部から横隔膜までの距離のグラフが領域406に表示される。この場合、オペレータが領域406の下部のスクロールバーを左右に操作すると、スクロールバーの位置に対応した動画像のフレームが領域402,405に表示される。すなわち、オペレータがこのスクロールバーを左端から右端まで移動させると、領域402,405には吸気から呼気までの一連の呼吸動作の動画像が表示される。 When the operator presses the “Lung area” button in the area 409, lung field recognition of the moving image is performed, and the movement amount is obtained by measuring the position of the diaphragm in the moving image, and the lung apex portion in each frame of the moving image A graph of the distance from the diaphragm to the diaphragm is displayed in area 406. In this case, when the operator operates the scroll bar at the bottom of the area 406 to the left and right, a frame of a moving image corresponding to the position of the scroll bar is displayed in the areas 402 and 405. That is, when the operator moves the scroll bar from the left end to the right end, moving images of a series of breathing motions from inspiration to expiration are displayed in areas 402 and 405.
図5は、図4に示した領域406を拡大したグラフであり、吸気から呼気までの一連の呼吸動作における肺尖部から横隔膜までの距離を示している。○(丸)印の特性は右の肺における前記距離であり、×(バツ)印は左の肺における前記距離である。 FIG. 5 is an enlarged graph of the region 406 shown in FIG. 4 and shows the distance from the lung apex to the diaphragm in a series of breathing operations from inspiration to expiration. The characteristic of the circle (circle) is the distance in the right lung, and the cross (cross) is the distance in the left lung.
図4に戻って、オペレータが領域409の「Isotropic data」ボタンを押下すると、CT画像のボクセル化が行われ、コロナル像、サジタル像及びレイサム画像が作成され、コロナル像が領域403に、サジタル像が領域404にそれぞれ表示される。尚、この段階では、領域403,405には、換気情報が重ね合わせ表示されていない。また、レイサム画像は画面には表示されない。この場合、オペレータが領域401の下部のスクロールバーを左右に操作すると、スクロールバーの位置に対応したサジタル像が領域404に表示される。また、オペレータが領域401の右部のスクロールバーを上下に操作すると、スクロールバーの位置に対応したコロナル像が領域403に表示される。また、オペレータが領域404の右部のスクロールバーを上下に操作すると、領域403,404においてスクロールバーの位置に対応したCT画像(アキシャル画像)が領域401に表示される。また、レイサム画像における肺尖部から横隔膜までの距離が計測される。 Returning to FIG. 4, when the operator presses the “Isotropic data” button in the area 409, the CT image is voxelized, and a coronal image, a sagittal image, and a latham image are created. Are displayed in the area 404 respectively. At this stage, ventilation information is not superimposed on the areas 403 and 405. In addition, the latham image is not displayed on the screen. In this case, when the operator operates the scroll bar at the bottom of the area 401 to the left or right, a sagittal image corresponding to the position of the scroll bar is displayed in the area 404. When the operator operates the scroll bar on the right side of the area 401 up and down, a coronal image corresponding to the position of the scroll bar is displayed in the area 403. When the operator operates the scroll bar on the right side of the region 404 up and down, CT images (axial images) corresponding to the positions of the scroll bars in the regions 403 and 404 are displayed in the region 401. In addition, the distance from the lung apex to the diaphragm in the Latham image is measured.
オペレータが領域409の「Registration」のボタンを押下すると、動画像の各フレームにおける肺尖部から横隔膜までの距離と、レイサム画像における肺尖部から横隔膜までの距離とを用いて呼吸レベルが一致する動画像のフレームが特定され、そのフレームが領域402,405に表示される。そして、該画像間で位置合わせが行なわれる。また、領域406のスクロールバーの位置は、その特定フレームに対応した位置となる。尚、この段階では、領域403,405には、換気情報が重ね合わせ表示されていない。 When the operator presses the “Registration” button in the region 409, the respiratory levels match using the distance from the apex to the diaphragm in each frame of the moving image and the distance from the apex to the diaphragm in the layham image. A frame of the moving image is specified, and the frame is displayed in the areas 402 and 405. Then, alignment is performed between the images. The position of the scroll bar in the area 406 is a position corresponding to the specific frame. At this stage, ventilation information is not superimposed on the areas 403 and 405.
オペレータが領域409の「Ventilation」ボタンを押下すると、胸部エリア毎の最大吸気/呼気時のピクセル差分値が換気情報として、領域403のコロナル像及び領域405の動画像の特定フレームに重ね合わせて表示される。また、胸部エリア毎のフレーム間ピクセル差分値が換気情報として、領域405の動画像に重ね合わせて表示される。この場合、オペレータが領域406の下部のスクロールバーを左右に操作すると、スクロールバーの位置に対応した動画像が領域405に表示されると共に、フレーム間ピクセル差分値の換気情報が重ね合わせて表示される。この領域405の動画像により、吸気及び呼気状態を把握することができる。 When the operator presses the “Ventilation” button in the region 409, the pixel difference value at the time of maximum inspiration / expiration for each chest area is displayed as ventilation information superimposed on a specific frame of the coronal image in the region 403 and the moving image in the region 405. Is done. Further, the inter-frame pixel difference value for each chest area is displayed as ventilation information superimposed on the moving image in the region 405. In this case, when the operator operates the scroll bar at the bottom of the area 406 left and right, a moving image corresponding to the position of the scroll bar is displayed in the area 405 and ventilation information of the inter-frame pixel difference value is displayed in an overlapping manner. The The inspiration and expiration state can be grasped from the moving image in this area 405.
図6は、図4に示した領域403を拡大した図であり、最大吸気/呼気時のピクセル差分値が換気情報としてコロナル像に重ね合わせて表示されている。図6に示すように、換気情報は、左右の胸部をそれぞれ8等分したエリア毎に表示され、ピクセル差分値が大きさに応じて濃淡で示され、ピクセル差分値が大きいほどエリアの色は濃くなっている。また、図7は、図4に示した領域405を拡大した図であり、最大吸気/呼気時のピクセル差分値が換気情報として動画像の特定フレームに重ね合わせて表示されている。図6と同様に、換気情報は、左右の胸部をそれぞれ8等分したエリア毎に表示され、ピクセル差分値が大きさに応じて濃淡で示され、ピクセル差分値が大きいほどエリアの色は濃くなっている。また、中央の山型の曲線は、中央の縦線を0とした各エリアのピクセル差分値を示しており、左側の曲線は左胸部の差分値であり、右側の曲線は右胸部の差分値である。胸部の下部に近いほど、ピクセル差分値が大きくなっているから、エリアの色は濃くなり、曲線は広がっている。 FIG. 6 is an enlarged view of the region 403 shown in FIG. 4, and the pixel difference value at the time of maximum inspiration / expiration is displayed superimposed on the coronal image as ventilation information. As shown in FIG. 6, the ventilation information is displayed for each area obtained by dividing the left and right chests into eight equal parts, and the pixel difference value is shown in shades according to the size. The larger the pixel difference value, the more the color of the area It is dark. FIG. 7 is an enlarged view of the area 405 shown in FIG. 4, and the pixel difference value at the time of the maximum inspiration / expiration is displayed as ventilation information superimposed on a specific frame of the moving image. As in FIG. 6, ventilation information is displayed for each area obtained by dividing the left and right breasts into eight equal parts, and the pixel difference value is shown in shades according to the size. The larger the pixel difference value, the darker the area color. It has become. The central mountain-shaped curve indicates the pixel difference value of each area with the central vertical line being 0, the left curve is the left chest difference value, and the right curve is the right chest difference value. It is. The closer to the lower part of the chest, the larger the pixel difference value, so the area becomes darker and the curve is wider.
図4及び図6を参照して、オペレータが領域403の下部のボタンを押下することにより、領域403において、オリジナルのコロナル像の表示、または、コロナル像及び最大吸気/呼気時のピクセル差分値(換気情報)の重ね合わせ表示を選択することができる。また、図4及び図7を参照して、領域405の下部のボタンを押下することにより、領域405において、動画像の表示、動画像の特定フレーム及び最大吸気/呼気時のピクセル差分値(換気情報)の重ね合わせ表示、または、動画像及びフレーム間ピクセル差分値の重ね合わせ表示を選択することができる。 4 and 6, when the operator presses a button at the bottom of the region 403, the display of the original coronal image or the pixel difference value at the time of coronal image and maximum inspiration / expiration ( The overlay display of ventilation information) can be selected. 4 and 7, when a button at the bottom of the area 405 is pressed, in the area 405, the display of the moving image, the specific frame of the moving image, and the pixel difference value during the maximum inspiration / expiration (ventilation) Information) or superimposed display of moving image and inter-frame pixel difference value can be selected.
ここで、図には示してないが、領域405において、胸部エリア毎のフレーム間ピクセル差分値が換気情報として動画像に重ね合わせて表示された場合は、ピクセル差分値がプラスのときは呼気フェーズの色、ピクセル差分値がマイナスのときは吸気フェーズの色に分け、さらに、ピクセル差分値の大きさに応じたその色の濃淡で示される。 Here, although not shown in the figure, in the area 405, when the inter-frame pixel difference value for each chest area is displayed as ventilation information superimposed on the moving image, when the pixel difference value is positive, the expiration phase When the pixel difference value and the pixel difference value are negative, the color is divided into the colors of the intake phase, and is further shown by the shading of the color according to the size of the pixel difference value.
また、領域409において、オペレータの操作により、画面のコントラスト等の設定や表示、ビットマップ等のデータ保存方式の設定やデータの保存等が行われる。 Further, in the area 409, setting and display of the screen contrast and the like, setting of a data storage method such as a bitmap, storage of data, and the like are performed by an operator's operation.
以上のように、本発明の実施の形態によるX線診断支援装置1によれば、統合部25が、分割された胸部エリア毎に、胸部の形態情報と換気情報である呼吸性動態情報とを重ね合わせて画面に表示するようにした。これにより、局所的な定量評価を実現し、読影や診断のために有効に利用できる評価結果を得ることが可能となる。 As described above, according to the X-ray diagnosis support apparatus 1 according to the embodiment of the present invention, the integration unit 25 obtains the chest shape information and the respiratory dynamics information that is ventilation information for each divided chest area. Overlaid on the screen. As a result, local quantitative evaluation can be realized, and an evaluation result that can be effectively used for interpretation and diagnosis can be obtained.
以上、実施の形態を挙げて本発明を説明したが、本発明は上記実施の形態に限定されるものではなく、本発明の精神及び意図を逸脱しない限り、種々変形が可能である。例えば、前記実施の形態では、胸部におけるピクセル差分値を換気情報として形態情報に重ね合わせたが、部位は胸部に限るものではなく、心臓等であってもよく、また、換気情報に限定するものではなく、他の動態情報、例えば、血管の拡散及び縮小に関する情報であってもよい。 The present invention has been described with reference to the embodiment. However, the present invention is not limited to the above embodiment, and various modifications can be made without departing from the spirit and intention of the present invention. For example, in the above embodiment, the pixel difference value in the chest is superimposed on the form information as ventilation information. However, the part is not limited to the chest but may be a heart or the like, and is limited to the ventilation information. Instead, it may be other dynamic information, for example, information related to blood vessel diffusion and contraction.
また、胸部X線動画像の画像サイズ、階調数、フレーム数等の条件、及びCT画像の画像サイズ、階調数、スライス厚等の条件は、前記実施の形態に示した条件に限定されるものではない。また、動画像処理部23により調整される画像サイズ、及びCT画像処理部24によりボクセル化されるスライス枚数は、前記実施の形態に示した画像サイズ及びスライス枚数に限定されるものではない。また、図4,6,7では、換気情報が、左右の胸部をそれぞれ8等分したエリア毎に表示されているが、8等分のエリアに限定されるものではない。 The conditions such as the image size, the number of gradations, and the number of frames of the chest X-ray moving image, and the conditions such as the image size, the number of gradations, and the slice thickness of the CT image are limited to the conditions described in the above embodiment. It is not something. The image size adjusted by the moving image processing unit 23 and the number of slices voxelized by the CT image processing unit 24 are not limited to the image size and the number of slices shown in the above embodiment. In FIGS. 4, 6, and 7, the ventilation information is displayed for each area obtained by dividing the left and right breasts into eight equal parts. However, the ventilation information is not limited to eight equal areas.
尚、X線診断支援装置1は、図1に示したように、CPU11、RAM12等の揮発性の記憶媒体、ROM13等の不揮発性の記憶媒体、マウス18やキーボード19、ポインティングデバイス等の入力装置、画像やデータを表示する表示器17、及び外部の装置と通信をするためのインタフェースI/F15を備えたコンピュータによって構成される。X線診断支援装置1に備えた動画像処理部23、CT画像処理部24及び統合部25の各機能は、これらの機能を記述したプログラムをCPU11に実行させることによりそれぞれ実現される。また、これらのプログラムは、磁気ディスク(フロッピィーディスク、ハードディスクHD14等)、光ディスク(CD−ROM、DVD等)、半導体メモリ等の記憶媒体に格納して頒布することもできる。
As shown in FIG. 1, the X-ray diagnosis support apparatus 1 includes a volatile storage medium such as a CPU 11 and a RAM 12, a non-volatile storage medium such as a ROM 13, an input device such as a mouse 18, a keyboard 19, and a pointing device. And a display 17 for displaying images and data, and a computer having an interface I / F 15 for communicating with an external device. Each function of the moving image processing unit 23, the CT image processing unit 24, and the integration unit 25 provided in the X-ray diagnosis support apparatus 1 is realized by causing the CPU 11 to execute a program describing these functions. These programs can also be stored and distributed in a storage medium such as a magnetic disk (floppy disk, hard disk HD14, etc.), an optical disk (CD-ROM, DVD, etc.), or a semiconductor memory.
Claims (6)
前記診断部位のX線静止画像が格納された静止画像格納部と、
前記動画像格納部から動画像を読み出し、該動画像に基づいて前記診断部位に関する機能を示す動態情報を生成する動画像処理部と、
前記静止画像格納部から静止画像を読み出し、該静止画像から前記診断部位に関する新たな画像を生成する静止画像処理部と、
前記動画像処理部により生成された動態情報と、前記動画像及び前記静止画像処理部により生成された新たな画像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示する統合部とを備えたことを特徴とするX線診断支援装置。A moving image storage unit in which an X-ray moving image of a series of operations of the diagnostic region is stored;
A still image storage unit storing an X-ray still image of the diagnostic region;
A moving image processing unit that reads out a moving image from the moving image storage unit and generates dynamic information indicating a function related to the diagnostic region based on the moving image;
A still image processing unit that reads out a still image from the still image storage unit and generates a new image related to the diagnostic region from the still image;
The dynamic information generated by the moving image processing unit and at least one form information of the moving image and the new image generated by the still image processing unit are integrated, and the integrated information is displayed on the screen. An X-ray diagnosis support apparatus, comprising:
前記動画像格納部に格納されたX線動画像を、胸部の後面から撮影し、吸気から呼気までの一連の呼吸動作を示す胸部X線動画像とし、
前記静止画像格納部に格納されたX線静止画像を、胸部のスライス面を示すCT画像とし、
前記動画像処理部は、胸部X線動画像の各フレームについて横隔膜の位置を計測し、該横隔膜の位置から最大吸気のフレーム及び最大呼気のフレームを特定し、予め設定された胸部のエリア毎に、前記最大吸気のフレームと最大呼気のフレームとの間のピクセル差分値を算出して換気情報とし、該換気情報を動態情報として生成し、
前記静止画像処理部は、CT画像に基づいて胸部の後面から見たレイサム画像、コロナル像、及びサジタル像のうちの少なくとも一つを新たな画像として生成することを特徴とするX線診断支援装置。The X-ray diagnosis support apparatus according to claim 1,
The X-ray moving image stored in the moving image storage unit is taken from the back of the chest, and is a chest X-ray moving image showing a series of breathing movements from inspiration to expiration,
The X-ray still image stored in the still image storage unit is a CT image showing a slice plane of the chest,
The moving image processing unit measures the position of the diaphragm for each frame of the chest X-ray moving image, specifies the frame of maximum inspiration and the frame of maximum exhalation from the position of the diaphragm, and sets each predetermined chest area , Calculating a pixel difference value between the maximum inspiration frame and the maximum expiration frame to obtain ventilation information, and generating the ventilation information as dynamic information,
The still image processing unit generates, as a new image, at least one of a ray-sum image, a coronal image, and a sagittal image viewed from the rear of the chest based on a CT image. .
前記動画像格納部に格納されたX線動画像を、胸部の後面から撮影し、吸気から呼気までの一連の呼吸動作を示す胸部X線動画像とし、
前記静止画像格納部に格納されたX線静止画像を、胸部のスライス面を示すCT画像とし、
前記動画像処理部は、胸部X線動画像の各フレームについて横隔膜の位置を計測し、予め設定された胸部のエリア毎に、胸部X線動画像の各フレームについて、時系列の前後におけるフレーム間のピクセル差分値を算出して換気情報とし、該換気情報を動態情報として生成し、
前記静止画像処理部は、CT画像に基づいて胸部の後面から見たレイサム画像、コロナル像、及びサジタル像のうちの少なくとも前記レイサム画像を含む新たな画像を生成し、該生成したレイサム画像について横隔膜の位置を計測し、
前記統合部は、静止画像処理部により計測された横隔膜の位置と、動画像処理部により各フレームについて計測された横隔膜の位置とを比較し、両位置が一致する動画像のフレームを特定し、動画像処理部により生成された動態情報と、前記動画像、前記特定した動画像のフレーム、静止画像処理部により生成されたレイサム画像、コロナル像、及びサジタル像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示することを特徴とするX線診断支援装置。The X-ray diagnosis support apparatus according to claim 1,
The X-ray moving image stored in the moving image storage unit is taken from the back of the chest, and is a chest X-ray moving image showing a series of breathing movements from inspiration to expiration,
The X-ray still image stored in the still image storage unit is a CT image showing a slice plane of the chest,
The moving image processing unit measures the position of the diaphragm for each frame of the chest X-ray moving image, and for each predetermined chest area, between each frame of the chest X-ray moving image before and after the time series The pixel difference value is calculated as ventilation information, the ventilation information is generated as dynamic information,
The still image processing unit generates a new image including at least the Latham image among a Latham image, a coronal image, and a sagittal image viewed from the back of the chest based on the CT image, and the diaphragm is generated with respect to the generated Latham image. Measure the position of
The integration unit compares the position of the diaphragm measured by the still image processing unit with the position of the diaphragm measured for each frame by the moving image processing unit, specifies a moving image frame in which both positions match, Dynamic information generated by the moving image processing unit, and at least one form information of the moving image, the frame of the identified moving image, a laysum image generated by the still image processing unit, a coronal image, and a sagittal image Are integrated, and the integrated information is displayed on a screen.
前記統合部は、特定した動画像のフレームと、静止画像処理部により生成されたレイサム画像とを、シフト及び回転させながら、各位置におけるピクセル差分値を算出し、該ピクセル差分値の絶対値の総和が最も小さいシフト及び回転位置を求め、該位置をマッチングした位置に決定し、該マッチング位置により、前記動態情報と形態情報とを統合して画面に表示することを特徴とするX線診断支援装置。The X-ray diagnosis support apparatus according to claim 3,
The integration unit calculates a pixel difference value at each position while shifting and rotating the frame of the identified moving image and the laysum image generated by the still image processing unit, and calculates an absolute value of the pixel difference value. X-ray diagnosis support characterized in that a shift and rotation position having the smallest sum is obtained, the position is determined as a matching position, and the dynamic information and form information are integrated and displayed on the screen by the matching position. apparatus.
前記動画像格納部から動画像を読み出し、該動画像に基づいて前記診断部位に関する機能を示す動態情報を生成する処理と、
前記静止画像格納部から静止画像を読み出し、該静止画像から前記診断部位に関する新たな画像を生成する処理と、
前記生成された動態情報と、動画像及び新たな画像のうちの少なくとも一つの形態情報とを統合し、該統合した情報を画面に表示する処理とを実行させるX線診断支援プログラム。A moving image storage unit storing an X-ray moving image of a series of operations of a diagnostic region; and a still image storage unit storing an X-ray still image of the diagnostic region, the X-ray moving image and the still image being stored A program executed by an X-ray diagnosis support apparatus that supports diagnosis by using a computer constituting the X-ray diagnosis support apparatus,
A process of reading out a moving image from the moving image storage unit and generating dynamic information indicating a function related to the diagnostic region based on the moving image;
A process of reading a still image from the still image storage unit and generating a new image related to the diagnostic region from the still image;
An X-ray diagnosis support program for executing a process of integrating the generated dynamic information and at least one form information of a moving image and a new image and displaying the integrated information on a screen.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007522243A JP4797173B2 (en) | 2005-06-21 | 2006-06-13 | X-ray diagnosis support apparatus, program, and recording medium |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005180758 | 2005-06-21 | ||
JP2005180758 | 2005-06-21 | ||
JP2007522243A JP4797173B2 (en) | 2005-06-21 | 2006-06-13 | X-ray diagnosis support apparatus, program, and recording medium |
PCT/JP2006/311847 WO2006137294A1 (en) | 2005-06-21 | 2006-06-13 | X-ray diagnosis supporting apparatus, program, and recording medium |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2006137294A1 JPWO2006137294A1 (en) | 2009-01-15 |
JP4797173B2 true JP4797173B2 (en) | 2011-10-19 |
Family
ID=37570323
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007522243A Active JP4797173B2 (en) | 2005-06-21 | 2006-06-13 | X-ray diagnosis support apparatus, program, and recording medium |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP4797173B2 (en) |
WO (1) | WO2006137294A1 (en) |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5896737B2 (en) * | 2008-04-03 | 2016-03-30 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Respirometer, Respirometer operating method, and Respiratory computer program |
JP5200656B2 (en) * | 2008-05-15 | 2013-06-05 | コニカミノルタエムジー株式会社 | Dynamic shooting system |
JP5157766B2 (en) * | 2008-09-05 | 2013-03-06 | コニカミノルタホールディングス株式会社 | Image display apparatus and program |
JP5192354B2 (en) * | 2008-11-04 | 2013-05-08 | 株式会社リガク | Region specifying method, volume measuring method using the same, region specifying device and region specifying program |
JP5190338B2 (en) * | 2008-12-04 | 2013-04-24 | 日立アロカメディカル株式会社 | X-ray image forming apparatus |
JP5157919B2 (en) * | 2009-01-06 | 2013-03-06 | コニカミノルタホールディングス株式会社 | Image display apparatus and program |
US9064302B2 (en) | 2010-08-27 | 2015-06-23 | Konica Minolta, Inc. | Diagnosis assistance system and computer readable storage medium |
WO2012026146A1 (en) * | 2010-08-27 | 2012-03-01 | コニカミノルタエムジー株式会社 | Thoracic diagnosis assistance system and program |
JP5976293B2 (en) | 2010-10-08 | 2016-08-23 | 東芝メディカルシステムズ株式会社 | Image processing apparatus and image processing method |
JP2012110400A (en) * | 2010-11-22 | 2012-06-14 | Konica Minolta Medical & Graphic Inc | Dynamic diagnosis support information generation system |
JP5672147B2 (en) * | 2011-05-24 | 2015-02-18 | コニカミノルタ株式会社 | Chest diagnosis support information generation system |
JP5238893B1 (en) * | 2012-02-22 | 2013-07-17 | 株式会社東芝 | X-ray CT apparatus, image display apparatus, and image display method |
US20150042677A1 (en) * | 2012-03-23 | 2015-02-12 | Konica Minolta, Inc. | Image-generating apparatus |
JP5500216B2 (en) * | 2012-07-13 | 2014-05-21 | コニカミノルタ株式会社 | Image display apparatus and program |
JP2012245395A (en) * | 2012-09-19 | 2012-12-13 | Konica Minolta Medical & Graphic Inc | Dynamic image-capturing system and diagnosis-assisting information-generating method |
CN105188541A (en) | 2013-05-16 | 2015-12-23 | 柯尼卡美能达株式会社 | Image processing device and program |
CN105246406A (en) | 2013-05-31 | 2016-01-13 | 柯尼卡美能达株式会社 | Image processing device and program |
JP5907185B2 (en) * | 2014-01-27 | 2016-04-26 | コニカミノルタ株式会社 | Dynamic image diagnosis support system |
JP6264361B2 (en) * | 2015-11-30 | 2018-01-24 | コニカミノルタ株式会社 | Dynamic image analyzer |
JP2017169830A (en) * | 2016-03-24 | 2017-09-28 | コニカミノルタ株式会社 | Dynamic analysis apparatus |
JP6790537B2 (en) * | 2016-07-15 | 2020-11-25 | コニカミノルタ株式会社 | Dynamic analyzer |
CN108471999B (en) * | 2016-12-23 | 2019-07-02 | 皇家飞利浦有限公司 | The system for being used to help the state of the lung of assessment object |
JP6436222B2 (en) * | 2017-12-21 | 2018-12-12 | コニカミノルタ株式会社 | Dynamic image analyzer |
JP2019171102A (en) * | 2019-05-30 | 2019-10-10 | コニカミノルタ株式会社 | Image processing device, image processing method, and program |
JP7345282B2 (en) | 2019-06-03 | 2023-09-15 | キヤノン株式会社 | Image processing device, image processing method and program |
CN110766735B (en) * | 2019-10-21 | 2020-06-26 | 北京推想科技有限公司 | Image matching method, device, equipment and storage medium |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002306483A (en) * | 2001-04-18 | 2002-10-22 | Toshiba Corp | Medical diagnosis image processing equipment and method thereof |
WO2004044612A2 (en) * | 2002-11-12 | 2004-05-27 | Accuray, Inc. | A method and apparatus for tracking an internal target region without an implanted fiducial |
JP2005020719A (en) * | 2003-06-06 | 2005-01-20 | Fuji Photo Film Co Ltd | Diagnostic image reading supporting method and apparatus, and program |
JP2005151099A (en) * | 2003-11-14 | 2005-06-09 | Canon Inc | Device and method for processing radiation image |
-
2006
- 2006-06-13 WO PCT/JP2006/311847 patent/WO2006137294A1/en active Application Filing
- 2006-06-13 JP JP2007522243A patent/JP4797173B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002306483A (en) * | 2001-04-18 | 2002-10-22 | Toshiba Corp | Medical diagnosis image processing equipment and method thereof |
WO2004044612A2 (en) * | 2002-11-12 | 2004-05-27 | Accuray, Inc. | A method and apparatus for tracking an internal target region without an implanted fiducial |
JP2005020719A (en) * | 2003-06-06 | 2005-01-20 | Fuji Photo Film Co Ltd | Diagnostic image reading supporting method and apparatus, and program |
JP2005151099A (en) * | 2003-11-14 | 2005-06-09 | Canon Inc | Device and method for processing radiation image |
Also Published As
Publication number | Publication date |
---|---|
WO2006137294A1 (en) | 2006-12-28 |
JPWO2006137294A1 (en) | 2009-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4797173B2 (en) | X-ray diagnosis support apparatus, program, and recording medium | |
JP5400326B2 (en) | Method for displaying tomosynthesis images | |
JP7462898B2 (en) | Diagnostic Support Program | |
US10540764B2 (en) | Medical image capturing apparatus and method | |
US8755565B2 (en) | X-ray CT apparatus and image display method of X-ray CT apparatus | |
JP5295562B2 (en) | Flexible 3D rotational angiography-computed tomography fusion method | |
JP2008537691A (en) | How to expand the field of imaging software in diagnostic workups | |
CN109419526A (en) | Method and system for locomotion evaluation and correction in the synthesis of digital breast tomography | |
JP6690774B2 (en) | Dynamic analysis system, program and dynamic analysis device | |
JP2018068400A (en) | Dynamic image processing device | |
JP6253085B2 (en) | X-ray moving image analysis apparatus, X-ray moving image analysis program, and X-ray moving image imaging apparatus | |
JP5839822B2 (en) | Image processing apparatus and X-ray CT apparatus | |
JP6662428B2 (en) | Dynamic analysis system | |
EP3545845B1 (en) | Method for processing computed tomography imaging data of a suspect`s respiratory system | |
US8391573B2 (en) | Method and apparatus for motion correcting medical images | |
JP2010154982A (en) | X-ray computer tomographic imaging apparatus and image processor | |
JP5157919B2 (en) | Image display apparatus and program | |
JP5500216B2 (en) | Image display apparatus and program | |
JP2017169830A (en) | Dynamic analysis apparatus | |
US20230329662A1 (en) | X-ray ct apparatus, image processing apparatus, and motion-corrected image reconstruction method | |
JP2019030478A (en) | Medical image diagnostic device and image processing method | |
US20220398720A1 (en) | Diagnostic support program | |
JPH07194583A (en) | X-ray diagnostic system | |
JP2004313513A (en) | X-ray ct apparatus | |
JP4503279B2 (en) | X-ray CT system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090609 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20090609 |
|
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: 20110705 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |