Patent Grant
9380944
1. Field of the Invention
The present invention relates to a photoacoustic image generation apparatus and an acoustic wave unit, and more particularly, to a photoacoustic image generation apparatus, which when a test object is irradiated with a laser beam having a plurality of wavelengths, generates a photoacoustic image on the basis of photoacoustic signals detected with respect to the respective wavelengths, and an acoustic wave unit.
2. Description of the Related Art
Hitherto, for example, as disclosed in JP2005-21380A and A High-Speed Photoacoustic Tomography System based on a Commercial Ultrasound and a Custom Transducer Array, Xueding Wang, Jonathan Cannata, Derek DeBusschere, Changhong Hu, J. Brian Fowlkes, and Paul Carson, Proc. SPIE Vol. 7564, 756424 (Feb. 23, 2010), a photoacoustic image forming apparatus that forms an image of the inside of a living body using a photoacoustic effect has been known. In the photoacoustic image forming apparatus, a living body is irradiated with pulsed light such as a pulse laser beam. Body tissues absorbing energy of the pulsed light expand in volume inside the living body irradiated with the pulsed light, and thus acoustic waves are generated. It is possible to detect the acoustic waves using an ultrasonic probe or the like, and to form a visible image of the inside of the living body on the basis of the detected signal (photoacoustic signal). In a photoacoustic image forming method, acoustic waves are generated in a specific light absorber, and thus it is possible to form an image of specific tissues in the living body, for example, blood vessels.
Incidentally, many of body tissues have an optical absorption property varying depending on a wavelength of light, and generally, the optical absorption property is unique for each tissue. For example,
Herein, with regard to a variable wavelength laser, JP2009-231483A discloses that a laser beam having a desired wavelength is obtained by disposing an etalon or a birefringent filter as a wavelength selection element within an optical resonator and adjusting the rotation angle thereof. In addition, JP2000-105464A discloses that an etalon as wavelength selection means is disposed within an optical resonator and that the etalon is scanned at a constant speed. JP2000-105464A discloses that laser oscillation is performed only when a transmission wavelength of the etalon is consistent with longitudinal mode oscillation of a laser beam and that the oscillation of the laser beam is performed in a pulsed manner when a scanning speed of the etalon is increased.
In JP2009-231483A, in order to switch and emit a laser beam having a plurality of wavelengths, it is necessary to adjust a rotation angle of the etalon or the birefringent filter at every laser emission. In photoacoustic imaging, for example, it is considered that when a test object is irradiated with pulse laser beams having a first wavelength and a second wavelength, the wavelength selection element is adjusted to irradiate the test object with the laser beam having the first wavelength, the detection of all photoacoustic signals of the laser beam having the first wavelength is terminated, and then the wavelength selection element is adjusted so as to emit the laser beam having the second wavelength, and the test object is irradiated with the laser beam having the second wavelength. In the photoacoustic imaging, an object having a movement such as a human is often selected. Therefore, when an object moves during switching from the first wavelength to the second wavelength, mismatching may occur between a photoacoustic signal at the time of irradiation with the laser beam having the first wavelength and a photoacoustic signal at the time of irradiation with the laser beam having the second wavelength.
In the photoacoustic imaging, in terms of the prevention of the above-mentioned mismatching, for example, it is considered that the irradiation with a laser beam may be performed by switching the first wavelength and the second wavelength for each pulse. In other words, for example, it is considered that the irradiation with the laser beam may be repeatedly performed in a predetermined wavelength sequence including the first wavelength and the second wavelength in this order. JP2000-105464A discloses a laser device that changes a wavelength of a laser beam with which the irradiation is performed for each pulse. However, in JP2000-105464A, since a laser is oscillated only when the transmission wavelength of the etalon is consistent with the longitudinal mode oscillation of the laser beam, a laser beam having only a specific wavelength sequence can be obtained, and a laser beam having any wavelength sequence cannot be obtained.
The invention is contrived in view of such situations, and an object thereof is to provide an acoustic wave unit capable of emitting a pulse laser beam in a desired wavelength sequence from a wavelength variable laser light source, and a photoacoustic image generation apparatus including the acoustic wave unit.
In order to achieve the above-described object, the invention provides a photoacoustic image generation apparatus including: a laser source unit that sequentially emits a plurality of pulse laser beams in a predetermined wavelength sequence having at least two different wavelengths, the laser source unit including a laser rod, an excitation light source that irradiates the laser rod with excitation light, an optical resonator that has a pair of mirrors facing each other with the laser rod interposed therebetween, a Q switch which is inserted into the optical resonator, and a birefringent filter which is inserted into the optical resonator and changes an oscillation wavelength of the optical resonator in association with rotational displacement of the birefringent filter; and an acoustic wave unit that generates a photoacoustic image, the acoustic wave unit including detection unit that detects a photoacoustic signal generated within an object when the object is irradiated with the pulse laser beam having each wavelength included in the predetermined wavelength sequence and generates pieces of photoacoustic data corresponding to the respective wavelengths, intensity ratio extraction unit that extracts a magnitude relation between relative signal intensities of the pieces of photoacoustic data corresponding to the respective wavelengths, photoacoustic image construction unit that generates the photoacoustic image on the basis of the extracted magnitude relation, and a trigger control circuit that causes the laser rod to be irradiated with excitation light from the excitation light source while rotating the birefringent filter at a predetermined rotation speed depending on the number of wavelengths included in the wavelength sequence, and after the irradiation with the excitation light, turns on the Q switch at a timing when a rotational displacement position of the birefringent filter is set to a position corresponding to the wavelength of the pulse laser beam to be emitted to cause the pulse laser beam to be emitted.
In the invention, the predetermined rotation speed may be determined on the basis of a change characteristic of the oscillation wavelength with respect to the rotational displacement position in the birefringent filter, the number of wavelengths included in the wavelength sequence, and the number of times of emission of the pulse laser beam per unit time.
When the number of times of a free spectral range repeated during one rotation is set to k[times/rotation], the number of wavelengths included in the wavelength sequence is set to n[pieces], and the number of times of emission of the pulse laser beam per unit time is set to m[times/second], the predetermined rotation speed of the birefringent filter may be determined as a value calculated by a relation of v=m/(k×n)[rotations/second].
In the invention, the trigger control circuit may continuously rotate the birefringent filter in a predetermined direction at the predetermined rotation speed.
The trigger control circuit may determine a timing at which the excitation light is irradiated and a timing at which the Q switch is turned on, on the basis of birefringent filter state information indicating the rotational displacement position of the birefringent filter.
When the birefringent filter state information is set to information indicating a position obtained by subtracting the amount of rotational displacement of the birefringent filter during a period of time required for the excitation of the laser rod from a position of the birefringent filter which corresponds to the wavelength of the pulse laser beam to be emitted, the trigger control circuit may cause the laser rod to be irradiated with excitation light.
The trigger control circuit may rotate the birefringent filter so that the amount of change in the birefringent filter state information during a predetermined period of time is set to the amount of change depending on the predetermined rotation speed.
Moreover, the laser source unit may further include driving unit that rotates the birefringent filter, rotational displacement detection unit that detects the rotational displacement of the birefringent filter, and a rotation control unit that controls the driving unit so that the amount of rotational displacement of the birefringent filter which is detected by the rotational displacement detection unit during a predetermined period of time is set to an amount depending on the predetermined rotation speed.
The rotational displacement detection unit may output a birefringent filter state signal indicating the rotational displacement position of the birefringent filter to the acoustic wave unit.
The trigger control circuit may output a birefringent filter control signal for controlling rotation of the birefringent filter to the driving unit.
The trigger control circuit may output a flash lamp trigger signal for controlling emission of the excitation light to the excitation light source on the basis of the birefringent filter state signal, and may cause the laser rod to be irradiated with excitation light from the excitation light source.
The trigger control circuit may output a Q switch trigger signal for turning on the Q switch when the birefringent filter state signal is set to signal indicating a position of the birefringent filter which transmits a wavelength of a pulse laser beam to be emitted, and may cause the pulse laser beam to be emitted from the laser source unit.
The acoustic wave unit may include an AD converting unit that performs sampling of the photoacoustic signal with a predetermined sampling period, and the trigger control circuit may output a sampling trigger signal indicating a sampling timing of the AD converting unit to the AD converting unit, in synchronization with an output timing of the Q switch trigger signal.
The acoustic wave unit may further include intensity information extraction unit that generates intensity information indicating signal intensity on the basis of the pieces of photoacoustic data corresponding to the respective wavelengths. The photoacoustic image construction unit may determine a gradation value of each pixel of the photoacoustic image on the basis of the intensity information and may determine a display color of each pixel on the basis of the extracted magnitude relation.
The predetermined wavelength sequence may include a first wavelength and a second wavelength. The acoustic wave unit may further include complex number creation unit that generates complex number data in which one of first photoacoustic data corresponding to a photoacoustic signal, detected when irradiation with the pulse laser beam having the first wavelength is performed, and second photoacoustic data corresponding to a photoacoustic signal, detected when irradiation with the pulse laser beam having the second wavelength is performed, is set to a real part and the other one is set to an imaginary part, and photoacoustic image reconstruction unit that generates a reconstructed image from the complex number data using a Fourier transform method. The intensity ratio extraction unit may extract phase information as the magnitude relation from the reconstructed image, and the intensity information extraction unit may extract the intensity information from the reconstructed image.
The detection unit may further detect reflected acoustic waves with respect to acoustic waves transmitted to the object to generate reflected acoustic wave data, and the acoustic wave unit may further include acoustic wave image generation unit that generates an acoustic wave image on the basis of the reflected acoustic wave data.
The invention also provides an acoustic wave unit including: detection unit that detects a photoacoustic signal generated within an object when the object is irradiated with a pulse laser beam having each wavelength included in a predetermined wavelength sequence including at least two different wavelengths, and generates pieces of photoacoustic data corresponding to the respective wavelengths; intensity ratio extraction unit that extracts a magnitude relation between relative signal intensities of the pieces of photoacoustic data corresponding to the respective wavelengths; photoacoustic image construction unit that generates a photoacoustic image on the basis of the extracted magnitude relation; and a trigger control circuit that causes the laser rod to be irradiated with excitation light from an excitation light source while rotating a birefringent filter which is inserted into an optical resonator that includes a pair of mirrors facing each other with a laser rod interposed therebetween and changes an oscillation wavelength of the optical resonator in association with rotational displacement of the birefringent filter, at a predetermined rotation speed depending on the number of wavelengths included in the wavelength sequence, and after the irradiation with the excitation light, turns on a Q switch inserted into the optical resonator at a timing when a rotational displacement position of the birefringent filter is set to a position corresponding to the wavelength of the pulse laser beam to be emitted, to cause the pulse laser beam to be emitted.
In a photoacoustic image generation apparatus and an acoustic wave unit of the invention, a birefringent filter which is inserted into an optical resonator and changes an oscillation wavelength in association with rotational displacement of the birefringent filter is rotated at a rotation speed depending on the number of wavelength sequences of a pulse laser beam to be emitted from a laser source unit, and a Q switch which is inserted into the optical resonator is turned on at a timing when a rotational displacement position of the birefringent filter is set to a position corresponding to the wavelength of the pulse laser beam to be emitted. As the rotation speed of the birefringent filter increases, the speed of the wavelength switching can be increased. On the contrary, as the rotation speed thereof decreases, the number of selectable oscillation wavelengths can be increased. The rotation speed of the birefringent filter of the present invention is controlled depending on the number of wavelengths included in a wavelength sequence. Based on such a configuration, the wavelengths of the pulse laser beam to be emitted from the laser source unit can be controlled to have any wavelength sequence by the acoustic wave unit. In addition, the acoustic wave unit of the present invention determines an emission timing of the pulse laser beam, and thus it is not necessary to acquire a signal such as a synchronization signal indicating laser emission from the laser source unit, in the start of sampling of a photoacoustic signal.
Hereinafter, embodiments of the invention will be described in detail with reference to the accompanying drawings. Meanwhile, in examples of the invention, ultrasonic waves are used as acoustic waves, but the acoustic waves may be acoustic waves having an audible frequency by selecting an appropriate frequency according to an object to be tested or measurement conditions.
For example, a wavelength of approximately 750 nm is considered as the first wavelength (center wavelength), and a wavelength of approximately 800 nm is considered as the second wavelength. Referring to
The pulse laser beam emitted from the laser source unit 13 is guided to a probe 11 using light guiding means such as an optical fiber, and is irradiated toward a test object from the probe 11. An irradiation position of the pulse laser beam is not particularly limited, and the pulse laser beam may be irradiated from any place other than the probe 11. Ultrasonic waves (acoustic waves) are generated within the test object by a light absorber absorbing energy of the irradiated pulse laser beam. The probe 11 includes an ultrasonic wave detector. The probe 11 includes, for example, a plurality of ultrasonic wave detector elements (ultrasonic wave vibrators) which are arranged one-dimensionally, and the acoustic waves (photoacoustic signal) from the inside of the test object are detected by the ultrasonic wave vibrators that are arranged one-dimensionally.
The ultrasonic wave unit 12 includes a reception circuit 21, AD conversion unit 22, a reception memory 23, complex number creation unit 24, photoacoustic image reconstruction unit 25, phase information extraction unit 26, intensity information extraction unit 27, detection and logarithmic transformation unit 28, photoacoustic image construction unit 29, a trigger control circuit 30, and control unit 31. The reception circuit 21 receives a photoacoustic signal detected by the probe 11. The AD conversion unit 22, which is detection unit, samples the photoacoustic signal received by the reception circuit 21 and generates photoacoustic data which is digital data. The AD conversion unit 22 samples the photoacoustic signal with a predetermined sampling period in synchronization with an AD clock signal.
The AD conversion unit 22 stores photoacoustic data in the reception memory 23. The AD conversion unit 22 stores, in the reception memory 23, photoacoustic data corresponding to the respective wavelengths of the pulse laser beam emitted from the laser source unit 13. In other words, the AD conversion unit 22 stores, in the reception memory 23, first photoacoustic data obtained by sampling a photoacoustic signal detected by the probe 11 when a test object is irradiated with a pulse laser beam having a first wavelength and second photoacoustic data obtained by sampling a photoacoustic signal detected by the probe 11 when the test object is irradiated with a second pulse laser beam.
The complex number creation unit 24 reads out the first photoacoustic data and the second photoacoustic data from the reception memory 23, and generates complex number data in which any one of the first photoacoustic data and the second photoacoustic data is set to a real part and the other one is set to an imaginary part. Hereinafter, a description will be given on the assumption that the complex number creation unit 24 generates the complex number data in which the first photoacoustic data is set to a real part and the second photoacoustic data is set to as an imaginary part.
The photoacoustic image reconstruction unit 25 inputs the complex number data from the complex number creation unit 24. The photoacoustic image reconstruction unit 25 performs image reconstruction from the input complex number data using a Fourier transform method (FTA method). A well-known method of the related art which is disclosed in, for example, a document “Photoacoustic Image Reconstruction-A Quantitative Analysis” Jonathan I. Sperl et al. SPIE-OSA, Vol. 6631 663103 can be applied to the image reconstruction using the Fourier transform method. The photoacoustic image reconstruction unit 25 inputs Fourier transformed data indicating the reconstructed image to the phase information extraction unit 26 and the intensity information extraction unit 27.
The phase information extraction unit 26 extracts a magnitude relation between relative signal intensities of pieces of photoacoustic data corresponding to the respective wavelengths. In this embodiment, the phase information extraction unit 26 sets the reconstructed image reconstructed by the photoacoustic image reconstruction unit 25 to input data. In addition, the phase information extraction unit generates, when the real part and the imaginary part are compared with each other, phase information indicating how relatively large either of the two parts is, using the input data which is the complex number data. For example, when the complex number data is expressed by X+iY, the phase information extraction unit 26 generates the relation of θ=tan−1(Y/X) as phase information. Meanwhile, when the relation of X=0 is satisfied, the relation of θ=90° is established. When first photoacoustic data (X) constituting the real part is equal to second photoacoustic data (Y) constituting the imaginary part, the phase information satisfies the relation of θ=45°. As the first photoacoustic data becomes relatively larger, the phase information becomes closer to the relation of θ=0°. As the second photoacoustic data becomes larger, the phase information becomes closer to the relation of θ=90°.
The intensity information extraction unit 27 generates intensity information indicating signal intensity on the basis of the pieces of photoacoustic data corresponding to the respective wavelengths. In this embodiment, the intensity information extraction unit 27 sets the reconstructed image reconstructed by the photoacoustic image reconstruction unit 25 to input data, and generates the intensity information from the input data which is complex number data. For example, when the complex number data is expressed by X+iY, the intensity information extraction unit 27 extracts (X2+Y2)1/2 as intensity information. The detection and logarithmic transformation unit 28 generates an envelope of data indicating the intensity information extracted by the intensity information extraction unit 27, and then widens a dynamic range by performing logarithmic transformation on the envelope.
The photoacoustic image construction unit 29 inputs the phase information from the phase information extraction unit 26, and inputs the intensity information after the detection and logarithmic transformation process from the detection and logarithmic transformation unit 28. The photoacoustic image construction unit 29 generates a photoacoustic image which is a distribution image of a light absorber, on the basis of the input phase information and intensity information. For example, the photoacoustic image construction unit 29 determines luminance (gradation value) of each pixel in the distribution image of the light absorber, on the basis of the input intensity information. In addition, for example, the photoacoustic image construction unit 29 determines color of each pixel (display color) in the distribution image of the light absorber, on the basis of the phase information. The photoacoustic image construction unit 29 determines color of each pixel on the basis of the input phase information, for example, using the range of phases 0° to 90° in a color map associated with a predetermined color.
Here, since the range of phases 0° to 45° is a range in which the first photoacoustic data is larger than the second photoacoustic data, a generation source of a photoacoustic signal is considered to be a vein through which blood flows, the blood mainly containing deoxygenated hemoglobin in which the amount of absorption of a wavelength of 756 nm is greater than that of a wavelength of 798 nm. On the other hand, since the range of phases 45° to 90° is a range in which the first photoacoustic data is smaller than the second photoacoustic data, the generation source of the photoacoustic signal is considered to be an artery through which blood flows, the blood mainly containing oxygenated hemoglobin in which the amount of absorption of a wavelength of 756 nm is less than that of a wavelength of 798 nm.
Consequently, as the color map, a color map is used in which color gradually changes so as to become blue at the phase of 0° and to become colorless (white) as the phase approaches 45° and in which color gradually changes so as to become red at the phase of 90° and to become white as the phase approaches 45°. In this case, in the photoacoustic image, a portion corresponding to the artery can be expressed by red, and a portion corresponding to the vein can be expressed by blue. Only color coding between the portion corresponding to the artery and the portion corresponding to the vein may be performed on the basis of the phase information by maintaining a constant gradation value, without using the intensity information. The image display unit 14 displays the photoacoustic image generated by the photoacoustic image construction unit 29 on a display screen.
Subsequently, a configuration of the laser source unit 13 will be described in detail.
The mirrors 53 and 54 face each other with the laser rod 51 interposed therebetween, and an optical resonator is constituted by the mirrors 53 and 54. The mirror 54 is assumed to be the output side. The Q switch 55 and the birefringent filter 56 are inserted into the optical resonator. An insertion loss within the optical resonator rapidly changes from a high loss (low Q) to a low loss (high Q) by the Q switch 55, and thus a pulse laser beam can be obtained. The birefringent filter 56 changes a transmission wavelength in association with rotational displacement and changes an oscillation wavelength of the optical resonator. The driving unit 57 rotates the birefringent filter 56. The rotational displacement detection unit 58 detects the rotational displacement of the birefringent filter 56. The rotational displacement detection unit 58 outputs BRF state information indicating the rotational displacement position of the birefringent filter 56 to the ultrasonic wave unit 12.
Referring back to
The trigger control circuit 30 outputs a BRF control signal for controlling the rotation of the birefringent filter 56. The driving unit 57 of the laser source unit 13 rotates the birefringent filter 56 in response to the BRF control signal. For example, the trigger control circuit 30 rotates the birefringent filter so that the amount of change in BRF state information during a predetermined period of time is set to the amount of change depending on a predetermined rotation speed, on the basis of the BRF control signal.
In addition to the above description, the trigger control circuit 30 outputs a flash lamp trigger signal for controlling the emission of the flash lamp 52 to the laser source unit 13, and causes the laser rod 51 to be irradiated with excitation light from the flash lamp 52. The trigger control circuit 30 outputs the flash lamp trigger signal on the basis of a BRF state signal. For example, when the BRF state information is set to information indicating the position obtained by subtracting the amount of rotational displacement of the birefringent filter 56 during a period of time required for the excitation of the laser rod 51 from the position of the birefringent filter 56 which corresponds to the wavelength of the pulse laser beam to be emitted, the trigger control circuit 30 outputs the flash lamp trigger signal and causes the laser rod 51 to be irradiated with excitation light.
After the irradiation with the excitation light, the trigger control circuit 30 outputs a Q switch trigger signal to the Q switch 55 at a timing when the rotational displacement position of the birefringent filter 56 is set to the position corresponding to the wavelength of the pulse laser beam to be emitted. In other words, the trigger control circuit 30 outputs the Q switch trigger signal when the BRF state information is set to information indicating the position of the birefringent filter 56 which transmits the wavelength of the pulse laser beam to be emitted. The Q switch 55 rapidly changes the insertion loss within the optical resonator from a high loss to a low loss (Q switch is turned on) in response to the Q switch trigger signal, and thus the pulse laser beam is emitted from the mirror 54 on the output side.
The trigger control circuit 30 outputs a sampling trigger signal (AD trigger signal) to the AD conversion unit 22 in accordance with a timing of the Q switch trigger signal, that is, the emission timing of the pulse laser beam. The AD conversion unit 22 starts the sampling of a photoacoustic signal on the basis of the sampling trigger signal.
First, the trigger control circuit 30 outputs the flash lamp trigger signal to the flash lamp 52 at time t1 in order to cause a pulse laser beam having a wavelength of 750 nm to be emitted from the laser source unit 13 ((b) of
Subsequently, the trigger control circuit 30 outputs the flash lamp trigger signal to the flash lamp 52 at time t3 in order to cause a pulse laser beam having a wavelength of 800 to be emitted from the laser source unit 13 ((b) of
Here, the time t1 when the flash lamp trigger signal is output is a time obtained by subtracting a time required for the excitation of the laser rod 51 from the time t2 when the rotational displacement position of the birefringent filter 56 is set to the position corresponding to the wavelength of 750 nm. In addition, the time t3 when the flash lamp trigger signal is output is a time obtained by subtracting a time required for the excitation of the laser rod 51 from the time t4 when the rotational displacement position of the birefringent filter 56 is set to the position corresponding to the wavelength of 800 nm. The rotational displacement positions of the birefringent filter 56 which correspond to the time t1 and the time t3 can be obtained from the rotational displacement positions of the birefringent filter 56 which correspond to the wavelengths of 750 nm and 800 nm, the rotation speed of the birefringent filter 56, and time required for the excitation of the laser rod 51.
Hereinafter, similarly, the trigger control circuit 30 outputs the flash lamp trigger signal to the flash lamp 52 at time t5, time t7, time t9, and time t11. In addition, the trigger control circuit outputs the Q switch trigger signal to the Q switch 55 at time t6, time t8, time t10, and time t12, and causes the pulse laser beam having a wavelength depending on the transmission wavelength of the birefringent filter 56 at each time to be emitted. The transmission wavelengths of the birefringent filter 56 at time t6 and time t10 are 750 nm, and the transmission wavelengths of the birefringent filter at time t8 and time t12 are 800 nm, and thus the laser source unit 13 sequentially and repeatedly emits the pulse laser beams having wavelengths of 750 nm and 800 nm in this order.
In the example of
When the photoacoustic signal is ready to be received, the trigger control circuit 30 outputs the flash lamp trigger signal to the laser source unit 13 at a predetermined timing in order to cause the pulse laser beam having a first wavelength (750 nm) constituting a wavelength sequence to be emitted (step A2). The flash lamp 52 of the laser source unit 13 is turned on in response to the flash lamp trigger signal, and thus the laser rod 51 starts to be excited (step A3). The trigger control circuit 30 turns on the flash lamp 52, for example, at a timing calculated back from a timing at which the rotational displacement position of the birefringent filter 56 is set to the position corresponding to the wavelength of 750 nm, on the basis of the BRF state information.
After the flash lamp 52 is turned on, the trigger control circuit 30 turns on the Q switch 55 at a timing when the rotational displacement position of the birefringent filter 56 is set to the position corresponding to the first wavelength (750 nm) constituting the wavelength sequence, on the basis of the BRF state information (step A4). The laser source unit 13 emits the pulse laser beam having a wavelength of 750 nm by the Q switch 55 being turned on.
The pulse laser beam having a wavelength of 750 nm which is emitted from the laser source unit 13 is guided to, for example, the probe 11, and a first partial region of the test object is irradiated with the pulse laser beam from the probe 11. A light absorber absorbs energy of the irradiated pulse laser beam within the test object, and thus a photoacoustic signal is generated. The probe 11 detects the photoacoustic signal generated within the test object. The photoacoustic signal detected by the probe 11 is received by the reception circuit 21.
The trigger control circuit 30 outputs the sampling trigger signal to the AD conversion unit 22 in accordance with a timing at which the Q switch trigger signal is output. The AD conversion unit 22 samples the photoacoustic signal received by the reception circuit 21 with a predetermined sampling period (step A5). The photoacoustic signal sampled by the AD conversion unit 22 is stored as first photoacoustic data in the reception memory 23.
The control unit 31 determines whether a remaining wavelength to be emitted is present or not, in other words, whether the pulse laser beams of all the predetermined wavelengths constituting the wavelength sequence have been emitted or not (step A6). When a remaining wavelength is present, the process returns to step A2 in order to emit the pulse laser beam having the next wavelength, and the flash lamp trigger signal is output to the laser source unit 13 from the trigger control circuit 30. In step A3, the flash lamp 52 is turned on in response to the flash lamp trigger signal, and in step A4, the trigger control circuit 30 turns on the Q switch 55 at a timing when the birefringent filter 56 is set to be at the rotational displacement position corresponding to the second wavelength (800 nm) constituting the wavelength sequence, to cause the pulse laser beam to be emitted.
The pulse laser beam having a wavelength of 800 nm which is emitted from the laser source unit 13 is guided to, for example, the probe 11, and the first partial region of the test object is irradiated with the pulse laser beam from the probe 11. The probe 11 detects a photoacoustic signal generated by the light absorber within the test object absorbing the pulse laser beam having a wavelength of 800 nm. The trigger control circuit 30 outputs the sampling trigger signal to the AD conversion unit 22 in accordance with the output of the Q switch trigger signal, and the AD conversion unit 22 samples the photoacoustic signal in step A5. The photoacoustic signal sampled by the AD conversion unit 22 is stored as second photoacoustic data in the reception memory 23. The photoacoustic image generation apparatus 10 performs step A1 to step A5 on the wavelengths constituting the wavelength sequence and irradiates the test object with the pulse laser beam having the wavelengths constituting the wavelength sequence, thereby detecting a photoacoustic signal from the test object.
When the control unit 31 determines in step A6 that a remaining wavelength is not present, the control unit determines whether all the partial regions have been selected (step A7). When the partial region to be selected remains, the process returns to step A2. The photoacoustic image generation apparatus 10 performs step A2 to step A6 on each partial region, sequentially irradiates each partial region with pulse laser beams having the wavelengths (750 nm and 800 nm) constituting the wavelength sequence, and stores the first photoacoustic data and the second photoacoustic data which correspond to each partial region, in the reception memory 23. When the irradiation with the pulse laser beam and the detection of the photoacoustic signal are performed on all the partial regions, photoacoustic data required to generate a photoacoustic image of one frame is gathered.
When the control unit 31 determines in step A7 that all the partial regions have been selected, the process proceeds to the generation of the photoacoustic image. The complex number creation unit 24 reads out the first photoacoustic data and the second photoacoustic data from the reception memory 23, and generates complex number data in which first photoacoustic image data is set to a real part and second photoacoustic image data is set to an imaginary part (step A8). The photoacoustic image reconstruction unit 25 performs image reconstruction from the complex number data generated in step A8, using a Fourier transform method (FTA method) (step A9).
The phase information extraction unit 26 extracts phase information from the reconstructed complex number data (reconstructed image) (step A10). For example, when the reconstructed complex number data is expressed by X+iY, the phase information extraction unit 26 extracts the relation of θ=tan−1(Y/X) as the phase information (but, when the relation of X=0 is satisfied, the relation of θ=90° is satisfied). The intensity information extraction unit 27 extracts intensity information from the reconstructed complex number data (step A11). For example, when the reconstructed complex number data is expressed by X+iY, the intensity information extraction unit 27 extracts (X2+Y2)1/2 as the intensity information.
The detection and logarithmic transformation unit 28 performs a detection and logarithmic transformation process on the intensity information extracted in step A11. The photoacoustic image construction unit 29 generates a photoacoustic image on the basis of the phase information extracted in step A10 and the performing of the detection and logarithmic transformation process, on the intensity information extracted in step A11 (step A12). For example, the photoacoustic image construction unit 29 generates the photoacoustic image by determining luminance (gradation value) of each pixel in a distribution image of a light absorber on the basis of the intensity information and by determining color of each pixel on the basis of the phase information. The generated photoacoustic image is displayed on the image display unit 14.
Here, the rotation speed of the birefringent filter 56 may be appropriately determined depending on the number of wavelengths included in a wavelength sequence of a pulse laser beam to be emitted. Hereinafter, a case where the wavelength sequence includes six wavelengths (720 nm, 740 nm, 760 nm, 780 nm, 800 nm, and 820 nm) will be described.
As illustrated in
Subsequently, the trigger control circuit 30 outputs the flash lamp trigger signal to the flash lamp 52 at time t23 in order to cause a pulse laser beam having a wavelength of 740 nm to be emitted from the laser source unit 13 ((b) of
Hereinafter, similarly, the trigger control circuit 30 outputs the flash lamp trigger signal to the flash lamp 52 at time t25, time t27, time t29, and time t31. In addition, the trigger control circuit outputs the Q switch trigger signal to the Q switch 55 at time t26, time t28, time t30, and time t32, and causes the pulse laser beam having a wavelength depending on the transmission wavelength of the birefringent filter 56 at each time to be emitted. The transmission wavelengths of the birefringent filter 56 at time t26, time t28, time t30, and time t32 are 760 nm, 780 nm, 800 nm, and 820 nm, respectively, and the laser source unit 13 emits six pulse laser beams having a wavelength increasing by 20 nm for ⅓ seconds in a range between 720 nm and 820 nm.
In an example of
It is preferable that the rotation speed of the birefringent filter 56 be set so that a pulse laser beam having wavelengths constituting a wavelength sequence can be emitted in one FSR. For example, when the number of times of the FSR repeated by the birefringent filter 56 during one rotation is set to k[times/rotation], the number of wavelengths included in the wavelength sequence is set to n[pieces], and the number of times of emission of the pulse laser beam per unit time is set to m[times/second], the rotation speed of the birefringent filter 56 can be set to a value determined by the relation of v=m/(k×n)[rotations/second]. In this case, m pulse lasers can be emitted per second while switching n wavelengths for each FSR (m Hz operation).
In this embodiment, the flash lamp 52 is turned on to excite the laser rod 51 while rotating the birefringent filter 56 at a predetermined rotation speed. After the excitation of the laser rod, the Q switch 55 is turned on at a timing when the rotational displacement position of the birefringent filter 56 is set to a position corresponding to a wavelength of a pulse laser beam to be emitted. As the rotation speed of the birefringent filter 56 decreases, for example, the number of oscillation wavelengths capable of being selected in one FSR of the birefringent filter 56 can be increased. On the other hand, when the number of wavelengths included in the wavelength sequence is two, the speed of the switching of the two wavelengths can be increased by increasing the rotation speed of the birefringent filter 56. In this manner, in this embodiment, it is possible to emit a pulse laser beam in a desired wavelength sequence from the laser source unit 13 by controlling the rotation speed of the birefringent filter 56. In this embodiment, the Q switch trigger signal is output from the ultrasonic wave unit 12, and thus it is not necessary to acquire information such as a synchronization signal indicating a laser emission timing from the laser source unit 13.
In this embodiment, complex number data is generated in which one of the first photoacoustic data and the second photoacoustic data which are obtained at two wavelengths is set to a real part and the other is set to an imaginary part, and a reconstructed image is generated from the complex number data using a Fourier transform method. In this case, it is possible to effectively perform the reconstruction as compared with a case where the first photoacoustic data and the second photoacoustic data are separately reconstructed. A pulse laser beam of a plurality of wavelengths is irradiated, and a photoacoustic signal (photoacoustic data) at the time of the irradiation with a pulse laser beam having each wavelength is used, and thus it is possible to perform functional imaging using optical absorption properties of the respective light absorbers being different from each other depending on wavelengths.
In addition, in this embodiment, for example, when a light irradiation region is divided into three partial regions, a first partial region is sequentially irradiated with a pulse laser beam having a first wavelength and a pulse laser beam having a second wavelength, and a second partial region is sequentially irradiated with the pulse laser beam having the first wavelength and the pulse laser beam having the second wavelength, and then a third partial region is sequentially irradiated with the pulse laser beam having the first wavelength and the pulse laser beam having the second wavelength. In this embodiment, any partial region is continuously irradiated with the pulse laser beam having the first wavelength and the pulse laser beam having the second wavelength, and then the irradiation moves to the next partial region. In this case, it is possible to shorten the time from the irradiation with the pulse laser beam having the first wavelength to the irradiation with the second wavelength at the same position, as compared with a case where the three partial regions are irradiated with the pulse laser beam having the first wavelength and are then irradiated with the pulse laser beam having the second wavelength. It is possible to suppress mismatching between the first photoacoustic data and the second photoacoustic data by shortening the time between the irradiation with the pulse laser beam having the first wavelength and the irradiation with the pulse laser beam having the second wavelength.
Subsequently, a second embodiment of the invention will be described.
In this embodiment, a probe 11 outputs (transmits) ultrasonic waves to a test object and detects (receives) reflected ultrasonic waves from the test object with respect to the transmitted ultrasonic waves, in addition to the detection of a photoacoustic signal. A trigger control circuit 30 transmits an ultrasonic wave transmission trigger signal for instructing the transmission of ultrasonic waves to the transmission control circuit 37 at the time of the generation of an ultrasonic image. When the transmission control circuit 37 receives the trigger signal, the transmission control circuit causes ultrasonic waves to be transmitted from the probe 11. The probe 11 detects reflected ultrasonic waves from the test object after the transmission of the ultrasonic waves.
The reflected ultrasonic waves detected by the probe 11 are input to AD conversion unit 22 through a reception circuit 21. The trigger control circuit 30 transmits a sampling trigger signal to the AD conversion unit 22 in accordance with the transmission timing of the ultrasonic waves, and starts to sample the reflected ultrasonic waves. The AD conversion unit 22 stores sampling data of the reflected ultrasonic waves (reflected ultrasonic data) in the reception memory 23.
The data separation unit 32 separates the reflected ultrasonic data stored in the reception memory 23 and first and second photoacoustic data from each other. The data separation unit 32 transmits the reflected ultrasonic data to the ultrasonic image reconstruction unit 33, and transmits the first and second photoacoustic data to complex number creation unit 24. The generation of the photoacoustic image on the basis of the first and second photoacoustic data is the same as that in the first embodiment. The data separation unit 32 inputs sampling data of the separated reflected ultrasonic waves to the ultrasonic image reconstruction unit 33.
The ultrasonic image reconstruction unit 33 generates pieces of data of lines of the ultrasonic image on the basis of reflected ultrasonic waves (sampling data thereof) which are detected by a plurality of ultrasonic vibrators of the probe 11. For example, the ultrasonic image reconstruction unit 33 adds data from 64 ultrasonic vibrators of the probe 11 on the basis of a delay time depending on the position of the ultrasonic vibrator to generate data for one line (delay addition method).
The detection and logarithmic transformation unit 34 obtains an envelope of the pieces of data of the lines which are output by the ultrasonic image reconstruction unit 33, and performs logarithmic transformation on the obtained envelope. The ultrasonic image construction unit 35 generates an ultrasonic image on the basis of the data of the lines on which the logarithmic transformation is performed. The ultrasonic image reconstruction unit 33, the detection and logarithmic transformation unit 34, and the ultrasonic image construction unit 35 constitute ultrasonic image generation unit that generates an ultrasonic image on the basis of reflected ultrasonic waves.
The image synthesis unit 36 synthesizes the photoacoustic image and the ultrasonic image. For example, the image synthesis unit 36 performs image synthesis by superimposing the photoacoustic image and the ultrasonic image on each other. At this time, it is preferable that the image synthesis unit 36 perform positioning so that corresponding points of the photoacoustic image and the ultrasonic image are set to be at the same position. The synthesized image is displayed on image display unit 14. It is also possible to display the photoacoustic image and the ultrasonic image on the image display unit 14 side by side without performing image synthesis, or to switch and display the photoacoustic image and the ultrasonic image.
When a photoacoustic signal is ready to be received, the trigger control circuit 30 outputs a flash lamp trigger signal in order to cause a pulse laser beam having a first wavelength (750 nm) constituting a wavelength sequence to be emitted (step B2). A flash lamp 52 is turned on in response to the flash lamp trigger signal, and thus a laser rod 51 starts to be excited (step B3).
After the flash lamp 52 is turned on, the trigger control circuit 30 turns on a Q switch 55 at a timing when a rotational displacement position of the birefringent filter 56 is set to the position corresponding to the first wavelength (750 nm) constituting the wavelength sequence, on the basis of the BRF control signal (step B4). The laser source unit 13 emits a pulse laser beam having a wavelength of 750 nm by the Q switch 55 being turned on.
The pulse laser beam having a wavelength of 750 nm which is emitted from the laser source unit 13 is guided to, for example, the probe 11, and a first partial region of the test object is irradiated with the pulse laser beam from the probe 11. A light absorber absorbs energy of the irradiated pulse laser beam within the test object, and thus a photoacoustic signal is generated. The probe 11 detects the photoacoustic signal generated within the test object. The trigger control circuit 30 outputs a sampling trigger signal to the AD conversion unit 22 in accordance with the output of a Q switch trigger signal. The AD conversion unit 22 receives the photoacoustic signal detected by the probe 11 through the reception circuit 21, and samples the photoacoustic signal with a predetermined sampling period (step B5). The photoacoustic signal sampled by the AD conversion unit 22 is stored as first photoacoustic data in the reception memory 23.
The control unit 31 determines whether a remaining wavelength is present, in other words, whether the pulse laser beams of all the wavelengths constituting the wavelength sequence have been emitted or not (step B6). When a remaining wavelength is present, the process returns to step B2 in order to emit the pulse laser beam having the next wavelength, and the flash lamp trigger signal is output to the laser source unit 13 from the trigger control circuit 30. In step B3, the flash lamp 52 is turned on in response to the flash lamp trigger signal, and in step B4, the trigger control circuit 30 turns on the Q switch 55 at a timing when the birefringent filter 56 is set to be at the rotational displacement position corresponding to the second wavelength (800 nm) constituting the wavelength sequence, to cause the pulse laser beam to be emitted.
The pulse laser beam having a wavelength of 800 nm which is emitted from the laser source unit 13 is guided to, for example, the probe 11, and the first partial region of the test object is irradiated with the pulse laser beam from the probe 11. The probe 11 detects the photoacoustic signal generated by the light absorber within the test object absorbing the pulse laser beam having a wavelength of 800 nm. The trigger control circuit 30 outputs the sampling trigger signal to the AD conversion unit 22 in accordance with the output of the Q switch trigger signal, and the AD conversion unit 22 samples the photoacoustic signal in step B5. The photoacoustic signal sampled by the AD conversion unit 22 is stored as second photoacoustic data in the reception memory 23. The photoacoustic image generation apparatus 10 performs step B1 to step B5 to each wavelengths constituting the wavelength sequence, and irradiates the test object with the pulse laser beam having each wavelengths constituting the wavelength sequence, thereby detecting a photoacoustic signal from the test object. The step B1 to step B5 may be the same as step A1 to step A5 of
When the control unit 31 determines in step B6 that a remaining wavelength is not present, the process proceeds to the transmission and reception of ultrasonic waves. The trigger control circuit 30 transmits the ultrasonic waves to the test object from the probe 11 through the transmission control circuit 37 (step B7). In step B7, the ultrasonic waves are transmitted to the same region as the partial region of the test object which is irradiated with the pulse laser beam. The probe 11 detects reflected ultrasonic waves with respect to the transmitted ultrasonic waves (step B8). The detected reflected ultrasonic waves are sampled in the AD conversion unit 22 through the reception circuit 21, and are stored as reflected ultrasonic data in the reception memory 23.
The control unit 31 determines whether all the partial regions have been selected (step B9). When the partial region to be selected remains, the process returns to step B2. The photoacoustic image generation apparatus 10a performs step B2 to step B6 on each partial region, sequentially irradiates each partial region with pulse laser beams having the wavelengths (750 nm and 800 nm) constituting the wavelength sequence, and stores the first photoacoustic data and the second photoacoustic data in the reception memory 23. In addition, step B7 and step B8 are performed to store the reflected ultrasonic data in the reception memory 23. When the irradiation with the pulse laser beam, the detection of the photoacoustic signal, and the transmission and reception of the ultrasonic waves are performed on all the partial regions, data required to generate a photoacoustic image and an ultrasonic image of one frame is gathered.
When the control unit 31 determines in step B9 that all the partial regions have been selected, the process proceeds to the generation of the photoacoustic image and the ultrasonic image. The data separation unit 32 separates the first and second photoacoustic data and the reflected ultrasonic data from each other. The data separation unit 32 transmits the separated first and second photoacoustic data to the complex number creation unit 24, and transmits the reflected ultrasonic data to the ultrasonic image reconstruction unit 33. The complex number creation unit 24 generates complex number data in which first photoacoustic image data is set to a real part and second photoacoustic image data is set to an imaginary part (step B10). The photoacoustic image reconstruction unit 25 performs image reconstruction from the complex number data generated in step B10, using a Fourier transform method (FTA method) (step B11).
The phase information extraction unit 26 extracts phase information from the reconstructed complex number data (step B12). The intensity information extraction unit 27 extracts intensity information from the reconstructed complex number data (step B13). The detection and logarithmic transformation unit 28 performs a detection and logarithmic transformation process on the intensity information extracted in step B13. The photoacoustic image construction unit 29 generates a photoacoustic image on the basis of the phase information extracted in step B12 and the performing of the detection and logarithmic transformation process, on the intensity information extracted in step B13 (step B14). Here, step B10 to step B14 may be the same as step A8 to step A12 of
The ultrasonic image reconstruction unit 33 generates pieces of data of lines of the ultrasonic image using, for example, a delay addition method. The detection and logarithmic transformation unit 34 obtains an envelope of the pieces of data of the lines which are output by the ultrasonic image reconstruction unit 33, and performs logarithmic transformation on the obtained envelope. The ultrasonic image construction unit 35 generates an ultrasonic image on the basis of the pieces of data of the lines on which the logarithmic transformation is performed (step B15). The image synthesis unit 36 synthesizes the photoacoustic image and the ultrasonic image and displays the synthesized image on the image display unit 14 (step B16).
In this embodiment, the photoacoustic image generation apparatus generates an ultrasonic image in addition to a photoacoustic image. It is possible to observe a portion not capable of being formed as an image in the photoacoustic image by referring to the ultrasonic image. Other effects are the same as those in the first embodiment.
Meanwhile, in the above-described embodiments, an example in which first photoacoustic data and second photoacoustic data are created as complex numbers has been described, but the first photoacoustic data and the second photoacoustic data may be separately reconstructed without being created as complex numbers. Furthermore, herein, a ratio between the first photoacoustic data and the second photoacoustic data is calculated by the creation of complex numbers and by using phase information, but the same effect is obtained even though the ratio is calculated from intensity information of both the pieces of data. In addition, the intensity information can be generated on the basis of signal intensity in a first reconstructed image and signal intensity in a second reconstructed image.
In the generation of a photoacoustic image, the number of wavelengths of a pulse laser beam with which a test object is to be irradiated is not limited two, and the test object may be irradiated with three or more pulse laser beams, and thus the photoacoustic image may be generated on the basis of pieces of photoacoustic data corresponding to the respective wavelengths. In this case, for example, the phase information extraction unit 26 may generate a magnitude relation between relative signal intensities of the pieces of photoacoustic data corresponding to the respective wavelengths, as phase information. In addition, the intensity information extraction unit 27 may generate the signal intensities in the pieces of photoacoustic data corresponding to the respective wavelengths, which are grouped into one, as intensity information.
In the above-described embodiments, a description has been made on the assumption that the trigger control circuit 30 monitors BRF state information and controls a rotation speed of the birefringent filter 56 to have a predetermined rotation speed on the basis of the BRF control signal, but is not limited thereto.
As described above, although the invention has been described on the basis of the preferred embodiments, the ultrasonic wave unit and the photoacoustic image generation apparatus of the invention are not limited to those in the above-described embodiments, and various corrections and modifications made to the configurations of the above-described embodiments may also be included in the scope of the invention.
| Number | Date | Country | Kind |
|---|---|---|---|
| 2011-166977 | Jul 2011 | JP | national |
| 2012-160426 | Jul 2012 | JP | national |
This application is a Continuation of PCT International Application No. PCT/JP2012/068933 filed on Jul. 26, 2012, which claims priority under 35 U.S.C §119(a) to Patent Application No. 2011-166977 filed in Japan on Jul. 29, 2011 and Patent Application No. 2012-160426 filed in Japan on Jul. 19, 2012, all of which are hereby expressly incorporated by reference into the present application.
| Number | Name | Date | Kind |
|---|---|---|---|
| 4272733 | Walling et al. | Jun 1981 | A |
| 4809283 | Harter | Feb 1989 | A |
| 5107509 | Esterowitz et al. | Apr 1992 | A |
| 5175736 | Woodward et al. | Dec 1992 | A |
| 6979292 | Kanayama et al. | Dec 2005 | B2 |
| 7777881 | Yamazoe | Aug 2010 | B2 |
| 20030072333 | Jacobowitz et al. | Apr 2003 | A1 |
| 20030118060 | Spuehler et al. | Jun 2003 | A1 |
| 20050187471 | Kanayama et al. | Aug 2005 | A1 |
| 20090187099 | Burcher | Jul 2009 | A1 |
| 20100049049 | Asao et al. | Feb 2010 | A1 |
| Number | Date | Country |
|---|---|---|
| 102095685 | Jun 2011 | CN |
| 2000-105464 | Apr 2000 | JP |
| 2006-324601 | Nov 2006 | JP |
| 2007-81065 | Mar 2007 | JP |
| 2010-046215 | Mar 2010 | JP |
| Entry |
|---|
| Extended European Search Report, dated Feb. 16, 2015, for European Application No. 12820356.9. |
| A High-Speed Photoacoustic Tomography System based on a Commercial Ultrasound and a Custom Transducer Array, Xueding Wang, Jonathan Cannata, Derek DeBusschere, Changhong Hu, J. Brian Fowlkes, and Paul Carson, Proc. SPIE vol. 7564, 756424 (Feb. 23, 2010). |
| International Search Report, mailed Jul. 26, 2012, issued in PCT/JP2012/068933. |
| Written Opinion of the International Search Authority, mailed Jul. 26, 2012, issued in PCT/JP2012/068933. |
| Notice of Reasons for Allowance dated Jul. 3, 2015 for Chinese Application No. 201280037135.4. |
| Number | Date | Country | |
|---|---|---|---|
| 20140148681 A1 | May 2014 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | PCT/JP2012/068933 | Jul 2012 | US |
| Child | 14167785 | US |
