JP5619942B2 - X-ray diagnostic apparatus and image processing apparatus - Google Patents

Description

  The present invention relates to an X-ray diagnostic apparatus for imaging a subject with pulsed X-rays and an image processing apparatus suitable for X-ray coronary angiography diagnosis.

  There is a technique of performing angiographic imaging immediately after intravascular intervention treatment, observing a temporal change of the contrast agent flowing into the myocardium, and assisting in determining the end of treatment. Since the publication (Gibson) that prognosis can be predicted with this technique in 2000, there have been many examples of evaluating medicinal efficacy and technique using this technique.

  Despite the evidence that this indicator is clinically useful, it has not yet been used daily in routine clinical practice. Causes include increased exposure and low quantitativeness.

  In an imaging method for examination of blood vessel morphology, imaging is generally performed as shown in FIG. That is, an image is obtained by irradiating X-rays in a pulse shape for 3 to 5 heartbeats for a heartbeat exercise and for about 5 seconds for an imaging time.

  On the other hand, when measuring myocardial perfusion is intended for examination, imaging over a long period of time is required. It is considered that continuous observation for at least 5 seconds, generally about 30 seconds, and at most about 60 seconds is preferable. Therefore, when imaging for the purpose of myocardial perfusion is added to imaging for the purpose of blood vessel morphology observation, an imaging method as shown in FIG. 7 is required.

  Here, as another plan, there is also an idea that imaging for the purpose of blood vessel morphology observation and imaging for the purpose of myocardial perfusion are separately performed and imaging is performed twice in total. However, in that case, the contrast medium to be administered is doubled, which is most undesirable.

  Further, in a conventional X-ray diagnostic apparatus for blood vessel observation, so-called X-ray angio apparatus, coronary angiography is performed, but quantitative measurement of myocardial blood flow using a contrast agent cannot be performed. It was necessary to measure myocardial blood flow with a nuclear medicine diagnostic device or an MRI device during the time period.

  When there are multiple stenosis in the coronary artery and there is ischemia due to any of these lesions, the myocardial blood flow cannot be measured in the catheter room, and the site of ischemia could not be easily identified.

When the presence of stenosis, embolism, or thrombus in peripheral blood vessels in coronary arteries was suspected, it was difficult to confirm or confirm recovery of myocardial blood flow after treatment with a thrombolytic agent.

Myocardial blood flow could not be measured and displayed using an X-ray diagnostic apparatus.
Relationship of TIMI Myocardial Perfusion Grade to Mortality After Administration of Thrombolytic Drugs, C. Michael Gibson, MS, MD; Christopher P. Cannon, MD; Sabina A. Murphy, MPH; Kathryn A. Ryan, BS; Rebecca Mesley, BS; Susan J. Marble, RN, MS; Carolyn H. McCabe, BS; Frans Van de Werf, MD, PhD; Eugene Braunwald, MD ;, Circulation, 101, 125-130, 2000 Arnoud WJ van 't Hof, MD; Aylee Liem, MD; Harry Suryapranata, MD; Jan CA Hoorntje, MD; Menko-Jan de Boer, MD; Felix Zijlstra, MD ;, Angiographic Assessment of Myocardial Reperfusion in Patients Treated With Primary Angioplasty for Acute Myocardial Infarction, Myocardial Blush Grade, Circulation, 97, 2302-2306, 1998)

An object of the present invention is to realize a reduction in exposure and a reduction in the amount of contrast medium in, for example, imaging having both a vascular morphology observation purpose and a myocardial perfusion purpose in which the imaging time is relatively long.
Another object of the present invention is to provide useful information regarding myocardial blood flow.

An X-ray diagnostic apparatus according to the present invention generates an image based on an X-ray source that generates pulsed X-rays, an X-ray detector that detects X-rays transmitted through a subject, and an output of the X-ray detector. An image generation unit; and a control unit that controls the X-ray source to repeatedly generate the pulse X-rays within the imaging period and to change the generation period of the pulse X-rays within the imaging period.
An image processing apparatus according to the present invention includes a storage unit that stores data of a plurality of images generated under a cardiac coronary angiography sequence, and a blood supply region to the myocardium based on the data of the plurality of images Generating a reference time density curve related to the reference area set to, and a plurality of time density curves related to a plurality of local areas set in the myocardial area, and And an arithmetic unit that calculates a plurality of indices each representing a correlation.

According to the present invention, it is possible to achieve a reduction in exposure and a reduction in the amount of contrast medium in imaging, for example, for both vascular morphology observation purposes and myocardial perfusion purposes for which the imaging time is relatively long.
Moreover, according to this invention, useful information regarding myocardial blood flow can be provided.

