JP2022135243A - Dynamic analysis apparatus and program - Google Patents

Abstract

To efficiently calculate a V/Q ratio.SOLUTION: A control unit of a diagnostic console acquires, from among a plurality of frame images obtained by radiographing the dynamic state of the chest of a subject, a plurality of frame images for ventilation analysis and a plurality of frame images for perfusion analysis, selects a reference frame from among the plurality of acquired frame images, divides a lung field area in the reference frame image into a plurality of small areas, and tracks, in a frame image other than the reference frame image, each of the small areas in the reference frame image to identify a position corresponding to the small area in the frame image other than the reference frame image, thereby setting a plurality of small areas corresponding between the acquired frame images. For each of the set small areas, a V/Q ratio is calculated.SELECTED DRAWING: Figure 3

Description

There is an application for exception to loss of novelty

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

Lung functions include ventilation and perfusion. Ventilation is the function of taking in oxygen as inspiration and sending it to the alveoli, and discharging carbon dioxide from the alveoli to the outside of the body as expiration. Blood flow is a function of sending blood from which carbon dioxide has been discharged from cells throughout the body to the alveoli, and sending oxygenated blood to the heart through gas exchange in the alveoli. In other words, it is important for gas exchange that ventilation and blood flow function in a well-balanced manner, and gas exchange will not go well if either function deteriorates.

The relationship between ventilation (V) and blood flow (Q) is important for pathophysiological understanding of respiratory diseases, and is clinically evaluated by the following VQ ratio.
VQ ratio = ventilation (V) / blood flow (Q)

For example, Patent Literature 1 describes calculating the VQ ratio for each small region from the ventilation analysis result and blood flow analysis result for each small region of the same frame image among a series of frame images showing the dynamics of the chest. It is

Japanese Unexamined Patent Application Publication No. 2017-202208

In the technique described in Patent Document 1, the VQ ratio is calculated from the ventilation analysis result and the blood flow analysis result for each same small region of the same frame image. That is, the VQ ratio is calculated using the ventilation analysis results and the blood flow analysis results of small regions located at the same pixel position.

However, in some cases, frame images during deep breathing or quiet breathing are used for ventilation analysis, and frame images during breath holding are used for blood flow analysis to calculate the VQ ratio. As shown in FIG. 11, the position where the same portion of the lung field (indicated by R in FIG. 11) is shown in the frame image for ventilation analysis (breathing image) and the frame image for blood flow analysis (breath-holding image). are different, an erroneous VQ ratio will be calculated if the VQ ratio is calculated using the ventilation analysis result and the blood flow analysis result of small regions located at the same pixel position.

In addition, in Patent Document 1, local matching and warping processing are performed to correct the positional deviation of the lung region between frame images, and then the division into small regions is performed, and VQ is performed between the small regions at the same pixel position. Calculating the ratio is also described, but there is a problem that the warping process takes a long time and is inefficient.

An object of the present invention is to enable efficient calculation of the VQ ratio.

In order to solve the above problems, the dynamic analysis device according to claim 1 is
Acquisition means for acquiring a plurality of frame images for ventilation analysis and a plurality of frame images for blood flow analysis from a plurality of frame images obtained by radiography of dynamics of the subject's chest;
a division means for selecting a reference frame image from among the plurality of frame images for ventilation analysis and blood flow analysis acquired by the acquisition means, and dividing a lung region in the reference frame image into a plurality of small regions;
tracking each of the plurality of small regions of the reference frame image in a frame image other than the reference frame image to specify a position corresponding to the small region in the other frame image; tracking means for setting a plurality of corresponding small areas between the frame images acquired by the acquisition means;
For each of the plurality of set small regions, an index value indicating a ventilation volume is calculated using the plurality of frame images for ventilation analysis, and a blood flow volume is calculated using the plurality of frame images for blood flow analysis. an analysis means for calculating an index value indicating
a ratio calculation means for calculating a ratio between an index value indicating the ventilation volume and an index value indicating the blood flow for each of the plurality of small regions calculated by the analysis means;
Prepare.

The invention according to claim 2 is the invention according to claim 1,
The acquiring means acquires a plurality of frame images obtained by imaging in a breathing state as a plurality of frame images for ventilation analysis, and acquires a plurality of frame images obtained by imaging in a breath-holding state as the plurality of frame images for ventilation analysis. Acquire multiple frame images for blood flow analysis.

The invention according to claim 3 is the invention according to claim 1 or 2,
The number of small regions set in the plurality of frame images for ventilation analysis and the plurality of frame images for blood flow analysis is the same.

The invention according to claim 4 is the invention according to any one of claims 1 to 3,
The dividing means divides the reference frame image so that the sizes of the small regions differ according to the positions of the plurality of small regions in the lung region.

The invention according to claim 5 is the invention according to claim 4,
The dividing means divides the reference frame image so that the size of the small region in the peripheral portion is larger than the size of the small region in the central portion of the lung field region.

The invention according to claim 6 is the invention according to claim 4,
The dividing means divides the reference frame image so that the size of the lower small region is smaller than the size of the upper small region of the lung field region.

The invention according to claim 7 is the invention according to any one of claims 1 to 6,
Division for determining whether or not the aspect ratio of each of the plurality of small areas set by the tracking means falls within a reference range, and dividing the small area determined not to fit or integrating it into adjacent small areas Equipped with integration means.

The invention according to claim 8 is the invention according to claim 7,
The dividing and integrating means uses a frame image captured in a respiratory state among the plurality of frame images acquired by the acquiring means to set the aspect ratio of each of the plurality of small regions set by the tracking means to the Determine whether or not it falls within the reference range.

The invention according to claim 9 is the invention according to claim 7 or 8,
The division/integration means adjusts the number of the plurality of small areas of all the frame images acquired by the acquisition means to be the same after the division or integration.

The invention according to claim 10 is the invention according to any one of claims 1 to 9,
determining whether or not the aspect ratio of each of the plurality of small areas set by the tracking means falls within a reference range; A deforming means for deforming is provided.

