JP2020168172A - Dynamic analysis apparatus, dynamic analysis system, and program - Google Patents

Abstract

To make it possible to grasp a state of a respiratory function by making it possible to grasp, in time series, relation between movement of a plurality of points on the ribs from a chest dynamic image.SOLUTION: A diagnostic console 3 includes a control unit 31 for: tracking a plurality of measurement points, which are set one by one on two ore more rib areas in an arbitrary frame image on a chest dynamic image, in frame images in the same section; acquiring a plurality of tracking points corresponding to the respective measurement points, in the time direction; and, for each of the acquired tracking points, calculating a movement amount and/or movement direction from the tracking point adjacent in the time direction to acquire time series data on movement amounts and/or movement directions of the plurality of measurement points.SELECTED DRAWING: Figure 3

Description

The present invention relates to a dynamic analysis device, a dynamic analysis system and a program.

Conventionally, a bone movement vector or movement amount map is created based on an X-ray motion image related to a bone shadow of the chest, and an image obtained by superimposing the calculated movement vector or movement amount map on one X-ray still image is generated. The technique is known (see, for example, Patent Document 1).

JP-A-2015-43894

By the way, the movement of a plurality of ribs contributes to the respiratory function, and grasping the relationship between the movements (amount and direction) of a plurality of ribs in chronological order is the state of the respiratory function (normal / abnormal, degree of disease). Etc.) is effective for diagnosis. However, in Patent Document 1, information on the relationship between the movements of a plurality of bones in a time series is missing. Therefore, it was insufficient to grasp the state of respiratory function.

An object of the present invention is to make it possible to grasp the state of respiratory function by making it possible to grasp the relationship between the movements of a plurality of ribs in a time series from a chest dynamic image.

In order to solve the above problems, the dynamic analysis device of the invention according to claim 1 is used.
A setting means for setting measurement points one by one on two or more rib regions of an arbitrary frame image of a chest dynamic image acquired by radiography of the chest dynamics.
The set plurality of measurement points are tracked in the frame image of the same section of the chest dynamic image, and a plurality of tracking points corresponding to each of the plurality of measurement points are acquired in the time direction, and for each acquired tracking point, A calculation means for calculating the movement amount and / or the movement direction from the tracking points adjacent to the time direction and acquiring the time-series data of the movement amount and / or the movement direction of the plurality of measurement points.
To be equipped.

The invention according to claim 2 is the invention according to claim 1.
The calculation means further calculates an index value indicating the synchronization of movements of the plurality of measurement points based on the time series data of the plurality of measurement points.

The invention according to claim 3 is the invention according to claim 1 or 2.
A display means for displaying the calculation result by the calculation means is provided.

The dynamic analysis system of the invention according to claim 4 is
An imaging device that acquires a chest dynamic image by radiographically photographing the chest dynamics,
The dynamic analysis apparatus according to any one of claims 1 to 3,
To be equipped.

The program of the invention according to claim 5
Computer,
A setting means for setting measurement points one by one on two or more rib regions of an arbitrary frame image of a chest dynamic image acquired by radiography of the chest dynamics.
The set plurality of measurement points are tracked in the frame image of the same section of the chest dynamic image, and a plurality of tracking points corresponding to each of the plurality of measurement points are acquired in the time direction, and for each acquired tracking point, A calculation means for calculating the movement amount and / or the movement direction from tracking points adjacent to the time direction and acquiring the time-series data of the movement amount and / or the movement direction of the plurality of measurement points.
To function as.

According to the present invention, since the relationship between the movements of a plurality of ribs can be grasped in time series from the chest dynamic image, it is possible to grasp the state of respiratory function.

