JP6852545B2 - Image display system and image processing equipment - Google Patents

Description

The present invention relates to an image display system and an image processing device.

Conventionally, the dynamic information (for example, the pixel difference value between the maximum inspiratory frame and the maximum expiratory frame) obtained by analyzing the dynamic image obtained by radiographing the dynamics of the diagnostic site is used as a radiographic apparatus. Has proposed a technique of superimposing and displaying an image taken by another different modality (for example, a coronal image taken by a CT device) (see, for example, Patent Document 1).

Japanese Patent No. 4797173

However, in the technique described in Patent Document 1, it is not possible to compare a functional information image (for example, lung scintigram) of a diagnostic site acquired by another modality with a dynamic image and perform comparative interpretation while matching the positions. Have difficulty. Further, it is more difficult to compare the mismatch of the analysis result image obtained by analyzing the functional information image of the diagnosis target site acquired by another modality and the dynamic image.

An object of the present invention is to make it possible to easily compare an analysis result image obtained from a dynamic image or a dynamic image with a functional information image acquired by another modality.

In order to solve the above problems, the image display system of the invention according to claim 1 is used.
A first lung field region identification means for identifying a lung field region from each of a plurality of frame images of a dynamic image acquired by dynamically photographing the chest of a subject with a radiography apparatus,
A second lung field region specifying means for identifying the lung field region from the function information image of the lung of the subject acquired by another modality different from the radiography apparatus, and
Alignment of the lung field region of the functional information image specified by the second lung field region specifying means and the lung field region of each of the plurality of frame images specified by the first lung field region specifying means. Deformation means to be performed and
An integrated means for integrating each of the plurality of frame images aligned by the deforming means and the functional information image, and
A display means for displaying an image integrated by the integration means, and a display means.
To be equipped.

The image display system of the invention according to claim 2 is
A first lung field region identification means for identifying a lung field region from each of a plurality of frame images of a dynamic image acquired by dynamically photographing the chest of a subject with a radiography apparatus,
An analysis means for generating an analysis result image showing the amount of change in signal value or the amount of change in morphology of a predetermined structure between the plurality of frame images, and an analysis means.
A second lung field region specifying means for identifying the lung field region from the function information image of the lung of the subject acquired by another modality different from the radiography apparatus, and
The lung of the analysis result image specified based on the lung field region of the functional information image specified by the second lung field region specifying means and the lung field region specified by the first lung field region specifying means. Deformation means for aligning with the field area,
An integrated means for integrating the analysis result image and the functional information image aligned by the deforming means, and
A display means for displaying an image integrated by the integration means, and a display means.
To be equipped.

The invention according to claim 3 is the invention according to claim 2.
A detection means for detecting an abnormal part based on the analysis result image and the function information image is provided.
The display means displays the abnormal portion detected by the detection means at a corresponding position of the image integrated by the integration means, the dynamic image, or the function information image.

The invention according to claim 4 is the invention according to any one of claims 1 to 3.
The functional information image is a pulmonary ventilation scintigram or a pulmonary blood flow scintigram.

The image processing apparatus of the invention according to claim 5 is
A first lung field region identification means for identifying a lung field region from each of a plurality of frame images of a dynamic image acquired by dynamically photographing the chest of a subject with a radiography apparatus,
A second lung field region specifying means for identifying the lung field region from the function information image of the lung of the subject acquired by another modality different from the radiography apparatus, and
Alignment of the lung field region of the functional information image specified by the second lung field region specifying means and the lung field region of each of the plurality of frame images specified by the first lung field region specifying means. Deformation means to be performed and
An integrated means for integrating each of the plurality of frame images aligned by the deforming means and the functional information image, and
To be equipped.

According to the present invention, it is possible to easily compare the analysis result image obtained from the dynamic image or the dynamic image with the functional information image acquired by another modality.

It is a figure which shows the whole structure of the image display system which concerns on this embodiment. It is a block diagram which shows the functional structure of the console of FIG. It is a flowchart which shows the image integration process A executed by the control part of FIG. 2 in 1st Embodiment. It is a figure for demonstrating the alignment of the frame image of a dynamic image, and a function information image. FIG. 5 is a flowchart showing an image integration process B executed by the control unit of FIG. 2 in the second embodiment. It is a figure which shows typically the calculation method of the movement amount at the representative point of a rib. It is a figure which shows an example of the analysis result image which shows the morphological change amount of the rib in each part of the right lung field. It is a figure which shows typically the generation of the integrated image in 2nd Embodiment.