The invention according to claim 11 is the invention according to any one of claims 1 to 10,
An adjusting means is provided for adjusting the size or shape of the small areas when the lung field area of the reference frame image is divided into the plurality of small areas by the dividing means.

The invention according to claim 12 is the invention according to claim 11,
The adjusting means determines that the calculated index value indicating the ventilation volume, the index value indicating the blood flow, or the ratio between the index value indicating the ventilation volume and the index value indicating the blood flow is an abnormal value. If so, make the above adjustments.

The invention according to claim 13 is the invention according to claim 11,
The adjustment means performs the adjustment in accordance with a user's operation.

The invention according to claim 14 is the invention according to any one of claims 1 to 12,
The number of divisions when the reference frame image is divided into the plurality of small areas by the division means can be selected from a plurality of predetermined numbers of divisions.

The invention according to claim 15 is the invention according to any one of claims 1 to 14,
The tracking means sets a marker in each of the plurality of small regions of the reference frame image, and tracks the set marker in a frame image other than the reference frame image, thereby position corresponding to the small area in the frame image of .

The invention according to claim 16 is the invention according to any one of claims 1 to 15,
An output means for outputting a calculation result by the ratio calculation means is provided.

The program of the invention according to claim 17,
the computer,
Acquisition means for acquiring a plurality of frame images for ventilation analysis and a plurality of frame images for blood flow analysis from a plurality of frame images obtained by radiography of dynamics of the subject's chest;
dividing means for selecting a reference frame image from among the plurality of frame images for ventilation analysis and blood flow analysis acquired by the acquisition means, and dividing the lung region in the reference frame image into a plurality of small regions;
tracking each of the plurality of small regions of the reference frame image in a frame image other than the reference frame image to specify a position corresponding to the small region in the other frame image; Tracking means for setting a plurality of corresponding small areas between the frame images acquired by the acquisition means;
For each of the plurality of set small regions, an index value indicating a ventilation volume is calculated using the plurality of frame images for ventilation analysis, and a blood flow volume is calculated using the plurality of frame images for blood flow analysis. analysis means for calculating an index value indicating
ratio calculation means for calculating a ratio between an index value indicating the ventilation volume and an index value indicating the blood flow for each of the plurality of small regions calculated by the analysis means;
function as

According to the present invention, it is possible to efficiently calculate the VQ ratio.

It is a figure showing the whole dynamics analysis system composition in an embodiment of the present invention. 4 is a flowchart showing imaging control processing executed by a control unit of the imaging console in FIG. 1; 2 is a flowchart showing VQ ratio calculation processing executed by a control unit of the diagnostic console of FIG. 1; FIG. 4 is a diagram for explaining a breathing period and a breath-holding period; FIG. 4 is a diagram for explaining the division number of lung fields; FIG. 4 is a diagram schematically showing how a lung region is divided into small regions; (a) is a diagram showing an example of dividing the lung field into equal widths according to the number of vertical and horizontal divisions, and (b) is a diagram showing an example of dividing the lower part of the lung field into finer intervals than the upper part. is. FIG. 4 is a flow chart showing a small area division and integration process executed in step S18 of FIG. 3; FIG. FIG. 4 is a diagram for explaining integration of small areas; FIG. 5 is a diagram showing an output example of a VQ ratio; FIG. 10 is a diagram showing changes in position within the lung region in a breathing image and a breath-holding image;

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 of the embodiment of the present invention will be described.
FIG. 1 shows the overall configuration of a dynamic analysis system 100 in this embodiment.
As shown in FIG. 1, the dynamic analysis system 100 includes an imaging apparatus 1 and an imaging console 2 connected by a communication cable or the like, and an imaging console 2 and a diagnostic console 3 connected via a LAN (Local Area Network) or the like. are connected via a communication network NT. 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 in accordance with DICOM.

<Structure of photographing device 1>
The imaging device 1 is imaging means for imaging the dynamics of the chest having a periodicity (cycle), such as morphological changes in expansion and contraction of the lungs associated with respiratory motion, heartbeats, and the like. Dynamic radiography involves irradiating the subject with pulsed radiation such as X-rays repeatedly at predetermined time intervals (pulse irradiation), or continuously irradiating the subject at a low dose rate without interruption (continuous irradiation). Acquisition of a plurality of images showing the dynamics of the subject. A series of images obtained by dynamic imaging are called dynamic images. Also, each of the plurality of images forming the dynamic image is called a frame image. Here, dynamic images include moving images, but do not include images obtained by capturing still images while displaying moving images.
In the following embodiments, an example of performing dynamic imaging by pulse irradiation will be described.

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 radiation irradiation conditions input from the imaging console 2 to perform radiation imaging. Radiation 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, and the like. is. The pulse rate is the number of radiation exposures per second and matches the frame rate described later. The pulse width is the radiation exposure time per radiation exposure. The pulse interval is the time from the start of one irradiation to the start of the next irradiation, and coincides with the frame interval described later.

The radiation detection unit 13 is composed of a semiconductor image sensor such as FPD. The FPD has, for example, a glass substrate or the like, and detects the 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 converts the detected radiation into an electrical signal. A plurality of detection elements (pixels) are arranged in a matrix. Each pixel includes a switching unit such as a TFT (Thin Film Transistor). FPDs include an indirect conversion type in which X-rays are converted into electric signals by a photoelectric conversion element via a scintillator, and a direct conversion type in which X-rays are directly converted into electric signals. Either type 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 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 conditions input from the imaging console 2 to switch the reading of the electric signal accumulated in each pixel. Then, image data is obtained by reading electrical signals accumulated in the radiation detection unit 13 . This image data is a frame image. Then, the reading control device 14 outputs the obtained frame image to the imaging console 2 . Image reading conditions include, for example, frame rate, frame interval, pixel size, 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 the acquisition operation of one frame image to the start of the acquisition operation of the next frame image, and coincides 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 to synchronize the radiation irradiation operation and the image reading operation.