(First embodiment)
Hereinafter, a first embodiment of the present invention will be described with reference to the drawings.
First, terms used in the following description are defined as follows.
Pulse rate; an index that regulates the frequency of irradiating a subject with pulse X-rays. The number of pulse X-ray irradiations per unit time (1 second or 1 minute) or 1 heart rate is expressed in units (times / second, times / The frame rate representing the number of images taken per unit time is substantially equivalent to the pulse rate. The pulse X-ray irradiation is not limited to a mode in which pulse X-rays are generated from an X-ray tube and irradiated on the subject as it is, but X-rays are continuously generated and pulse X-rays are emitted by an X-ray shutter or the like. The form which produces | generates and irradiates a subject is also included.
Perfusion; refers to the flow of blood to the heart muscle.
R wave; represents the peak wave of the electrocardiogram waveform.
RR interval; represents the time interval between R and R waves.
mA: A tube current flowing between the electrodes of the X-ray tube and used as an index representing the height of the pulse X-ray.
Pulse width: represents the duration of a pulse X-ray, and its unit is msec.
mAs: tube current time product obtained by multiplying the tube current and the pulse width, and used as an index representing the intensity of X-rays.
Protocol: A set of parameters in which various conditions are stored together
Blush: Gibson et al. Semi-quantified the amount of contrast medium flowing into the myocardium.
Control Image; Reference Image FIG. 1 shows an X-ray imaging apparatus according to this embodiment. The X-ray imaging apparatus has a gantry 100. The gantry 100 has a C-arm 7. The C arm 7 is rotatably supported by the support mechanism 6. An X-ray generator 2 is attached to one end of the C arm 7. The X-ray generator 2 includes an X-ray tube 21 and an X-ray collimator 22. The high voltage generator 1 generates a high voltage (tube voltage) applied between the electrodes of the X-ray tube 22 and generates a filament current to be supplied to the filament of the X-ray tube 22. The high voltage controller 20 controls the tube voltage and / or the filament current generated by the high voltage generator 1 according to the control of the system controller 8.

  An X-ray detector 5 is attached to the other end of the C arm 7. The X-ray detection unit 5 faces the X-ray tube 22 of the X-ray generation unit 2 across the subject 3 placed on the bed 4. The X-ray detection unit 5 is typically a solid flat detector in which a plurality of detection elements (pixels) that directly or indirectly convert incident X-rays into electric charges are two-dimensionally arranged. The X-ray detector 5 typically repeats a one-cycle detection operation consisting of charge accumulation, charge readout, and reset at a constant period under the control of the system controller 8.

  The image calculation / storage unit 10 has a function of generating image data based on an output from the X-ray detection unit 5, a function of storing image data, and a function of processing image data. The system control unit 8 outputs an injection start signal output from the injector 15 when the injection of the contrast agent from the injector 15 to the subject 3 is started, and is output from the injector 15 when the injection of the contrast agent into the subject is completed. Based on the injection end signal and the electrocardiogram (ECG) of the subject 3 measured by the electrocardiograph 16, the main function is to control the imaging operation as described later. An operation unit 9 is connected to the system control unit 8. The operation unit 9 is provided with a user interface 14 having a hand switch 12, a display, a touch panel, and the like.

  FIG. 2 shows an imaging operation by the system control unit 8. Prior to imaging, a catheter is inserted into the subject and advanced to, for example, the heart coronary artery entrance, and a contrast examination is started at that position.

  With the start of inspection, generation of pulse X-rays is started. The imaging period from the start of generation of pulse X-rays to the end of generation of pulse X-rays is divided into three sections (A) in order due to differences in imaging control by the system control unit 8, more specifically, generation control of pulse X-rays. , (B), (C).

  The section (A) is a period in which the contrast medium is not yet administered, and the time length is about 1 to 3 heartbeats. In the section (A), pulse X-rays are repeatedly generated at an indefinite period under electrocardiographic synchronization. For example, in the section (A), the system control unit 8 detects an R wave of the electrocardiogram (ECG), identifies the end diastole from the RR interval, and generates a pulse X-ray only at that time. That is, the pulse rate in the section (A) is 1 time / heartbeat, and pulse X-rays are repeatedly generated at an indefinite period.

  The detection operation of the X-ray detection unit 5 includes, for example, 30 cycles / second in synchronization with the clock over the entire imaging period including the section (A), the section (A), the section (B), and the section (C). (30 fps) is repeated at a constant cycle that is the inverse of the frame rate.

  Note that a collection method may be employed as the electrocardiogram synchronization and at a constant cycle (for example, the reciprocal of 30 fps). Specifically, as shown in FIG. 8, for example, if an R wave is detected during collection at a fixed period (for example, the reciprocal of 30 fps), the rhythm that has been irradiated at a fixed period until then is reset and fixed again. Collect at a period (for example, the reciprocal of 30 fps). According to this technique, the frame immediately after the R wave is synchronized with the electrocardiogram at any heartbeat, and thereafter, it is collected at a high speed, for example, 30 fps.

  Further, in a section in which collection is performed at a constant cycle (for example, the reciprocal of 30 fps), as shown in FIG. 9, ECG synchronization may be performed and collection may be performed at a certain cycle (for example, 1/10 of the RR interval). Specifically, for example, when an R wave is detected while collecting at a fixed period (for example, the reciprocal of 30 fps), the rhythm that has been irradiated at a fixed period is reset, and the RR interval is calculated to calculate the RR interval. 1/10 of the time is calculated, and the next heartbeat cycle is collected at the calculated cycle (constant cycle). According to this method, the frame immediately after the R wave is synchronized with the electrocardiogram at any heartbeat, and thereafter, it is collected at a high speed.

  Image data corresponding to the generation period of the pulse X-ray is stored in the image calculation storage unit 10. However, image data corresponding to a period in which no pulse X-ray is generated is not stored in the image calculation storage unit 10. The image data is stored together with a code for distinguishing the section (A) before contrast from the sections (B) and (C) described later.

  Injection of a contrast agent is started by a trigger operation of the injector 16. A contrast medium injection start signal is supplied from the injector 16 to the system controller 8. The system control unit 8 ends the section (A) and starts the section (B) when a predetermined delay time Δt1 has elapsed from the start of contrast agent injection. The delay time Δt1 is set to, for example, about 1 second from the time when the contrast medium injection signal is output from the injector 16 until the contrast medium exits from the catheter tip.

  In the section (B), the pulse X-ray is repeatedly generated at a constant period. For example, in the section (B), the system control unit 8 repeats at a constant cycle of the reciprocal number of the same pulse rate (30 times / second) as the frame rate of 30 cycles / second (30 fps) of the detection operation in synchronization with the clock. Pulse X-rays are generated. The time width (duration time) of the pulse X-ray is set to be equivalent to the charge accumulation period length in the detection operation of the X-ray detector 5. As a result, the pulse X-ray is generated in synchronization with the detection cycle of the X-ray detector 5. In the section (B), image data is generated at a frame rate of 30 fps, and is all stored in the image calculation storage unit 10. Image data is stored together with a code that distinguishes the section (B) in the contrast period from the sections (A) and (C).

  The injector 16 supplies a contrast medium injection end signal to the system control unit 8 when the injection of a predetermined amount of contrast medium is completed. The system control unit 8 ends the section (B) and starts the section (C) when a predetermined delay time Δt2 has elapsed from the end of contrast agent injection. The delay time Δt2 is longer than the delay time Δt1, and is set to, for example, a time (about 2 seconds) that is twice the time required for the contrast agent to exit from the catheter tip after the contrast agent injection signal is output from the injector 16. . The section (B) is, for example, about 3 to 5 heartbeats or about 5 seconds.

  In the section (C), similarly to the section (A), pulse X-rays are repeatedly generated at an indefinite period under the synchronization of the electrocardiogram. Even in the section (C), the system control unit 8 detects the R wave of the electrocardiogram (ECG), identifies the end diastole from the RR interval, and generates a pulse X-ray only at that time. The pulse rate in the section (C) is 1 time / heartbeat, and pulse X-rays are repeatedly generated at an indefinite period. The image data is stored together with a code for distinguishing the post-contrast section (C) from the sections (A) and (B).

  In the above description, in the sections (A) and (C), pulse X-rays are generated at an indefinite period due to electrocardiogram synchronization. However, as shown in FIG. 3, pulse X-rays are repeatedly generated at a constant period. You may do. Typically, the pulse X-rays in the sections (A) and (C) are repeatedly generated at a constant period (1 second) longer than the period (1/30 seconds) in the section (B). In this case, the pulse rate in the section (B) is 1 time / second.

  According to the pulse rate control described above, exposure can be reduced as compared with the case where the high pulse rate is continued over the entire period. In addition, a relatively high time resolution can be ensured during the period in which the contrast medium is distributed to the site of interest. During the other periods, images of the end diastole required for diagnosis can be collected even with relatively low time resolution.

  The end of the section (C) is determined by the system control unit 8 in response to the image processing result of the image calculation storage unit 10. For example, a reference image (control image) generated from a plurality of images collected in the section (A) is subtracted from each of the images collected in the section (C), and the density of the whole or local region of the obtained difference image The total or average concentration (measured value) is measured in real time. The control image is generated and stored by the image calculation storage unit 10 as an average image of a plurality of images collected in the section (A). Random noise can be suppressed by the averaging process.

  As illustrated in FIG. 4, the measured value of the difference image changes from the initial value L0 with the passage of the contrast agent to the imaging region or the region of interest as time passes (lower in FIG. 4). Along with this, it returns to the initial value L0. Therefore, a threshold value approximating the initial density value L0 is set to a value of a · L0 in advance, for example, 0.9 <a <1, and at time t4 when the measured value reaches the threshold value. The system control unit 8 outputs an X-ray irradiation stop signal to the high voltage control unit 20. The high voltage control unit 20 receives the X-ray irradiation stop signal from the system control unit 8 and ends the application of the tube voltage from the high voltage generation unit 1 to the X-ray tube 21. In addition, when the section (C) has passed, for example, 30 seconds, the system control unit 8 outputs an X-ray irradiation stop signal and forcibly X-ray irradiation even if the measured value does not become the threshold value or less. Quit.

  The end of the section (C) may be determined by the slope of the concentration time curve indicated by the broken line in FIG. The section (C) ends at time t3 when the gradient of the density time curve repeatedly calculated in the image calculation storage unit 10 reaches a threshold value approximated horizontally.

  The image calculation storage unit 10 has various image processing functions. For example, an image necessary for myocardial perfusion is extracted from the image obtained in the section (B). For this, an image closest to the time of the end diastole is identified for each heartbeat from the images collected in the section (B). Since the image is collected at 30 fps in the section (B), the time shift error is within 33 msec at the maximum. This is an amount that is sufficiently problematic for myocardial perfusion measurement. A composite image is generated by combining each of the plurality of images specified in the section (B) and an image selected from the plurality of images obtained in the section (C). A plurality of difference images are generated by subtracting the control images in the section (A) from the plurality of synthesized images. By arranging and displaying a plurality of difference images according to the passage of time in the section (B), it is possible to grasp the state of the contrast agent flowing out.

  In the above description, the section is switched automatically, that is, the pulse rate (frame rate) is changed by the start / end of the injection of the contrast agent. However, according to the pulse rate switching instruction of the operation unit 9 by the operator, The pulse rate may be manually switched at an arbitrary time.

  In the above description, the three sections are used. However, the first section (A) may be omitted and two sections may be used. In this case, the control frame uses the first frame of the section (B).

  After obtaining the difference image, quantitative evaluation may be performed. In this case, a region of interest (ROI) is set on the image and the value is read.

  When strictness is required, it is preferable to consider the time difference between the vascular region of interest and the myocardial region of interest. In other words, since the blood flow flows through a thick blood vessel and then obtains a small blood vessel and reaches the myocardium, there is a time difference. Actually, there is a time difference of about 300 msec to 1 sec. This time difference is even longer in patients with lesions. Therefore, when a vascular region of interest is set for a thick blood vessel, as a result, the data of the region of interest can be accumulated from time 0 to (T−τ), and the myocardial region of interest for this can take data at time T. preferable. When it is desired to do such strict things, it is preferable to take two pieces of data for one heartbeat in the previous image acquisition, that is, X-ray pulse irradiation is performed twice for one heartbeat.

  In the above-described quantitative measurement, when it is desired to further reduce the noise, it is preferable to perform an averaging process on neighboring frames. That is, as shown in FIG. 5, at the moment when the heart motion is small, a plurality of, for example, two pulse X-rays are generated as one body, and the measurement is performed from an averaged image of the two images. Thereby, noise is reduced as compared with the case of only one image. However, since the elements of the heart motion are added, it is essential to output a plurality of pulses in a time when there is almost no heart motion. Alternatively, there is a means for correcting motion.

  It is okay to obtain the myocardial perfusion image and quantitative value even if it is spatially rough. Therefore, for example, even if the raw image is an image of 1024 × 1024 matrix, the myocardial perfusion image may be a 256 × 256 matrix. In view of this, in the heartbeat synchronization collection period, the image is collected by adding pixels so that the X-ray absorption amount per pixel may be small instead of making the image coarse. For example, when adding 4 pixels, the X-ray dose is reduced to 1/4.

  In order to calculate myocardial perfusion, the X-ray conditions must not change between frames. If the X-ray conditions change, the density of the image also changes, so it makes no sense what is being measured. Therefore, at least one image per heartbeat in the proposed collection method is taken with the same kV, mA, and msec.

  Here, as a further contrivance, only kV may be fixed and mAs may be varied. In this case, the mAs value for each frame is recorded, and when creating a myocardial perfusion image, the image is corrected based on the recorded mAs value. The reason why this is true may be referred to a simple X-ray physical equation, and the fluctuation in mAs is simply proportional to the fluctuation in concentration.

In addition, various modifications are possible as follows.
The image density value was compared with a threshold value to determine the end of the final section (C). The desired perfusion calculation value was calculated in real time, and the end determination was made when sufficient data was obtained for the calculation. May be. In addition, the imaging may be terminated when sufficient data is obtained in the section (C), and the imaging may be terminated when sufficient data is obtained to draw an approximate curve. At that time, if sufficient data is obtained in the section (C), it may be notified that the imaging may be ended, and an operator's end instruction may be waited. The notification is performed by sounding a message display, a beep, or the like.

The perfusion measurement value is calculated in real time during the imaging period in the image calculation storage unit 10 and displayed on the display unit 11.
A section where pulse X-rays are generated at a constant period and a section where pulse X-rays are generated in synchronization with electrocardiogram may be manually switched by an operator. For example, the imaging button provided in the operation unit 9 has two stages and is switched by pressing. Although it has a sub button and only the main button generates pulse X-rays at a fixed period, it generates pulse X-rays in sync with the heartbeat while pressing the main and sub simultaneously. The sub button may be in the vicinity of the main button or may be separated (a foot switch and a hand switch). When the operator releases the image pickup button in the last section, the operation ends.

  In the heartbeat synchronization, when one irradiation (one pulse) is performed for one heartbeat, one of the end systole, the end diastole, and the middle diastole is selected for one irradiation by the operator.

In the case of 2 irradiations (2 pulses) per heartbeat, for example, the first irradiation is after a predetermined time from the R wave, and the second irradiation is after a predetermined time different from the R wave. The term “after” may be managed by an absolute time from the R wave, or may be managed by a relative time based on the RR interval. For example, the first irradiation is any one of the end systole, the end diastole, and the middle diastole. The second irradiation has a timing that is a predetermined time before the first irradiation. The predetermined time is the time required for the contrast medium to flow from measurement point A to measurement point B. Three irradiations (3 pulses) or more per heartbeat, 3 pulses are adjacent, and the obtained multiple images are averaged. Used for. One irradiation (one pulse) per multiple heartbeats. When the time of heartbeat synchronous imaging in the section (B) exceeds a predetermined time, the system automatically shifts to 1 irradiation (1 pulse) per K heartbeat.

  The X-ray energy is the same in the heartbeat-synchronized irradiation. Although the same kV and mAs are basically used for all the pulse X-rays, the kV is the same, but the mAs may be different. At that time, the mAs value is recorded and corrected by a subsequent process.

  During playback, the image is processed and displayed for easy viewing. Considering the generation of image data at a fixed interval from the data of a plurality of images collected under control with varying pulse rates, both sections (A, C) and (B) are used as a common clock. On the other hand, the generation timing of the pulse X-ray is synchronized. As a result, the data with the same cardiac phase can always be generated when generating data.

  In the sections (A, C), images may be collected at the same cycle as the section (B). When an image sequence based on electrocardiogram synchronization is generated, an image at the closest time is specified and generated. An image group with a fixed period and an image group with a heartbeat interval are reproduced separately. When playing back, it is displayed in real time.

(Second Embodiment)
The second embodiment will be described below with reference to the drawings. First, terms used in the following description are defined as follows.
Myocardial perfusion: Myocardial dense staining, myocardial perfusion, myocardial brush (Blush), which are technically different, express the phenomenon of blood flowing into and out of capillaries.
Micro circulation; blood flow in capillaries.
X-ray angio apparatus; mainly used for blood vessel imaging in an X-ray diagnostic apparatus.
X-ray image: an image representing the intensity distribution of X-rays transmitted through a subject, also referred to as an X-ray image.
X-ray moving image: A data set of a series of X-ray images taken repeatedly over time t by a two-dimensional detector.
Coronary artery imaging: Imaging of an X-ray image of the coronary artery lumen emphasized by a contrast medium.
Coronary artery; coronary artery
Pixel; pixel in the collected X-ray image
Fusion image: An image obtained by combining an X-ray image and an X-ray image by overlay.
Natural logarithm; Ln
Cardiac phase: The R wave detected by the electrocardiogram signal is used as a mark, the time from the R wave to the next R wave is normalized by 100%, and the time of the current frame is expressed by%. For example, the cardiac phase at the end systole is around 25%.

ROI: Region Of Interest
Catheter room; almost synonymous with catheter room, angiography room, and vascular interventional room.

  FIG. 10 shows the appearance of the X-ray diagnostic apparatus according to the present embodiment, here, the X-ray angiography apparatus, and FIG. 11 shows a functional block diagram. The X-ray imaging apparatus has a gantry 100. The gantry 100 has a C-arm 7. The C arm 7 is rotatably supported by the mechanism control unit 6. An X-ray generator 2 is attached to one end of the C arm 7. The X-ray generation unit 2 includes an X-ray tube 20 and an X-ray restrictor 21. The high voltage generator 1 generates a high voltage (tube voltage) applied between the electrodes of the X-ray tube 20 and generates a filament current to be supplied to the filament of the X-ray tube 20. The high voltage controller 17 controls the tube voltage and / or the filament current generated in the high voltage generator 1 according to the control of the system controller 8.

  An X-ray detector 5 is attached to the other end of the C arm 7. The X-ray detection unit 5 includes a detector 18 and an image data generation unit 19 that processes the output of the detector 18 to generate image data. The detector 18 faces the X-ray tube 20 of the X-ray generation unit 2 with the subject 3 placed on the bed 4 interposed therebetween. The detector 18 is typically a solid flat detector in which a plurality of detection elements (pixels) that directly or indirectly convert incident X-rays into electric charges are two-dimensionally arranged. The X-ray detection unit 5 repeats a one-cycle detection operation including charge accumulation, charge read-out, and reset under a certain period under the control of the system control unit 8. An electrocardiogram monitor terminal 15 is attached to the subject 3. The electrocardiogram monitor receiving unit 16 receives a signal from the electrocardiogram monitor terminal 15 and generates electrocardiogram data of the subject 3. An operation unit 9 is connected to the system control unit 8. The operation unit 9 is provided with a user interface 14 having a hand switch 12, a display, a touch panel, and the like.

  The arithmetic processing unit 23 includes an image collection / storage unit 10 and an image processing unit 11. The image acquisition and storage unit 10 stores the image data output from the X-ray detection unit 5 in association with cardiac phase data. Based on the data of a plurality of images generated under the cardiac coronary angiography imaging sequence, the image processing unit 11 performs a time density curve (reference density region) related to a reference region set on the coronary artery as a blood supply region to the myocardium. A plurality of indices each representing the correlation of a plurality of time density curves with respect to a plurality of local regions set in the myocardial region with respect to a reference time concentration curve), a map of the index is created, and the map is defined as a coronary artery image. A composite image is generated by combining. Note that the time density curve may target the density of the original image, or may target the density (approximate value of the contrast agent density) on the difference image between the mask image before contrast and the contrast image after contrast. Here, a description will be given as a time density curve related to the density on the difference image.