Hereinafter, preferred embodiments according to the present invention will be described with reference to the drawings. The present invention is not limited to the illustrated examples.

<First Embodiment>
(Configuration of image display system 100)
First, the configuration of the first embodiment according to the present invention will be described.
FIG. 1 shows an example of the overall configuration of the image display system 100 according to the first embodiment.
As shown in FIG. 1, the image display system 100 includes a radiography apparatus 1, a SPECT (Single Photon Emission Computed Tomography) apparatus 2, and a console 3. The console 3 can be communicatively connected to the radiography apparatus 1 and the SPECT apparatus via a communication network such as a LAN (Local Area Network).

(Radiation imaging device 1)
The radiography apparatus 1 is composed of a radiation source, an exposure switch, a radiation detection unit (FPD: Flat Panel Detector), a communication unit, and the like, and performs dynamic imaging on the chest of the subject to acquire a dynamic image of the chest. It is a device that transmits to the console 3.
Dynamic imaging means that the subject is repeatedly irradiated with radiation such as X-rays in the form of pulses at predetermined time intervals (pulse irradiation), or the subject is continuously irradiated at a low dose rate without interruption (continuous irradiation). By doing so, it means to acquire a plurality of images. In the dynamic imaging, for example, the dynamics of the chest having a periodicity (cycle) such as morphological changes of lung expansion and contraction associated with respiratory movements and heartbeats are photographed. A series of images obtained by this continuous shooting is called a dynamic image. Further, each of the plurality of images constituting the dynamic image is called a frame image.

(SPECT device 2)
The SPECT device 2 is a device that images the radiation emitted to the outside of the body when a radioisotope (RI = radioisotope) is inhaled or intravenously injected into the body and transmits it to the console 3. In the present embodiment, the SPECT device 2 externally when the subject inhales xenon gas (Xe), krypton gas (Kr), technetium gas (TC), or the like as a functional information image showing the functional state of lung ventilation. Technetium / macroaggregated albumin (Tc-MAA) from the veins of the arm is used as a pulmonary ventilation scintigram showing the distribution of gas in the lung field, which is an image of the emitted radiation, and as a functional information image showing the functional state of pulmonary blood flow. Obtain a pulmonary blood flow scintigram showing the distribution of blood flow in the lung field, which is an image of the radiation emitted to the outside when injected. Each pixel of the pulmonary ventilation scintigram and the pulmonary blood flow scintigram indicates an RI value.

(Console 3)
The console 3 controls the dynamic image acquisition operation by the radiography apparatus 1 and the image acquisition operation by the SPECT apparatus 2. Further, the console 3 is an image processing device that integrates a dynamic image or an analysis result image thereof and a functional information image acquired by the SPECT device 2 to generate and display an integrated image.

FIG. 2 shows an example of the functional configuration of the console 3. As shown in FIG. 2, the console 3 includes a control unit 31, a storage unit 32, an operation unit 33, a display unit 34, a communication unit 35, and the like, and each unit is connected by a bus 36.

The control unit 31 includes a CPU (Central Processing Unit) and a RAM (Random Access Memory).
) Etc. 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, expands them in the RAM, and operates each unit of the console 3 according to the expanded program. In addition, the operations of the radiography apparatus 1 and the SPECT apparatus 2 are centrally controlled. Further, the control unit 31 executes various processes including the image integration process A, which will be described later, according to the developed program, and performs the first lung field region specifying means, the second lung field region specifying means, the deforming means, and the integration. Act as a 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 required for executing processing by various programs and programs executed by the control unit 31. For example, the storage unit 32 stores a program for executing the image integration process A shown in FIG. 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.

Further, the storage unit 32 stores the dynamic image transmitted from the radiography apparatus 1 and the functional information image transmitted from the SPECT apparatus 2 in association with the patient information, the examination information, and the like of the subject. The patient information includes the patient ID, the patient's name, age, gender, and the like. The test information includes a test ID, a test date and time, a test site (here, the chest) or a type of test target (for example, pulmonary ventilation, pulmonary blood flow, etc.).

The operation unit 33 includes 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 (Liquid Crystal Display) or a CRT (Cathode Ray Tube), and displays input instructions, data, and the like from the operation unit 33 according to instructions of display signals input from the control unit 31. To do.

