WO2012026146A1 - 胸部診断支援システム及びプログラム - Google Patents
胸部診断支援システム及びプログラム Download PDFInfo
- Publication number
- WO2012026146A1 WO2012026146A1 PCT/JP2011/055524 JP2011055524W WO2012026146A1 WO 2012026146 A1 WO2012026146 A1 WO 2012026146A1 JP 2011055524 W JP2011055524 W JP 2011055524W WO 2012026146 A1 WO2012026146 A1 WO 2012026146A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- value
- image
- frame
- frame images
- imaging
- Prior art date
Links
- 238000003745 diagnosis Methods 0.000 title claims abstract description 42
- 210000000115 thoracic cavity Anatomy 0.000 title abstract 2
- 238000003384 imaging method Methods 0.000 claims abstract description 160
- 210000004072 lung Anatomy 0.000 claims abstract description 126
- 210000000038 chest Anatomy 0.000 claims abstract description 40
- 230000003434 inspiratory effect Effects 0.000 claims abstract description 37
- 238000004458 analytical method Methods 0.000 claims description 64
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 30
- 230000000241 respiratory effect Effects 0.000 claims description 20
- 238000000605 extraction Methods 0.000 claims description 4
- 210000000481 breast Anatomy 0.000 claims 1
- 238000009423 ventilation Methods 0.000 abstract description 13
- 230000033001 locomotion Effects 0.000 abstract description 6
- 239000000284 extract Substances 0.000 abstract description 3
- 230000005855 radiation Effects 0.000 description 171
- 238000000034 method Methods 0.000 description 67
- 238000012545 processing Methods 0.000 description 50
- 238000012937 correction Methods 0.000 description 45
- 238000001514 detection method Methods 0.000 description 43
- 238000003860 storage Methods 0.000 description 40
- 238000004891 communication Methods 0.000 description 30
- 230000002950 deficient Effects 0.000 description 24
- 238000006243 chemical reaction Methods 0.000 description 21
- 201000010099 disease Diseases 0.000 description 20
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 20
- 238000010191 image analysis Methods 0.000 description 16
- 230000000284 resting effect Effects 0.000 description 15
- 208000011623 Obstructive Lung disease Diseases 0.000 description 12
- 208000019693 Lung disease Diseases 0.000 description 10
- 230000002159 abnormal effect Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 230000000414 obstructive effect Effects 0.000 description 10
- 239000004065 semiconductor Substances 0.000 description 9
- 238000009826 distribution Methods 0.000 description 7
- 230000005856 abnormality Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 5
- 238000012935 Averaging Methods 0.000 description 4
- 230000002238 attenuated effect Effects 0.000 description 4
- 230000001678 irradiating effect Effects 0.000 description 4
- 230000004199 lung function Effects 0.000 description 4
- 230000001575 pathological effect Effects 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 239000000758 substrate Substances 0.000 description 4
- 238000012790 confirmation Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000002601 radiography Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 238000004040 coloring Methods 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 239000010408 film Substances 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 230000004660 morphological change Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 239000010409 thin film Substances 0.000 description 2
- YDLQKLWVKKFPII-UHFFFAOYSA-N timiperone Chemical compound C1=CC(F)=CC=C1C(=O)CCCN1CCC(N2C(NC3=CC=CC=C32)=S)CC1 YDLQKLWVKKFPII-UHFFFAOYSA-N 0.000 description 2
- 229950000809 timiperone Drugs 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Classifications
-
- 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
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/113—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb occurring during breathing
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4035—Arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
-
- 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/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
-
- 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
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4233—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using matrix detectors
-
- 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/56—Details of data transmission or power supply, e.g. use of slip rings
- A61B6/563—Details of data transmission or power supply, e.g. use of slip rings involving image data transmission via a network
-
- 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/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/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- 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
Definitions
- the present invention relates to a chest diagnosis support system and program.
- a dynamic image of the chest is taken using a semiconductor image sensor such as FPD (flat panel detector), Attempts have been made to apply it to diagnosis.
- FPD flat panel detector
- pulsed radiation is continuously irradiated from the radiation source in accordance with the reading / erasing timing of the semiconductor image sensor. Take multiple shots per second to capture the dynamics of the chest.
- a doctor can observe a series of movements of the chest accompanying respiratory motion, heart beats, and the like.
- Patent Document 1 a difference image indicating a difference in pixel values between a plurality of frame images constituting a dynamic image is generated, and the maximum value (average value, intermediate value, etc.) of each pixel value is generated from the generated difference image. ) Is displayed.
- An object of the present invention is to facilitate understanding of a pathological condition related to lung field ventilation and to provide information useful for diagnosis to a doctor.
- it is to provide diagnostic information that enables accurate diagnosis even by a doctor who has little practical experience with a stethoscope.
- the chest diagnosis support system Photographing means for photographing the dynamics of the chest including at least one respiratory cycle, and generating a plurality of continuous frame images; Extracting means for extracting a lung field region from each of a plurality of frame images generated by the imaging means; A region dividing unit that divides the lung field region extracted by the extracting unit into a plurality of small regions and associates the small regions across the plurality of frame images; Analysis is performed for each small region associated between the plurality of frame images to calculate a feature amount of inspiration and a feature amount of exhalation for each of the small regions, and the calculated feature amount of inspiration and the feature of exhalation An analysis means for calculating a ratio value with respect to a quantity and creating a histogram of the calculated ratio value; Display means for displaying the histogram created by the analyzing means; Is provided.
- the feature amount of the inspiration is the maximum absolute value of the inter-frame difference value obtained by taking the difference of the pixel signal values between temporally adjacent frame images of the frame image group corresponding to the inspiratory period
- the feature quantity of expiration is preferably the maximum absolute value of the inter-frame difference value obtained by taking the difference of pixel signal values between temporally adjacent frame images of the frame image group corresponding to the expiration period.
- the feature quantity of inspiration is the number of frame images corresponding to the inspiratory period in one respiratory cycle
- the feature quantity of exhalation is the number of frame images corresponding to the exhalation period in one respiratory cycle. preferable.
- the feature quantity of inspiration is a maximum value of pixel signal values in a plurality of frame images for one breath cycle
- the feature quantity of breath is a pixel signal value in a plurality of frame images for one breath cycle.
- a minimum value is preferred.
- the region dividing means divides a lung field region of one frame image among the plurality of frame images into a plurality of small regions, and the lung field regions of the other frame images are the same as the small regions of the one frame image. It is preferable to divide into small regions at pixel positions and associate the small regions at the same pixel position between the plurality of frame images with each other.
- the region dividing unit divides a lung field region of one frame image among the plurality of frame images into a plurality of small regions, and each of the other frame images has a matching degree with each of the small regions in the one frame image. It is preferable that a high area is extracted, a lung field area of the other frame image is divided into the extracted areas, and the small areas having a high degree of matching between the plurality of frame images are associated with each other.
- the photographing means performs dynamic photographing without using a grid.
- the program Computer Extraction means for extracting a lung field region from each of a plurality of frame images showing the dynamics of the chest including at least one respiratory cycle; A region dividing unit for dividing the lung field region extracted by the extracting unit into a plurality of small regions and associating the small regions between the plurality of frame images; Analysis is performed for each small region associated between the plurality of frame images to calculate a feature amount of inspiration and a feature amount of exhalation for each of the small regions, and the calculated feature amount of inspiration and the feature of exhalation An analysis means for calculating a ratio value with respect to a quantity and creating a histogram of the calculated ratio value; Display means for displaying a histogram created by the analyzing means; To function as.
- the analysis result (pathologic information) relating to lung field ventilation is provided in a visible manner, so that it is possible to grasp the pathologic condition independent of the degree of experience of a doctor, and early treatment based on early detection of a lesioned part Is possible.
- FIG. 1 is a diagram illustrating an overall configuration of a chest diagnosis support system according to a first embodiment of the present invention.
- FIG. 1 shows the overall configuration of a chest diagnosis support system 100 in the present embodiment.
- an imaging device 1 and an imaging console 2 are connected by a communication cable or the like, and the imaging console 2 and the diagnostic console 3 are connected to a LAN (Local Area Network). Etc., and connected via a communication network NT.
- LAN Local Area Network
- Each device constituting the chest diagnosis support system 100 conforms to the DICOM (Digital Image and Communications in Medicine) standard, and communication between the devices is performed in accordance with DICOM.
- DICOM Digital Image and Communications in Medicine
- the imaging apparatus 1 is an apparatus that images the dynamics of the chest with periodicity (cycle), such as pulmonary expansion and contraction morphological changes, heart pulsation, and the like accompanying respiratory motion. Dynamic imaging is performed by continuously irradiating a human chest with radiation such as X-rays to acquire a plurality of images (that is, continuous imaging). A series of images obtained by this continuous shooting is called a dynamic image. Each of the plurality of images constituting the dynamic image is called a frame image. As shown in FIG. 1, the imaging apparatus 1 includes a radiation source 11, a radiation irradiation control device 12, a radiation detection unit 13, a reading control device 14, and the like.
- the radiation source 11 is disposed at a position facing the radiation detection unit 13 across the subject M, and irradiates the subject M with radiation (X-rays) according to the control of the radiation irradiation control device 12.
- the radiation irradiation control device 12 is connected to the imaging console 2 and controls the radiation source 11 based on the radiation irradiation conditions input from the imaging console 2 to perform radiation imaging.
- the radiation irradiation conditions input from the imaging console 2 are, for example, pulse rate, pulse width, pulse interval during continuous irradiation, number of imaging frames per imaging, X-ray tube current value, X-ray tube voltage value, Filter type and the like.
- the pulse rate is the number of times of radiation irradiation per second, and matches the frame rate described later.
- the pulse width is a radiation irradiation time per one irradiation.
- the pulse interval is the time from the start of one radiation irradiation to the start of the next radiation irradiation in continuous imaging, and coincides with a frame interval described later.
- the radiation detection unit 13 is configured by a semiconductor image sensor such as an FPD.
- the FPD has, for example, a glass substrate or the like, detects radiation that has been irradiated from the radiation source 11 and transmitted through at least the subject M at a predetermined position on the substrate according to its intensity, and detects the detected radiation as an electrical signal.
- a plurality of pixels to be converted and stored are arranged in a matrix.
- Each pixel is configured by a switching unit such as a TFT (Thin Film Transistor).
- the reading control device 14 is connected to the imaging console 2.
- the reading control device 14 controls the switching unit of each pixel of the radiation detection unit 13 based on the image reading condition input from the imaging console 2 to switch the reading of the electrical signal accumulated in each pixel.
- the image data is acquired by reading the electrical signal accumulated in the radiation detection unit 13.
- This image data is a frame image.
- the read control unit 14 outputs the frame image acquired in imaging console 2.
- the image reading conditions are, for example, a frame rate, a frame interval, a pixel size, an image size (matrix size), and the like.
- the frame rate is the number of frame images acquired per second and matches the pulse rate.
- the frame interval is the time from the start of one frame image acquisition operation to the start of the next frame image acquisition operation in continuous shooting, and coincides with the pulse interval.
- the radiation irradiation control device 12 and the reading control device 14 are connected to each other, and exchange synchronization signals with each other to synchronize the radiation irradiation operation and the image reading operation.
- the time of calibration for acquiring a plurality of dark images for calculating an offset correction coefficient used for offset correction which will be described later, it is not synchronized with the radiation irradiation operation, and reset, accumulation, and data reading are performed without radiation.
- the reading operation of a series of reset images is performed, but may be performed at any timing before or after a series of dynamic imaging.
- the imaging console 2 outputs radiation irradiation conditions and image reading conditions to the imaging apparatus 1 to control radiation imaging and radiographic image reading operations by the imaging apparatus 1, and also captures dynamic images acquired by the imaging apparatus 1. Displayed for confirmation of whether the image is suitable for confirmation of positioning or diagnosis.
- the imaging console 2 includes a control unit 21, a storage unit 22, an operation unit 23, a display unit 24, and a communication unit 25, and each unit is connected by a bus 26.
- the control unit 21 includes a CPU (Central Processing Unit), a RAM (Random Access Memory), and the like.
- the CPU of the control unit 21 reads the system program and various processing programs stored in the storage unit 22 in accordance with the operation of the operation unit 23, expands them in the RAM, and performs shooting control processing described later according to the expanded programs.
- Various processes including the beginning are executed to centrally control the operation of each part of the imaging console 2 and the radiation irradiation operation and the reading operation of the imaging apparatus 1.
- the storage unit 22 is configured by a nonvolatile semiconductor memory, a hard disk, or the like.
- the storage unit 22 stores various programs executed by the control unit 21 and data such as parameters necessary for execution of processing by the programs or processing results.
- the storage unit 22 stores a shooting control processing program for executing the shooting control process shown in FIG.
- the storage unit 22 stores radiation irradiation conditions and image reading conditions in association with the examination target region.
- Various programs are stored in the form of readable program code, and the control unit 21 sequentially executes operations according to the program code.
- the operation unit 23 includes a keyboard having a cursor key, numeric input keys, various function keys, and the like, and a pointing device such as a mouse.
- the control unit 23 controls an instruction signal input by key operation or mouse operation on the keyboard.
- the operation unit 23 may include a touch panel on the display screen of the display unit 24. In this case, the operation unit 23 outputs an instruction signal input via the touch panel to the control unit 21.
- the display unit 24 includes a monitor such as an LCD (Liquid Crystal Display) or a CRT (Cathode Ray Tube), and displays input instructions, data, and the like from the operation unit 23 in accordance with display signal instructions input from the control unit 21. To do.
- a monitor such as an LCD (Liquid Crystal Display) or a CRT (Cathode Ray Tube)
- LCD Liquid Crystal Display
- CRT Cathode Ray Tube
- the communication unit 25 includes a LAN adapter, a modem, a TA (Terminal Adapter), and the like, and controls data transmission / reception with each device connected to the communication network NT.
- the diagnostic console 3 is a moving image processing apparatus for acquiring a dynamic image from the imaging console 2 and displaying diagnostic support information such as the acquired dynamic image and a histogram to be described later for a diagnostic interpretation by a doctor.
- the diagnostic console 3 includes a control unit 31, a storage unit 32, an operation unit 33, a display unit 34, and a communication unit 35, and each unit is connected by a bus 36.
- the control unit 31 includes a CPU, a RAM, and the like.
- the CPU of the control unit 31 reads out the system program and various processing programs stored in the storage unit 32 in accordance with the operation of the operation unit 33 and expands them in the RAM, and performs image analysis described later according to the expanded programs.
- Various processes including the process are executed to centrally control the operation of each part of the diagnostic console 3.
- the control unit 31 realizes functions as an extraction unit, a region division unit, and an analysis unit by executing an image analysis process described later.
- the storage unit 32 is configured by a nonvolatile semiconductor memory, a hard disk, or the like.
- the storage unit 32 stores data such as parameters necessary for execution of processing or data such as processing results by various programs and programs including an image analysis processing program for executing image analysis processing by the control unit 31. These various programs are stored in the form of readable program codes, and the control unit 31 sequentially executes operations according to the program codes.
- the operation unit 33 includes a keyboard having cursor keys, numeric input keys, various function keys, and the like, and a pointing device such as a mouse.
- the control unit 33 controls an instruction signal input by key operation or mouse operation on the keyboard.
- the operation unit 33 may include a touch panel on the display screen of the display unit 34, and in this case, an instruction signal input via the touch panel is output to the control unit 31.
- the display unit 34 is configured by a monitor such as an LCD or a CRT, and displays an input instruction, data, or the like from the operation unit 33 in accordance with an instruction of a display signal input from the control unit 31.
- the communication unit 35 includes a LAN adapter, a modem, a TA, and the like, and controls data transmission / reception with each device connected to the communication network NT.
- FIG. 2 shows photographing control processing executed in the control unit 21 of the photographing console 2.
- the photographing control process is executed in cooperation with the photographing control processing program stored in the control unit 21 and the storage unit 22.
- the operation section 23 of the imaging console 2 is operated by the imaging technician, and patient information (patient name, height, weight, age, sex, etc.) of the imaging target (subject M) is input (step S1).
- the radiation irradiation conditions are read from the storage unit 22 and set in the radiation irradiation control device 12, and the image reading conditions are read from the storage unit 22 and set in the reading control device 14 (step S2).
- a radiation irradiation instruction by the operation of the operation unit 23 is waited (step S3).
- the imaging operator prepares for imaging such as positioning of the subject M in the imaging device 1 and also instructs the subject (subject M) to take comfort in order to capture the dynamics of resting breathing to encourage resting breathing. .
- the shooting is ready to enter the radiation irradiation instruction by operating the operation unit 23.
- a photographing start instruction is output to the radiation irradiation control device 12 and the reading control device 14, and dynamic photographing is started (step S4). That is, radiation is emitted from the radiation source 11 at a pulse interval set in the radiation irradiation control device 12, and a frame image is acquired by the radiation detection unit 13.
- the control unit 21 outputs a photographing end instruction to the radiation irradiation control device 12 and the reading control device 14, and the photographing operation is stopped.
- the number of frames to be captured is the number of frames that can be captured for at least one respiratory cycle.
- the frame images acquired by shooting are sequentially input to the shooting console 2, and correction processing is performed on each frame image (step S5).
- three correction processes are performed: an offset correction process, a gain correction process, and a defective pixel correction process.
- an offset correction process is performed on each acquired frame image, and an offset value caused by a dark current superimposed on each acquired frame image is removed.
- the offset correction process for example, a process of subtracting an offset correction coefficient stored in advance from each pixel value (density value; hereinafter referred to as a signal value) of each acquired frame image is performed.
- the offset correction coefficient is an image obtained by averaging a plurality of frame images acquired in advance when radiation is not irradiated.
- a gain correction process is performed, and pixel-to-pixel variations caused by individual differences between detection elements corresponding to the pixels of each frame image and gain unevenness of the readout amplifier are removed.
- the gain correction process for example, a process of multiplying each frame image after the offset correction by a gain correction coefficient stored in advance is performed.
- the gain correction coefficient is an image obtained by averaging a plurality of offset-corrected frame images acquired when the radiation detector 13 is uniformly irradiated with radiation, and an output signal value expected under the radiation irradiation conditions at this time.
- the coefficient is calculated and stored in advance so that the corrected signal value of each pixel becomes uniform.
- defective pixel correction processing is performed to remove pixels whose sensitivity is non-linear compared to surrounding pixels and missing pixels that have no sensitivity.
- the defective pixel correction processing for example, in accordance with the defective pixel position information map stored in advance, in each defective pixel registered in the defective pixel position information map, the signal value of the defective pixel is averaged of the signal values of the pixels that are not defective in the vicinity thereof. Processing to replace with a value is performed.
- the defective pixel position information map a plurality of defective pixels are recognized in advance from the offset-corrected and gain-corrected frame images acquired when the radiation detecting unit 13 is uniformly irradiated with radiation, and the positions of the defective pixels are detected. Is a registered map.
- the offset correction coefficient, gain correction coefficient, and defective pixel position information map optimum values are stored in advance according to collection modes such as binning and dynamic range, and optimum values corresponding to the respective collection modes are stored. Is read out.
- each frame image after the correction process is stored in the storage unit 22 in association with the number indicating the shooting order (step S6) and displayed on the display unit 24 (step S7).
- the signal value of each pixel of each frame image may be stored after being subjected to logarithmic conversion processing for converting from a true number to a logarithm.
- the imaging engineer confirms the positioning and the like based on the displayed dynamic image, and determines whether an image suitable for diagnosis is acquired by imaging (imaging OK) or re-imaging is necessary (imaging NG). Then, the operation unit 23 is operated to input a determination result. Note that the frame images acquired by shooting may be input together after the completion of all shooting.
- step S8 When a determination result indicating that the shooting is OK is input by a predetermined operation of the operation unit 23 (step S8; YES), an identification ID for identifying the dynamic image or each of a series of frame images acquired by the dynamic shooting is displayed.
- Information such as patient information, examination target region, radiation irradiation condition, image reading condition, imaging order number, imaging date and time are attached (for example, written in the header area of the image data in DICOM format), and the communication unit 25 is To the diagnostic console 3 (step S9). Then, this process ends.
- step S8 when a determination result indicating photographing NG is input by a predetermined operation of the operation unit 23 (step S8; NO), a series of frame images stored in the storage unit 22 is deleted (step S10), and this processing is performed. finish. In this case, re-photographing is executed.
- the diagnostic console 3 receives a series of frame images of the dynamic image from the imaging console 2 via the communication unit 35, it cooperates with the control unit 31 and the image analysis processing program stored in the storage unit 32. As a result, the image analysis processing shown in FIG. 3 is executed.
- a lung field region is extracted from each frame image (step S21). Any method may be used to extract the lung field region. For example, a threshold value is obtained by discriminant analysis from a histogram of signal values (density values) of each pixel of an arbitrary frame image (here, the first (first) frame image in the shooting order) in a series of frame images. A region having a signal higher than this threshold is first extracted as a lung field region candidate. Next, edge detection is performed in the vicinity of the boundary of the first extracted lung field region candidate, and the boundary of the lung field region can be extracted by extracting along the boundary the point where the edge is maximum in a small region near the boundary. it can.
- the lung field region of each frame image is divided into a plurality of small regions, and the small regions of each frame image are associated with each other (step S22).
- the position of the pixel in each small area is stored in the RAM of the control unit 31.
- T time phases
- An image D1 in FIG. 5 is a frame image of a resting expiratory position (timing at which the diaphragm position is highest during rest breathing).
- An image D2 in FIG. 5 is a frame image of a rest inspiratory position (timing at which the diaphragm position is lowest during rest breathing).
- the images D1 and D2 in FIG. 5 are images taken at the timing with the largest shape difference in one respiration cycle.
- the positional displacement is slight even in the lower region of the lung field region where the positional displacement is the largest (A11 in the image D2 has the same pixel position as A1 in the image D1).
- A2 in the image D2 indicates a region in which the same portion in the lung field as the image A1 in the image D1 is drawn).
- one frame image is set as a reference image from a series of frame images.
- the lung field region from which the reference image is extracted is divided into a plurality of small regions (for example, a 2 mm ⁇ 2 mm rectangular region) (see FIG. 5).
- the lung field regions of the other frame images are divided into small regions (regions of signal values output from the same image sensor of the radiation detection unit 13) at the same pixel positions as the small regions of the reference image.
- the small areas at the same pixel position between the frame images are associated with each other. In this process, it is possible to perform division and association of the frame image into small areas at high speed.
- the reference image is preferably a frame image at a resting expiratory position.
- the position of the diaphragm is the highest during rest breathing, i.e., the area of the lung field region is the smallest, so when the small region of the reference image is associated with another frame image, the small region is This is because it is not associated with an area outside the lung field of another frame image.
- An image of the resting breath position can be acquired by extracting an image having the highest diaphragm position from a series of frame images.
- the diaphragm reference position B1 shown in FIG. 6 is defined in advance as an average position in the vertical direction of the diaphragm curve C (shown by a dotted line in FIG. 6). C (the lower end of the lung field region) is extracted, the average position in the vertical direction is obtained, and the obtained position is specified as the reference position B1 of the diaphragm.
- the captured dynamic image is an image during deep breathing, as shown in FIG. 7, the pixel position indicating the same portion of the lung field is greatly shifted.
- the image D3 in FIG. 7 is a frame image at the maximum expiratory level during deep breathing
- the image D4 in FIG. 7 is a frame image at the maximum inspiratory position during deep breathing (A31 in the image D4 is A3 in the image D3).
- A4 in the image D4 indicates a region in which the same portion of the lung field as the A3 in the image D3 is drawn).
- a lung field region extracted from the first (first) frame image in the imaging order is divided into small regions made of, for example, a 2 mm square.
- P1 is a frame image having the first shooting order
- P2 is a frame image adjacent to the frame image (a frame image having a shooting order is adjacent (that is, a temporally adjacent frame image; the same applies hereinafter)).
- a search area for each small area of P1 is set.
- the search area of P2 has the same center point (x, y), and the vertical and horizontal widths are smaller than those of the small area of P1, where the coordinates of the center point in each small area of P1 are (x, y). Set to be large (for example, 1.5 times).
- a corresponding position on P2 for each small region of P1 is calculated by obtaining a region having the highest matching degree in the search range of P2.
- the matching degree the least square method or the cross correlation coefficient is used as an index.
- the lung field area of P2 is divided at the corresponding position of each small area of P1.
- P2 is newly regarded as P1
- the next frame image whose shooting order is P2 is regarded as new P2
- the corresponding position of P2 in each small region of P1 is calculated.
- P1 is a frame image with the first shooting order
- P2 is a frame image that is adjacent (temporally adjacent) to the shooting order, and between adjacent frame images calculated by the local matching process.
- a shift vector from P1 to P2 is calculated for each small region.
- the calculated shift vector is fitted with a polynomial, and the shift vector of each pixel in each small region is calculated using this polynomial.
- warping processing is performed based on the calculated shift vector of each pixel, and the position of each pixel in each small area of P2 is shifted to the position of the corresponding pixel in the frame image of P1.
- the warping process P2 is newly regarded as P1, and the next frame image whose shooting order is P2 is regarded as a new P2, and the above process is performed.
- the positions of the small areas of all the frame images can be made to substantially coincide with the frame image (reference image in this process) whose shooting order is 1. It becomes possible.
- the correspondence of the position of the small area between the frame images is stored in the RAM of the control unit 31.
- an analysis is performed for each small region associated in step S22 of the series of frame images, and an inspiration feature value and an expiration feature value are calculated (step S23).
- a representative value (maximum absolute value) of inter-frame difference values (differential values) for each small region in each of the expiration period and the inspiration period is calculated.
- the inter-frame difference value is a value indicating the signal change amount at the timing when the frame image is captured. When breathing in and out through breathing, the density of the lungs changes according to the flow of the breath, and the X-ray transmission amount (that is, the output signal value of the pixel) changes accordingly. Therefore, the signal variation can be regarded as a value indicating the air velocity at that time.
- the representative value is not limited to the maximum absolute value, and may be a median value, an average value, or a mode value.
- step S24 the value of the ratio of the inhalation feature quantity and the expiration feature quantity (inspiration feature quantity / exhalation feature quantity) for each small region is calculated (step S24).
- “the maximum value of the inter-frame difference value in the inspiration period / the maximum value of the inter-frame difference value in the expiration period” is calculated.
- index value here, mean, standard deviation
- the “inspiration feature amount / expiration feature amount” value obtained for each small area is stored in advance in the storage unit 32, and the “inspiration feature amount / expiration feature amount” value is displayed.
- Is converted into a display parameter value based on the parameter value conversion table and an image displaying each small region of a reference image (for example, a frame image of a resting breath position) is created using the converted parameter value (Ste S26).
- the conversion table includes, for example, threshold values (threshold values for each category) and hues that define the range of the size of the feature amount of each category when the feature amount is classified into each category of normal / abnormal (severity 1 to n).
- the conversion table of the parameter values at the time of display it is preferable to associate the threshold value of each category with the hue in order to enhance the discrimination ability with respect to the size of the feature amount.
- several (for example, 5 to 6) hues are associated with the threshold values of the respective categories, and further, intermediate hues are assigned (graded) to the feature value between them.
- High-performance display can be achieved.
- the image colored based on the display parameter value may be displayed as an overlay on the frame image of the reference image.
- the created histogram, the created still image, and the like are displayed in parallel on the display unit 34 (step S27), and the image analysis process ends.
- the area of the histogram is colored according to the same standard as that of each small area of the lung field area of the reference image, based on the conversion table of the above-described “inspiration feature quantity / expiration feature quantity” value and the parameter value at the time of display. Is displayed.
- FIG. 8A to 8C show examples of display screens displayed on the display unit 34 in step S27.
- FIG. 8A is a display screen displaying an analysis result obtained by analyzing a dynamic image of a lung field of a normal person.
- FIG. 8B is a display screen displaying an analysis result obtained by analyzing a dynamic image of a lung field of COPD (obstructive disease).
- FIG. 8C is a display screen displaying an analysis result obtained by analyzing a dynamic image of a lung field of a mixed disease.
- step S27 “inspiratory feature quantity / exhalation feature quantity” of each small region (2 mm square size) in the lung field region extracted from the half-cut size frame image. ”Value histogram 34a, still image 34b displaying a list of“ inspiration feature quantity / expiration feature quantity ”of each small region, hues displayed in histogram 34a and still image 34b, and“ inspiration feature quantity ”
- a display 34c indicating a relationship with the value of “/ expired feature value” and an index value 34d indicating a tendency of “inspired feature value / exhaled feature value” in the entire lung field are displayed. It is possible to grasp the distribution characteristics of “inhalation feature quantity / expiration feature quantity”.
- the histogram 34a is displayed with the horizontal axis color-coded in five hues according to the magnitude of the value of “inhalation feature quantity / expiration feature quantity”.
- each small region has the same standard as the color coding of the histogram according to the value of “inspiration feature amount / expiration feature amount”. Therefore, it is possible for a doctor to easily grasp a local abnormal part (occlusion part, restriction part) in the lung field.
- the index value 34d indicating the tendency of “inhalation feature quantity / expiration feature quantity” in the entire lung field is calculated and the average value and standard deviation are calculated and displayed together on the screen.
- the tendency of “inhalation feature quantity / expiration feature quantity” in can be provided as a numerical value to a doctor.
- the maximum value (absolute value) of the inter-frame difference value in the expiration period is the exhalation feature amount
- the maximum value (absolute value) of the inter-frame difference value in the inspiration phase is the inspiration feature amount
- the lungs of a normal person it is known that the average value of “inhalation feature quantity / expiration feature quantity” in the entire lung field is about 0.9 to 1.2, and the standard deviation is about 0.10 to 0.22. Therefore, when the display screen shown in FIG. 8A is displayed in step S27, the doctor can easily grasp that the photographed lung field is normal.
- the average value of “inhalation feature value / expiration feature value” of the entire lung field does not fall within 0.9 to 1.2 (larger than normal subjects) It is known that the standard deviation does not fall within 0.10 to 0.22 (larger than that of a normal person). Therefore, when the display screen shown in FIG. 8B is displayed in step S27, the doctor can easily grasp that the photographed lung field is COPD.
- the value of “inspiratory feature value / exhalation feature value” of the entire lung field increases both in the number of data of 0.66 or less and the number of data of 1.5 or more. It has been known. Therefore, when the display screen shown in FIG.
- the doctor can easily grasp that the photographed lung field is a mixed disease.
- the chest diagnosis support system 100 can provide doctors with useful diagnostic information that can identify pathologies such as COPD (obstructive pulmonary disease) and mixed lung disease.
- the feature quantity of exhalation and the feature quantity of inspiration may be used.
- the feature quantity of exhalation may be the number of frame images corresponding to the expiration period in one breath cycle (expiration time)
- the feature quantity of inspiration may be the number of frame images corresponding to the inspiration period in one breath cycle (inspiration period).
- the inspiratory time and the expiratory time should be approximately the same length or slightly longer. Therefore, by looking at the value of “the number of frame images corresponding to the expiration period / the number of frame images corresponding to the inspiration period”, the doctor can grasp whether or not there is a suspicion of lung disease.
- the region where “the number of frame images during the expiration period / the number of frame images during the inspiration period”> 1.5 may be an obstructive portion that is difficult to ventilate during expiration and delays the discharge of the captured air. Recognize.
- the relationship of “maximum value of inter-frame difference value in inspiratory period / maximum value of inter-frame difference value in expiratory period ⁇ expiratory time (number of expiratory frame images) / inspiratory time (number of inspiratory frame images)” is satisfied.
- Pulmonary diseases and mixed lung diseases can be identified.
- the signal value (average signal value) of the pixels in each small region is calculated, and the minimum value and the maximum value of the signal value in one respiratory cycle are determined for each small region.
- the minimum value may be the feature value of exhalation in the region
- the maximum value may be the feature value of inspiration. It can be considered that the difference between the maximum value and the minimum value of the signal value is large at a normal location, and the difference between the two is extremely small at a location where there is an abnormality. This is because the alveoli move poorly at abnormal locations, and the density change of the alveoli becomes small.
- the doctor refers to the histogram of “maximum value of signal value / minimum value of signal value” and confirms the average value and standard deviation to determine whether the lung field is normal or disease. can do. For example, when the average value of “maximum value of signal value / minimum value of signal value” of the entire lung field is greater than 1 and the standard deviation is small, it can be determined that the lung function is normal. On the other hand, if the average value of “maximum value of signal value / minimum value of signal value” of the entire lung field is a value close to 1 and the standard deviation is large, it can be determined that there is a disease in the lung function.
- the count number (number of blocks (small areas) is the peak of the histogram.
- the value of “amount / exhalation feature amount”, the peak count number (the number of blocks), or the ratio of the count number that the “inspiration feature amount / exhalation feature amount” value is not less than a predetermined value and not more than a predetermined value It may be used as Or it may create a new index value by combining these plurality of index values.
- inspiratory feature value / expiratory feature value is converted into a logarithm, and a histogram of the above-described histogram is created, an index value is calculated, and a colorized image based on a conversion table with a display parameter is converted to the logarithmically converted value. You may make it.
- the control unit 31 extracts and extracts a lung field region from each of a plurality of frame images obtained by imaging the dynamics of the chest.
- the lung field region is divided into a plurality of small regions, and each small region is associated between a plurality of frame images.
- analysis is performed for each corresponding small region of the plurality of frame images to calculate an inspiratory feature value and an expiratory feature value for each small region, and a ratio between the calculated inspiratory feature value and the exhaled feature value. value is calculated and a histogram of the calculated ratio value. Then, the created histogram is displayed on the display unit 34.
- the doctor can provide the doctor with a histogram, which is information useful for diagnosis, showing how the ratio of the inspiratory feature quantity to the expiratory feature quantity is distributed in the lung field. It becomes possible to easily grasp the pathological condition related to the ventilation of the field.
- the maximum absolute value of the inter-frame difference value between temporally adjacent frame images of the frame image group corresponding to the inspiration period is calculated, and the frame corresponding to the expiration period is calculated as the feature value of expiration.
- the number of frame images corresponding to the inspiratory period in one breath cycle is calculated as the feature quantity of inspiration, and the number of frame images corresponding to the exhalation period in one breath cycle is calculated as the feature quantity of breath.
- the maximum value of the pixel signal values in a plurality of frame images for one breath cycle is calculated as the inspiratory feature quantity, and the minimum value of the pixel signal values in the plurality of frame images for one breath cycle is used as the exhalation feature quantity.
- the doctor determines whether the captured lung field is normal or abnormal by displaying the distribution of the ratio of the calculated inhalation feature quantity and exhalation feature quantity as a histogram. Can provide useful information.
- the lung field area of one frame image among the plurality of frame images is divided into a plurality of small areas, and the lung field areas of the other frame images are
- the frame image is divided into small areas at the same pixel position as each small area of the frame image, and each small area at the same pixel position between a plurality of frame images is associated with each other, thereby dividing the frame image into small areas of the lung field and the frame image It is possible to associate small areas between them at high speed. As a result, the analysis process can be speeded up.
- the lung field region of one frame image among the plurality of frame images is divided into a plurality of small regions, and one frame image in each of the other frame images Region with a high matching degree is extracted, the lung field region of the other frame image is divided into the extracted regions, and the small regions with a high matching degree between a plurality of frame images are associated with each other. It becomes possible to perform analysis with high accuracy.
- the irradiation dose per still image shooting and the total irradiation dose related to a series of dynamic imaging is calculated so that the two generally match.
- the dose reaching the radiation detector (detector) is attenuated. Therefore, when dynamic imaging is performed under the radiation irradiation condition described in the known document 2, each frame image is displayed. The dose to reach is greatly insufficient, and the S / N ratio decreases. Therefore, when dynamic analysis is performed using a series of frame images, accurate analysis cannot be performed, and re-imaging is performed by changing the irradiation dose. As a result, the total exposure dose of the subject may increase. There is also.
- a radiation imaging system 200 that performs imaging without a grid during dynamic imaging will be described.
- FIG. 9 shows the overall configuration of the radiation imaging system 200 in the present embodiment.
- the radiation imaging system 200 is a system that functions as a chest diagnosis support system.
- the imaging device 4 and the imaging console 5 are connected by a communication cable or the like, and the imaging console 5 and the diagnostic console 6 are connected to a LAN (Local Area Network) or the like. And connected via a communication network NT2.
- Each device constituting the radiation imaging system 200 conforms to the DICOM (Digital Image and Communications in Medicine) standard, and communication between the devices is performed according to DICOM.
- DICOM Digital Image and Communications in Medicine
- the photographing device 4 is a device that performs dynamic photographing or still image photographing of the subject M by irradiating the subject M (the chest of the human body) with radiation.
- Dynamic imaging refers to acquiring a plurality of images (that is, continuous imaging) by continuously irradiating the subject M with radiation such as X-rays in a pulsed manner.
- the dynamics of the subject M having periodicity (cycle) such as morphological changes of lung expansion and contraction accompanying the respiratory motion, heart pulsation, and the like are captured.
- a series of images obtained by this continuous shooting is called a dynamic image.
- Each of the plurality of images constituting the dynamic image is called a frame image.
- Still image photography is used for diagnosis based on the density resolution of the imaging region as in the conventional film method and CR method.
- the subject M is irradiated with radiation such as X-rays once to obtain a single sheet. This refers to acquiring a still image.
- the imaging device 4 includes a radiation source 41, a radiation irradiation control device 42, a radiation detection unit 43, a reading control device 44, a holding unit 45, a grid 46, and the like.
- the radiation source 41 is a radiation generator capable of single shot and continuous shot (pulse irradiation). That is, the radiation generator is compatible with both still image capturing and dynamic image capturing.
- the radiation source 41 is disposed at a position facing the radiation detection unit 43 across the subject M, and irradiates the subject M with radiation (X-rays) according to the control of the radiation irradiation control device 42.
- the radiation irradiation control device 42 is connected to the imaging console 5 and controls the radiation source 41 based on the radiation irradiation conditions input from the imaging console 5 to perform radiation imaging.
- the radiation irradiation conditions input from the imaging console 5 are, for example, pulse rate, pulse width, pulse interval at the time of continuous irradiation, number of imaging frames per imaging, X-ray tube current value, X-ray tube voltage value, Filter type and the like.
- the pulse rate is the number of times of radiation irradiation per second, and matches the frame rate described later.
- the pulse width is a radiation irradiation time per one irradiation.
- the pulse interval is the time from the start of one radiation irradiation to the start of the next radiation irradiation in continuous imaging, and coincides with a frame interval described later.
- the radiation detection unit 43 is configured by an FPD or the like that supports dynamic shooting and still image shooting.
- the FPD has, for example, a glass substrate, detects radiation that has been irradiated from the radiation source 41 and transmitted through at least the subject M at a predetermined position on the substrate according to the intensity, and detects the detected radiation as an electrical signal.
- a plurality of pixels to be converted and stored are arranged in a matrix.
- Each pixel includes a switching unit such as a TFT (Thin Film Transistor).
- the radiation detection unit 43 is held by a holding unit 45 so as to face the radiation source 41 with the subject M interposed therebetween.
- the reading control device 44 is connected to the photographing console 5.
- the reading control device 44 controls the switching unit of each pixel of the radiation detection unit 43 based on the image reading condition input from the imaging console 5 to switch the reading of the electrical signal accumulated in each pixel.
- the image data (still image or frame image) is acquired by reading the electrical signal stored in the radiation detection unit 43. Then, the read control unit 44 outputs the image data acquired in imaging console 5.
- the image reading conditions include, for example, a frame rate, a frame interval, a pixel size (binning size), an image size (matrix size), and the like.
- the frame rate is the number of frame images acquired per second and matches the pulse rate.
- the frame interval is the time from the start of one frame image acquisition operation to the start of the next frame image acquisition operation in continuous shooting, and coincides with the pulse interval.
- the radiation irradiation control device 42 and the reading control device 44 are connected to each other, and exchange synchronization signals to synchronize the radiation irradiation operation and the image reading operation.
- the reset, accumulation, and data are not synchronized with the radiation irradiation operation and are not irradiated with radiation.
- a series of image reading operations from reading to reset is performed, but may be performed at any timing before or after a series of dynamic imaging.
- the holding unit 45 has a detector holding unit 451 and holds the radiation detection unit 43 so as to face the radiation source 41 and the subject M during imaging.
- the holding unit 45 includes a grid mounting unit 452 for mounting a grid 46 for removing scattered radiation on the subject side (radiation source 41 side) of the radiation detection unit 43. That is, the holding part 45 is configured so that the grid 46 can be attached and detached.
- the grid mounting part 452 is provided with a grid mounting detection MS (microswitch) 453 for detecting whether or not the grid 46 is mounted.
- the holding unit 45 detects the grid mounting detection.
- the detection signal of MS453 is output to the reading control device 44.
- the holding unit 45 is provided with a subject detection sensor 454 for detecting whether or not the subject M exists at a predetermined distance as shown in FIG.
- the detection signal 454 is output to the imaging console 5 via the reading control device 44.
- the imaging console 5 outputs radiation irradiation conditions and image reading conditions to the imaging device 4 to control radiation imaging and radiation image reading operations by the imaging device 4, and still images or dynamic images acquired by the imaging device 4. For example, a preview image that has been subjected to thinning binning, a processed image that has been subjected to gradation processing, and the like are appropriately created to confirm whether the image is suitable for positioning confirmation or diagnosis by a photographer. indicate.
- the imaging console 5 includes a control unit 51, a storage unit 52, an operation unit 53, a display unit 54, and a communication unit 55, and each unit is connected by a bus 56.
- the control unit 51 includes a CPU (Central Processing Unit), a RAM (Random Access Memory), and the like.
- the CPU of the control unit 51 reads the system program and various processing programs stored in the storage unit 52 in accordance with the operation of the operation unit 53, expands them in the RAM, and performs shooting control processing described later according to the expanded programs.
- Various processes including the beginning are executed to centrally control the operation of each part of the imaging console 5 and the radiation irradiation operation and reading operation of the imaging apparatus 4.
- the storage unit 52 is configured by a nonvolatile semiconductor memory, a hard disk, or the like.
- the storage unit 52 stores various programs executed by the control unit 51 and parameters such as parameters necessary for execution of processing by the program or data such as processing results.
- the storage unit 52 stores a shooting control processing program for executing the shooting control processing shown in FIGS. 12A to 12B.
- Various programs are stored in the form of readable program codes, and the control unit 51 sequentially executes operations according to the program codes.
- the storage unit 52 stores radiation irradiation conditions and image reading conditions for dynamic imaging and still image imaging.
- the operation unit 53 includes a keyboard having cursor keys, numeric input keys, various function keys, and the like, and a pointing device such as a mouse, and controls an instruction signal input by key operation or mouse operation on the keyboard. To 51. Further, the operation unit 53 may include a touch panel on the display screen of the display unit 54. In this case, the operation unit 53 outputs an instruction signal input via the touch panel to the control unit 51.
- the display unit 54 is configured by a monitor such as an LCD (Liquid Crystal Display) or a CRT (Cathode Ray Tube), and displays input instructions, data, and the like from the operation unit 53 in accordance with display signal instructions input from the control unit 51. To do.
- a monitor such as an LCD (Liquid Crystal Display) or a CRT (Cathode Ray Tube)
- LCD Liquid Crystal Display
- CRT Cathode Ray Tube
- the communication unit 55 includes a LAN adapter, a modem, a TA (Terminal Adapter), and the like, and controls data transmission / reception with each device connected to the communication network NT2.
- the diagnostic console 6 acquires a still image or a series of frame images of a dynamic image from the imaging console 5 and displays diagnostic support information such as the acquired image and a histogram to be described later for a diagnostic interpretation by a doctor.
- the diagnostic console 6 includes a control unit 61, a storage unit 62, an operation unit 63, a display unit 64, and a communication unit 65, and each unit is connected by a bus 66.
- the control unit 61 includes a CPU, a RAM, and the like.
- the CPU of the control unit 61 reads out the system program stored in the storage unit 62 and various processing programs in accordance with the operation of the operation unit 63 and develops them in the RAM.
- Various processes including the process are executed to centrally control the operation of each part of the diagnostic console 6.
- the storage unit 62 is configured by a nonvolatile semiconductor memory, a hard disk, or the like.
- the storage unit 62 stores various types of programs including an image analysis processing program for executing image analysis processing by the control unit 61 and data such as parameters necessary for execution of processing or processing results. These various programs are stored in the form of readable program codes, and the control unit 61 sequentially executes operations according to the program codes.
- the operation unit 63 includes a keyboard having cursor keys, numeric input keys, various function keys, and the like, and a pointing device such as a mouse, and controls an instruction signal input by key operation or mouse operation on the keyboard. To 61.
- the operation unit 63 may include a touch panel on the display screen of the display unit 64. In this case, the operation unit 63 outputs an instruction signal input via the touch panel to the control unit 61.
- the display unit 64 includes a monitor such as an LCD or a CRT, and displays an input instruction, data, and the like from the operation unit 63 in accordance with an instruction of a display signal input from the control unit 61.
- the communication unit 65 includes a LAN adapter, a modem, a TA, and the like, and controls data transmission / reception with each device connected to the communication network NT2.
- FIG. 12A to FIG. 12B show photographing control processing executed in the control unit 51 of the photographing console 5.
- the imaging control process is executed in cooperation with the imaging control processing program stored in the control unit 51 and the storage unit 52.
- the operation section 53 of the imaging console 5 is operated by the imaging operator, and patient information (patient name, height, weight, age, sex, etc.) and imaging type (dynamic imaging, still image) of the imaging object (subject M) (Shooting distinction) is input (step S31).
- patient information patient name, height, weight, age, sex, etc.
- imaging type dynamic imaging, still image
- data transmitted from another device (not shown) connected to the communication network NT2 and received via the communication unit 55 can also be used as similar input information.
- step S32 it is determined whether the input shooting type is dynamic shooting or still image shooting. If it is determined that the input photographing type is dynamic photographing (step S32; YES), the radiation irradiation conditions for dynamic photographing are read from the storage unit 52 and set in the radiation irradiation control device 42, and dynamics are also obtained. Image reading conditions for photographing are read from the storage unit 52 and set in the reading control device 44 (step S33). When it is determined that the input shooting type is still image shooting (step S32; NO), radiation irradiation conditions for still image shooting are read from the storage unit 52 and set in the radiation irradiation control device 42. The image reading conditions for still image shooting are read from the storage unit 52 and set in the reading control device 44 (step S37).
- radiation irradiation conditions are set on the premise that imaging is performed without using the grid 46 for dynamic imaging (without the grid), and the grid 46 is installed for capturing a still image
- grid Radiation irradiation conditions are set on the premise of imaging (with and without).
- the radiation irradiation conditions for dynamic imaging are as follows: a dose that reaches the radiation detection unit 43 without using the grid 46 for each frame image is used, and the radiation detection unit uses the grid 46 for still image shooting. Radiation irradiation conditions are set so as to be the same as the dose reaching 43. That is, in dynamic imaging, radiation irradiation conditions are set so that the irradiation dose per sheet is lower than that in still image imaging.
- the radiation transmitted through the body is scattered in the body tissue.
- noise is generated in the radiation image. Therefore, in photographing a still image in which the absolute output value (signal value) of each pixel is important, which is used for detection of a lesioned part or observation of a lesioned part, a surface on the subject side of the radiation detection unit 43, that is, a subject is used. It is preferable that the grid 46 be provided on the surface on the side where the transmitted radiation is incident to perform imaging.
- the amount of radiation reaching the radiation detection unit 43 is attenuated by the grid 46 (for example, it is attenuated to about 1/2 in a grid with an exposure multiple of 2). It is necessary to irradiate.
- shooting using a grid has been performed as in still image shooting.
- the inventors of the present invention have used chest dynamic images obtained by dynamic imaging mainly for dynamic analysis of ventilation function and blood flow, etc. It has been found that even if a dynamic image taken without arranging is used, substantially the same result can be obtained as when a dynamic image taken with a grid is used. In other words, not only the presence / absence of the grid but also the same analysis result can be obtained as long as the dose reaching the radiation detector is the same (see FIG. 13, which will be described later in detail).
- the still image shooting is performed with the grid 46 (with a grid), and the dynamic imaging without the grid 46 (without the grid). ) Take a picture. Then, by taking a picture under a radiation irradiation condition in which the dose reaching the radiation detection unit 43 is substantially the same at the time of still image photography using the grid 46 and the photography per frame of dynamic photography without the grid 46.
- the radiation dose to the radiation detection unit 43 is configured to reduce the exposure dose of the subject M while maintaining substantially the same dose as before.
- the radiation exposure conditions can be the same as the total irradiation dose for one still image and a series of dynamic images, and still images can be taken with a grid, and dynamic images can be taken without a grid.
- the S / N ratio of each frame image may be improved as compared with the conventional method to improve the analysis accuracy.
- Step S34 it is determined whether or not the grid 46 is not attached to the grid attachment unit 452 based on the output from the grid attachment detection MS 453.
- step S34 If it is determined that the grid 46 is not attached to the grid attachment part 452 (step S34; YES), a radiation irradiation instruction by the operation of the operation part 53 is waited (step S35).
- the imaging operator instructs the subject (subject M) to take comfort in order to capture the dynamics of resting breathing, and prompts resting breathing.
- the operation unit 53 is operated to input a radiation irradiation instruction.
- a photographing start instruction is output to the radiation irradiation control device 42 and the reading control device 44, and dynamic photographing is started (step S36). That is, radiation is emitted from the radiation source 41 at a pulse interval set in the radiation irradiation control device 42, and a frame image is acquired by the radiation detection unit 43.
- the control unit 51 outputs a photographing end instruction to the radiation irradiation control device 42 and the reading control device 44, and the photographing operation is stopped.
- the number of frames to be captured is the number of frames that can be captured for at least one respiratory cycle.
- step S38 the control unit 51 performs control so that still image shooting is not performed without the grid 46 being attached.
- step S38 If it is determined that the grid 46 is mounted on the grid mounting portion 452 (step S38; YES), a radiation irradiation instruction by the operation of the operation portion 53 is waited (step S39).
- the shooter instructs the subject to stop after inhaling.
- the operation unit 53 is operated to input a radiation irradiation instruction.
- a photographing start instruction is output to the radiation irradiation control device 42 and the reading control device 44, and still image photographing is performed (step S40).
- step S41 When dynamic shooting or still image shooting is completed, images (each frame image or still image) acquired by shooting are sequentially input to the shooting console 5 and correction processing is performed (step S41).
- an offset correction process is performed on the acquired image, and the offset value due to the dark current superimposed on the acquired image is removed.
- the offset correction process for example, a process of subtracting a previously stored offset correction coefficient from each pixel value (density value; hereinafter referred to as a signal value) of the acquired image is performed.
- the offset correction coefficient is an image obtained by averaging a plurality of images acquired in advance when radiation is not irradiated.
- gain correction processing is performed, and pixel-to-pixel variations caused by individual differences between detection elements corresponding to the pixels of the acquired image and gain unevenness of the readout amplifier are removed.
- the gain correction process for example, a process of multiplying the image after offset correction by a gain correction coefficient stored in advance is performed.
- the gain correction coefficient is an image obtained by averaging a plurality of offset-corrected images acquired when the radiation detection unit 43 is uniformly irradiated with radiation, and an output signal value expected under the radiation irradiation conditions at this time.
- the coefficient is calculated and stored in advance so that the signal value of each pixel after correction is uniform.
- defective pixel correction processing is performed to remove pixels whose sensitivity is non-linear compared to surrounding pixels and missing pixels that have no sensitivity.
- the signal value of the defective pixel is averaged of the signal values of the pixels that are not defective in the vicinity thereof. Processing to replace with a value is performed.
- the defective pixel position information map a plurality of defective pixels are recognized in advance from the offset corrected and gain corrected frame images acquired when the radiation detecting unit 43 is uniformly irradiated with radiation, and the positions of the defective pixels are detected. Is a registered map.
- the offset correction coefficient, gain correction coefficient, and defective pixel position information map optimum values are stored in advance according to collection modes such as binning and dynamic range, and optimum values corresponding to the respective collection modes are stored. Is read out. These correction processes can be omitted to shorten the processing time. Moreover, it is preferable to perform a thinning process, a gradation process, etc. with respect to the acquired image as needed.
- each frame image is stored in the association storage unit 52 with a number indicating the shooting order.
- it may be stored after performing logarithmic conversion processing for converting the signal value of each pixel from a true number to a logarithm.
- the imaging operator confirms the positioning and the like from the displayed dynamic image, and determines whether an image suitable for diagnosis is acquired by imaging (imaging OK) or re-imaging is necessary (imaging NG).
- the operation unit 53 is operated to input a determination result. Note that the frame images acquired by shooting may be input together after the completion of all shooting.
- each of a still image acquired by still image shooting or a series of frame images acquired by dynamic shooting is displayed.
- Information such as identification ID for identifying an image, patient information, examination site, radiation irradiation condition, image reading condition, number indicating imaging order, imaging date and time, information indicating presence / absence of grid at imaging (grid presence / absence information), etc. Is attached (for example, written in the header area of the image data in DICOM format) and transmitted to the diagnostic console 6 via the communication unit 55 (step S45). Then, this process ends.
- step S44 when a determination result indicating photographing NG is input by a predetermined operation of the operation unit 53 (step S44; NO), a series of frame images stored in the storage unit 52 is deleted (step S46), and this processing is performed. finish. In this case, re-photographing is executed.
- a lung field region is extracted from each frame image (step S51).
- the lung field region of each frame image is divided into a plurality of small regions, and the small regions of each frame image are associated with each other (step S52). Since the specific processing of steps S51 to S52 is the same as that described in steps S21 to S22 of FIG. 3 in the first embodiment, the description is used.
- an analysis is performed for each small region associated in step S52 of the series of frame images, and an inspiration feature amount and an expiration feature amount are calculated (step S53).
- a representative value (maximum absolute value) of inter-frame difference values (differential values) for each small region in each of the expiration period and the inspiration period is calculated.
- the inter-frame difference value is a value indicating the signal change amount at the timing when the frame image is captured.
- the signal change amount can be regarded as a value indicating the air velocity at the timing.
- the representative value is not limited to the maximum absolute value, and may be a median value, an average value, or a mode value.
- signal values (average signal values) of pixels in each small region of each frame image are calculated.
- filter processing in the time axis direction is performed on each small region associated with each frame image.
- This filtering process is a process for removing a high-frequency signal change such as blood flow and extracting a time change of a signal value due to ventilation.
- Low-pass filtering is performed at a cutoff frequency of 0.7 Hz for the respiratory image group and at a cutoff frequency of 0.5 Hz for the deep respiratory image group.
- the cut-off frequency of the low-pass filter is preferably optimized for each captured dynamic image, rather than a fixed value.
- the position of the diaphragm of a series of frame images is analyzed, and in the case of rest ventilation, a frame that is a resting expiratory position and a resting inspiratory position is detected, and the frame of the resting expiratory position and the next resting inspiratory position are detected.
- the time of the inspiration period is obtained from the number of frames between the frames, and a low-pass filter is performed with a value obtained by multiplying the reciprocal by a predetermined coefficient as a cutoff frequency.
- the cut-off frequency that is automatically set is preferably limited to 0.2 to 1.0 Hz.
- step S31 vital signs such as a one-minute breathing rate and pulse rate measured at rest may be input as patient information, and the cutoff frequency may be calculated from these values.
- the low-pass filter may be implemented using a value obtained by converting a respiration rate per minute input as patient information into a respiration rate per second and multiplying the respiration rate by a predetermined coefficient.
- the input one-minute pulse rate may be converted into a one-second pulse rate, and the low-pass filter may be implemented using the average value of the one-second respiratory rate and the one-second heart rate as a cutoff frequency. .
- an inter-frame difference process for calculating a difference between the signal values of the small areas between the frame images adjacent in the shooting order.
- N + 1 ⁇ N difference values are calculated for frame images of frame numbers N and N + 1 (N is 1, 2, 3,...)
- N is 1, 2, 3,
- the maximum value of the inter-frame difference value in the expiration period absolute maximum value
- the maximum value of the inter-frame difference value in the inspiration period maximum absolute value
- the maximum value of the inter-frame difference value corresponds to the maximum differential value.
- the period in which the sign of the inter-frame difference value in each small region is positive is the inspiration period
- the negative period is the expiration period.
- step S54 the value of the ratio of the inhalation feature quantity and the exhalation feature quantity (inspiration feature quantity / exhalation feature quantity) for each small region is calculated (step S54).
- the maximum value of the inter-frame difference value in the inspiration period / the maximum value of the inter-frame difference value in the expiration period (referred to as the maximum flow velocity ratio) is calculated.
- a histogram of the value of “inhalation feature value / expiration feature value” for each calculated small area is created, and the tendency of “inspiration feature value / expiration feature value” in the entire lung field is shown.
- An index value (here, average value, standard deviation) is calculated (step S55).
- the count number on the vertical axis of the histogram is preferably normalized by dividing by the number of all small regions in the lung field.
- the value of “inhalation feature quantity / expiration feature quantity” obtained for each small area is stored in advance in the storage unit 62 and the “inspiration feature quantity / expiration feature quantity” value is displayed.
- Is converted into a display parameter value based on the parameter value conversion table and an image displaying each small region of a reference image (for example, a frame image of a resting breath position) is created using the converted parameter value (Ste S56).
- the conversion table includes, for example, threshold values (threshold values for each category) and hues that define the range of the size of the feature amount of each category when the feature amount is classified into each category of normal / abnormal (severity 1 to n).
- the conversion table of the parameter values at the time of display it is preferable to associate the threshold value of each category with the hue in order to enhance the discrimination ability with respect to the size of the feature amount.
- several (for example, 5 to 6) hues are associated with the threshold values of the respective categories, and an intermediate hue is assigned (gradation) to the feature value between them.
- the image colored based on the display parameter value may be displayed as an overlay on the frame image of the reference image.
- the storage unit 62 stores a conversion table corresponding to a dynamic image photographed with a grid and a conversion table corresponding to a dynamic image photographed without a grid, and in step S56 and the next step S57. Is determined based on grid presence / absence information attached to a series of frame images, whether the image is captured with or without the grid, and coloring is performed using a conversion table according to the determination result.
- the created histogram, the created still image, and the like are displayed in parallel on the display unit 64 (step S57), and the image analysis process ends.
- the area of the histogram is colored according to the same standard as that of each small area of the lung field area of the reference image, based on the conversion table of the above-described “inspiration feature quantity / expiration feature quantity” value and the parameter value at the time of display. Is displayed.
- FIG. 13 is a diagram showing a comparison of dynamic analysis results of dynamic images obtained by photographing a normal lung field with and without a grid.
- the analysis result shows an image colored based on the “maximum value of the inter-difference value” (maximum expiratory airflow velocity), an image colored based on the maximum flow velocity ratio, and a histogram of the maximum flow velocity ratio.
- FIG. 13 shows the analysis result of a dynamic image shot under the following shooting conditions.
- the conversion table of the maximum inspiratory airflow velocity, the maximum expiratory airflow velocity, the maximum flow velocity ratio and the color is the same (in this case, a conversion with a grid) for comparison. Table).
- the maximum inspiratory airflow velocity, the maximum expiratory airflow velocity, and the maximum flow velocity ratio are almost equal.
- FIG. 13 the shape of the histogram of the dynamic image photographed without the grid is wider than that with the grid. For this reason, for example, when the areas in the lung field and the areas of the histogram are classified into normal and abnormal categories 1 to n according to the maximum flow rate ratio, the same maximum value is used if the same threshold value (conversion table) is used with and without the grid.
- the degree of difference between the result of analyzing a dynamic image taken with a grid and the result of analyzing a dynamic image taken without a grid depends on the characteristics of the photographing system and the analysis content. Therefore, by analyzing multiple images of the same subject taken with and without the grid, the threshold used for the dynamic image with the grid and the threshold used for the dynamic image without the grid using the results are analyzed. Is preferably calculated in advance.
- the imaging device 4 since the imaging device 4 controls to capture a dynamic image without a grid, it may be sufficient to store a threshold for a dynamic image captured without a grid.
- the diagnostic console 6 is connected to an imaging system that can capture a dynamic image with a grid. In this case, if the shooting conditions for the presence / absence of the grid are different, the determination is erroneous. Therefore, in the present embodiment, grid presence / absence information is attached to the additional information of each frame image constituting the dynamic image, and the control unit 61 responds to the presence / absence of the grid at the time of shooting based on the grid presence / absence information. Analysis is performed with an analysis algorithm using the threshold value.
- FIG. 14A to 14C show examples of display screens displayed on the display unit 64 in step S57.
- FIG. 14A is a display screen displaying an analysis result obtained by analyzing a dynamic image of a lung field of a normal person photographed with a grid.
- FIG. 14B is a display screen displaying an analysis result obtained by analyzing a dynamic image of a lung field of COPD (obstructive disease).
- FIG. 14C is a display screen displaying an analysis result obtained by analyzing a dynamic image of a lung field of a mixed disease.
- step S57 “inhalation feature quantity / expiration feature quantity of each small region (2 mm square size) in the lung field region extracted from the half-cut frame image. ”Value histogram 64a, still image 64b displaying a list of“ inspiratory feature value / expiratory feature value ”of each small region, hue displayed in histogram 64a and still image 64b, and“ inspiratory feature value ”.
- a display 64c indicating the relationship with the value of “/ expiratory feature value” and an index value 64d indicating a tendency of “inspiratory feature value / expiratory feature value” in the entire lung field are displayed. Also, as shown in FIGS.
- the histogram 64a is displayed with the horizontal axis color-coded in six hues according to the value of “inhalation feature quantity / expiration feature quantity”.
- the still image 64b indicating “inspiration feature amount / expiration feature amount” of each small region each small region has the same standard as the color coding of the histogram according to the value of “inspiration feature amount / expiration feature amount”. Therefore, it is possible for a doctor to easily grasp a local abnormal part (occlusion part, restriction part) in the lung field.
- the average value and standard deviation are calculated and displayed together on the screen as an index value 64d indicating the tendency of “inhalation feature quantity / expiration feature quantity” in the entire lung field, thereby displaying the entire lung field.
- the tendency of “inhalation feature quantity / expiration feature quantity” in can be provided as a numerical value to a doctor.
- the maximum value (absolute value) of the inter-frame difference value in the expiration period is the exhalation feature amount
- the maximum value (absolute value) of the inter-frame difference value in the inspiration phase is the inspiration feature amount
- the lungs of a normal person when analysis was performed using dynamic images taken with a grid, the average value of “inspiratory feature value / expiratory feature value” of the entire lung field was 0.9 to 1.2, and the standard deviation was 0. It is known to be about 10 to 0.22. Therefore, when the display screen shown in FIG. 14A is displayed in step S57, the doctor can easily grasp that the photographed lung field is normal.
- the average value of “inhalation feature value / expiration feature value” of the entire lung field is 0.9. It is known that the standard deviation is not within ⁇ 1.2 (larger than that of a normal person) and the standard deviation is not within 0.10 to 0.22 (larger than that of a normal person). Therefore, when the display screen shown in FIG. 14B is displayed in step S57, the doctor can easily grasp that the captured lung field is COPD.
- the data of “inspiratory feature value / expiratory feature value” for the entire lung field is 0.66 or less. It is known that both the number and the number of data of 1.5 or more increase. Therefore, when the display screen shown in FIG. 14C is displayed in step S57, the doctor can easily grasp that the photographed lung field is a mixed disease.
- an uneven ventilation state such as COPD (obstructive pulmonary disease), mixed lung disease or the like, or It is possible to provide useful diagnostic information to the doctor so that the severity can be specified.
- the normality / abnormality of the subject M may be determined using an index value or the like indicating a tendency of “inhalation feature quantity / expiration feature quantity”. Even in this case, it is preferable to change the threshold value for determining normality / abnormality based on the grid presence / absence information. For example, when determining normal / abnormal based on the average value of “inhalation feature quantity / expiration feature quantity” for the entire lung field, the average value of 0.9 to 1.2 is assumed to be normal when there is a grid. In the absence, it is preferable to judge that an average value of 0.8 to 1.3 is normal.
- the feature quantity of exhalation and the feature quantity of inspiration may be used.
- the feature quantity of exhalation may be the number of frame images corresponding to the expiration period in one breath cycle (expiration time)
- the feature quantity of inspiration may be the number of frame images corresponding to the inspiration period in one breath cycle (inspiration period).
- the inspiratory time and the expiratory time should be approximately the same length or slightly longer. Therefore, by looking at the value of “the number of frame images corresponding to the expiration period / the number of frame images corresponding to the inspiration period”, the doctor can grasp whether or not there is a suspicion of lung disease.
- the region where “the number of frame images during the expiration period / the number of frame images during the inspiration period”> 1.5 may be an obstructive portion that is difficult to ventilate during expiration and delays the discharge of the captured air. Recognize.
- the relationship of “maximum value of inter-frame difference value in inspiratory period / maximum value of inter-frame difference value in expiratory period ⁇ expiratory time (number of expiratory frame images) / inspiratory time (number of inspiratory frame images)” is satisfied.
- Pulmonary diseases and mixed lung diseases can be identified.
- the signal value (average signal value) of the pixels in each small region is calculated, and the minimum value and the maximum value of the signal value in one respiratory cycle are determined for each small region.
- the minimum value may be the feature value of exhalation in the region
- the maximum value may be the feature value of inspiration. It can be considered that the difference between the maximum value and the minimum value of the signal value is large at a normal location, and the difference between the two is extremely small at a location where there is an abnormality. This is because the alveoli move poorly at abnormal locations, and the density change of the alveoli becomes small.
- the doctor refers to the histogram of “maximum value of signal value / minimum value of signal value” and confirms the average value and standard deviation to determine whether the lung field is normal or disease. can do. For example, when the average value of “maximum value of signal value / minimum value of signal value” of the entire lung field is greater than 1 and the standard deviation is small, it can be determined that the lung function is normal. On the other hand, if the average value of “maximum value of signal value / minimum value of signal value” of the entire lung field is a value close to 1 and the standard deviation is large, it can be determined that there is a disease in the lung function.
- the count number (number of blocks (small areas) is the peak of the histogram.
- the grid 46 is used when capturing a still image, and imaging is performed without using the grid 46 when capturing dynamic images.
- the dose that reaches the radiation detection unit 43 increases with respect to the irradiation dose to the same subject as the dose is not attenuated by the grid, so that high-definition video analysis can be performed.
- the arrival dose to the radiation detection unit 43 is constant, it is possible to reduce the irradiation dose to the subject (subject exposure dose) compared to the conventional case where a grid is used during dynamic imaging. .
- Each frame image obtained by dynamic imaging is accompanied by grid presence / absence information, and the control unit 61 of the diagnostic console 6 changes the analysis algorithm based on the grid presence / absence information. For example, a threshold value used for analysis is set. Since the change is made, it is possible to provide an accurate analysis result regardless of the presence or absence of the grid during dynamic imaging.
- imaging without a grid it is preferable to perform imaging by separating the subject M from the radiation detection unit 43 in order to reduce the influence of scattered radiation and perform better dynamic observation and dynamic analysis.
- the influence of scattered radiation can be reduced by the air gap. Enlarged shooting is performed by the separation, but there is no effect on the dynamic analysis result.
- imaging conditions for phase contrast imaging may be used.
- the subject detection sensor 454 may detect whether the subject is separated or in close contact, and control may be performed to prohibit dynamic imaging when the subject is in close contact without a grid.
- the maximum flow rate ratio is displayed as a histogram as an analysis result of “inhalation feature quantity / expiration feature quantity”, but the present invention is not limited to this.
- FIG. 15 shows the entire lung field calculated from dynamic images, with the X axis being the average value of “the maximum value of the interframe difference value in the inspiratory period / the maximum value of the interframe difference value in the expiration period” and the Y axis being the standard deviation.
- FIG. 6 is a graph in which an average value of “inhalation feature value / expiration feature value” and an index value of standard deviation are plotted.
- a value obtained by linearly combining two index values (for example, an average value and a standard deviation of “inspiratory feature value / exhalation feature value”) is used as a new index value, and the two indices as shown in FIG.
- a graph with one of the values on the X axis and the other on the Y axis shows thresholds th1 to th4 for classifying new index values (index values obtained by linearly combining two index values) according to severity.
- a new index value calculated from the dynamic image may be plotted.
- the first principal component calculated by principal component analysis from a large number of measured value data with respect to the average value and the standard deviation can be used as the index value.
- the first principal component is calculated by calculating the maximum eigenvalue of the covariance matrix for N ⁇ M data from a large number (N) of measured value data for a plurality (M) of index values,
- the calculated principal component may be used as an index value.
- an index value indicating a trend may be calculated for each of the inspiration feature value or the expiration feature value.
- inspiratory feature value / expiratory feature value is converted into a logarithm, and a histogram of the above-described histogram is created, an index value is calculated, and a colorized image based on a conversion table with a display parameter is converted to the logarithmically converted value. You may make it.
- an imaging device capable of both dynamic imaging and still image imaging has been described.
- an imaging device capable of only dynamic imaging can also perform imaging without using a grid. It is effective to reduce the exposure amount of the subject or improve the image quality.
- a hard disk, a semiconductor nonvolatile memory, or the like is used as a computer-readable medium of the program according to the present invention, but the present invention is not limited to this example.
- a portable recording medium such as a CD-ROM can be applied.
- a carrier wave is also applied as a medium for providing program data according to the present invention via a communication line.
- each device constituting the chest diagnosis support system 100 and the radiation imaging system 200 can be changed as appropriate without departing from the spirit of the present invention.
- It may be used to support diagnosis of chest dynamics in the medical field.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Physiology (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Pulmonology (AREA)
- Cardiology (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Human Computer Interaction (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
胸部診断支援システムは、
少なくとも1サイクルの呼吸周期を含む胸部の動態を撮影し、連続する複数のフレーム画像を生成する撮影手段と、
前記撮影手段により生成された複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段と、
前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段と、
前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段と、
前記解析手段により作成されたヒストグラムを表示する表示手段と、
を備える。
プログラムは、
コンピュータを、
少なくとも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線等の放射線を連続照射して複数の画像を取得(即ち、連続撮影)することにより行う。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
撮影装置1は、図1に示すように、放射線源11、放射線照射制御装置12、放射線検出部13、読取制御装置14等を備えて構成されている。
放射線照射制御装置12は、撮影用コンソール2に接続されており、撮影用コンソール2から入力された放射線照射条件に基づいて放射線源11を制御して放射線撮影を行う。撮影用コンソール2から入力される放射線照射条件は、例えば、連続照射時のパルスレート、パルス幅、パルス間隔、1撮影あたりの撮影フレーム数、X線管電流の値、X線管電圧の値、フィルタ種等である。パルスレートは、1秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射1回当たりの放射線照射時間である。パルス間隔は、連続撮影において、1回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。
撮影用コンソール2は、放射線照射条件や画像読取条件を撮影装置1に出力して撮影装置1による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置1により取得された動態画像を撮影技師によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
撮影用コンソール2は、図1に示すように、制御部21、記憶部22、操作部23、表示部24、通信部25を備えて構成され、各部はバス26により接続されている。
診断用コンソール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に示す画像解析処理が実行される。
まず、各フレーム画像から肺野領域が抽出される(ステップS21)。
肺野領域の抽出方法は何れの方法であってもよい。例えば、一連のフレーム画像中の任意のフレーム画像(ここでは、撮影順が一番(最初)のフレーム画像とする。)の各画素の信号値(濃度値)のヒストグラムから判別分析によって閾値を求め、この閾値より高信号の領域を肺野領域候補として1次抽出する。次いで、1次抽出された肺野領域候補の境界付近でエッジ検出を行い、境界付近の小領域でエッジが最大となる点を境界に沿って抽出すれば肺野領域の境界を抽出することができる。
しかし、安静呼吸時に撮影された画像においては、上記の位置ずれは小さく、後述する解析結果が狂ってしまうほどの位置ずれはおきない。図5の画像D1は、安静呼気位(安静呼吸時において横隔膜の位置が最も高くなったタイミング)のフレーム画像である。図5の画像D2は、安静吸気位(安静呼吸時において横隔膜の位置が最も低くなったタイミング)のフレーム画像である。即ち、図5の画像D1とD2とは、呼吸1サイクルにおいて最も形状の差の大きいタイミングで撮影された画像である。しかし、図5の画像D1、D2間においては、最も位置ずれの大きい肺野領域の下部領域においても位置ずれはわずかであることがわかる(画像D2のA11は画像D1のA1と同じ画素位置を示し、画像D2のA2は画像D1のA1と肺野における同一部分を描画した領域を示している)。
基準画像としては、安静呼気位のフレーム画像とすることが好ましい。安静呼気位では、安静呼吸時において横隔膜の位置が最も高くなる、即ち、肺野領域の面積が最も小さくなるので、基準画像の小領域を他のフレーム画像に対応付けたときに、小領域が他のフレーム画像の肺野外の領域に対応付けられることがないからである。
安静呼気位の画像は、一連のフレーム画像の中から横隔膜の位置が最も高い位置にある画像を抽出することで取得することができる。横隔膜の位置は、例えば、図6に示す横隔膜の基準位置B1を横隔膜のカーブC(図6に点線で示す)の垂直方向の平均位置と予め定義しておき、肺野領域Rから横隔膜のカーブC(肺野領域の下端)を抽出し、その垂直方向の平均位置を求め、求めた位置を横隔膜の基準位置B1として特定する。
次いで、撮影順が1番のフレーム画像をP1、これと隣接するフレーム画像(撮影順が隣接するフレーム画像(即ち、時間的に隣接するフレーム画像。以下同様。))をP2とし、P2に、P1の各小領域の探索領域を設定する。ここで、P2の探索領域は、P1における各小領域における中心点の座標を(x,y)とすると、同一の中心点(x,y)をもち、P1の小領域よりも縦横の幅が大きくなるように設定する(例えば、1.5倍)。そして、P1の各領域毎に、P2の探索範囲で最もマッチング度合いが高くなる領域を求めることで、P1の各小領域に対するP2上での対応位置を算出する。マッチング度合いとしては、最小二乗法や相互相関係数を指標に用いる。そして、P2の肺野領域をP1の各小領域の対応位置で分割する。
次いで、P2を、新たにP1とみなし、撮影順がP2の次のフレーム画像を新たなP2とみなして、P1の各小領域におけるP2の対応位置を算出する。以上の処理を繰り返すことで、各フレーム画像の各小領域が隣接するフレーム画像のどの位置に対応するかが求まる。求めた処理結果は、制御部31のRAMに記憶される。
このとき、例えば、数個(例えば、5~6個)の色相を上記各カテゴリの閾値と対応させ、さらにその間の特徴量の値に対しては中間の色相を割り当てる(グラデーションさせる)ことで識別能の高い表示とすることができる。
表示用パラメータ値に基づいて色付けされた画像は、基準画像のフレーム画像にオーバーレイ表示してもよい。
図8Aは、正常者の肺野の動態画像を解析した解析結果を表示した表示画面である。図8Bは、COPD(閉塞性疾患)の肺野の動態画像を解析した解析結果を表示した表示画面である。図8Cは、混合性疾患の肺野の動態画像を解析した解析結果を表示した表示画面である。
一方、COPD(閉塞性疾患)の肺野では、肺野全体の「吸気の特徴量/呼気の特徴量」の平均値は0.9~1.2に収まらず(正常者に比べて大きく)、標準偏差も、0.10~0.22に収まらない(正常者に比べて大きい)ことが知られている。よって、ステップS27で図8Bに示す表示画面が表示された場合、医師は、撮影した肺野がCOPDであることを容易に把握することができる。
一方、混合性肺疾患の肺野では、肺野全体の「吸気の特徴量/呼気の特徴量」の値が0.66以下のデータ数、及び1.5以上のデータ数がともに増加することが知られている。よって、ステップS27で図8Cに示す表示画面が表示された場合、医師は、撮影した肺野が混合性疾患であることを容易に把握することができる。
このように、胸部診断支援システム100においては、COPD(閉塞性肺疾患)、混合性肺疾患等の病態を特定できるような、有用な診断情報を医師に提供することが可能となる。
例えば、呼気の特徴量を呼吸1サイクルにおける呼気期に相当するフレーム画像数(呼気時間)、吸気の特徴量を呼吸1サイクルにおける吸気期に相当するフレーム画像数(吸気期間)としてもよい。ここで、肺の換気機能が正常である場合、吸気時間と呼気時間はほぼ同じ長さ、もしくは若干呼気時間が長くなる筈である。よって、「呼気期に相当するフレーム画像数/吸気期に相当するフレーム画像数」の値を見れば、医師は肺疾患の疑いがあるか否かを把握することができる。特に、「呼気期のフレーム画像数/吸気期のフレーム画像数」>1.5の領域は、呼気で換気がしにくく、取り込んだ空気を排出するのが遅れてしまう閉塞性部分であることがわかる。なお、「吸気期のフレーム間差分値の最大値/呼気期のフレーム間差分値の最大値≒呼気時間(呼気のフレーム画像数)/吸気時間(吸気のフレーム画像数)」の関係が成り立つため、医師は、呼気の特徴量を呼気期のフレーム間差分値の最大値、吸気の特徴量を吸気期のフレーム間差分値の最大値とした場合と同様の判断基準により、正常、COPD(閉塞性肺疾患)、混合性肺疾患の識別を行うことが可能である。
以下、図面を参照して本発明の第2の実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
また、近年、特開2005-312775号公報(公知文献2)に記載のように、FPD(flat panel detector)等の放射線検出器を利用して胸部の動態撮影を行い、診断に応用するシステムが提案されているが、このような動態撮影用のシステムにおいても、静止画撮影と同様に、グリッドを使用して撮影を行う構成となっている。
従って、動態撮影においては、グリッド無しで撮影を行うことができる。グリッド無しでの撮影において、検出器に到達する線量を従来並みに保つ場合、被写体への照射線量はグリッドによる減衰分を予め加算しておく必要が無く、従来手法に比べ、被写体の被曝線量を半減することが可能となる。また、被写体の被曝線量(照射線量)を一定のままとすると、各フレーム画像あたりの検出器到達線量を増加することができるので、より精度を有する解析を行うことができる。
まず、構成を説明する。
図9に、本実施の形態における放射線撮影システム200の全体構成を示す。この放射線撮影システム200は、胸部診断支援システムとして機能するシステムである。
図9に示すように、放射線撮影システム200は、撮影装置4と、撮影用コンソール5とが通信ケーブル等により接続され、撮影用コンソール5と、診断用コンソール6とがLAN(Local Area Network)等の通信ネットワークNT2を介して接続されて構成されている。放射線撮影システム200を構成する各装置は、DICOM(Digital Image and Communications in Medicine)規格に準じており、各装置間の通信は、DICOMに則って行われる。
撮影装置4は、被写体M(人体の胸部)に放射線を照射して被写体Mの動態撮影又は静止画撮影を行う装置である。
動態撮影とは、被写体Mに対し、X線等の放射線をパルス的に連続照射して複数の画像を取得(即ち、連続撮影)することをいう。動態撮影では、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性(サイクル)を持つ被写体Mの動態を撮影する。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
静止画撮影とは、従来のフィルム方式やCR方式と同様に撮影部位の濃度分解能に基づく診断に使用されるもので、被写体Mに対し、X線等の放射線を1回照射して一枚の静止画像を取得することをいう。
放射線照射制御装置42は、撮影用コンソール5に接続されており、撮影用コンソール5から入力された放射線照射条件に基づいて放射線源41を制御して放射線撮影を行う。撮影用コンソール5から入力される放射線照射条件は、例えば、連続照射時のパルスレート、パルス幅、パルス間隔、1撮影あたりの撮影フレーム数、X線管電流の値、X線管電圧の値、フィルタ種等である。パルスレートは、1秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射1回当たりの放射線照射時間である。パルス間隔は、連続撮影において、1回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。
放射線検出部43は、図10に示すように、被写体Mを挟んで放射線源41と対向するように保持部45によって保持されている。
撮影用コンソール5は、放射線照射条件や画像読取条件を撮影装置4に出力して撮影装置4による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置4により取得された静止画像又は動態画像に基づく画像、例えば間引きビニング処理されたプレビュ画像や諧調処理等を施した処理済画像等を適宜作成し、撮影実施者によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
撮影用コンソール5は、図9に示すように、制御部51、記憶部52、操作部53、表示部54、通信部55を備えて構成され、各部はバス56により接続されている。
また、記憶部52は、動態撮影用、静止画撮影用のそれぞれの放射線照射条件及び画像読取条件を記憶している。
診断用コンソール6は、撮影用コンソール5から静止画像、又は動態画像の一連のフレーム画像を取得し、取得した画像や後述するヒストグラム等の診断支援情報を表示して医師が読影診断するためのコンピュータ装置である。
診断用コンソール6は、図9に示すように、制御部61、記憶部62、操作部63、表示部64、通信部65を備えて構成され、各部はバス66により接続されている。
次に、上記放射線撮影システム200における動作について説明する。
まず、撮影装置4、撮影用コンソール5による撮影動作について説明する。
図12A~図12Bに、撮影用コンソール5の制御部51において実行される撮影制御処理を示す。撮影制御処理は、制御部51と記憶部52に記憶されている撮影制御処理プログラムとの協働により実行される。
制御部51は、このステップS38によって、グリッド46が装着されずに静止画撮影が行われることのないように制御している。
次に、診断用コンソール6における動作について説明する。
診断用コンソール6においては、通信部65を介して撮影用コンソール5から静止画像が受信され、操作部63により当該画像の表示指示が入力されると、表示部64に受信された静止画像が表示され、医師の診断に供される。
一方、通信部65を介して撮影用コンソール5から動態画像の一連のフレーム画像が受信され、操作部63により動態解析が指示されると、制御部61と記憶部62に記憶されている画像解析処理プログラムとの協働により画像解析処理が実行される。
第2の実施の形態における画像解析処理は、第1の実施の形態で図3で示した画像解析処理と大まかな処理ステップは同じであるが詳細が異なるので、以下に図3及び図13以降の図を参照して説明する。
次いで、各フレーム画像間で対応づけられた各小領域に対して、時間軸方向のフィルタ処理が施される。このフィルタ処理は、血流等の高周波数な信号変化を除去し、換気による信号値の時間変化を抽出するための処理であり、例えば、小領域毎の信号値の時間変化に対して、安静呼吸画像群ではカットオフ周波数0.7Hz、深呼吸画像群ではカットオフ周波数0.5Hzでローパスフィルタリングを行う。ここで、ローパスフィルタのカットオフ周波数は、固定値とするよりも、撮影された動態画像毎に最適化することが好ましい。例えば、前述のように一連のフレーム画像の横隔膜の位置を解析し、安静換気の場合は安静呼気位および安静吸気位となるフレームを検出し、安静呼気位のフレームとその次の安静吸気位のフレームの間のフレーム数から吸気期の時間を求め、その逆数に所定の係数を乗じた値をカットオフ周波数としたローパスフィルタを実施する。このとき安静換気の場合、自動的に設定されるカットオフ周波数は、0.2~1.0Hzの間に制限することが好ましい。また、ステップS31において、別途測定された安静時における1分間の呼吸数、脈拍数等のバイタルサインを患者情報として入力しておき、これらの値からカットオフ周波数が算出されるようにしてもよい。例えば、患者情報として入力された1分間の呼吸数を1秒間の呼吸数に変換し、その呼吸数に所定の係数を乗じた値をカットオフ周波数としてローパスフィルタを実施してもよい。さらには、入力された1分間の脈拍数を1秒間の脈拍数に変換し、1秒間の呼吸数と1秒間の心拍数の平均値をカットオフ周波数としてローパスフィルタを実施するようにしてもよい。
フィルタ処理に次いで、撮影順が隣接するフレーム画像間で各小領域の信号値の差分を算出するフレーム間差分処理が行われる。ここでは、各小領域毎に、フレーム番号NとN+1(Nは1、2、3・・・)のフレーム画像について、N+1-Nの差分値が算出される。そして、呼気期のフレーム間差分値の最大値(絶対値の最大値)が呼気の特徴量として取得され、吸気期のフレーム間差分値の最大値(絶対値の最大値)が吸気の特徴量として取得される。フレーム間差分値の最大値(絶対値の最大値)は、最大微分値に相当する。ここで、各小領域内のフレーム間差分値の符号が正の期間が吸気期、負の期間が呼気期である。
このとき、例えば、数個(例えば、5~6個)の色相を上記各カテゴリの閾値と対応させ、その間の特徴量の値に対しては中間の色相を割り当てる(グラデーションさせる)ことで識別能の高い表示とすることができる。
表示用パラメータ値に基づいて色付けされた画像は、基準画像のフレーム画像にオーバーレイ表示してもよい。
図13は、或る正常の肺野をグリッド有りとグリッド無しで撮影した動態画像の動態解析結果の比較を示す図である。図13においては、グリッド有りとグリッド無しで撮影した動態画像の各小領域を「吸気期のフレーム間差分値の最大値」(最大吸気気流速度)に基づいて色付けした画像、「呼気期のフレーム間差分値の最大値」(最大呼気気流速度)に基づいて色付けした画像、最大流速比に基づいて色付けした画像、及び最大流速比のヒストグラム、を解析結果として示している。
検出器サイズ40×30cm、検出器画素サイズ194μm、グリッドピッチ80line/cm、グリッド比12:1、管球~検出器間距離2m、フレーム数75フレーム(約10秒間撮影)、トータル被曝線量(検出器への到達線量を一定とした場合の被写体被曝線量)グリッド有り0.24mGy、グリッド無し0.14mGy
また、最大吸気気流速度、最大呼気気流速度、最大流速比の大きさと色(図13においては濃度で示す)との変換テーブルは、両者を比較するために同一(ここでは、グリッド有り用の変換テーブル)を用いている。
また、同一患者の過去の画像の解析結果と現在の解析結果とを比較するような場合も想定される。互いにグリッド有無情報が付帯されている場合、一方の解析結果を他方の解析結果に合わせるべく、何れか一方の解析結果の値を補正して使用すること(例えば、両者の変換テーブルの閾値の差分を解析結果としての特徴量に加算又は減算する等)が好ましい。
図14Aは、グリッド有りで撮影された正常者の肺野の動態画像を解析した解析結果を表示した表示画面である。図14Bは、COPD(閉塞性疾患)の肺野の動態画像を解析した解析結果を表示した表示画面である。図14Cは、混合性疾患の肺野の動態画像を解析した解析結果を表示した表示画面である。
尚、本表示例に於いては、縦軸は、2mm角の小ブロックの出現個数/全ブロック数((430mm÷2mm)×(354mm÷2mm)=38055個)の割合としているが、図8に示すように、縦軸を2mm角の小ブロックの出現個数自体としても良い。動態撮影に使用する検出器のサイズが複数存在する場合には、本表示例のような全体の領域個数で正規化した割合表示によって同一症例に対する解析結果表示の統一性が達成され、解析結果に対する医師等の誤認を招き難いので好ましい。
一方、COPD(閉塞性疾患)の肺野では、グリッド有りで動態撮影した動態画像で解析を行った場合、肺野全体の「吸気の特徴量/呼気の特徴量」の平均値は0.9~1.2に収まらず(正常者に比べて大きく)、標準偏差も、0.10~0.22に収まらない(正常者に比べて大きい)ことが知られている。よって、ステップS57で図14Bに示す表示画面が表示された場合、医師は、撮影した肺野がCOPDであることを容易に把握することができる。
一方、混合性肺疾患の肺野では、グリッド有りで動態撮影した動態画像で解析を行った場合、肺野全体の「吸気の特徴量/呼気の特徴量」の値が0.66以下のデータ数、及び1.5以上のデータ数がともに増加することが知られている。よって、ステップS57で図14Cに示す表示画面が表示された場合、医師は、撮影した肺野が混合性疾患であることを容易に把握することができる。
このように、放射線撮影システム200においては、「吸気の特徴量/呼気の特徴量」の傾向を示す指標値から、COPD(閉塞性肺疾患)、混合性肺疾患等の換気不均等の病態もしくはその重症度を特定できるような、有用な診断情報を医師に提供することが可能となる。
例えば、呼気の特徴量を呼吸1サイクルにおける呼気期に相当するフレーム画像数(呼気時間)、吸気の特徴量を呼吸1サイクルにおける吸気期に相当するフレーム画像数(吸気期間)としてもよい。ここで、肺の換気機能が正常である場合、吸気時間と呼気時間はほぼ同じ長さ、もしくは若干呼気時間が長くなる筈である。よって、「呼気期に相当するフレーム画像数/吸気期に相当するフレーム画像数」の値を見れば、医師は肺疾患の疑いがあるか否かを把握することができる。特に、「呼気期のフレーム画像数/吸気期のフレーム画像数」>1.5の領域は、呼気で換気がしにくく、取り込んだ空気を排出するのが遅れてしまう閉塞性部分であることがわかる。なお、「吸気期のフレーム間差分値の最大値/呼気期のフレーム間差分値の最大値≒呼気時間(呼気のフレーム画像数)/吸気時間(吸気のフレーム画像数)」の関係が成り立つため、医師は、呼気の特徴量を呼気期のフレーム間差分値の最大値、吸気の特徴量を吸気期のフレーム間差分値の最大値とした場合と同様の判断基準により、正常、COPD(閉塞性肺疾患)、混合性肺疾患の識別を行うことが可能である。
200 放射線撮影システム
1、4 撮影装置
11、41 放射線源
12、42 放射線照射制御装置
13、43 放射線検出部
14、44 読取制御装置
2、5 撮影用コンソール
21、51 制御部
22、52 記憶部
23、53 操作部
24、54 表示部
25、55 通信部
26、56 バス
3、6 診断用コンソール
31、61 制御部
32、62 記憶部
33、63 操作部
34、64 表示部
35、65 通信部
36、66 バス
Claims (8)
- 少なくとも1サイクルの呼吸周期を含む胸部の動態を撮影し、連続する複数のフレーム画像を生成する撮影手段と、
前記撮影手段により生成された複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段と、
前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段と、
前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段と、
前記解析手段により作成されたヒストグラムを表示する表示手段と、
を備える胸部診断支援システム。 - 前記吸気の特徴量が、吸気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値であり、前記呼気の特徴量が、呼気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値である請求項1に記載の胸部診断支援システム。
- 前記吸気の特徴量が、呼吸1サイクルにおける吸気期に相当するフレーム画像の数であり、前記呼気の特徴量が、前記呼吸1サイクルにおける呼気期に相当するフレーム画像の数である請求項1に記載の胸部診断支援システム。
- 前記吸気の特徴量が、呼吸1サイクル分の複数のフレーム画像における画素信号値の最大値であり、前記呼気の特徴量が、前記呼吸1サイクル分の複数のフレーム画像における画素信号値の最小値である請求項1に記載の胸部診断支援システム。
- 前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像の肺野領域を前記一のフレーム画像の各小領域と同じ画素位置の小領域に分割し、前記複数のフレーム画像間の前記同じ画素位置の各小領域を互いに対応付ける請求項1~4の何れか一項に記載の胸部診断支援システム。
- 前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像のそれぞれにおいて前記一のフレーム画像における各小領域とマッチング度合いの高い領域を抽出し、前記他のフレーム画像の肺野領域を前記抽出された各領域に分割し、前記複数のフレーム画像間の前記マッチング度合いの高い各小領域を互いに対応付ける請求項1~4の何れか一項に記載の胸部診断支援システム。
- 前記撮影手段は、グリッドを使用せずに動態撮影を行う請求項1~6の何れか一項に記載の胸部診断支援システム。
- コンピュータを、
少なくとも1サイクルの呼吸周期を含む胸部の動態を示す複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段、
前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段、
前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段、
前記解析手段により作成されたヒストグラムを表示する表示手段、
として機能させるためのプログラム。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012530550A JP5765343B2 (ja) | 2010-08-27 | 2011-03-09 | 胸部診断支援システム及びプログラム |
CN201180041386.5A CN103079466B (zh) | 2010-08-27 | 2011-03-09 | 胸部诊断辅助*** |
US13/819,195 US9044194B2 (en) | 2010-08-27 | 2011-03-09 | Thoracic diagnosis assistance system and computer readable storage medium |
US14/556,410 US9237877B2 (en) | 2010-08-27 | 2014-12-01 | Thoracic diagnosis assistance system and computer readable storage medium |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010190241 | 2010-08-27 | ||
JP2010-190241 | 2010-08-27 | ||
JP2011-013538 | 2011-01-26 | ||
JP2011013538 | 2011-01-26 |
Related Child Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/819,195 A-371-Of-International US9044194B2 (en) | 2010-08-27 | 2011-03-09 | Thoracic diagnosis assistance system and computer readable storage medium |
US14/556,410 Continuation US9237877B2 (en) | 2010-08-27 | 2014-12-01 | Thoracic diagnosis assistance system and computer readable storage medium |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012026146A1 true WO2012026146A1 (ja) | 2012-03-01 |
Family
ID=45723166
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/055524 WO2012026146A1 (ja) | 2010-08-27 | 2011-03-09 | 胸部診断支援システム及びプログラム |
Country Status (4)
Country | Link |
---|---|
US (2) | US9044194B2 (ja) |
JP (1) | JP5765343B2 (ja) |
CN (1) | CN103079466B (ja) |
WO (1) | WO2012026146A1 (ja) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013140749A1 (ja) * | 2012-03-23 | 2013-09-26 | 富士フイルム株式会社 | 画像解析装置および方法並びにプログラム |
US20130331725A1 (en) * | 2012-06-11 | 2013-12-12 | Konica Minolta, Inc. | Thoracic diagnosis assistance information generation method, thoracic diagnosis assistance system, and dynamic state image processing apparatus |
JP2014000351A (ja) * | 2012-06-21 | 2014-01-09 | Canon Inc | 画像診断支援装置 |
WO2014162911A1 (ja) * | 2013-04-05 | 2014-10-09 | 東芝メディカルシステムズ株式会社 | 医用画像処理装置及び医用画像処理方法 |
CN104470427A (zh) * | 2012-06-08 | 2015-03-25 | 皇家飞利浦有限公司 | 用于监测患者的肺功能的方法和*** |
JP2017012444A (ja) * | 2015-06-30 | 2017-01-19 | キヤノン株式会社 | 画像処理装置および画像処理方法、画像処理プログラム |
JP2018007949A (ja) * | 2016-07-15 | 2018-01-18 | コニカミノルタ株式会社 | 動態解析装置 |
JP2018148965A (ja) * | 2017-03-10 | 2018-09-27 | コニカミノルタ株式会社 | 動態解析システム |
WO2019202695A1 (ja) * | 2018-04-19 | 2019-10-24 | 株式会社島津製作所 | X線撮影システム |
JP2020501818A (ja) * | 2016-12-23 | 2020-01-23 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 被験者の肺状態の評価を支援するシステム |
JP2020048685A (ja) * | 2018-09-25 | 2020-04-02 | キヤノンメディカルシステムズ株式会社 | 乳がん診断支援装置、乳がん診断支援システム、および乳がん診断支援方法 |
JP2020081279A (ja) * | 2018-11-22 | 2020-06-04 | コニカミノルタ株式会社 | 計測装置及びプログラム |
Families Citing this family (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5408399B1 (ja) * | 2012-03-23 | 2014-02-05 | コニカミノルタ株式会社 | 画像生成装置 |
JP6191160B2 (ja) * | 2012-07-12 | 2017-09-06 | ノーリツプレシジョン株式会社 | 画像処理プログラムおよび画像処理装置 |
JP6283875B2 (ja) * | 2013-09-05 | 2018-02-28 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置、x線診断装置およびx線コンピュータ断層撮影装置 |
JP5968928B2 (ja) * | 2014-02-14 | 2016-08-10 | 富士フイルム株式会社 | 医用画像表示制御装置およびその作動方法並びに医用画像表示制御プログラム |
JP6318739B2 (ja) * | 2014-03-17 | 2018-05-09 | コニカミノルタ株式会社 | 画像処理装置、およびプログラム |
ES2841414T3 (es) | 2014-04-15 | 2021-07-08 | 4Dmedical Ltd | Método de obtención de imágenes |
JP6072102B2 (ja) * | 2015-01-30 | 2017-02-01 | キヤノン株式会社 | 放射線撮影システム及び放射線撮影方法 |
JP6540807B2 (ja) * | 2015-07-24 | 2019-07-10 | コニカミノルタ株式会社 | 撮影用コンソール |
US11723617B2 (en) | 2016-02-03 | 2023-08-15 | 4DMedical Limited | Method and system for imaging |
US11278256B2 (en) | 2016-03-04 | 2022-03-22 | 4DMedical Limited | Method and system for imaging |
JP6740699B2 (ja) * | 2016-05-10 | 2020-08-19 | コニカミノルタ株式会社 | 画像解析システム |
JP6686733B2 (ja) * | 2016-06-23 | 2020-04-22 | コニカミノルタ株式会社 | 動態解析システム |
JP2018000281A (ja) * | 2016-06-28 | 2018-01-11 | コニカミノルタ株式会社 | 動態解析システム |
US10417745B2 (en) * | 2016-06-28 | 2019-09-17 | Raytheon Company | Continuous motion scene based non-uniformity correction |
EP3936048A1 (en) * | 2016-07-19 | 2022-01-12 | Paramevia Pte. Ltd. | Diagnosis support program |
JP2018122093A (ja) * | 2017-02-02 | 2018-08-09 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置、x線ct装置、及び医用画像処理方法 |
JP6878038B2 (ja) * | 2017-02-15 | 2021-05-26 | キヤノン株式会社 | 放射線撮影装置、放射線撮影システム、画像処理装置、画像処理装置の作動方法、及びプログラム |
CN111837196A (zh) * | 2017-08-16 | 2020-10-27 | 皇家飞利浦有限公司 | 动态空气支气管图的超声确定及相关联的设备、***和方法 |
JP7071095B2 (ja) * | 2017-11-15 | 2022-05-18 | キヤノン株式会社 | 放射線撮像装置 |
JP6950483B2 (ja) * | 2017-11-20 | 2021-10-13 | コニカミノルタ株式会社 | 動態撮影システム |
EP3494885A1 (en) | 2017-12-07 | 2019-06-12 | Koninklijke Philips N.V. | Apparatus for presentation of dark field x-ray image information |
JP6927020B2 (ja) * | 2017-12-22 | 2021-08-25 | コニカミノルタ株式会社 | 動態画像処理方法、動態画像処理装置及びプログラム |
JP7047574B2 (ja) * | 2018-04-26 | 2022-04-05 | コニカミノルタ株式会社 | 動態画像解析装置、動態画像解析システム、動態画像解析プログラム及び動態画像解析方法 |
CN111918611B (zh) * | 2018-05-16 | 2023-12-29 | 松下控股株式会社 | 胸部x线图像的异常显示控制方法、记录介质及装置 |
CN109102568A (zh) * | 2018-05-28 | 2018-12-28 | 杭州阿特瑞科技有限公司 | 基于区域分解的血管血流模拟方法及相关装置 |
WO2020044840A1 (ja) * | 2018-08-31 | 2020-03-05 | 富士フイルム株式会社 | 領域分割装置、方法およびプログラム、類似度決定装置、方法およびプログラム、並びに特徴量導出装置、方法およびプログラム |
EP3656308A1 (en) * | 2018-11-20 | 2020-05-27 | Koninklijke Philips N.V. | Apparatus for alveoli based visualization in dark-field or phase-contrast x-ray imaging |
CN109949233B (zh) * | 2019-02-18 | 2022-12-13 | 深圳蓝影医学科技股份有限公司 | 滤除x射线图像中散射线的方法、***、设备及存储介质 |
JP7344769B2 (ja) * | 2019-11-22 | 2023-09-14 | キヤノン株式会社 | 放射線検出装置及び出力方法 |
CN113576488B (zh) * | 2021-06-18 | 2023-06-23 | 深圳技术大学 | 基于心率的肺影像组学的确定方法及装置、设备及介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63168152A (ja) * | 1986-12-29 | 1988-07-12 | 富士電機株式会社 | 医用画像処理方法 |
JP2006020716A (ja) * | 2004-07-06 | 2006-01-26 | Toshiba Corp | X線循環器診断システム |
JP2006043200A (ja) * | 2004-08-05 | 2006-02-16 | Japan Science & Technology Agency | 脳内出血/くも膜下出血診断支援システム |
WO2006137294A1 (ja) * | 2005-06-21 | 2006-12-28 | National University Corporation Kanazawa University | X線診断支援装置、プログラム及び記録媒体 |
JP2009273671A (ja) * | 2008-05-15 | 2009-11-26 | Konica Minolta Medical & Graphic Inc | 動態撮影システム |
JP2010069099A (ja) * | 2008-09-19 | 2010-04-02 | Toshiba Corp | 画像処理装置及びx線コンピュータ断層撮影装置 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1098587A (ja) | 1996-09-20 | 1998-04-14 | Fuji Photo Film Co Ltd | 放射線画像読取装置 |
US6850597B2 (en) | 1998-09-11 | 2005-02-01 | Canon Kabushiki Kaisha | X-ray image photographing apparatus and grid device |
JP3477166B2 (ja) * | 2000-12-07 | 2003-12-10 | 学校法人慶應義塾 | 監視装置 |
JP3639825B2 (ja) * | 2002-04-03 | 2005-04-20 | キヤノン株式会社 | 動画像表示方法、プログラム、コンピュータ可読記憶媒体、及び動画像表示装置 |
JP3697233B2 (ja) * | 2002-04-03 | 2005-09-21 | キヤノン株式会社 | 放射線画像処理方法及び放射線画像処理装置 |
JP2004177251A (ja) | 2002-11-27 | 2004-06-24 | Canon Inc | 放射線画像撮影装置 |
JP4404291B2 (ja) | 2003-04-08 | 2010-01-27 | キヤノン株式会社 | 画像処理装置及び方法及びシステム |
JP4612802B2 (ja) | 2004-04-30 | 2011-01-12 | キヤノン株式会社 | 放射線画像取得装置および放射線画像取得方法、プログラム |
US7570732B2 (en) * | 2005-11-09 | 2009-08-04 | Dexela Limited | Methods and apparatus for obtaining low-dose imaging |
WO2009078297A1 (ja) * | 2007-12-19 | 2009-06-25 | Konica Minolta Medical & Graphic, Inc. | 動態画像処理システム |
JP2009153677A (ja) * | 2007-12-26 | 2009-07-16 | Konica Minolta Medical & Graphic Inc | 動態画像処理システム |
JP5136478B2 (ja) * | 2009-03-17 | 2013-02-06 | 株式会社島津製作所 | 放射線撮影装置 |
-
2011
- 2011-03-09 CN CN201180041386.5A patent/CN103079466B/zh active Active
- 2011-03-09 JP JP2012530550A patent/JP5765343B2/ja active Active
- 2011-03-09 US US13/819,195 patent/US9044194B2/en active Active
- 2011-03-09 WO PCT/JP2011/055524 patent/WO2012026146A1/ja active Application Filing
-
2014
- 2014-12-01 US US14/556,410 patent/US9237877B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63168152A (ja) * | 1986-12-29 | 1988-07-12 | 富士電機株式会社 | 医用画像処理方法 |
JP2006020716A (ja) * | 2004-07-06 | 2006-01-26 | Toshiba Corp | X線循環器診断システム |
JP2006043200A (ja) * | 2004-08-05 | 2006-02-16 | Japan Science & Technology Agency | 脳内出血/くも膜下出血診断支援システム |
WO2006137294A1 (ja) * | 2005-06-21 | 2006-12-28 | National University Corporation Kanazawa University | X線診断支援装置、プログラム及び記録媒体 |
JP2009273671A (ja) * | 2008-05-15 | 2009-11-26 | Konica Minolta Medical & Graphic Inc | 動態撮影システム |
JP2010069099A (ja) * | 2008-09-19 | 2010-04-02 | Toshiba Corp | 画像処理装置及びx線コンピュータ断層撮影装置 |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9117287B2 (en) | 2012-03-23 | 2015-08-25 | Fujifilm Corporation | Image analysis apparatus, method, and program |
JP2013192912A (ja) * | 2012-03-23 | 2013-09-30 | Fujifilm Corp | 画像解析装置および方法並びにプログラム |
WO2013140749A1 (ja) * | 2012-03-23 | 2013-09-26 | 富士フイルム株式会社 | 画像解析装置および方法並びにプログラム |
US10342457B2 (en) | 2012-06-08 | 2019-07-09 | Koninjklijke Philips N.V. | Method and system for monitoring the lung function of a patient |
CN104470427A (zh) * | 2012-06-08 | 2015-03-25 | 皇家飞利浦有限公司 | 用于监测患者的肺功能的方法和*** |
EP2674912A1 (en) | 2012-06-11 | 2013-12-18 | Konica Minolta, Inc. | Thoracic diagnosis assistance information generation method, thoracic diagnosis assistance system, and dynamic state image processing apparatus |
CN103479375A (zh) * | 2012-06-11 | 2014-01-01 | 柯尼卡美能达株式会社 | 胸部诊断辅助信息生成方法及***以及动态图像处理装置 |
US20130331725A1 (en) * | 2012-06-11 | 2013-12-12 | Konica Minolta, Inc. | Thoracic diagnosis assistance information generation method, thoracic diagnosis assistance system, and dynamic state image processing apparatus |
JP2013255540A (ja) * | 2012-06-11 | 2013-12-26 | Konica Minolta Inc | 胸部診断支援情報生成方法、胸部診断支援システム及び動態画像処理装置 |
JP2014000351A (ja) * | 2012-06-21 | 2014-01-09 | Canon Inc | 画像診断支援装置 |
US9730666B2 (en) | 2013-04-05 | 2017-08-15 | Toshiba Medical Systems Corporation | Medical image processing apparatus and medical image processing method |
JP2014210171A (ja) * | 2013-04-05 | 2014-11-13 | 東芝メディカルシステムズ株式会社 | 医用画像処理装置及び医用画像処理方法 |
WO2014162911A1 (ja) * | 2013-04-05 | 2014-10-09 | 東芝メディカルシステムズ株式会社 | 医用画像処理装置及び医用画像処理方法 |
JP2017012444A (ja) * | 2015-06-30 | 2017-01-19 | キヤノン株式会社 | 画像処理装置および画像処理方法、画像処理プログラム |
JP2018007949A (ja) * | 2016-07-15 | 2018-01-18 | コニカミノルタ株式会社 | 動態解析装置 |
JP2020501818A (ja) * | 2016-12-23 | 2020-01-23 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 被験者の肺状態の評価を支援するシステム |
JP2018148965A (ja) * | 2017-03-10 | 2018-09-27 | コニカミノルタ株式会社 | 動態解析システム |
WO2019202695A1 (ja) * | 2018-04-19 | 2019-10-24 | 株式会社島津製作所 | X線撮影システム |
JP2020048685A (ja) * | 2018-09-25 | 2020-04-02 | キヤノンメディカルシステムズ株式会社 | 乳がん診断支援装置、乳がん診断支援システム、および乳がん診断支援方法 |
JP7223539B2 (ja) | 2018-09-25 | 2023-02-16 | キヤノンメディカルシステムズ株式会社 | 乳がん診断支援装置、乳がん診断支援システム、および乳がん診断支援方法 |
JP2020081279A (ja) * | 2018-11-22 | 2020-06-04 | コニカミノルタ株式会社 | 計測装置及びプログラム |
JP7135775B2 (ja) | 2018-11-22 | 2022-09-13 | コニカミノルタ株式会社 | 計測装置及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
US9237877B2 (en) | 2016-01-19 |
US20130156158A1 (en) | 2013-06-20 |
CN103079466B (zh) | 2015-11-25 |
US9044194B2 (en) | 2015-06-02 |
JPWO2012026146A1 (ja) | 2013-10-28 |
JP5765343B2 (ja) | 2015-08-19 |
CN103079466A (zh) | 2013-05-01 |
US20150092911A1 (en) | 2015-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5765343B2 (ja) | 胸部診断支援システム及びプログラム | |
JP5672147B2 (ja) | 胸部診断支援情報生成システム | |
JP5962237B2 (ja) | 胸部診断支援情報生成方法 | |
JP5874636B2 (ja) | 診断支援システム及びプログラム | |
JP6597548B2 (ja) | 動態解析システム | |
JP6772873B2 (ja) | 動態解析装置及び動態解析システム | |
JP6418091B2 (ja) | 胸部画像表示システム及び画像処理装置 | |
JP5556413B2 (ja) | 動態画像処理装置及びプログラム | |
JP6743662B2 (ja) | 動態画像処理システム | |
JP2012110400A (ja) | 動態診断支援情報生成システム | |
JP2009153678A (ja) | 動態画像処理システム | |
JP2009153677A (ja) | 動態画像処理システム | |
JP2017176202A (ja) | 動態解析システム | |
US20180018772A1 (en) | Dynamic analysis apparatus | |
JP6662428B2 (ja) | 動態解析システム | |
JP5617577B2 (ja) | 胸部診断支援情報生成方法 | |
JP2017169830A (ja) | 動態解析装置 | |
JP2019122449A (ja) | 動態画像解析装置 | |
JP5625799B2 (ja) | 動態診断支援情報生成システム | |
JP2020044445A (ja) | 動態解析システム、プログラム及び動態解析装置 | |
JP2012115582A (ja) | 胸部診断支援システム | |
JP2014147844A (ja) | 動態診断支援情報生成システム | |
JP2019092612A (ja) | 動態撮影システム | |
JP2019010135A (ja) | 動態解析装置及び動態解析システム | |
JP2022135243A (ja) | 動態解析装置及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201180041386.5 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11819622 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2012530550 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13819195 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 11819622 Country of ref document: EP Kind code of ref document: A1 |