In this embodiment, three types are provided as indices.
First index K ₁ as the subject a contrast agent inflow period until the contrast medium injected ends from contrast injection start, input function a reference time density curve of a coronary artery, and the output function of time density curve of the local myocardial region At this time, it is calculated as “an index K ₁ representing the blood inflow state” regarding the local myocardium. Second indicator K ₂ are, as a target the contrast agent outflow period to completion taken from contrast injection ends, the input function a reference time density curve of a coronary artery, when the output function of time density curve of the local myocardial region , Calculated as “index K ₂ representing blood outflow state” for the local myocardium. The third indicator _{K 3} is calculated from the index _{K 1} and the index _{K 2} Prefecture. Note that the time density curve at this time is a signal component proportional only to the amount of contrast agent when the contrast agent administered to the human body is distributed to or passes through the region of interest (coronary artery and myocardium in this embodiment). Define it.

  12 and 13 illustrate coronary angiographic images (hereinafter referred to as CAG images) obtained by general X-ray angiography. In the CAG image, the X-ray absorbed dose of the contrast agent is large, so that the contrast can be identified so that the shape of the coronary artery can be distinguished from other tissues while passing through the coronary artery (see FIG. 12). As shown in FIG. 13, the catheter 101 is inserted up to the coronary artery 102, and a contrast medium is continuously injected from the catheter 101 at that position for a certain period of time. X-ray imaging is performed at least over a period from before the start of contrast agent injection to the lapse of a predetermined time after the contrast agent injection ends.

  When the contrast medium is injected into the coronary artery via the catheter 101, an X-ray image is collected together with an electrocardiogram by the X-ray diagnostic apparatus of FIG. After the image collection, the reference region (myocardial blood supply region) 103 on the coronary artery and a plurality of myocardial local regions 104 are set on the myocardium on the CAG image by the operator via the operation unit 9 (FIG. 32A, FIG. 33).

  The plurality of myocardial local regions 104 are set on the myocardial region and typically have a plurality of pixels. The density corresponding to the contrast agent amount in the myocardial local region 104 is typically calculated as a pixel average value. However, the myocardial local region 104 may have a single pixel. The myocardial blood supply region 103 typically has a rectangular shape having a width substantially equal to or slightly smaller than the blood vessel, is set in an arbitrary direction along the blood vessel, and includes a plurality of pixels. The concentration corresponding to the amount of contrast medium in the myocardial blood supply region 103 is typically calculated as a pixel average value. The myocardial blood supply region 103 may be any portion of the flow path between the injector and the myocardial region of interest, more specifically, any portion of the catheter, or between the catheter outlet (equivalent to the coronary artery inlet) and the myocardial region of interest. Set to any part.

  The image processing unit 11 generates a time density curve (TDC) regarding the myocardial blood supply region 103. Similarly, the image processing unit 11 generates a plurality of time density curves related to the plurality of myocardial local supply regions 104. The details of this generation process are shown in FIG.