The communication unit 35 includes a LAN adapter and the like, and controls data transmission / reception with external devices such as a radiography apparatus 1 and a SPECT apparatus 2 connected to a communication network N such as a LAN.

(Operation of image display system 100)
Next, the operation in the image display system 100 will be described.
FIG. 3 shows the flow of the image integration process A executed by the console 3 in response to an instruction from the operation unit 33. The image integration process A is executed in collaboration with the program stored in the control unit 31 and the storage unit 32. Hereinafter, the image integration process A will be described with reference to FIG.

First, the control unit 31 acquires a dynamic image of the processing target (step S1). For example, the control unit 31 displays an image search screen or the like on the display unit 34, and from the dynamic images stored in the storage unit 32, the control unit 31 displays a dynamic image that matches the search condition input by the operation unit 33 from the storage unit 32. Read and get.

Next, the control unit 31 identifies the lung field region from each frame image of the acquired dynamic image (step S2).
The method for specifying the lung field is not particularly limited, and a known method can be applied. For example, as disclosed in Japanese Patent No. 2987633, in the X-ray image, the image density of the left and right lung portions in the lung field region is higher than that in the periphery. Therefore, a density histogram of an arbitrary frame image may be created, an image portion of a high density region corresponding to the lung field region may be determined from the shape and area of the density histogram, and the image portion may be specified as the lung field. Alternatively, as disclosed in Japanese Patent Application Laid-Open No. 2003-6661, the lung field region is specified by performing template matching on an arbitrary frame image using a standard template that defines the contour of the lung field region. You can also do it.

Next, the control unit 31 reads out the functional information image (pulmonary ventilation scintigram and / or pulmonary blood flow scintigram) corresponding to the acquired dynamic image from the storage unit 32 and acquires it (step S3).
For example, the control unit 31 reads out from the storage unit 32 a functional information image relating to the same patient as the dynamic image and having the same shooting date (the latest shooting date) as a functional information image corresponding to the dynamic image.

Next, the control unit 31 identifies the lung field region from the functional information image (step S4).
For example, the control unit 31 binarizes the functional information image with a predetermined threshold value, and specifies a region having an RI value equal to or higher than the threshold value as a lung field region.

Next, the control unit 31 deforms the lung field region of the functional information image and aligns it with the lung field region of each frame image of the dynamic image (step S5).
For example, as shown in FIG. 4, the height Y1 and width X1 of the lung field region of each frame image of the dynamic image and the height Y2 and width X2 of the lung field region of the functional information image are obtained, and the height Y2 and width X2 of each frame image are obtained. The function information image is deformed so that Y1 and Y2 and X1 and X2 are equal.

Next, the control unit 31 integrates (combines) the function information image with each frame image of the dynamic image to generate an integrated image (step S6). For example, alpha blending is performed for each corresponding pixel of each frame image and the function information image to generate an integrated image.
Then, the control unit 31 displays the integrated image on the display unit 34 (step S7), and ends the image integration process A. The integrated image may be displayed as a moving image, the frame images may be displayed side by side, or the frame images may be sequentially switched and displayed according to the operation of the operation unit 33.

In the first embodiment, a lung function information image (pulmonary ventilation scintigram or pulmonary blood flow scintigram) is integrated and displayed in each frame image of the dynamic image. Therefore, the user can easily compare the dynamic image with the functional information image acquired by another modality (here, SPECT device 2) different from the radiography apparatus.

<Second embodiment>
Next, a second embodiment of the present invention will be described.
In the second embodiment, an example in which the analysis result image generated based on the dynamic image and the functional information image are integrated and displayed will be described.

The configuration in the second embodiment is the same as that described in the first embodiment except that the program for executing the image integration process B is stored in the storage unit 32 of the console 3. The image integration process B executed on the console 3 will be described below with reference to the image integration process B.

FIG. 5 is a flowchart showing an image integration process B executed by the console 3 in response to an instruction from the operation unit 33 in the second embodiment. The image integration process B is executed in collaboration with the program stored in the control unit 31 and the storage unit 32.

First, the control unit 31 acquires a dynamic image of the processing target (step S11). Since the process of step S11 is the same as the process of step S1 of FIG. 3, the description is incorporated.

Next, the control unit 31 identifies the lung field region from each frame image of the acquired dynamic image (step S12). Since the process of step S12 is the same as the process of step S2 of FIG. 3, the description is incorporated.

Next, the control unit 31 analyzes the dynamic image and generates an analysis result image (step S13).

