The present invention relates to a biological information imaging apparatus and a biological information imaging method.
Many imaging apparatuses using X-rays, ultrasonic waves, MRI (Magnetic Resonance Imaging), or the like are adopted as apparatuses acquiring biological information in a medical field, and are used for diagnosis.
With respect to an ultrasonic diagnostic imaging apparatus among these, Japanese patent application laid-open No. S62-189054 describes a method of diagnosing the heart, etc., using difference images of ultrasonic tomographic images that are taken at a high frame rate. The apparatus of this Japanese patent application laid-open No. S62-189054 once stores in a memory a plurality of time series tomographic images that are acquired by means of a probe by transmitting and receiving ultrasonic waves to and from a sample or object to be inspected, generates a difference image by subtracting the corresponding pixels of tomographic images for two frames sequentially stored, and displays the difference image thus generated on a screen.
In an ultrasonic diagnostic imaging apparatus, it is difficult to acquire a difference image in an entire region of a range (a part to be diagnosed) in which a diagnose diagnosis is intended to be made, due to the principles of diagnosis by means of ultrasonic waves. That is, in an ultrasonic diagnostic imaging apparatus, first of all, an operator sets up an ultrasonic transmission/reception direction towards a part the condition of which is to be diagnosed, such as for example a moving tissue part in a sample such as a heart, a blood vessel, etc., or a part having a blood flow, and places an ultrasonic probe in abutment with the sample. Then, ultrasonic wave pulses are transmitted from the ultrasonic probe to the interior of the sample. The ultrasonic wave propagated into the sample is reflected on a boundary where the acoustic impedance in the sample varies. The reflected wave (echo) is received by the ultrasonic probe. The echo thus received by this probe is subjected to various processing such as amplification, phasing, detection, etc., by means of a reception circuit, and then, it is digitally converted by an analog/digital converter, and stored into one of line memories in a buffer memory circuit. As a result, data for one line of ultrasonic scanning is stored in the buffer memory circuit.
Subsequently, a CPU controls a transmission circuit so that an ultrasonic wave is caused to be transmitted in a direction different from that of a first transmission/reception wave, e.g., a direction contiguous to the first transmission/reception wave direction. An echo is taken in by making its direction of reception identical to that of this transmission/reception. In synchronization with this second transmission/reception, a first echo signal is read from the buffer memory circuit, and transmitted and stored in a first image memory. On the other hand, a second echo signal is stored into the other line memory of the buffer memory circuit. Thus, the transmission and reception of an ultrasonic wave is performed by changing the direction thereof in a sequential manner, and when the scanning of one image is completed, the transmission/reception direction is returned to the initial direction, and the transmission/reception of an ultrasonic wave are carried out in a repeated manner. Then, a plurality of tomographic images are taken in a manner like a first image, a second image, a third image , . . . , an n-th image.
Thus, in an ultrasonic diagnostic imaging apparatus, a plurality of line images acquired in line sequential order have to be arranged in order to form one tomographic image. That is, a tomographic image obtained by the ultrasonic diagnostic imaging apparatus does not correctly represent the shape of a certain momentary image, but becomes an image containing a shape distortion resulting from the difference in acquisition time of the line images. For example, in cases where a tomographic image of moving tissue such as a heart is acquired, there arises a shift or difference in phase of a motion on the right and the left of the tomographic image. Accordingly, in conventional ultrasonic images, a comparison in image shape between the images of different frames has not been able to be correctly made.
In view of the aforementioned problems, the present invention has for its object to provide a biological information imaging apparatus and a biological information imaging method which are capable of performing comparison between images of different frames in a correct manner.
The present invention in a first aspect provides a biological information imaging apparatus comprising: a measurement unit that measures an acoustic wave emitted from a tissue in a living body at a time when the tissue absorbs light irradiated onto the living body; a storage unit that stores first information obtained in a first measurement by the measurement unit; and an information processing unit that generate an image of information of the tissue in the living body by using second information obtained in a second measurement by the measurement unit and the first information stored in the storage unit.
The present invention in a second aspect provides a biological information imaging method comprising: a first measurement step to measures an acoustic wave emitted from a tissue in a living body at a time when the tissue absorbs light irradiated onto the living body; a storage step to store first information obtained in the first measurement step; a second measurement step to measures an acoustic wave emitted from the tissue in the living body at a time when the tissue absorbs light irradiated onto the living body; and a processing step to generate an image of information of the tissue in the living body by using second information obtained in the second measurement step and the first information stored in the storage step.
According to the present invention, a comparison between images of different frames can be performed in a correct manner. By doing so, it becomes possible to correctly catch a change of a moving tissue, such as a motion of visceral including a heart, a change in thickness or size of a blood vessel due to a blood flow, etc.
Further features of the present invention will become apparent from the following description of exemplary embodiments with reference to the attached drawings.
The present invention provides a biological information imaging apparatus that is constructed in the following manner.
The biological information imaging apparatus of the present invention is a biological information imaging apparatus using a photoacoustic tomography (PAT). This biological information imaging apparatus is provided with a measurement unit that measures an acoustic wave emitted from a tissue in a living body at a time when the tissue absorbs light irradiated onto the living body. The measurement unit comprises, for example, a light source, and an acoustic wave detector that detects acoustic waves (elastic waves, ultrasonic waves) generated from a light absorber (optical absorber) which has absorbed a part of the energy of light irradiated to the living body, and transduces them into electric signals. In addition, the biological information imaging apparatus is provided with a storage unit (memory) that stores first information acquired in a first measurement by the measurement unit. It is desirable that the storage unit have a capacity capable of storing the information acquired in at least one measurement (e.g., a reconstructed image for one frame). Also, the biological information imaging apparatus is provided with an information processing unit, and this information processing unit generates an image of the information of the tissue in the living body by using second information acquired in a second measurement by the measurement unit and the first information stored in the storage unit. Specifically, the second measurement is carried out at a timing at which the shape of the tissue is different from that at the timing of the first measurement, and the information processing unit generates a difference image representing a difference between the first information and the second information. Here, “the shape of the tissue is different” means that at least either of the shape, position and size of the tissue differs. The tissue, which is a target of imaging, is a blood vessel, for example.
PAT is a technique that irradiates a living body with pulsed light generated from a light source, detects acoustic waves emitted from a living tissue which absorbed the energy of the light propagated and diffused in the living body, and performs analytic processing of their signals. With this technology, because of the use of light, it can be found how the light is absorbed in the living tissue, and at the same time, because of performing detection by acoustic waves (elastic waves, typically ultrasonic waves), an optical property distribution in the living body can be obtained at a resolution substantially equivalent to that in an ultrasonic diagnostic imaging apparatus.
In PAT, a sound pressure P of the acoustic waves obtained from a light absorber in the living body due to the light absorption is given by the following formula (1).
P=Γ·μ
a·Φ (1)
Here Γ is a Grüneisen parameter which is an elastic property value, and which is obtained by dividing the product of a volume expansion coefficient β and the square of a speed of sound c by a specific heat CP. In addition, μa is an optical absorption coefficient of the absorber, and Φ is a quantity of light in a local region (a quantity of localized light irradiated to the absorber).
Since it is known that Γ will take an almost constant value if the tissue is decided, the product of μa and Φ, i.e., an optical energy absorption density, can be estimated by measuring the sound pressure P, which is the magnitude of the acoustic waves, in a time sharing manner. In addition, by measuring this at multipoints, a spatial distribution of the absorber in the living body can be acquired.
The basic principle of image picking-up by the PAT technique is briefly as follows.
As stated in the above-mentioned problems, one frame period is required for the conventional ultrasonic diagnostic imaging apparatus to form one image. That is, in general, in cases where a video of 30 frames per second is displayed, images at time points different by 1/30 second from each other will be displayed at the first line and the last line, respectively, in one frame.
On the other hand, in cases where the PAT technique is used, ultrasonic waves can be simultaneously generated from a desired diagnostic region by irradiating it with light. For example, if the diagnostic region is a square area with each side of 10 centimeter, the arrival time differences of the light are about a maximum of 0.3 picoseconds. This is equivalent to an error of 0.5 micrometers when it is assumed that the propagation speed of ultrasonic waves in the living body is 1,500 meters per second. In view of the fact that the resolution required in the PAT technique is on the order of from submillimeters to several millimeters, such an error can be completely disregarded. That is, in the PAT technique, it is possible to acquire an image (image shape) at substantially the same time.
The light source 50 is a unit for irradiating a living body 51 with pulsed light 53. An optical waveguide such as an optical fiber for guiding the light to the living body 51 may be provided as needed. The light source 50 irradiates light of a specific wave length that can be absorbed by a specific component existing in the living body. As the light source, there is provided at least one pulsed light source that can generate pulsed light on the order of from several to hundreds nanoseconds. A laser is desirable as the light source, but it is also possible to use a light emitting diode or the like instead of the laser. As a laser, there can be used various kinds of laser such as a solid-state laser, a gas laser, a dye laser, a semiconductor laser, and so on. In addition, although in this embodiment, an example of a single light source is shown, a plurality of light sources can be used. In the case of the plurality of light sources, in order to raise the irradiation intensity of the light to be irradiated on the living body, there can be used two or more light sources which oscillate at the same wave length, or in order to measure differences in the optical property value distribution according to the wave lengths, two or more light sources having different oscillation wavelengths can be used. Here, note that if a dye of which the oscillating wave length is convertible and OPO (Optical Parametric Oscillators) can be used as the light source, it will also become possible to measure the difference in the optical property value distributions depending upon the wave lengths. It is preferable that the wave lengths to be used be in a range of from 700 to 1,100 nm in which there is little absorption in the living body. However, in cases where an optical property value distribution of the living body tissue relatively near a living body surface thereof, it is also possible to use a wave length range, such as for example a range of from 400 to 1,600 nm, wider than the above-mentioned wavelength range.
The acoustic wave detector 55 is a unit for detecting an acoustic wave 54 emitted from the light absorber 52 which absorbed a part of the energy of light irradiated to the living body 51, and transducing (changing) it into an electric signal. The electric signal outputted from the acoustic wave detector 55 is inputted into the signal processing part 56. As the acoustic wave detector 55, there can be used a transducer using a piezo-electric phenomenon, a transducer using the resonance of light, a transducer using a change of capacitance, and so on. As the light absorber 52, it corresponds, for example, to a blood vessel, a tumor, or like other tissues similar to these. Here, note that this embodiment is suitable in particular for imaging a light absorber in motion.
The signal processing part 56 is a function to analyze the electric signal outputted from the acoustic wave detector 55, and to obtain optical property value distribution information on the living body. In addition, the signal processing part 56 has an image memory (frame memory) that stores image information (image data of tomographic images or three dimensional images) for at least one screen. Moreover, the signal processing part 56 also has a function to compare a first image acquired in a first state at one timing and stored in this image memory with a second image acquired in a second state at another timing, and to calculate a difference between them. The signal processing part 56 can also take any construction and arrangement as long as functions, such as inputting of electric signals, storage and operation of data (including image information), outputting of the results of the operation, and so on, can be achieved. For example, the signal processing part 56 can be comprised of a computer with an analysis program. Here, note that the optical property value distribution information, the difference image information, etc., which are the results of operations of the signal processing part 56, are outputted to a display device (not shown).
The image picking-up by the PAT technique in this embodiment is performed in the following flow. Here, note that a mechanism in which acoustic waves are emitted from the light absorber by irradiation of light is the same as that in the flow of the basic principle as mentioned above.
According to the above-mentioned flow, after performing the first measurement and the second measurement at different timings and reconstructing the images in the two different states, it is possible to extract only those parts which have changed by acquiring the difference between these two images. Here, note that the first measurement and the second measurement are performed at timings at which tissues (a blood vessel, a lesion part, etc.) to be measured are mutually different in shape from each other.
According to this embodiment of the present invention as described above, a comparison between frames becomes able to be correctly made, so it is possible to extract only those parts which are changing by acquiring the difference between the frames. That is, in the conventional ultrasonic diagnostic apparatus which acquires an image in line sequential order, it is impossible to make a comparison between frames in a correct manner as shown in the aforementioned problems, but in contrast, with the construction of this embodiment, photoacoustic tomography is used upon imaging. Since this method performs imaging by irradiating light without scanning in a line sequential fashion unlike the conventional ultrasonic diagnostic imaging apparatus, it becomes possible to acquire simultaneous information on different places inside the living body. As a result, it becomes possible from the difference images between the frames to grasp how the state of the tissue changes in a correct manner. By making use of this, it will become possible to take a vivid or clear image of a moving tissue such as for example a blood vessel figure, a lesion part that has changed over time, etc.
A first example of the present invention will be described below. The first example is an example in which image reconstruction is carried out from ultrasonic waves detected in the entire circumference of the sample. Here, reference will be made to a simulation that was performed under the condition that a situation at the time of actual measurement is simulated while using drawings. Although a three-dimensional space has to be taken into consideration in actual measurement or diagnosis, computations on a two-dimensional plane are shown in this computation process.
It is known that, as for an acoustic wave generated due to an photoacoustic effect, an acoustic wave emitted from an absorber, if its shape is spherical, can be analytically calculated through calculation, and that the acoustic wave thus calculated is in match with that obtained through experiments. In a calculation limited on a two-dimensional plane, if the shape of an absorber is circular, an acoustic wave emitted from the absorber can be calculated in a similar manner. That is, it is known that the profile of the ultrasonic wave exhibits a change of “N” shape in which the change of its amplitude once rises from zero to a positive pressure, then falls to a negative pressure, and thereafter returns to zero. It is also known that the width of this N shape is proportional to the size of the light absorber, and that the height (amplitude) of the N shape is proportional to the product of an amount of light arrived and an optical absorption coefficient, as shown by the above-mentioned formula (1).
The acoustic waves generated in the light absorbers spread in an isotropic manner, propagate through the interior of the living body, and arrive at the ultrasonic wave receiving elements. The arrival acoustic waves, which have been measured at a plurality of points, can be imaged by using a well-known analysis method. In this example, image reconstruction was performed using a well-known circular back-projection method (CBP method), and a first image was obtained from acoustic wave signals detected at sixty points in the sample of
Even if the image of
Accordingly, in this example, a difference between the image of
A second example of the present invention will be described below. The second example is an example in which image reconstruction is carried out from ultrasonic waves detected at one side of a sample. That is, in the first example, the ultrasonic wave receiving elements are arranged so as to surround the two light absorbers, but in contrast, in the second example, ultrasonic wave receiving elements are arranged only at one side with respect to a light absorber. Here, reference will be made to a simulation that was performed under the condition that a situation at the time of actual measurement is simulated while using drawings. Although a three-dimensional space has to be taken into consideration in actual measurement or diagnosis, computations on a two-dimensional plane are shown in this computation process.
In such arrangements, image reconstruction was carried out by using the CBP method, as in the first example. The first image, which is the result reconstructed from ultrasonic signals detected in the sample arrangement of
Thus, in both
In contrast to these,
A third example of the present invention will be described below. In the first and second examples, those parts which have changed in their state are extracted by generating a difference image between image data after performing image reconstruction. That is, difference calculation is performed in real space. In contrast to this, in the third example, a difference image is generated from data that is obtained by transforming acoustic wave signals into corresponding signals in a frequency domain, taking a difference between signal data in the frequency domain, and performing inverse transformation of the difference to a space domain.
Hereafter, the result of a simulation of the third example will be explained by using
A blood vessel can be imaged by using the methods according to the above-mentioned examples. For example, the method of the first example can be used in cases where an ultrasonic probe 121 is arranged around a finger for measuring ultrasonic waves, as shown in
Here, note that, in either of the above cases, it is desirable to use a matching agent or the like which is arranged between the ultrasonic probe and the living body for taking the matching of acoustic impedance therebetween.
The blood vessel is momentarily changing its size in accordance with the flow of blood therein. That is, the blood vessel contracts and expands its size to cause blood to flow therein in synchronization with the heart beat. Thus, the biological information imaging apparatus monitors the heart beat, and performs light irradiation in synchronization with the heart beat, whereby it becomes possible for the apparatus to acquire the state at that timing as an image. Specifically, when blood pressure is substantially systolic blood pressure (first state), a first measurement is carried out, and when the blood pressure is substantially diastolic blood pressure (second state), a second measurement is carried out. Thereby, an image of the blood vessel in a high blood pressure state and an image of the blood vessel in a low blood pressure state are obtained. Then, by performing a comparison between these images, it becomes possible to represent only those parts in which the diameter of the blood vessel has changed following the blood pressure change.
Now, in cases where a part in which the heart beat is measured and a part which is to be diagnosed are separated from each other, a synchronized image can not necessarily be obtained. In this case, if a measurement is made in a rest state as the first state and a measurement is then made in a high heart rate as the second state, a blood vessel figure can be imaged in a correct manner. That is, the timing at which the size of a blood vessel in the first state coincides with that in the second state in which the heart rate is different from that in the first state is very rare, and hence, if a difference is taken between such states, it will become possible to image the blood vessel. In cases where such a difference can not be taken well, the probability at which the blood vessel can be imaged well will be increased if a measurement is performed again in a high heart rate state.
Here, note that, in order to create these two states, a state in which the measurement is carried out after one takes a rest such as sleeping on a bed or sitting on a chair is made the first state, and a state in which the measurement is carried out after one is put into a high heart rate state by doing exercise spontaneously or by being given a stimulus or stress from the outside is made the second state. Thus, the blood vessel can be imaged from a difference image between the two states.
Alternatively, a device with a function to cause such a high heart rate state (e.g., a stress applying function) may also be built into the above-mentioned biological information imaging apparatus. In addition, the biological information imaging apparatus may have a function to measure the heart rate, so that it can automatically determine whether one is in a high heart rate state, and decide the timing at which a second measurement is made.
In addition, a compare may be made between different blood pressure states. For example, an early morning hypotensive state is assumed to be the first state, and a measurement at this time is assumed as the first measurement. Subsequently, a high blood pressure state in an active time zone is assumed to be the second state, and a measurement at this time is assumed as the second measurement. By acquiring a difference between these two states, it becomes possible to image a blood vessel. Also, the biological information imaging apparatus may have a function to measure the blood pressure, so that it can automatically determine whether one is in a low blood pressure state or in a high blood pressure state.
Moreover, it is known that a difference between tissues can be observed by using two kinds of light of different wave lengths. For example, since the degree of oxidation of hemoglobin changes with the wave length of light, it is known that an artery and a vein can be imaged by performing measurements while changing the wave length of light.
When imaging the above-mentioned blood vessels, the change in the size of the blood vessel in the artery is large, whereas in the vein, the size of its blood vessel is little changed even due to the flow of blood therein. That is, it can be considered that the vein is a reference tissue of which the state does not change. By making use of this, at the time of acquiring a difference between the two states, it becomes possible to obtain a difference image with substantially no occurrence of a positional shift by performing registration in such a manner that the reference tissues, i.e., veins, of which the states do not change, are in register with each other.
In order to perform imaging according to a vein(s), imaging in the first state and imaging in the second state are carried out in a wavelength band (700 through 800 nm) in which a lot of hemoglobin uncombined with oxygen can be absorbed. Subsequently, imaging in the first state and imaging in the second state are carried out in a wavelength band (800 through 1,000 nm) in which a lot of oxygenated hemoglobin can be absorbed. The reconstruction of images thus obtained is carried out in such a manner that a reconstructed image in the wavelength band in which a lot of hemoglobin uncombined with oxygen can be absorbed is taken to be a venous image, and registration is performed in such a manner that a venous image in the first state and a venous image in the second state are in match with each other. Thereafter, in that position thus registered, a difference between the reconstructed images at the wave length at which the oxygenated hemoglobin is absorbed is acquired.
It becomes possible to image the artery in a more accurate manner by obtaining a difference image with such a procedure.
While the present invention has been described with reference to exemplary embodiments, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all such modifications and equivalent structures and functions.
Number | Date | Country | Kind |
---|---|---|---|
2008-211993 | Aug 2008 | JP | national |
This application is a continuation of application Ser. No. 14/481,984, filed Sep. 10, 2014, which is a divisional of application Ser. No. 12/544,834, filed Aug. 20, 2009, claims benefit of those filing dates under 35 U.S.C. § 120, and claims the benefit of Japanese Patent Application No. 2008-211993, filed on Aug. 20, 2008; the entire contents of all mentioned earlier applications are incorporated herein by reference.
Number | Date | Country | |
---|---|---|---|
Parent | 12544834 | Aug 2009 | US |
Child | 14481984 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 14481984 | Sep 2014 | US |
Child | 15912849 | US |