It is a figure which shows the whole structure of the dynamic analysis system in embodiment of this invention. It is a flowchart which shows the shooting control processing executed by the control part of the shooting console of FIG. It is a flowchart which shows the image analysis processing which is executed by the control part of the diagnostic console of FIG. It is a figure which shows an example of setting of a measurement point. It is a figure for demonstrating the calculation method of the movement amount and movement direction. (A) is a graph showing the diaphragm position calculated from each frame image of the chest dynamic image in time series, and (b) is a graph showing the differential value of the diaphragm position in (a) in time series. It is a figure which shows the synchronism of the movement direction of the upper, middle, and lower posterior ribs of a respiratory normal person, a COPD mild patient, and a COPD severe patient. It is an image diagram showing the movement direction of the upper, middle, and lower posterior ribs of a respiratory normal person, a COPD mild patient, and a COPD severe patient with arrows. It is a figure which shows the measurement point set on the upper, middle, and lower posterior ribs. In (a), when measurement points are set on the upper, middle, and lower posterior ribs in the chest dynamic image of a normal breathing person, the movement directions of the measurement points calculated by the image analysis process of FIG. 3 are arranged in time series. The graphs shown side by side, (b) is a diagram showing the cosine similarity in the movement directions of the upper and middle posterior ribs, the middle and lower posterior ribs, and the upper and lower posterior ribs of the inspiratory phase shown in (a). is there. In (a), when measurement points are set on the upper, middle, and lower posterior ribs in the chest dynamic image of a patient with mild COPD, the movement directions of the measurement points calculated by the image analysis process of FIG. 3 are arranged in time series. A side-by-side graph, (b) is a diagram showing the cosine similarity in the movement directions of the upper and middle posterior ribs, the middle and lower posterior ribs, and the upper and lower posterior ribs of the inspiratory phase shown in (a). is there. In (a), when measurement points are set on the upper, middle, and lower posterior ribs in the chest dynamic image of a severely ill COPD patient, the movement directions of the measurement points calculated by the image analysis process of FIG. 3 are arranged in time series. A side-by-side graph, (b) is a diagram showing the cosine similarity in the movement directions of the upper and middle posterior ribs, the middle and lower posterior ribs, and the upper and lower posterior ribs of the inspiratory phase shown in (a). is there.

Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. However, the scope of the invention is not limited to the illustrated examples.

[Configuration of dynamic analysis system 100]
First, the configuration will be described.
FIG. 1 shows the overall configuration of the dynamic analysis system 100 according to the present embodiment.
As shown in FIG. 1, in the dynamic analysis system 100, the imaging device 1 and the 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) or the like. It is configured to be connected via the communication network NT of. Each device constituting the dynamic analysis system 100 conforms to the DICOM (Digital Image and Communications in Medicine) standard, and communication between the devices is performed according to DICOM.

[Structure of imaging device 1]
The imaging device 1 is an imaging means for photographing the dynamics of the chest having a periodicity (cycle), such as morphological changes of lung expansion and contraction associated with respiratory movements and heartbeat. In dynamic photography, radiation such as X-rays is pulsed and repeatedly irradiated at predetermined time intervals (pulse irradiation), or low dose rate is continuously irradiated (continuous irradiation). By doing so, it means to acquire a plurality of images showing the dynamics. A series of images obtained by dynamic photography is called a dynamic image. Further, each of the plurality of images constituting the dynamic image is called a frame image. In the following embodiment, a case where dynamic imaging of the chest is performed by pulse irradiation will be described as an example.

The radiation source 11 is arranged at a position facing the radiation detection unit 13 with the subject M interposed therebetween, and irradiates the subject M with radiation (X-rays) under 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 irradiation conditions input from the imaging console 2 to perform radiation imaging. Irradiation conditions input from the imaging console 2 include, for example, pulse rate, pulse width, pulse interval, number of imaging frames per imaging, X-ray tube current value, X-ray tube voltage value, additional filter type, etc. Is. The pulse rate is the number of irradiations per second, which is consistent with the frame rate described later. The pulse width is the irradiation time per irradiation. The pulse interval is the time from the start of one irradiation to the start of the next irradiation, and is consistent with the frame interval described later.

The radiation detection unit 13 is composed of a semiconductor image sensor such as an FPD. The FPD has, for example, a glass substrate or the like, detects radiation emitted from the radiation source 11 at a predetermined position on the substrate and transmitted through at least the subject M according to its intensity, and detects the detected radiation as an electric signal. A plurality of detection elements (pixels) that are converted into and accumulated in a matrix are arranged in a matrix. Each pixel is configured to include a switching unit such as a TFT (Thin Film Transistor). The FPD has an indirect conversion type that converts X-rays into an electric signal by a photoelectric conversion element via a scintillator and a direct conversion type that directly converts X-rays into an electric signal, and either of them may be used.
The radiation detection unit 13 is provided so as to face the radiation source 11 with the subject M interposed therebetween.