For example, the control unit 31 calculates the amount of change in the signal value for each pixel or each block of a plurality of pixels in the lung field region of the dynamic image, and generates an analysis result image showing the amount of change in the signal value for each pixel or each block. ..
For example, when the type of inspection target is lung ventilation, first, the control unit 31 applies a time change of the pixel signal value (density value) for each pixel in the lung field region of the dynamic image to a low-pass filter in the time direction (for example,). Filter at a cutoff frequency of 0.85 Hz). Alternatively, the lung field region of each frame image of the dynamic image is divided into blocks of a plurality of pixels, and representative values (for example, mean value, median value, maximum value, etc.) of pixel signal values are calculated for each block of each frame image. Then, the time change of the calculated representative value of each block is filtered by the low-pass filter in the time direction. Thereby, the high frequency component due to the pulmonary blood flow and the like can be removed to obtain the time change of the signal component of the pulmonary ventilation. It should be noted that the above processing may be performed after the lung field region is aligned between the frame images by a warping process or the like. Then, for each pixel or block of each frame image, the difference value of the pixel signal value from the corresponding pixel or block of the reference frame image (for example, the frame image of the maximum expiratory position or the frame image of the maximum intake position) is calculated. , An analysis result image is generated in which the pixel signal value of each pixel indicates the difference value (that is, the amount of signal change) of the pixel signal value from the reference frame image.
For example, when the type of inspection target is pulmonary blood flow, first, the control unit 31 applies a time change of the pixel signal value (density value) for each pixel in the lung field region of the dynamic image to a high-pass filter in the time direction (for example,). , Cutoff frequency 0.80 Hz). Alternatively, the lung field region of each frame image of the dynamic image is divided into blocks of a plurality of pixels, and representative values (for example, mean value, median value, maximum value, etc.) of pixel signal values are calculated for each block of each frame image. Then, the time change of the calculated representative value of each block is filtered by the high-pass filter in the time direction. As a result, it is possible to remove the low frequency component due to pulmonary ventilation or the like and obtain the time change of the signal component of the pulmonary blood flow. It should be noted that the above processing may be performed after the lung field region is aligned between the frame images by a warping process or the like. Then, for each pixel or block of each frame image, the difference value of the pixel signal value from the corresponding pixel or block of the reference frame image (for example, the frame image of the maximum expiratory position or the frame image of the maximum intake position) is calculated. , An analysis result image is generated in which the pixel signal value of each pixel indicates the difference value (that is, the amount of signal change) of the pixel signal value from the reference frame image.

Further, for example, the control unit 31 calculates the amount of morphological change of a predetermined structure in the lung field for each block of a plurality of pixels in the lung field region of the dynamic image, and obtains an analysis result image showing the amount of morphological change for each block. Generate.
For example, when the type of examination target is lung ventilation, first, as shown in FIG. 6A, the control unit 31 extracts the contour of the ribs and represents the representative points of each rib (for example, the mediastinal side, the center of the lung field, Calculate the moving distance between each frame image on the outer thoracic side). The method for extracting the contour of the ribs is not particularly limited, and a known method can be applied. For example, as disclosed in Japanese Patent Application Laid-Open No. 5-176919, a number of contour lines are defined in the vertical direction (head-foot direction) in the lung field region of an arbitrary frame image, and the contour lines are defined in advance. The contour of the rib region is estimated by applying the defined model function. Then, a plurality of image regions to be processed are defined in the estimated contour portion, and the gradient and the orientation corresponding to the gradient are obtained for each pixel of each image region by the Sobel operator. The largest gradient and its orientation in each pixel are defined as the gradient and orientation of the image area. When the gradient and orientation of each image area are plotted in the coordinate space as coordinates, a set is formed for each area such as the upper edge, lower edge of the rib, the inside of the rib edge, and other than the rib in the coordinate space. The image region classified into the peripheral edges such as the upper edge and the lower edge of the ribs can be extracted as the contour of the rib region.
Next, as shown in the right lung field of FIG. 6B, the control unit 31 divides the lung field into upper, middle, and lower parts, and further divides each into the mediastinal side, the center of the lung field, and the outer thoracic side, and is divided. The average amount of movement of the representative points of the ribs included in each region is calculated as the amount of morphological change in that region, and the pixel signal value of each pixel in the lung field indicates the amount of morphological change in the region to which the pixel belongs. Generate. The analysis result image is an image colored according to the pixel signal value.