The pixel value at time t in the coordinates (x, y) of the CAGX line image in FIG. 32 is CAG (x, y, t), and the pixel value in FIG. 34 is ln (CAG (x, y, t)). Accordingly, if the attenuation coefficient per unit volume of the tissue or contrast agent is μ and the distribution thickness of the tissue or the preparation agent in the X-ray transmission direction of FIG. 21 is L, the pixel of FIG. 14 is represented by ln (CAG (x, y, t)) There is a relationship of ∝μL. The total amount of contrast agent distributed in the coronary artery or myocardial image recognized on the two-dimensional image shown in FIG. 13 is proportional to μL. That the pixel values of the image T0 at time t = t ₀ with no contrast agent _{CAG (x, y, t 0} ) When,
ln (CAG (x, y, t)) - ln (CAG (x, y, t 0))
Is proportional to the total amount of contrast medium distributed in the coronary artery or myocardial part image recognized on the X-ray image of FIG. 13 and is relatively coronary in the series of images indicated by images T1 to T5 shown in FIG. And the total amount of contrast medium distributed in the myocardial part image and the like.

  That is, since the CAG image is a transmission image from one direction, an X-ray hits substantially perpendicularly to the myocardium in the portion to be examined in the vicinity of the center of the image, usually in the center of the image, so that the thickness effect does not occur much. However, since the X-ray strikes the myocardium in a direction close to the tangent at the peripheral portion outside the center of the image, the thickness effect appears strongly. That is, the integral distance of the contrast agent concentration is longer at the peripheral portion than at the central portion. Accordingly, the reliability of the contrast agent concentration measured at the peripheral portion of the heart is lower than the reliability of the contrast agent concentration measured at the central portion. However, as described above, the myocardial examination target portion is usually positioned near the center of the image, so that an index can be obtained with high reliability for that portion. In addition, you may exclude the parameter | index regarding the peripheral part of the heart with low reliability which was separated from the image center which is not a test object from a display object.

  The blood supply region 103 in FIG. 14 is drawn in a circle for easy understanding on the illustration, but its shape passes through the catheter 101 and flows into the coronary artery as shown in FIG. 13 when mixed with arterial blood. For example, a shape close to a rectangle is set so as to surround the coronary artery so that the TDC can be measured. A local region 104 on the myocardium is a dominant region of the coronary artery into which the contrast medium is injected by the catheter 101 (an area of the myocardium that is active by supplying blood that has passed through the coronary artery). For the entire myocardial region suspected of the presence or absence of microcirculation in the previous capillary blood vessels, for example, N × N pixels (N = 1, 3, 5,... Odd number) indicated by the region of interest 104 are sampled with N = 5, for example. For example, after the smoothing process of 5 × 5 size, the calculation result described below is recorded with the center position of the pixel as the representative position.