The reading control device 14 is connected to the photographing 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 photographing console 2 to switch the reading of the electric signal stored in each pixel. Then, the image data is acquired by reading the electric signal stored in the radiation detection unit 13. This image data is a frame image. Then, the reading control device 14 outputs the acquired frame image to the shooting 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, which is consistent with 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, and is consistent with the pulse interval.

Here, 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.

[Configuration of shooting console 2]
The imaging console 2 outputs radiation irradiation conditions and image reading conditions to the imaging device 1 to control radiation imaging and radiation image reading operations by the imaging device 1, and a camera operator captures a dynamic image acquired by the imaging device 1. It is displayed for confirmation of positioning by the photographer and confirmation of whether or not the image is suitable for diagnosis.
As shown in FIG. 1, the photographing 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) and a RAM (Random Access Memory).
) Etc. The CPU of the control unit 21 reads out the system program and various processing programs stored in the storage unit 22 in response to the operation of the operation unit 23, develops them in the RAM, and performs the shooting control processing described later according to the expanded programs. Various processes including the above are executed to centrally control the operation of each part of the photographing console 2 and the irradiation operation and reading operation of the photographing device 1.

The storage unit 22 is composed of a non-volatile semiconductor memory, a hard disk, or the like. The storage unit 22 stores data such as parameters or processing results necessary for executing processing by various programs and programs executed by the control unit 21. For example, the storage unit 22 stores a program for executing the photographing control process shown in FIG. In addition, the storage unit 22 stores the irradiation conditions and the image reading conditions corresponding to the imaging site (here, the chest). Various programs are stored in the form of a readable program code, and the control unit 21 sequentially executes an operation according to the program code.

The operation unit 23 is configured to include a keyboard equipped with cursor keys, number input keys, various function keys, and a pointing device such as a mouse, and controls instruction signals input by key operations on the keyboard or mouse operations. Output to 21. Further, the operation unit 23 may include a touch panel on the display screen of the display unit 24, and 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 is composed of 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 according to instructions of display signals input from the control unit 21. To do.

The communication unit 25 includes a LAN (Local Area Network) adapter, a modem, a TA (Terminal Adapter), and the like, and controls data transmission / reception with each device connected to the communication network NT.

[Configuration of diagnostic console 3]
The diagnostic console 3 is a dynamic analysis device for acquiring a dynamic image from the imaging console 2 and displaying the acquired dynamic image and the analysis result of the dynamic image to support the diagnosis of a doctor.
As shown in FIG. 1, 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 is composed of a CPU, 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 response to the operation of the operation unit 33, develops them in the RAM, and performs image analysis described later according to the expanded programs. Various processes including the process are executed, and the operation of each part of the diagnostic console 3 is centrally controlled. The control unit 31 functions as a setting means and a calculation means.

The storage unit 32 is composed of a non-volatile semiconductor memory, a hard disk, or the like. The storage unit 32 stores data such as parameters or processing results necessary for executing the processing by various programs and programs including a program for executing the image analysis processing in the control unit 31. These various programs are stored in the form of a readable program code, and the control unit 31 sequentially executes an operation according to the program code.

The operation unit 33 is configured to include a keyboard equipped with cursor keys, number input keys, various function keys, and a pointing device such as a mouse, and controls an instruction signal input by key operation on the keyboard or mouse operation. Output to 31. Further, the operation unit 33 may include a touch panel on the display screen of the display unit 34, and in this case, the operation unit 33 outputs an instruction signal input via the touch panel to the control unit 31.

The display unit 34 is composed of a monitor such as an LCD or a CRT, and performs various displays according to 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.

[Operation of dynamic analysis system 100]
Next, the operation in the dynamic analysis system 100 will be described.