<Structure 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 radiographic image reading operations by the imaging device 1, and transmits dynamic images acquired by the imaging device 1 to the imaging technician. The image is displayed for confirmation of positioning by the person performing the imaging, and confirmation of whether or not the image is suitable for diagnosis.
As shown in FIG. 1, the imaging console 2 includes a control section 21, a storage section 22, an operation section 23, a display section 24, and a communication section 25, which are connected by a bus 26. FIG.

The control unit 21 is configured by 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 executes shooting control processing described later according to the expanded programs. , and centrally controls the operation of each part of the imaging console 2 and the radiation irradiation operation and reading operation of the imaging apparatus 1 .

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

The operation unit 23 includes a keyboard having cursor keys, numeric input keys, various function keys, etc., and a pointing device such as a mouse. 21. Further, the operation unit 23 may have a touch panel on the display screen of the display unit 24 , and in this case, outputs an instruction signal input through the touch panel to the control unit 21 .

The display unit 24 is configured by a monitor such as an LCD (Liquid Crystal Display) or a CRT (Cathode Ray Tube), and displays input instructions, data, etc. from the operation unit 23 according to instructions of display signals input from the control unit 21. do.

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

<Configuration of Diagnosis Console 3>
The diagnostic console 3 is a dynamic analysis device that acquires a dynamic image from the imaging console 2 and displays the acquired dynamic image and analysis results of the dynamic image to assist the doctor's diagnosis. In this embodiment, the diagnostic console 3 calculates the VQ ratio based on the dynamic chest image and displays the result.
As shown in FIG. 1, the diagnosis console 3 includes a control section 31, a storage section 32, an operation section 33, a display section 34, and a communication section 35, which are connected by a bus .

The control unit 31 is composed of a CPU, a RAM, and the like. The CPU of the control unit 31 reads the system program and various processing programs stored in the storage unit 32 according to the operation of the operation unit 33, expands them in the RAM, and according to the expanded programs, VQ ratio described later Various processing including calculation processing is executed, and the operation of each section of the diagnostic console 3 is centrally controlled. The control unit 31 functions as acquisition means, division means, tracking means, analysis means, ratio calculation means, division/integration means, transformation means, and adjustment means.

The storage unit 32 is configured by a nonvolatile semiconductor memory, hard disk, or the like. The storage unit 32 stores various programs such as a program for executing the VQ ratio calculation processing in the control unit 31, parameters required for execution of processing by the program, or data such as processing results. 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, etc., and a pointing device such as a mouse. 31. Further, the operation unit 33 may have a touch panel on the display screen of the display unit 34 , and in this case, outputs an instruction signal input through the touch panel to the control unit 31 .

The display unit 34 is configured by a monitor such as an LCD or a CRT, and performs various displays according to instructions of display signals input from the control unit 31 . The display unit 34 functions as output means.

The communication unit 35 includes a LAN adapter, modem, TA, etc., and controls data transmission/reception with each device connected to the communication network NT.

<Operation of dynamic analysis system 100>
Next, the operation of the dynamic analysis system 100 will be described.

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