[Observation]
As described above, the observed value ln (CAG (x, y, t)) - a ln (CAG (x, y, t 0)). Here, the reference TDC of the contrast agent relative concentration per unit blood vessel volume for the myocardial blood supply region 103 is Ca (t), and similarly, the contrast agent relative concentration per myocardial unit volume is Cmyo (t). Then, the relationship between the observed value and Ca and Cmyo is expressed by the following equation. Since Ca represents a difference value before and after contrast agent injection, the initial value Ca (0) is zero.

Here, (x _a , y _a ) is a pixel in the region of interest 103 set upstream of the coronary artery in FIG. 13, and represents an average value or a total amount in the region of interest. Either may be used, but in the following operation, only one selected is used. Do not mix. t represents time, and the unit is seconds. a means an artery.

Here, (x _myo , y _myo ) is a pixel in the local region of interest 104 set in the myocardium shown in FIG. 13, and represents an average value or a total amount in the region of interest 104. The plurality of regions of interest 104 have the same shape and the same shape obtained by dividing the perfusion calculation range 105 set on the CAG image via the operation unit 9 shown in FIGS. 15A and 15B as shown in FIGS. 32A and 33. Multiple local regions of size. The subscript myo means myocardium (coronary artery).

[model]
Next, (theory 1) described later is introduced as a model for the inflow and outflow of the contrast medium into the myocardium. As a result, the following relationship is obtained.

[Adaptation of observation value to model]
Images T <b> 1 to T <b> 5 are CAG images corresponding to a specific cardiac phase 106 among a series of CAG images captured in synchronization with the electrocardiogram. The collection time is t1, t2,. For the set of electrocardiogram synchronization images T1 to T5, the image processing unit 11 performs processing described by the equations (1), (2), and (3) derived from [Observation] and [Model].

The image processing unit 11 calculates an index K ₁ ′ reflecting the actual state of blood inflow into the local region of the myocardium that is proportional to the myocardial blood flow. First, the contrast agent is a contrast agent injection time to the point of complete injection from the time it is started injected from the catheter is set as a calculation target period for the index K _1. As shown in FIG. 16, from a plurality of X-ray images collected during the contrast agent inflow period, a time concentration curve Ca (t) relating to the blood supply region 103 to the myocardium and a plurality of times relating to the plurality of myocardial local regions 103 are obtained. A density curve Cmyo (t) is generated.

Next, in the image processing unit 11, as shown in FIG. 17, the time integration (contrast medium inflow amount) ∫Ca (t) passing through the blood vessel supply region 103 is taken as the horizontal axis, and the myocardium is taken along the vertical axis. A value for each time is plotted as the contrast medium existing amount Cmyo (t) in the local region 104. That is, a discrete distribution of the time variation of the blood uptake amount in the myocardial region with respect to the blood supply amount to the entire myocardium is generated. An index (first index) K ₁ ′ is calculated as the slope of the straight line obtained by the straight line fitting process for this discrete distribution. This index K ₁ ′ has significance as a quantification of how much the myocardial region follows the blood supply to the whole myocardium and accepts blood. When the index K ₁ ′ deviates from the normal range to the lower value side, the blood inflow in the myocardial region does not follow the blood inflow of the coronary artery, that is, the injured region in the myocardial region there is a possibility. When the index K ₁ ′ deviates from the normal range to the high value side, the myocardial local region may have an outflow disorder. The image processing unit 11 calculates the index K ₁ ′ by the same process for all the local regions 104 and creates a map of the index K ₁ ′. The obtained map is superimposed (fused) on the CAG image to generate a composite image, which can be displayed on the display unit 24. Thereby, the myocardial function can be determined together with the positional relationship with the myocardial tissue.