(Operation of shooting device 1 and shooting console 2)
First, a shooting operation by the shooting device 1 and the shooting console 2 will be described.
FIG. 2 shows a shooting control process executed by the control unit 21 of the shooting console 2. The photographing control process is executed in collaboration with the control unit 21 and the program stored in the storage unit 22.

First, the operator operates the operation unit 23 of the imaging console 2, and the patient information (patient's name, height, weight, age, gender, etc.) and the imaging site (here, the chest) of the imaging target (subject M) are operated. Is input (step S1).

Next, the radiation irradiation condition is read from the storage unit 22 and set in the radiation irradiation control device 12, and the image reading condition is read from the storage unit 22 and set in the reading control device 14 (step S2).

Next, the instruction of radiation irradiation by the operation of the operation unit 23 is waited for (step S3). Here, the photographer arranges the subject M between the radiation source 11 and the radiation detection unit 13 for positioning. Further, in the present embodiment, since the image is taken under the breathing state, the subject (subject M) is instructed to make it easier and the patient is encouraged to breathe at rest. Alternatively, deep breathing may be instructed. When the shooting preparation is completed, the operation unit 23 is operated to input the irradiation instruction.

When the radiation irradiation instruction is input by the operation unit 23 (step S3; YES), the imaging start instruction is output to the radiation irradiation control device 12 and the reading control device 14, and the dynamic imaging is started (step S4). That is, radiation is irradiated by the radiation source 11 at the pulse interval set in the radiation irradiation control device 12, and a frame image is acquired by the radiation detection unit 13.

When the imaging of a predetermined number of frames is completed, the control unit 21 outputs an instruction to end the imaging to the radiation irradiation control device 12 and the reading control device 14, and the imaging operation is stopped. The number of frames captured is the number of frames that can be captured in at least one breath cycle.

The frame images acquired by shooting are sequentially input to the shooting console 2, stored in the storage unit 22 in association with the number indicating the shooting order (frame number) (step S5), and displayed on the display unit 24. (Step S6). The photographer confirms the positioning and the like from the displayed dynamic image, and determines whether the image suitable for the diagnosis is acquired by the image (shooting OK) or the re-shooting is necessary (shooting NG). Then, the operation unit 23 is operated to input the determination result.

When a determination result indicating that shooting is OK is input by a predetermined operation of the operation unit 23 (step S7; YES), an identification ID for identifying the dynamic image and each of the series of frame images acquired in the dynamic shooting are used. , Patient information, imaging site, irradiation condition, image reading condition, number indicating imaging order (frame number), etc. are attached (for example, written in the header area of image data in DICOM format), and the communication unit 25 is It is transmitted to the diagnostic console 3 via (step S8). Then, this process ends. On the other hand, when a determination result indicating shooting NG is input by a predetermined operation of the operation unit 23 (step S7; NO), a series of frame images stored in the storage unit 22 are deleted (step S9), and this process is performed. finish. In this case, re-shooting is required.

(Operation of diagnostic console 3)
Next, the operation in the diagnostic console 3 will be described.
In the diagnostic console 3, when a series of frame images of dynamic images are received from the photographing console 2 via the communication unit 35, the figure is shown in collaboration with the program stored in the control unit 31 and the storage unit 32. The image analysis process shown in 3 is executed.

Hereinafter, the flow of the image analysis process will be described with reference to FIG.
First, 1 is set in the variable n (step S11).

Next, a plurality of measurement points are set on the rib region of the nth frame image (step S12).
Here, measurement points are set one by one on two or more rib regions of the nth frame image.
For example, as shown in FIG. 4, the first frame image is displayed on the display unit 34, and a plurality of points designated one by one on the plurality of rib regions by the operation of the operation unit 33 by the user are set as measurement points. Set.
Alternatively, the control unit 31 may automatically set a plurality of points on the rib region as measurement points. For example, the rib region is extracted from the first frame image, and the representative points on the extracted rib region (for example, the 0th rib, the △ rib, the most curved point of the □ rib, etc.) are automatically extracted. The measurement point may be set. The rib region can be extracted by a known image processing technique, for example, as in the rib extraction method using the model function and the Sobel operator described in JP-A-5-176919.
In the present embodiment, the case where the measurement point is set in the first frame image will be described as an example, but the measurement point can be set in any frame image, not limited to the first image.

Next, n + 1 is set in n (step S13), each measurement point is tracked in the nth frame image, and the tracking points corresponding to each measurement point are extracted (step S14).
For example, in the nth frame image, the pixel extracted as the tracking point in the n-1st frame image (in the case of n = 2, the pixel set as the measurement point in the n-1st frame image) is the center. A search area is set in an area consisting of M × N pixels (M and N are positive integers, for example, 3 × 3), and the pixel value is extracted as a tracking point in the n-1th frame image in the search area. The pixel closest to the pixel value of the pixel is extracted as a tracking point.

Next, the movement amount and / or the movement direction from the tracking point extracted in the n-1st frame image is calculated for each tracking point corresponding to each measurement point extracted in the nth frame image (step). S15).
For example, as shown in FIG. 5, when the tracking point moves from At to At + 1, the distance from At to At + 1 is the amount of movement, and the angle between the line connecting At and At + 1 and the x-axis direction. θ is calculated as the moving direction. As a result, the movement amount and / or the movement direction between the frame images of each measurement point is calculated.