Next, the control unit 31 reads out the functional information image (pulmonary ventilation scintigram and / or pulmonary blood flow scintigram) corresponding to the acquired dynamic image from the storage unit 32 and acquires it (step S14).
For example, the control unit 31 reads out from the storage unit 32 a functional information image relating to the same patient as the dynamic image and having the same shooting date (the latest shooting date) as a functional information image corresponding to the dynamic image.

Next, the control unit 31 identifies the lung field region from the functional information image (step S15).
For example, the control unit 31 binarizes the functional information image with a predetermined threshold value, and specifies a region having an RI value equal to or higher than the threshold value as a lung field region.

Next, the control unit 31 deforms the lung field region of the functional information image and aligns it with the lung field region of each frame image of the analysis result image (step S16).
For example, as shown in FIG. 4, the height Y1 and width X1 of the lung field region of each frame image of the analysis result image and the height Y2 and width X2 of the lung field region of the functional information image are obtained, and each frame image is obtained. The function information image is deformed so that Y1 and Y2 and X1 and X2 are equal to each other. The position (coordinates) of the lung field region of each frame image of the analysis result image can be specified by the lung field region identified from each frame image of the dynamic image in step S12. In addition, the pulmonary ventilation scintigram is aligned with the pulmonary ventilation analysis result image, and the pulmonary blood flow scintigram is aligned with the pulmonary blood flow analysis result image.

Next, the control unit 31 integrates (synthesizes) the function information image with each frame image of the analysis result image (step S17). For example, alpha blending is performed for each frame image of the analysis result image and each corresponding pixel of the function information image to generate an integrated image. In the present embodiment, the pulmonary ventilation scintigram is aligned with the pulmonary ventilation analysis result image, and the pulmonary blood flow scintigram is aligned with the pulmonary blood flow analysis result image.
Then, the integrated image is displayed on the display unit 34 (step S18), and the image integration process B is terminated. The integrated image may be displayed as a moving image, the frame images may be displayed side by side, or the frame images may be sequentially switched and displayed according to the operation of the operation unit 33.

G1 of FIG. 7 is an image surrounding the lung field region in an arbitrary frame image of the dynamic image acquired by the radiographing apparatus 1. G2 in FIG. 7 is an analysis result image of the frame image of G1 (here, the amount of change in the signal value of lung ventilation). G3 of FIG. 7 is a lung ventilation scintigram taken by SPECT apparatus 2. In addition, G2 and G3 show the outline of the lung field region for reference. FIG. 7 is a diagram showing an example of an integrated image in which G2 and G3 are integrated.
As shown in FIG. 7, in the second embodiment, the pulmonary ventilation scintigram (or pulmonary blood flow scintigram) is integrated and displayed in each frame image of the analysis result image of the dynamic image. Therefore, the user can input the analysis result image of the dynamic image and the functional information image (pulmonary ventilation scintigram, pulmonary blood flow scintigram) acquired by another modality different from the radiography apparatus (SPECT apparatus 2 in this case). It becomes possible to easily compare. For example, it is possible to easily recognize a mismatch between the analysis result image and the feature amount represented by the function information image.

In addition, for example, patients with COPD may have a longer expiratory time due to airflow limitation, which may result in a difference in the amount of rib movement. Therefore, the diagnostic performance can be improved by displaying the analysis result image of the amount of morphological change of the ribs in the dynamic image and the lung ventilation scintigram in an integrated manner.

As a detection means, the control unit 31 may detect an abnormal portion based on the analysis result image and the function information image, and display the detection result on the display unit 34.

For example, in a region where the signal value change amount of the analysis result image showing the signal value change amount related to lung ventilation is larger than a predetermined threshold value (that is, a region where the alveolar density change is large), the RI value of lung ventilation scintillography is predetermined. If there is a region smaller than the threshold value of, the region is detected as an abnormal site because there is a mismatch that the alveolar density change is large but the ventilation volume is small. Then, for example, the region of the abnormal portion of the integrated image, the dynamic image, or the functional information image is annotated and displayed on the display unit 34. This makes it possible for the user to more easily recognize the abnormal part.

Further, for example, in a region where the signal value change amount of the analysis result image showing the morphological change amount of a predetermined structure (for example, rib) is larger than the predetermined threshold value, the RI value of the lung ventilation scintigraphy is smaller than the predetermined threshold value. If there is, it is detected as an abnormal part because there is a mismatch that the rib movement is large but the ventilation volume is small in that area. Then, for example, the region of the abnormal portion of the integrated image, the dynamic image, or the functional information image is annotated and displayed on the display unit 34. This makes it possible for the user to more easily recognize the abnormal part.