First, the imaging operator operates the operation unit 23 of the imaging console 2, and the patient information (patient's name, height, weight, age, sex, etc.) of the imaging target (subject M) and examination information (examination ID, imaging region) are displayed. (here, chest), breathing category (deep breathing/quiet breathing/breath holding/mixing, etc.), etc.) are input (step S1).

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

Next, a radiation irradiation instruction by operating the operation unit 23 is awaited (step S3). Here, the imaging operator positions the subject M by placing it between the radiation source 11 and the radiation detection unit 13 . When imaging preparations are completed, the operation unit 23 is operated to input a radiation irradiation instruction.

When a radiation irradiation instruction is input by the operation unit 23 (step S3; YES), an imaging start instruction is output to the radiation irradiation control device 12 and the reading control device 14, and dynamic imaging is started (step S4). That is, radiation is emitted from the radiation source 11 at pulse intervals set in the radiation exposure control device 12 , and frame images are acquired by the radiation detection unit 13 . Note that, for example, when the photographing apparatus 1 includes an audio output unit and a display unit, and an instruction to start photographing is output, a respiration-inducing voice such as "breathe in", "breathe out", "hold your breath", etc. or display. After the imaging start instruction is given, radiation may be emitted from the radiation source 11 at a timing at which predetermined respiration guidance (for example, “breathe in”) is performed.

When imaging of a predetermined number of frames is completed, the control unit 21 outputs an imaging end instruction 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 for at least one respiratory cycle.

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

When a determination result indicating OK for photographing is input by a predetermined operation of the operation unit 23 (step S7; YES), an identification ID for identifying the dynamic image and , patient information, examination information (examination ID, imaging site, irradiation conditions, image reading conditions, breathing classification (deep breathing/quiet breathing/breath holding/mixed breathing, etc.)), number indicating the order of imaging (frame number), etc. It is appended (for example, written in the header area of the image data in DICOM format) and transmitted to the diagnostic console 3 via the communication section 25 (step S8). Then, the process ends. Note that when respiratory guidance is output or displayed by voice, information on respiratory guidance that was output when the frame image was captured may be attached to each frame image as supplementary information. In the diagnostic console 3 , the dynamic image transmitted from the imaging console 2 is stored in the storage section 32 by the control section 31 in association with the incidental information.
On the other hand, when a judgment result indicating that the shooting is NG is input by a predetermined operation of the operation unit 23 (step S7; NO), the series of frame images stored in the storage unit 22 are deleted (step S9), and the present process is terminated. finish. In this case, re-shooting is required.

When calculating the VQ ratio, it is necessary to image chest dynamics including breathing state (deep breathing or quiet breathing) for ventilation analysis, and chest dynamics including breath holding state or breathing state for blood flow analysis. . In the examination for calculating the VQ ratio, the breathing state and the breath-holding state may be separately photographed to obtain two dynamic images, or the breathing state and the breath-holding state may be continuously photographed at one time, A single dynamic image may be acquired. Alternatively, dynamic imaging may be performed once in the respiratory state for both ventilation analysis and blood flow analysis.

(Operation of Diagnosis Console 3)
Next, the operation of the diagnostic console 3 will be described.
In the diagnosis console 3, for example, the frame image of the dynamic image acquired by the examination is received from the imaging console 2 via the communication unit 35, and when the operation unit 33 instructs the calculation of the VQ ratio, FIG. is executed. The VQ ratio calculation process is executed by cooperation between the control unit 31 and a program stored in the storage unit 32 .

In the VQ ratio calculation process, first, a plurality of frame images for ventilation analysis and a plurality of frame images for blood flow analysis are acquired from the received dynamic images (step S11).

Here, if the respiratory state and the breath-holding state are captured separately and two dynamic images are received as images captured by the examination, multiple frame images of the dynamic images captured in the respiratory state are used for ventilation analysis. A plurality of frame images are acquired, and a plurality of frame images of dynamic images taken in a breath-holding state are acquired as a plurality of frame images for blood flow analysis.

When the breathing state and the breath-holding state are continuously photographed at one time and one dynamic image is received as an image photographed by the examination, a plurality of frame images of the breathing interval (the interval photographed in the breathing state) are obtained. A plurality of frame images are acquired for ventilation analysis, and a plurality of frame images of a breath-hold interval (a interval captured while holding breath) are acquired as a plurality of frame images for blood flow analysis. For example, as shown in FIG. 4, the diaphragm is extracted from each frame image of the dynamic image to acquire the temporal change in the position of the diaphragm. Images are acquired as a plurality of frame images for ventilation analysis, and a plurality of frame images in a section (breath-hold section) in which the displacement of the diaphragm position is equal to or less than a predetermined threshold are acquired as frame images for blood flow analysis. The diaphragm position can be extracted using a known method, as described in Japanese Patent Application Laid-Open No. 2018-148964, for example. In addition, if the incidental information of each frame image of the dynamic image includes respiratory induction information at the time of capturing the frame image, multiple frame images for ventilation analysis and blood A plurality of frame images for flow analysis may be acquired. Alternatively, for example, a graph showing changes in the position of the diaphragm over time may be displayed on the display unit 34, and the user may use the operation unit 33 to specify a frame image interval for ventilation analysis and a frame image interval for blood flow analysis. good.

When imaging is performed once in the respiratory state for both ventilation analysis and blood flow analysis, and one dynamic image is received as an image taken by examination, the received dynamic image is used for ventilation analysis and blood flow analysis. are acquired as multiple frame images for

Frame images for one respiratory cycle may be obtained as the frame images used for the ventilation analysis, and frame images for one heartbeat cycle may be obtained as the frame images for the blood flow analysis. As the frame image for blood flow analysis, it is preferable to use a frame image captured in a breath-holding state in terms of analysis accuracy and calculation accuracy of the VQ ratio.

Next, a lung region is extracted from each of the acquired plurality of frame images for ventilation analysis and blood flow analysis (step S12).
A known method such as a method using edge detection described in Japanese Patent Application Laid-Open No. 2018-148964 can be used to extract the lung region.

Next, a reference frame image is selected from among the acquired plurality of frame images for ventilation analysis and blood flow analysis (step S13).
The reference frame image is preferably, for example, the frame image at the maximum exhalation position with the smallest lung area. As for the frame image of the maximum expiratory position, for example, the frame image having the smallest area (the number of pixels) of the lung region of each frame image can be extracted as the frame image of the maximum expiratory position.

Next, the division number of the lung region for calculating the VQ ratio is set (step S14).
In step S14, for example, as shown in FIG. 5, the vertical division number M and the horizontal division number N of the lung region are set (M and N are positive integers). For example, the number of divisions selected by the user by operating the operation unit 33 may be set from a plurality of different combinations of M and N prepared in advance, or the numerical values of M and N input by the user may be set. It may be set as a number. Also, a predetermined combination of M and N may be set as the number of divisions.

Next, the lung region of the reference frame image is divided into a plurality of small regions based on the number of divisions set in step S14 (step S15).
For example, the lung region is divided into a plurality of small regions by dividing a rectangular region circumscribing the lung region into the number of divisions M×N. FIG. 6 shows how the lung region is divided into small regions. As shown in FIG. 6, it is preferable to exclude small regions where calculation of the VQ ratio is unnecessary, such as regions outside the lung field, heart region and bronchus region where gas exchange does not occur. In addition, in FIG. 6, the drawing of the inside of the lung region is omitted in order to clearly show how the image is divided into small regions.

As shown in FIG. 7A, the lung field region may be divided into equal widths according to the number of divisions, or may be divided into different widths depending on the position of the lung field region. For example, the upper and lower portions of the lung region may be divided by different widths. For example, since the movement of the lower lung field is large, as shown in FIG. Good (e.g. subdivide the lower third, etc.). In addition, blood vessels are dense in the vicinity of the center of the lung, but the blood vessels in the peripheral part are too thin and blood flow cannot be seen in some cases, making analysis meaningless. Therefore, the lung field may be divided so that the peripheral (outer) small region is larger in size than the central (inner) small region.
As for the size of the small region, a small size (for example, about 1 cm×1 cm) is preferable for detailed evaluation, and a large size (for example, about 3 cm×3 cm) is preferable when priority is given to processing speed.

A unique ID is then assigned to each small region of the reference frame image (step S16).
Each small area is assigned a unique ID according to a predetermined rule. For example, an ID such as 1-5 is assigned to the fifth small area from the left in the first row.

Each small region of the reference frame image is then tracked in another frame image (step S17).
As a tracking method, for example, a method of tracking density information (luminance information) of each small region of the reference frame image in another frame image by optical flow or local matching can be used. For example, for each small area of the reference frame image, a marker is set at the vertex or center of the small area to obtain the density information around the marker, and the pixel position corresponding to the pixel position of the small area is obtained in another frame image. search for an area with density information similar to the density information around the marker, and specify the position with similar density information as the position to which the small area has moved (i.e., the position corresponding to the small area). . With this method, the same part (for example, a structure (blood vessel shadow, bronchial shadow, etc.)) inside the lung field can be tracked, so that small regions can be precisely aligned between frame images.

Alternatively, small regions may be simply tracked based on the lung field contours. For example, trace the lung field contour in another frame image, divide (equally) between the top end and the bottom end of the lung field contour by the division number M, and divide between the right end and the left end of the lung field contour by the division number N split by . Then, the position corresponding to each small region of the reference frame image in the other frame image (for example, the position corresponding to the upper right small region of the reference frame image is the upper right small region of the other frame image) is tracked. is identified as the position of With this tracking method, small areas can be tracked at high speed.

By tracking each small area of the reference frame image in the other frame image, the other frame image is divided into the same number of divisions as the reference frame image, and the corresponding plurality of frame images are obtained among the plurality of frame images acquired in step S11. can be set. In addition, the same number of small regions can be set in the plurality of frame images for ventilation analysis and the plurality of frame images for blood flow analysis.

Next, a small area division and integration process is executed (step S18).
FIG. 8 is a flow chart showing the flow of the small area segmentation and integration process executed in step S18. The small area segmentation and integration process of FIG.

First, a small area is selected and the aspect ratio (horizontal to vertical ratio) is checked (step S181).
Here, if either the vertical or horizontal dimension of a small region becomes extremely small due to movement and deformation of the lung field due to respiration, the aspect ratio becomes extremely large or approaches zero. In such a case, it is considered that the small region has collapsed vertically or horizontally due to the movement and deformation of the lung field due to respiration. is large, and the analysis accuracy and the calculation accuracy of the VQ ratio deteriorate. Also, if either the length or the width is extremely large, the aspect ratio becomes extremely large or becomes infinitely close to 0, but this is also preferable in terms of analysis accuracy and VQ ratio calculation accuracy. do not have. Therefore, in step S181, the aspect ratio of each small area is checked.

The aspect ratio of the small area may be checked for each frame image acquired in step S11, but since the small area is deformed by breathing, it may be checked for frame images captured in a breathing state. Specifically, the frame image of the maximum inspiratory position (the resting inspiratory position in the case of quiet breathing) where the change in the shape of the lung field region from the reference frame image is the largest (for example, the maximum lung field area or the position of the diaphragm is the largest). It suffices to check in the lower frame image). Therefore, in the present embodiment, the frame image at the maximum inspiratory position (the resting inspiratory position in the case of quiet breathing) is used as a check image for checking.