Next, the calculated movement amount and / or movement direction of each measurement point is recorded in the RAM as time-series data in association with the frame number (step S16).

Steps S13 to S16 are repeatedly executed until n reaches the number of frame images of the dynamic image, and when n reaches the number of frame images (step S17; YES), the calculation result of the movement amount and / or the movement direction of the plurality of measurement points. (Time series data) is graphed on the same coordinate plane and displayed on the display unit 34 (step S18).

Further, based on the time-series data of the movement amount and / or the movement direction of the plurality of measurement points, an index value indicating the synchronization of the movement of the measurement points is calculated for each set of the two measurement points (step S19). An index value indicating the synchronization of the calculated movements of the measurement points is displayed on the display unit 34 (step S20), and the image analysis process ends.
The index value indicating the synchronization of the movements of the two measurement points is an index value indicating how similar the movements (movement amount and movement direction) of the two measurement points are, for example, of the two measurement points. The mutual correlation coefficient and cosine similarity of time series data can be mentioned.

Here, the graph display and the calculation of the index value may be displayed and calculated separately for the expiratory phase and the inspiratory phase.
The respiratory phase can be specified, for example, by the following procedures (1) to (5). In the following description, the upper left coordinate of the image is set to (0, 0), and the coordinate value increases toward the right side and the lower side of the image.
(1) First, the position of the diaphragm is calculated from each frame image of the dynamic image (see FIG. 6A).
For example, the lower edge portion of the lung field region of each frame image is extracted as the diaphragm boundary portion, a reference point is set at a certain x-coordinate position of the diaphragm boundary portion, and the y-coordinate of the set reference point is set as the diaphragm. It can be obtained as a position. The lung field region may be extracted using any known method. For example, a threshold value is obtained by discriminant analysis from a histogram of pixel values of each pixel, and a region having a signal higher than this threshold value is first-order extracted as a lung field region candidate. Next, edge detection is performed near the boundary of the primary extracted lung field region candidate, and the boundary of the lung field region can be extracted by extracting the point where the edge is maximum in the small region near the boundary along the boundary. it can.
(2) Calculate the highest and lowest positions of the diaphragm.
(3) Calculate the differential value of (1) (see FIG. 6 (b)).
(4) Regarding the expiratory phase, the first inflection point (point with a differential value of 0) where the differential value changes from “+” to “-” after the highest frame image is set as the “expiratory phase start”, and the expiratory phase starts. Later, the inflection point where "-" → "+" first changes is called "expiratory phase end". In consideration of the influence of noise, "i-frame (i is a positive integer) derivative value 0 continues" may be added to the start / end condition.
(5) Regarding the intake phase, the first inflection point where the differential value becomes “-” → “+” after the lowest frame image is set as “intake phase start”, and the first “+” → after the start of the intake phase. The inflection point that becomes "-" is defined as "end of intake phase". In consideration of the influence of noise, "the i-frame derivative value 0 continues" may be added to the start / end condition.