Modification A)
In the above description, a series of ECG-synchronized CAG images T1 to T5 are extracted from continuous images. In order to reduce exposure to the patient, X-rays are generated only in a specific ECG phase at the time of imaging. It is also possible to collect only an image set necessary for (theory 1) described later by controlling the apparatus.
Modification B)
Further, according to (Theory 2) described later, K obtained by linear approximation by the graph plot method shown in FIG. 19 from a series of TDC curves (FIG. 18) in the outflow period in which the contrast medium reaches the myocardial region and is excreted further. ₁ and can be measured back diffusion (the amount of the contrast medium returns from the myocardium into a blood vessel) from the myocardium of K ₂ by blood and contrast agent may be as before display with reference to FIG. Similarly, at time t _n, to reduce patient exposure, only the set of images required for Theory 1 is acquired by controlling the apparatus of FIG. 10 to generate X-rays only at specific ECG phases during imaging. You can also.

Modification C)
Further, by calculating and displaying the ratio of K ₁ and K ₂ (K ₁ / K ₂ ), “C. Michael Gibson, MS, MD; Christopher P. Cannon, MD; Sabina A. Murphy, MPH; Kathryn A Ryan, BS; Rebecca Mesley, BS; Susan J. Marble, RN, MS; Carolyn H. McCabe, BS; Frans Van de Werf, MD, PhD; Eugene Braunwald, MD; for the TIMI (Thrombolysis In Myocardial Infarction) Study Group, Relationship of TIMI Myocardial Perfusion Grade to Mortality After Administration of Thrombolytic Drugs, Circulation. 2000; 101: 125-130. ”Can be quantified (automated). That is, it is possible to automatically classify according to the classification table of FIG. 28, or to plot and express on a graph representing the relationship between the indices K ₁ and K ₂ as shown in FIGS. 29A and 29B.

Modification D)
According to the present embodiment, two types of ECG-synchronized CAG images in the same patient and the same direction are collected. That is, as shown in FIGS. 20A and 20B, a myocardial blood flow image 205 (= Prest (x, y)) at rest and a drug after administration of a drug having an effect of increasing myocardial blood flow, such as adenosine, for example. A myocardial blood flow image 206 (= Pstress (x, y)) at the time of loading is collected. By defining (myocardial blood flow image 206) / (myocardial blood flow image 205) = Pstress (x, y) / Prest (x, y) as a comparison image between resting and loading, myocardial blood flow images 205 and 206 are Since the myocardial thickness in the imaging direction is the same at the same position (x, y) collected in the same direction, it is approximately canceled by division, and a blood flow increase ratio (myocardial flow reserve) per unit myocardial volume can be obtained. .

Similarly, by measuring and dividing an index K ′ _before treatment such as stenosis and an index K ′ _after treatment, a blood flow increase ratio (recovery ratio by treatment) per unit myocardial volume can be obtained.