Next, it is determined whether or not the aspect ratio of the small area is within the reference range (step S182).
That is, it is determined whether or not the aspect ratio of the selected small area in the check image is within the reference range.

If it is determined that the aspect ratio of the small area is within the reference range (step S182; YES), the process proceeds to step S185.

If it is determined that the aspect ratio of the small area does not fall within the reference range (step S182; NO), the small area is divided or integrated (step S183).
For example, if the size of the small area of the check image is smaller than a predetermined threshold, it is merged with adjacent small areas. For example, in the example shown in FIG. 9, since the small areas with ID=2 and 4 are smaller than the predetermined threshold, they are integrated into the adjacent small areas with ID=1 and 3, respectively. If the size of the small area of the check image is greater than or equal to a predetermined threshold, the small area is divided into a plurality of small areas. After integration or division, IDs are provisionally reassigned according to predetermined rules. For example, in the case of consolidation, the smaller ID number is assigned. In the case of division, a branch number is assigned. For example, when a small area with ID=1 is divided into two, IDs 1-1, 1-2, etc. are assigned.

When the subregions are divided or integrated, the subregions having the same ID in other frame images are similarly integrated or divided, and the IDs are provisionally reassigned (step S184). goes to step S185. Through step S184, adjustment is made so that the number of small areas of all acquired frame images is the same.

In step S185, it is determined whether or not the aspect ratios have been checked for all small areas (step S185).
If it is determined that the aspect ratios of all the small areas have not been checked yet (step S185; NO), the process returns to step S181, and the small areas that have not been checked are selected. is repeatedly executed.
If it is determined that all small areas have been checked (step S185; YES), the process proceeds to step S19 in FIG.

In this embodiment, the case where the aspect ratio is checked for all small areas is described as an example, but the aspect ratio may be checked for only some of the small areas. For example, since expansion and contraction due to respiration is greater in the lower lung field than in the upper lung field, the aspect ratio may be checked only for a small area of the lower ⅓ of the lung field area. Thereby, the processing speed can be increased.

In step S19 of FIG. 3, in each acquired frame image, the small areas are transformed so that the aspect ratios of the set small areas fall within the reference range (step S19).
For example, in each acquired frame image, it is determined whether or not the aspect ratio of each of the set small regions falls within the reference range. The vertices are adjusted and the small area is deformed so that it fits inside.

Then, reassignment of IDs to each small area is performed (step S20). For example, an ID is assigned to each small area of the reference frame image according to the same rule as in step S16, and the same ID is assigned to corresponding small areas of other frame images.

Next, ventilation analysis is performed using a plurality of frame images for ventilation analysis, and blood flow analysis is performed using a plurality of frame images for blood flow analysis (step S21).

In the ventilation analysis, an index value indicating the ventilation volume is calculated for each pixel of each frame image for ventilation analysis. For example, the signal value of each pixel in each frame image for ventilation analysis and the signal value of the corresponding pixel in the analysis reference frame image used as a reference in the analysis (absolute value of the difference value; ratio may also be used). , is calculated as an index value indicating the ventilation volume for each pixel of each frame image. The corresponding pixels between the frame images referred to here are, for example, pixels arranged in the same position (coordinates within the small area) in the small area to which the same ID is assigned. For the analysis reference frame image, for example, the time change of the signal value is obtained for each corresponding pixel between the frame images, and the frame image with the lowest signal value (that is, the frame image when the region is contracted the most) is the pixel of that pixel. It can be an analysis reference frame image.