As described above, in the above image analysis process, a plurality of measurement points set one by one on two or more rib regions of an arbitrary frame image on the chest dynamic image are tracked in the frame image of the same section, and a plurality of measurement points are tracked. A plurality of tracking points corresponding to each of the measurement points of the above are acquired in the time direction, and for each acquired tracking point, the movement amount and / or the movement direction from the tracking points adjacent to the time direction is calculated, and the plurality of measurement points are obtained. Acquire time-series data of the movement amount and / or the movement direction of.
Therefore, users such as doctors can grasp the relationship between the movements of multiple ribs in the chest dynamic image in chronological order, and can grasp the state of respiratory function (normal, abnormal, degree of disease, etc.). Become. Therefore, it is possible to reduce the burden on the patient due to the respiratory examination. In addition, the time series data can be used for confirming the effect of the treatment effect, and can support the formulation of the treatment plan.

In addition, the movement amount and / or movement direction calculation result (time series data) of a plurality of measurement points can be graphed on the same coordinate plane and displayed on the display unit 34, or a plurality of measurement points can be displayed based on the time series data. By calculating an index value indicating the synchronism of the movements of the above and displaying the calculation result on the display unit 34, the user can more easily grasp the state of the respiratory function.

Here, an example of the calculation result (graph, index value indicating synchronism) displayed by the image analysis process will be described with reference to a specific example.
FIG. 7 is a diagram showing the synchronization relationship of the movement directions of the upper, middle, and lower posterior ribs (see FIG. 9) of a respiratory normal person, a COPD mild patient, and a COPD severe patient in the inspiratory phase, and FIG. 8 shows the movement. It is a figure which shows the direction schematically by an arrow. As shown in FIGS. 7 and 8, the synchrony of the movement directions of the upper, middle, and lower posterior ribs in the inspiratory phase is different in the respiratory normal person, the COPD mild patient, and the COPD severe patient. For normal breathing individuals, the upper and middle posterior ribs and the lower posterior ribs tend to move in different directions during the inspiratory / expiratory phase. For mildly ill patients with COPD, the middle posterior rib and the lower posterior rib tend to move synchronously in the same direction during the inspiratory phase. For critically ill patients with COPD, the upper, middle, and lower posterior ribs tend to move synchronously in the inspiratory phase.

FIG. 10A shows an image of FIG. 3 when measurement points are set at three points above, in, and below the posterior ribs as shown in FIG. 9 with respect to a chest dynamic image of a certain respiratory normal person. It is a graph which showed the moving direction of each measurement point calculated by an analysis process in chronological order. FIG. 10B is a diagram showing the cosine similarity in the movement directions of the upper and middle posterior ribs, the middle and lower posterior ribs, and the upper and lower posterior ribs of the inspiratory phase shown in FIG. 10A. As shown in FIGS. 10A and 10B, when the chest dynamic image of the respiratory normal person is analyzed by the above image analysis processing, the synchronization (similarity) of the movement directions of the upper and middle posterior ribs is obtained in the inspiratory phase. The results show that the high and upper and lower, middle and lower posterior ribs are less synchronous (similar) in the direction of movement.

FIG. 11A shows an image of FIG. 3 when measurement points are set at three points above, in, and below the posterior ribs as shown in FIG. 9 with respect to a chest dynamic image of a certain COPD mild patient. It is a graph which showed the moving direction of each measurement point calculated by an analysis process in chronological order. FIG. 11B is a diagram showing the cosine similarity in the movement directions of the upper and middle posterior ribs, the middle and lower posterior ribs, and the upper and lower posterior ribs shown in FIG. 11A. As shown in FIGS. 11A and 11B, when the chest dynamic image of a patient with mild COPD was analyzed by the above image analysis processing, the synchronization (similarity) of the movement directions of the middle and lower posterior ribs was observed in the inspiratory phase. The results show that it is high, and the synchronization (similarity) of the movement directions of the posterior ribs above and below and above and inside is lower.