Modification E)
In addition, according to the embodiment, as illustrated in FIG. 27, Cmyo (t) = ln (CAG (x, y, t) −ln ( CAG (x, y, t ₀ )), the image of Cmyo (t) -Cmyo (t ₀ ) shows the blood flow to the myocardium, and this differential moving image can be displayed during the same examination. Only possible, where t ₀ = immediately before contrast agent injection.

(Theory 1)
The contrast medium relative concentration Cmyo (t) in the myocardium at time t, the contrast medium relative concentration time curve in the coronary artery blood is represented as Ca (t), the myocardial inflow blood flow is represented as KK1, and the myocardial outflow blood flow is represented as KK2. In FIG. 23, the mass balance of the contrast agent in the myocardium and coronary blood can be expressed by equation (1).

The blood flow in the myocardial part or a parameter correlated with the blood flow can be calculated by the Patlak plot method using the phenomenon that the contrast agent described below flows into the myocardial part. Conditions under which contrast medium begins to flow into the myocardium (0 ≦ t ≦ MTT) (MTT = mean tangent time = average transit time from contrast medium entering the myocardium and exiting the vein, usually 5-10 sec for myocardium In this case, the outflow amount k2C _myo (t) is extremely small, and the condition of C _myo (t) ≈0 and K ₁ Ca (t) >> K ₂ C _myo (t) is satisfied. ) Can be simplified as follows.

K ₁ represents a transition constant (sec ⁻¹ ) or (ml / min / g) at which the contrast medium moves from the blood to the myocardial stromal region, and is proportional to the blood flow.
Rutland MD.A single injection technique for subtraction of blood background in 131I-hippuran renograms.Br J Radiol 1979; 52: 134-137.
Patlak CS, Blasberg RG, Fenstermacher JD. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data.J Cereb Blood Flow Metab 1983; 3: 1-7.
(Theory 2)
The contrast medium relative concentration in the myocardium at time t is expressed as Cmyo (t), and the contrast medium relative concentration in coronary artery blood is expressed as Ca (t). In FIG. 10, the mass balance of the contrast agent in the myocardial part and coronary artery blood can be expressed by equation (1).

  The blood flow in the myocardial part or a parameter correlated with the blood flow can be calculated by a method using the clearance of the contrast agent flowing out from the myocardium described below.

Eq. (1) is integrated to obtain (2).

When Y (t) is plotted on the Y axis and X (t) is plotted on the X axis, the slope of the straight line is K ₂ , and the intercept of the Y axis is K ₁ . K ₁ is a transition constant (sec ⁻¹ ) or (ml / min / g) at which the contrast medium moves from the blood to the myocardial stromal region, and is proportional to the blood flow. K ₂ represents a transition constant (sec ⁻¹ ) in which the contrast agent moves from the myocardial stromal region to the blood (refer to the literature with FIG. 31 below).

Yokoi T, Iida H, Itoh H, Kanno I. A new graphic plot analysis for cerebral blood flow and partition coefficient with iodine-123-iodoamphetamine and dynamic SPECT validation studies using oxygen-15-water and PET.J Nucl Med 1993; 34 : 498-505.
(Actual processing flow)
Hereinafter, an example of an actual image processing flow will be described. The example shown here is an example, and other processing methods based on the above-described theory may be used.
The overall specific processing flow of this embodiment will be described with reference to FIG. The X-ray image is a moving image showing an angiographic subject (S11). For example, a moving image has a duration of 3 to 60 seconds. It is a moving image having 10 to 30 images per second. The images are taken in the same direction and the same X-ray conditions. It is assumed that there is no patient movement or bed movement during imaging.

  The image processing unit 11 extracts one image corresponding to the specific heart phase for each heartbeat (S14), and differentially processes the image from the image of the concentric phase before contrast agent injection (S12). Specifically, when there are a plurality of concentric phases before contrast agent injection, the averaged image is subjected to differential processing.

  The image processing unit 11 performs motion compensation processing (S13). A motion compensation amount is detected by pattern matching processing based on an image pattern (blood vessel, heart wall, cate). In a frame having no characteristic pattern, the movement amount obtained at the same phase at another time is substituted. The motion compensation has a table for moving the entire image or knowing the ROI movement position. Motion compensation also compensates for some body movements and respiratory movements of the patient.

  The specific cardiac phase is designated in advance by the operator via the operation unit 9. For example, it is specified at the end diastole. The frame at the closest time is extracted from the multi-frames that have been pulse-imaged. Different cardiac phases may be used in the concentration time curve Ca of the blood supply region 103 and the concentration time curve Cmyo of the myocardial local region 104. Cmyo (t + T) is used for Ca (t). Here, T is a delay value required for the contrast medium to flow from the blood supply region 103 to the myocardial local region 104. Specifically, it is often about 1 to 10 frames (1/30 to 1/3 sec). This delay time is generally called TIMI frame count (TFC) or corrected TIMI frame count (CTFC).

In the image processing unit 11 calculates an index K ₁ as slope representing the correlation of Cmyo (t) to the concentration-time integral Σca (t) (S15). A regression line is calculated by fitting as illustrated in FIG. The regression is calculated only in a predetermined time interval, and this predetermined interval is determined in conjunction with the contrast agent injector. That is, if the contrast medium injection start time from the injector is J1 and the injection end time is J2, the start time t1 of the predetermined section is, for example, t1 = J1 + 1 sec, and the end of the predetermined section is, for example, J2 + 2 sec. Further, the predetermined section may be determined from the value of the image. That is, the time when the TDC rises is defined as t1, and the time when the TDC becomes flat is determined as the end of the predetermined section. Calculate the correlation value with the regression line. It is preferable to exclude heartbeats including arrhythmias from the calculation. Calculating one K ₁ for one of the moving image. However for the purpose of speeding up, or a method of calculating the needed K ₁ for each frame. As described above, the average value of the surrounding N × M pixels (may be Median) is used as the value of the pixel, not for each pixel. The display unit 24 displays a blood flow image (S17). To view a map of the index _{K 1.} Map of K ₁ may be displayed by converting the density values to the color values. Collected displays the original motion picture, may be displayed overlapping the maps of K ₁ immediately after that, may be displayed side by side map of the original image and the K ₁ to the left and right.

  Corresponding to Modification A described above, this embodiment may integrate image processing into the X-ray diagnostic apparatus (FIGS. 10 and 22), or an image processing apparatus that is independent from the X-ray diagnostic apparatus. (FIGS. 30 and 23) may be provided. The original image (blood vessel image, captured image) is superimposed and displayed on the blood flow image (index image, map image) (S16). A black and white original image (blood vessel image) is superimposed and displayed on a color blood flow image obtained by applying the index to the color table as it is. The original image is a moving image, and a blood flow color still image may be superimposed on the blood vessel moving image. Actually, the blood flow color still image is superimposed on the blood vessel still image. In order to create a still image from a moving image, the darkest one frame is selected from the moving images, or the smallest pixel value of the moving image is detected for each pixel, and an image is obtained. The original image is an image before the difference or an image after the difference. It is characterized by displaying a graph (S21). The graphs of FIGS. 16 to 19 are displayed for a specific ROI (x, y). In the graphs of FIGS. 17 and 18, the slope, intercept, and correlation value of the fitting function (regression line) are displayed. The blood supply area 103 is designated by the operator. A graphical user interface 14 for designating the blood supply region 103 on the image is provided. One pulse X-ray is irradiated at a specific heart phase of one heartbeat (FIG. 23).

The image processing unit 11 calculates an index _{K 1} and _{K 2} (Figure 24). The ratio of K _{2 to} index K ₁ (K ₂ ) / K ₁ is calculated (FIG. 25A). The image processing unit 11 generates a map of the ratio of _{K 1} and _{K 2.}

Further, the image processing unit 11, as shown in FIG. 25B, the K ₁ on the horizontal axis, by arranging the K ₂ on the vertical axis, plots points corresponding to K ₁ and K ₂ for each local region. In the distribution, each myocardial local region can be classified into four types representing myocardial functions such as “normal”, “inflow disorder”, “outflow disorder”, and “outflow inflow disorder”. That classifies index K ₁ the local myocardial region and distinguishing high and low with respect to the first threshold value TH1, the combination of the distinction between high and low for the second threshold TH2 of the index K _2. Thereby, the myocardial function can be determined for each region.

Further, the image processing unit 11 assigns four codes (1 to 4) to the four sections as shown in FIG. 32B, applies the codes to the myocardial local region as shown in FIG. 32A, and displays them over the original image. Alternatively, as shown in FIG. 33, the four sections may be distinguished from each other by different display modes, for example, brightness or hue, and may be displayed superimposed on the original image. By this display, the state of myocardial dysfunction can be spatially grasped together with the distinction.
The division may be classified as shown in FIGS. 28 and 29B.

Further, the image processing unit 11 is equipped with a comparative analysis process (FIG. 26). For example the time A and time B in the same operation (C, D, E ,,,) are displayed side by side on each calculated dominate the index K ₁ of. Or, display by switching over time. By calculating the ratio of K ₁ time A and time B (C, D, E ,,, ), to the image display ratio. The original image is displayed superimposed on the color image of the ratio. For example, by displaying a color image of the ratio before and after treatment superimposed on the original image on which blood vessels are depicted, it becomes possible to visualize and know which region of the myocardium has recovered blood flow by treatment.

The features of the present embodiment described above are as follows.
As a difference from a perfusion calculation method using a contrast agent emitted from an X-ray CT apparatus, image processing is greatly different between CT and X-ray, and processing specific to X-ray is required. Since X-rays are moving images, they require special measures such as extraction of one image per heartbeat, background difference, motion compensation, thickness direction processing, and blood vessel ROI setting. It was.

As a difference from the known method of looking at the slope of the time density curve with X-rays, the idea itself described in [Theory 1, 2] is different. Not in the conventional idea to draw a graph in which the horizontal axis inflow in the index K _1. As a result, subsequent analysis and result images are different. The processing method is completely different from the processing for obtaining the slope of the time density curve.

  The literature (or FIG. 28) MBG (TMP) method was not meaningful for image figures. In the present embodiment, the numerical value of the image has a meaning that it is a flow velocity into the myocardium, which is clearly more useful than MBG. Since normalization is performed with the amount of inflow, it is not affected by the amount of contrast medium to be injected, and quantitativeness can be obtained even in a clinical setting.

  Since it is not visual, it is quantitative. Since the prognosis can be predicted, this value can be used as a reference to determine whether or not the intervention treatment can be completed. Since comparisons can be made before and after treatment (treatment), it is possible to know how much the flow rate has increased, that is, how much the blood supply to the myocardium has improved. By automating the process, it can be easily used in clinical practice.

  Note that the present invention is not limited to the above-described embodiment as it is, and can be embodied by modifying the constituent elements without departing from the scope of the invention in the implementation stage. In addition, various inventions can be formed by appropriately combining a plurality of components disclosed in the embodiment. For example, some components may be deleted from all the components shown in the embodiment. Furthermore, constituent elements over different embodiments may be appropriately combined.

The figure which shows the structure of the X-ray diagnostic apparatus which concerns on 1st Embodiment of this invention. The figure which shows an example of the imaging operation by the system control part of FIG. The figure which shows the other imaging operation by the system control part of FIG. The figure which shows the time density curve for supplementing the imaging completion determination process by the system control part of FIG. FIG. 10 is a diagram illustrating still another imaging operation by the system control unit in FIG. 1. The figure which shows the conventional imaging operation. The figure which shows the conventional imaging operation. The figure which shows the other imaging operation by the system control part of FIG. The figure which shows the other imaging operation by the system control part of FIG. The figure which shows the structure of the X-ray diagnostic apparatus by 2nd Embodiment of this invention. FIG. 11 is a detailed view of a specific part of FIG. 10. The figure which shows an example of a coronary arteriography image. The figure which shows ROI set on a coronary arteriography image in this embodiment. The figure which shows an image collection sequence in this embodiment. The figure which shows the calculation range of perfusion in this embodiment. The figure which shows the calculation range of the perfusion by which setting operation was carried out on the CAG image in this embodiment. The figure which shows the time density curve in the contrast agent inflow period produced | generated by the image processing part of FIG. It shows an index K ₁ calculated by the image processing unit in FIG. 11. The figure which shows the time density curve in the imaging period produced | generated by the image process part of FIG. It shows an index K ₂ calculated by the image processing unit in FIG. 11. The figure which shows the calculation range of the perfusion at the time of rest and load in this embodiment. The figure which shows the comparison image of the rest image and load image which were produced | generated by the image process part of FIG. In this embodiment, the supplementary figure of the thickness effect of a transparent image. The flowchart which shows the process sequence of this embodiment. The flowchart which shows the other process sequence of this embodiment. The flowchart which shows the other process sequence of this embodiment. The flowchart which shows the other process sequence of this embodiment. The flowchart which shows the other process sequence of this embodiment. The figure which shows distribution of 2 types of parameter | index of this embodiment. The flowchart which shows the other process sequence of this embodiment. The flowchart which shows the other process sequence of this embodiment. The figure which shows the classification table for the classification | category of the local region using the 2 types of parameter | index by the image processing part of FIG. The figure which shows distribution of 2 types of parameter | indexes produced | generated by the image processing part of FIG. The figure which shows the division on the parameter | index distribution corresponding to the classification table of FIG. 1 is a diagram showing an image processing apparatus that is independent of an X-ray diagnostic apparatus in the present embodiment. In this embodiment, illustration of the index _K 1, _{K 2.} It shows the myocardial function map generated from the two indices K _1, K ₂ by the image processing unit in FIG. 11. The figure which shows the map code of FIG. 32A. The figure which shows the other myocardial function map of FIG. 32A.

  DESCRIPTION OF SYMBOLS 1 ... High voltage generation part, 2 ... X-ray generation part, 3 ... Subject, 4 ... Bed, 5 ... X-ray detection part, 6 ... Mechanism control part, 7 ... C arm, 8 ... System control part, 9 ... Operation , 10 ... Image calculation storage unit, 11 ... Display unit, 12 ... Hand switch, 13 ... Switch holder, 14 ... User interface, 15 ... Injector, 20 ... High voltage control unit.

Claims (18)

A storage unit for storing data of a plurality of images representing an intensity distribution of X-rays transmitted through an angiographic subject;
Based on the data of the plurality of images, a reference time concentration curve related to a reference region set in a blood supply region to the myocardium and a plurality of time concentration curves related to a plurality of local regions set in the myocardial region are generated. A processing unit that calculates a plurality of indices representing the correlation of each of the plurality of time density curves with respect to the reference time density curve, and generates a map relating to the index;
A display unit that displays a map related to the index superimposed on the image,
The image processing apparatus, wherein the processing unit performs motion compensation on the image data to generate the time density curve. The said process part performs the difference process of the data of the said image in which angiography was carried out, and the data of the image before contrast agent injection | pouring, and produces | generates the said time density | concentration curve. Image processing device. Wherein the display unit, the image processing apparatus according to claim 1, wherein a is intended to overlay the map for the index before Kiga image.   The image processing according to claim 1, wherein the processing unit calculates, for each of the local regions, a correlation of the time density curve with respect to the reference time density curve in the contrast agent inflow period into the myocardium as the index. apparatus.   The image processing according to claim 1, wherein the processing unit calculates, for each of the local regions, a correlation of the time density curve with respect to the reference time density curve in a contrast medium outflow period from the myocardium as the index. apparatus.   The processing unit includes a first index representing a correlation of the time concentration curve with the reference time concentration curve in the contrast agent inflow period into the myocardium, and the reference time concentration curve in the contrast agent outflow period from the myocardium. The image processing apparatus according to claim 1, wherein a second index representing the correlation of the time density curve is calculated for each local region.   The image processing apparatus according to claim 6, wherein the processing unit calculates a ratio between the first index and the second index for each local region, and generates a map of the ratio.   The image processing apparatus according to claim 7, further comprising a display unit that displays the map of the ratio so as to overlap the image.   The image processing apparatus according to claim 1, further comprising a display unit configured to display the reference time density curve as a graph together with the time density curve.   The said processing part classifies the said local area | region into the combination of high / low distinction with respect to the 1st threshold value of the said 1st parameter | index, and high / low distinction with respect to the 2nd threshold value of the said 2nd parameter | index. Image processing apparatus.   The image processing apparatus according to claim 10, wherein the processing unit generates a map related to the classified combinations.   The processing unit calculates a correlation of the reference time concentration curve in the contrast agent inflow period to the myocardium collected after the treatment with respect to the reference time concentration curve in the contrast agent inflow period to the myocardium collected before the treatment. The image processing apparatus according to claim 1, wherein:   The image processing according to claim 1, wherein the processing unit calculates the index obtained from the image data collected before the treatment and the index obtained from the image data collected after the treatment. apparatus.   The image processing apparatus according to claim 13, wherein the index is obtained from data of the image during a contrast agent inflow period into the myocardium. The image processing apparatus according to claim 13, wherein the index is obtained from data of the image in a contrast medium outflow period from the myocardium.   2. The processing unit according to claim 1, wherein the processing unit calculates the index obtained from the data of the image collected before drug administration and the index obtained from the data of the image collected after drug administration. Image processing device.   The image processing apparatus according to claim 6, wherein the processing unit generates a map related to the first index or a map related to the second index. An X-ray tube that generates X-rays;
An X-ray detector for detecting X-rays transmitted through the subject;
An image generator for generating image data representing an intensity distribution of transmitted X-rays based on an output of the X-ray detector;
A storage unit for storing data of a plurality of images representing an intensity distribution of X-rays transmitted through an angiographic subject;
Based on the data of the plurality of images, a reference time concentration curve related to a reference region set in a blood supply region to the myocardium and a plurality of time concentration curves related to a plurality of local regions set in the myocardial region are generated. A processing unit that calculates a plurality of indices representing the correlation of each of the plurality of time density curves with respect to the reference time density curve, and generates a map relating to the index;
A display unit that displays a map related to the index superimposed on the image,
The X-ray diagnostic apparatus characterized in that the processing unit performs motion compensation on the image data to generate the time density curve.

Images