In blood flow analysis, an index value indicating blood flow is calculated for each pixel of each frame image for blood flow analysis. For example, the difference value (absolute value of the difference value; ratio may be used) between the signal value of each pixel of each frame image for blood flow analysis and the signal value of the corresponding pixel of the analysis reference frame image that serves as a reference in the analysis. , is calculated as an index value indicating the blood flow volume for each pixel of each frame image. For the analysis reference frame image, for example, the time change of the signal value is obtained for each corresponding pixel between the frame images, and the frame image with the highest signal value (that is, the frame image with the lowest blood flow) is used as the analysis reference for that pixel. It can be a frame image. Note that when the frame images for blood flow analysis are images captured in a respiratory state, temporal changes in the signal value for each corresponding pixel between the frame images are filtered by a high-pass filter (for example, a cutoff frequency of 0.7 Hz) in the time direction. ), and then calculate the analysis reference frame image and the difference value. As a result, low-frequency signal changes due to ventilation can be removed, and time changes in signal values due to blood flow can be extracted.

In the ventilation analysis and blood flow analysis, first, the representative value (e.g., average value) of the signal values of the pixels in each small region is calculated, and the signal value of each small region is used as the calculated representative value. The difference value between the signal value of each small area of each frame image for analysis (for blood flow analysis) and the signal value of the corresponding small area of the analysis reference frame image used as the reference in the analysis is calculated as the small area of each frame image. It may also be calculated as an index value indicating the ventilation volume (blood flow volume) for each period.

Next, a representative value of the index values indicating the ventilation volume and a representative value of the index values indicating the blood flow are calculated for each small region (step S22).
For example, in each small region, the maximum index value indicating the ventilation volume calculated in step S21 is obtained for each pixel in the small region, and the representative value (for example, the average value, median value, maximum value, minimum value, etc.) is calculated as an index value indicating the ventilation volume of each small region. Further, in each small region, the maximum value of the index values indicating the blood flow calculated in step S21 is acquired for each pixel in the small region, and the representative value (for example, the average value, median value, maximum value, minimum value, etc.) is calculated as an index value indicating the blood flow rate of each small region.
In step S21, as an index value indicating the ventilation volume (blood flow) for each small region of each frame image, the signal value of each small region of each frame image corresponds to the analysis reference frame image serving as a reference in the analysis. When the difference value from the representative value of the signal values of the small regions is calculated, in step S22, the maximum value of the difference values is calculated as the index value indicating the ventilation volume (blood flow volume) of each small region.

Next, the VQ ratio is calculated between small regions with the same ID (step S23).
Specifically, VQ representing the ratio of V and Q is obtained by dividing the representative value (V) of the index values indicating the ventilation volume between the small regions with the same ID by the representative value (Q) of the index values indicating the blood flow volume. A ratio is calculated.
Note that the VQ ratio may be calculated by multiplying each of the V value and the Q value by a coefficient. For example, the value of V may be multiplied by a coefficient α and the value of Q may be multiplied by a coefficient β to calculate V×α/Q×β as the VQ ratio. In addition, the value verified by the clinical experiment is used for the coefficient.
Here, in a small region with sufficient ventilation and sufficient blood flow, the VQ ratio is close to 1, and it can be determined that this small region is normal.
A small area with adequate ventilation but insufficient blood flow will have a very large VQ ratio, and this small area can be judged to have an abnormality such as a pulmonary embolism.
In a small region with insufficient ventilation and blood flow, the VQ ratio approaches 1, and it can be judged that this small region has a normal balance between ventilation and blood flow.
In a small area with insufficient ventilation and sufficient blood flow, the VQ ratio is close to 0, and this small area can be judged to have an abnormality such as acute atelectasis.

Next, the calculated VQ ratio for each small region is mapped onto the image and output (step S24), and the VQ ratio calculation process ends.
In step S24, for example, each small region of one frame image (for example, the reference frame image) is displayed on the display unit 34 with a color corresponding to the calculation result of the VQ ratio calculated in the small region. . This enables the user to easily grasp the VQ ratio of each small region, whether or not there is a small region with a poor VQ ratio, and if so, which region it is.
In addition, the lung region of one frame image (for example, reference frame image) is divided into a plurality of ROIs (regions of interest), and for each ROI, the representative value (average value, etc.) of the VQ ratio of the small region in the ROI is calculated. It is also possible to calculate and display the calculation result (numerical value). Alternatively, as shown in FIG. 10(a), for each ROI, the sum (total) of the VQ ratios of the small regions in the ROI is calculated, and the total VQ ratio for each ROI is The ratio may be calculated and the calculation result (numerical value) may be displayed in association with the ROI on the image. This allows the user to easily grasp the VQ ratio for each ROI.
Further, as shown in FIG. 10B, when an ROI is designated by the operation unit 33 on a frame image (for example, a reference frame image) displayed on the display unit 34, a small area within the designated ROI is displayed. A representative value (average value, etc.) of the VQ ratio may be calculated and the calculation result (numerical value) may be displayed. This allows the user to easily check the VQ ratio of any ROI.

Note that the control unit 31, for example, sets each of the ventilation analysis result, the blood flow analysis result, and the VQ ratio to a predetermined threshold for determining an abnormal value (each of the ventilation analysis result, the blood flow analysis result, and the VQ ratio If the result of the comparison is an abnormal value, the size (the shape is constant) or the shape of the small region when dividing the lung region of the reference frame image into small regions (the size is constant) may be automatically adjusted. The shape is adjusted, for example, by adjusting the vertices of the small region so that the aspect ratio changes while the area of the small region remains constant.
For example, when a structure such as a bone is included in a small area, the VQ ratio becomes abnormally large, resulting in numerical values influenced by artifacts. In order to avoid this, it is conceivable to make the small area large enough to eliminate the influence of bones and the like. Therefore, the control unit 31 adjusts the size or shape of the small region when an abnormal value is found in the ventilation analysis result, the blood flow analysis result, or the VQ ratio. The size or shape of all small regions may be adjusted collectively, or small regions with abnormal values may be individually adjusted. The adjustment is made little by little, and after the adjustment, the process returns to step S15 in FIG. value).