Although the first and second embodiments of the present invention have been described above, the description contents in the above-described embodiments and modifications are suitable examples of the image display system according to the present invention, and are not limited thereto. Absent.

For example, in the above embodiment, the case where the SPECT device 2 is another modality different from the radiography device 1 and the functional information image is the pulmonary ventilation scintigram or the pulmonary blood flow scintigram has been described as an example. Not limited to. For example, lung ventilation images obtained by imaging with an MRI (Magnetic Resonance Imaging) device using enzyme gas as a contrast medium (for example, Tae Iwasawa "Current Status and Future Prospects of MR Ventilation Images", INNER VISION (28/10) ) Refer to 2013) may be used as a functional information image. In addition, for example, xenone inhalation double-energy CT image (for example, Kohei Aoki "Evaluation of lung morphology and pulmonary ventilation blood flow ratio before and after lung resection using double-energy CT after xenon inhalation and contrast medium administration", Saitama Medical University Magazine Vol. 42, No. 1, see August 2015) may be used as the functional information image.

Further, in the second embodiment described above, the analysis result image and the function information image are integrated and displayed, but the analysis result image and the function information image are integrated on each frame image of the dynamic image. It may be displayed as.
Further, the dynamic analysis image and the functional information image may be integrated and displayed on the schema prepared in advance.

Further, in the above description, an example in which a hard disk, a semiconductor non-volatile memory, or the like is used as a computer-readable medium for the program according to the present invention has been disclosed, 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 data of the program according to the present invention via a communication line.

In addition, the detailed configuration and detailed operation of each device constituting the image display system can be appropriately changed without departing from the spirit of the invention.

100 Image display system 1 Radiation imaging device 2 SPECT device 3 Console 31 Control unit 32 Storage unit 33 Operation unit 34 Display unit 35 Communication unit 36 Bus

Claims (5)

A first lung field region identification means for identifying a lung field region from each of a plurality of frame images of a dynamic image acquired by dynamically photographing the chest of a subject with a radiography apparatus,
A second lung field region specifying means for identifying the lung field region from the function information image of the lung of the subject acquired by another modality different from the radiography apparatus, and
Alignment of the lung field region of the functional information image specified by the second lung field region specifying means and the lung field region of each of the plurality of frame images specified by the first lung field region specifying means. Deformation means to be performed and
An integrated means for integrating each of the plurality of frame images aligned by the deforming means and the functional information image, and
A display means for displaying an image integrated by the integration means, and a display means.
Image display system with. A first lung field region identification means for identifying a lung field region from each of a plurality of frame images of a dynamic image acquired by dynamically photographing the chest of a subject with a radiography apparatus,
An analysis means for generating an analysis result image showing the amount of change in signal value or the amount of change in morphology of a predetermined structure between the plurality of frame images, and an analysis means.
A second lung field region specifying means for identifying the lung field region from the function information image of the lung of the subject acquired by another modality different from the radiography apparatus, and
The lung of the analysis result image specified based on the lung field region of the functional information image specified by the second lung field region specifying means and the lung field region specified by the first lung field region specifying means. Deformation means for aligning with the field area,
An integrated means for integrating the analysis result image and the functional information image aligned by the deforming means, and
A display means for displaying an image integrated by the integration means, and a display means.
Image display system with. A detection means for detecting an abnormal part based on the analysis result image and the function information image is provided.
The image display system according to claim 2, wherein the display means displays an abnormal portion detected by the detection means at a corresponding position of an image integrated by the integration means, the dynamic image, or the function information image. The image display system according to any one of claims 1 to 3, wherein the functional information image is pulmonary ventilation scintigraphy or pulmonary blood flow scintigraphy. A first lung field region identification means for identifying a lung field region from each of a plurality of frame images of a dynamic image acquired by dynamically photographing the chest of a subject with a radiography apparatus,
A second lung field region specifying means for identifying the lung field region from the function information image of the lung of the subject acquired by another modality different from the radiography apparatus, and
Alignment of the lung field region of the functional information image specified by the second lung field region specifying means and the lung field region of each of the plurality of frame images specified by the first lung field region specifying means. Deformation means to be performed and
An integrated means for integrating each of the plurality of frame images aligned by the deforming means and the functional information image, and
An image processing device comprising.

Images