FIG. 12A shows an image of FIG. 3 when measurement points are set at three points above, in, and below the posterior ribs as shown in FIG. 9 with respect to a chest dynamic image of a certain severely ill COPD patient. It is a graph which showed the moving direction of each measurement point calculated by an analysis process in chronological order. FIG. 12 (b) is a diagram showing the cosine similarity in the movement directions of the upper and middle posterior ribs, the middle and lower posterior ribs, and the upper and lower posterior ribs shown in FIG. 12 (a). As shown in FIGS. 12 (a) and 12 (b), when the chest dynamic image of the severely ill COPD patient was analyzed by the above image analysis processing, all the synchronizations of the movement directions of the upper, middle, and lower posterior ribs in the inspiratory phase ( The result is that it is highly similar.

That is, a graph showing the moving directions of the measurement points (upper, middle, lower) of the posterior ribs obtained by analyzing the chest dynamic image of the subject by the above image analysis process, and a plurality of measurements. By observing the index value indicating the synchrony of the movement direction of the points, the user can grasp whether the subject's breathing is normal, COPD is mild, COPD is severe, and the like.

Although the embodiment of the present invention has been described above, the description in the present embodiment is an example of a suitable dynamic analysis system according to the present invention, and is not limited thereto.
For example, in the above embodiment, a plurality of measurement points are tracked from the frame image of the entire section of the chest dynamic image to acquire time series data, but the present invention is not limited to this, and a part of the chest dynamic image is obtained. Time series data may be acquired by tracking a plurality of measurement points from the frame image of the section. In this case, tracking of a plurality of measurement points is performed in the same section of the chest dynamic image, and time series data of the same section is acquired. In the present invention, tracking and time-series data may be acquired in the same section at a plurality of set measurement points, and in addition, tracking and time-series data may be performed in other sections at some measurement points. It does not prevent you from getting.

Further, for example, in the above description, an example in which a hard disk, a non-volatile memory of a semiconductor, or the like is used as a computer-readable medium of the program according to the present invention has been disclosed, but the present invention is not limited to this example. CD-ROM as another computer-readable medium
It is possible to apply a portable recording medium such as. A carrier wave is also applied as a medium for providing data of a program according to the present invention via a communication line.

In addition, the detailed configuration and detailed operation of each device constituting the dynamic analysis system 100 can be appropriately changed without departing from the spirit of the present invention.

100 Dynamic analysis system 1 Imaging device 11 Radiation source 12 Radiation irradiation control device 13 Radiation detection unit 14 Reading control device 2 Imaging console 21 Control unit 22 Storage unit 23 Operation unit 24 Display unit 25 Communication unit 26 Bus 3 Diagnostic console 31 Control Unit 32 Storage unit 33 Operation unit 34 Display unit 35 Communication unit 36 Bus

Claims (5)

A setting means for setting measurement points one by one on two or more rib regions of an arbitrary frame image of a chest dynamic image acquired by radiography of the chest dynamics.
The set plurality of measurement points are tracked in the frame image of the same section of the chest dynamic image, and a plurality of tracking points corresponding to each of the plurality of measurement points are acquired in the time direction, and for each acquired tracking point, A calculation means for calculating the movement amount and / or the movement direction from the tracking points adjacent to the time direction and acquiring the time-series data of the movement amount and / or the movement direction of the plurality of measurement points.
A dynamic analysis device including. The dynamic analysis device according to claim 1, wherein the calculation means further calculates an index value indicating synchronization of movements of the plurality of measurement points based on time-series data of the plurality of measurement points. The dynamic analysis apparatus according to claim 1 or 2, further comprising a display means for displaying the calculation result by the calculation means. An imaging device that acquires a chest dynamic image by radiographically photographing the chest dynamics,
The dynamic analysis apparatus according to any one of claims 1 to 3,
A dynamic analysis system equipped with. Computer,
A setting means for setting measurement points one by one on two or more rib regions of an arbitrary frame image of a chest dynamic image acquired by radiography of the chest dynamics.
The set plurality of measurement points are tracked in the frame image of the same section of the chest dynamic image, and a plurality of tracking points corresponding to each of the plurality of measurement points are acquired in the time direction, and for each acquired tracking point, A calculation means for calculating the movement amount and / or the movement direction from tracking points adjacent to the time direction and acquiring the time-series data of the movement amount and / or the movement direction of the plurality of measurement points.
A program to function as.

Images