Further, the control unit 31 may adjust the size or shape of the small area in response to an instruction to adjust the size or shape of the small area by the operation of the operation unit 33 by the user. As a result, for example, when the user wants to find a small disease, by setting the size of the small region small, the size of the small region is adjusted to be small, and the small disease can be found. When it is desired to reduce artifacts, it is possible to reduce the influence of artifacts by setting the size of the small region to be large. The setting (adjustment) of the size of the small area by the user may be performed instead of setting the number of divisions in step S14, or the adjustment may be received after the VQ ratio calculation process is completed, and the reference frame size may be adjusted based on the adjusted size. The processing of steps S16 to S24 may be performed after segmenting the lung region of the image. Further, the adjustment of the shape of the small region by the user is performed by displaying a plurality of divided small regions on the display unit 34 after the lung region of the reference frame image is divided in step S3, and adjusting the shape of the small region by the user. may be accepted, or the adjustment may be accepted after the VQ ratio calculation process is completed, and the processes of steps S17 to S24 may be performed based on the adjusted shape.

As described above, according to the diagnostic console 3, the control unit 31 selects a plurality of frame images for ventilation analysis and blood flow from a plurality of frame images obtained by radiographic imaging of the chest dynamics of the subject. A plurality of frame images for analysis are acquired, a reference frame image is selected from among the acquired frame images, the lung region in the reference frame image is divided into a plurality of small regions, and a plurality of For each small region, by tracking the small region in a frame image other than the reference frame image and specifying the position corresponding to the small region in the other frame image, set a small region of Then, for each of the plurality of set small regions, an index value indicating the ventilation volume is calculated using the plurality of frame images for ventilation analysis, and the blood flow volume is indicated using the plurality of frame images for blood flow analysis. An index value is calculated, and a ratio (VQ ratio) between the calculated index value indicating the ventilation volume and the index value indicating the blood flow is calculated for each of the plurality of small regions.
Therefore, by tracing the small region for VQ ratio calculation of the reference frame image in other frame images for ventilation analysis and blood flow analysis, corresponding small regions are set between frame images to calculate the VQ ratio. , it is possible to efficiently calculate the VQ ratio without performing warping processing.

Further, for example, the control unit 31 acquires a plurality of frame images obtained by imaging in a breathing state as a plurality of frame images for ventilation analysis, and obtains a plurality of frame images obtained by imaging in a breath-holding state. Images are acquired as multiple frame images for blood flow analysis. Therefore, it is possible to accurately perform ventilation analysis and blood flow analysis to calculate the VQ ratio.

Further, for example, the control unit 31 sets the number of small regions to be the same for the plurality of frame images for ventilation analysis and the plurality of frame images for blood flow analysis. It is often possible to calculate the VQ ratio.

Further, for example, the control unit 31 can set small regions with sizes corresponding to the characteristics of the lung regions by dividing the reference frame image so that the sizes of the small regions differ according to the positions in the lung regions. can be done.
For example, by dividing the reference frame image so that the size of the small region in the obliterated portion is larger than the size of the small region in the center of the lung region, the processing time for small regions where accuracy is not guaranteed can be shortened. , and the VQ ratio can be calculated efficiently.
Also, for example, by dividing the reference frame image so that the size of the lower small region is smaller than the size of the upper small region of the lung field region, it is possible to perform division according to the movement of the lung field. becomes.

Also, for example, the control unit 31 determines whether or not the aspect ratio of each of the plurality of small regions falls within the reference range, and divides the small regions determined to be out of range or integrates them into adjacent small regions. This makes it possible to accurately calculate the VQ ratio while suppressing the influence of artifacts.
Further, for example, the control unit 31 determines whether or not the aspect ratio of each of the plurality of small regions falls within the reference range using the frame image captured in the respiratory state among the plurality of acquired frame images. By doing so, it is possible to efficiently determine whether or not the aspect ratio of each of the plurality of small regions falls within the reference range.

Further, for example, after dividing or integrating the small regions, the control unit 31 adjusts the number of the plurality of small regions in all the acquired frame images to be the same. It is possible to prevent small areas from becoming incompatible between frame images.

Further, for example, the control unit 31 determines whether or not the aspect ratio of each of the plurality of small regions falls within the reference range, and adjusts the small regions determined to not fit within the reference range so that their aspect ratios fall within the reference range. By deforming, it is possible to suppress the influence of artifacts and calculate the VQ ratio with high accuracy.

Further, the control unit 31 determines that the calculated index value indicating the ventilation volume, the index value indicating the blood flow, or the ratio between the index value indicating the ventilation volume and the index value indicating the blood flow becomes an abnormal value. When it is determined that the lung region of the reference frame image is divided into a plurality of small regions, the VQ ratio can be calculated with high accuracy by adjusting the size or shape of the small regions.

Further, the control unit 31 adjusts the size or shape of the small regions when dividing the lung field region of the reference frame image into a plurality of small regions according to the user's operation. can be adjusted according to the user's needs, such as reducing the size of the small area and increasing the size of the small area to suppress the influence of artifacts.

Further, by allowing the user to select the number of divisions when dividing the reference frame image into the plurality of small areas from a plurality of predetermined numbers of divisions, the user can easily set the number of divisions. Become.

In addition, the control unit 31 sets a marker in each of the plurality of small regions of the reference frame image, and tracks the set marker in the other frame image other than the reference frame image, so that other frames can be displayed. By specifying the position corresponding to the small area in the image, it is possible to associate the position of the small area between the frame images with high accuracy.

Note that the description in this embodiment is an example of a suitable dynamic analysis device and program according to the present invention, and is not limited to this.

For example, in the above embodiment, the case where the VQ ratio is calculated from a plurality of frame images acquired by dynamic photography has been described as an example. It may be an image. For example, as a plurality of frame images for ventilation analysis, two still images obtained by photographing the maximum expiratory position and the maximum inspiratory position may be used, and as frame images for blood flow analysis, diastole and systole may be used. may be used.

Further, in the above embodiment, if the aspect ratio does not fit within the reference range and it is necessary to divide or integrate the small areas, the division or integration is performed in steps S183 and S184 of FIG. Although it was decided to reassign the ID, the information of the ID of the small region to be integrated, the ID of the small region to be divided and the number of divisions are stored in RAM or the like, and the ventilation analysis, blood flow analysis, and calculation of the VQ ratio Sometimes, a small region with an ID specified to be integrated is treated as one small region for calculation, or a small region with an ID specified to be divided is divided based on the number of divisions into one small region. Calculation may be performed by treating it as a region.

In the above embodiment, the output means for outputting the calculation result of the VQ ratio is the display unit, and the display unit outputs the calculation result of the VQ ratio. It may be output to an external device by the unit 35, or may be printed out by a printing device or the like.

Further, for example, in the above description, an example using a hard disk, a semiconductor non-volatile memory, or the like is disclosed as a computer-readable medium for the program according to the present invention, but the present invention is not limited to this example. As another computer-readable medium, 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.

In addition, the detailed configuration and detailed operation of each device that constitutes dynamic analysis system 100 can be changed as appropriate without departing from the scope 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 Diagnosis console 31 Control Unit 32 Storage unit 33 Operation unit 34 Display unit 35 Communication unit 36 Bus

Claims (17)

Acquisition means for acquiring a plurality of frame images for ventilation analysis and a plurality of frame images for blood flow analysis from a plurality of frame images obtained by radiography of dynamics of the subject's chest;
a division means for selecting a reference frame image from among the plurality of frame images for ventilation analysis and blood flow analysis acquired by the acquisition means, and dividing a lung region in the reference frame image into a plurality of small regions;
tracking each of the plurality of small regions of the reference frame image in a frame image other than the reference frame image to specify a position corresponding to the small region in the other frame image; tracking means for setting a plurality of corresponding small areas between the frame images acquired by the acquisition means;
For each of the plurality of set small regions, an index value indicating a ventilation volume is calculated using the plurality of frame images for ventilation analysis, and a blood flow volume is calculated using the plurality of frame images for blood flow analysis. an analysis means for calculating an index value indicating
a ratio calculation means for calculating a ratio between an index value indicating the ventilation volume and an index value indicating the blood flow for each of the plurality of small regions calculated by the analysis means;
A dynamic analysis device comprising a The acquiring means acquires a plurality of frame images obtained by imaging in a breathing state as a plurality of frame images for ventilation analysis, and acquires a plurality of frame images obtained by imaging in a breath-holding state as the plurality of frame images for ventilation analysis. The dynamics analysis device according to claim 1, which is acquired as a plurality of frame images for blood flow analysis. 3. The dynamic analysis apparatus according to claim 1, wherein the plurality of frame images for ventilation analysis and the plurality of frame images for blood flow analysis have the same number of small regions. 4. The reference frame image according to any one of claims 1 to 3, wherein the dividing means divides the reference frame image so that the size of each of the plurality of small regions differs according to the position of each of the plurality of small regions in the lung region. kinetic analysis equipment. 5. The dynamic analysis apparatus according to claim 4, wherein the dividing means divides the reference frame image so that the size of the small area of the peripheral part is larger than the size of the small area of the central part of the lung field area. 5. The dynamic analysis apparatus according to claim 4, wherein the dividing means divides the reference frame image so that the size of the lower small region is smaller than the size of the upper small region of the lung field region. Division for determining whether or not the aspect ratio of each of the plurality of small areas set by the tracking means falls within a reference range, and dividing the small area determined not to fit or integrating it into adjacent small areas The dynamic analysis device according to any one of claims 1 to 6, comprising integration means. The dividing and integrating means uses a frame image captured in a respiratory state among the plurality of frame images acquired by the acquiring means to set the aspect ratio of each of the plurality of small regions set by the tracking means to the 8. The dynamic analysis device according to claim 7, which determines whether or not it falls within the reference range. 9. The dynamic state according to claim 7 or 8, wherein said dividing and integrating means adjusts the numbers of said plurality of small regions of all frame images acquired by said acquiring means to be the same after said dividing or said integrating. analysis equipment. determining whether or not the aspect ratio of each of the plurality of small areas set by the tracking means falls within a reference range; The dynamic analysis device according to any one of claims 1 to 9, comprising deformation means for deformation. 11. The method according to any one of claims 1 to 10, further comprising adjusting means for adjusting the size or shape of the small regions when the lung field region of the reference frame image is divided into the plurality of small regions by the dividing means. Dynamic analysis device. The adjusting means determines that the calculated index value indicating the ventilation volume, the index value indicating the blood flow, or the ratio between the index value indicating the ventilation volume and the index value indicating the blood flow is an abnormal value. 12. The dynamic analysis device according to claim 11, wherein the adjustment is performed when the 12. The dynamic analysis apparatus according to claim 11, wherein said adjustment means performs said adjustment according to user's operation. 13. The division number according to any one of claims 1 to 12, wherein the number of divisions when the reference frame image is divided into the plurality of small areas by the division means can be selected from among a plurality of predetermined division numbers. kinetic analysis equipment. The tracking means sets a marker in each of the plurality of small regions of the reference frame image, and tracks the set marker in a frame image other than the reference frame image, thereby 15. The dynamics analysis device according to any one of claims 1 to 14, which specifies a position corresponding to the small region in the frame image of . The dynamic analysis apparatus according to any one of claims 1 to 15, further comprising an output means for outputting a result of calculation by said ratio calculation means. the computer,
Acquisition means for acquiring a plurality of frame images for ventilation analysis and a plurality of frame images for blood flow analysis from a plurality of frame images obtained by radiography of dynamics of the subject's chest;
dividing means for selecting a reference frame image from among the plurality of frame images for ventilation analysis and blood flow analysis acquired by the acquisition means, and dividing the lung region in the reference frame image into a plurality of small regions;
tracking each of the plurality of small regions of the reference frame image in a frame image other than the reference frame image to specify a position corresponding to the small region in the other frame image; Tracking means for setting a plurality of corresponding small areas between the frame images acquired by the acquisition means;
For each of the plurality of set small regions, an index value indicating a ventilation volume is calculated using the plurality of frame images for ventilation analysis, and a blood flow volume is calculated using the plurality of frame images for blood flow analysis. analysis means for calculating an index value indicating
ratio calculation means for calculating a ratio between an index value indicating the ventilation volume and an index value indicating the blood flow for each of the plurality of small regions calculated by the analysis means;
A program to function as

Images