WO2013098942A1 - 情報信号生成方法 - Google Patents
情報信号生成方法 Download PDFInfo
- Publication number
- WO2013098942A1 WO2013098942A1 PCT/JP2011/080223 JP2011080223W WO2013098942A1 WO 2013098942 A1 WO2013098942 A1 WO 2013098942A1 JP 2011080223 W JP2011080223 W JP 2011080223W WO 2013098942 A1 WO2013098942 A1 WO 2013098942A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- signal
- information
- extended
- extrapolated
- peak
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000001131 transforming effect Effects 0.000 claims abstract description 32
- 230000003287 optical effect Effects 0.000 claims description 116
- 238000001228 spectrum Methods 0.000 claims description 99
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 53
- 238000005259 measurement Methods 0.000 claims description 42
- 238000012545 processing Methods 0.000 claims description 40
- 238000003325 tomography Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 3
- 239000002131 composite material Substances 0.000 abstract description 6
- 238000012014 optical coherence tomography Methods 0.000 description 28
- 239000013307 optical fiber Substances 0.000 description 25
- 230000003595 spectral effect Effects 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 238000007493 shaping process Methods 0.000 description 7
- 239000010410 layer Substances 0.000 description 6
- 230000035945 sensitivity Effects 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000013213 extrapolation Methods 0.000 description 3
- 230000003252 repetitive effect Effects 0.000 description 3
- 230000002194 synthesizing effect Effects 0.000 description 3
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 239000002356 single layer Substances 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02083—Interferometers characterised by particular signal processing and presentation
- G01B9/02084—Processing in the Fourier or frequency domain when not imaged in the frequency domain
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/0209—Low-coherence interferometers
- G01B9/02091—Tomographic interferometers, e.g. based on optical coherence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
Definitions
- the present invention relates to a signal generation method for generating an information signal using Fourier transform, such as a signal generation method for acquiring an optical spectrum interference signal and generating an information signal related to a tomography to be measured by Fourier transforming the optical spectrum interference signal. Regarding the method.
- optical coherence tomographic imaging apparatus Fourier Domain-Optical Coherence Tomography: FD-OCT
- FD-OCT Fourier Domain-Optical Coherence Tomography
- FD-OCT divides the light source output into two or more, and irradiates the object to be measured as one reference light and other light as measurement light. Scattered light or reflected light returns from the object to be measured, and an optical spectrum interference signal with the aforementioned reference light is acquired.
- the optical spectrum interference signal is observed on the wave number space axis, and a signal oscillated on the wave number space axis is obtained according to the optical path length difference between the reference light path and the measurement light path. Therefore, a tomographic information signal indicating a peak according to the optical path length difference is obtained by Fourier transforming the acquired optical spectrum interference signal.
- the intensity of the optical spectrum interference signal is proportional to the product of the reference light and the intensity of the return light from the measurement object, even if the return light from the measurement object is attenuated by absorption, scattering, or transmission, it is highly sensitive. An information signal can be obtained.
- the tomographic information signal obtained by Fourier transforming the optical spectrum interference signal is a convolution operation of a Fourier transform signal of a sine wave having a frequency corresponding to the optical path length difference and a shape obtained by Fourier transforming the spectrum shape.
- the wider the spectral band for example, depending on the wavelength band of the light source
- the higher the resolution in the depth direction the ability to decompose and display the layer structure
- the spectrum band is finite and the spectrum shape has a certain shape. Therefore, the shape of the tomographic information signal reflects the shape obtained by Fourier transform of the spectrum shape.
- the image is displayed as if the tomographic layer which is originally a single layer is divided into a plurality of layers, or the layers are displayed in a single layer without being decomposed, or the bottom of the tomographic information signal having a high strength is low in strength. Degradation of the tomographic information signal such that the tomographic information signal is buried occurs.
- Such signal deterioration is not limited to optical tomographic measurement, but is a common problem in signal processing and signal analysis processing involving Fourier transform.
- a super-spatially-variant-apodization as described in Patent Document 1 is performed a plurality of times of extrapolation, and decomposition is performed beyond the diffraction limit.
- a method has been proposed. In this method, an extrapolated signal is generated by extracting a waveform in a wide band around the peak of the Fourier-transformed signal and performing inverse Fourier transform. Then, the resolution is increased by extrapolating the extrapolated signal to the original waveform. By repeating this signal processing, signals are acquired with higher resolution.
- Patent Document 2 erases adjacent peaks, and there is a fear that information that should originally be erased.
- An information signal generation method provided by the present invention includes: In the signal generation method for generating the information signal F (b) by Fourier-transforming the signal f (a), Detecting the maximum peak intensity of the information signal F (b); Calculating the amplitude, phase and frequency of the signal f (a) corresponding to the maximum value of the peak intensity; Generating a signal in which the signal f (a) corresponding to the maximum value of the peak intensity is extended on the a-axis based on the amplitude, phase and frequency information; a1 ⁇ a2 and extracting a signal in a region smaller than a1 from the extended signal and a signal in a region larger than a2 in the extended signal to generate an extrapolated signal; Combining the extrapolated signal and the signals f (a) from the a1 to a2 to generate a combined signal; Further Fourier transforming the combined signal; It is characterized by having.
- an extrapolated signal is generated based on the signal f (a) corresponding to the maximum value of the peak intensity of the information signal F (b), and a synthesized signal is generated using the extrapolated signal. Is Fourier transformed.
- an extrapolated signal is generated using the amplitude, phase, and frequency indicated by the peak, and a composite signal is generated, thereby generating side noise associated with the peak. Can be suppressed. As a result, the information signal having a small peak intensity value buried in the side noise remains without being erased, and a more detailed information signal can be obtained.
- the schematic diagram for demonstrating the process of the signal generation method of this invention The schematic diagram for demonstrating the process of the signal generation method of this invention Schematic diagram for explaining repetitive processing in the present invention Flowchart for explaining repetitive processing in the present invention Schematic diagram for explaining the effect of the window function Graph for explaining the effect of the present invention Schematic diagram of an optical coherence tomography apparatus to which the present invention is applied Schematic diagram of an optical coherence tomography apparatus to which the present invention is applied
- the information signal generation method for generating the information signal F (b) by performing Fourier transform on the signal f (a) of the present invention includes the following steps.
- Detecting the maximum peak intensity of the information signal F (b). Calculating the amplitude, phase and frequency of the signal f (a) corresponding to the maximum value of the peak intensity; Generating a signal f ′ (a) obtained by extending the signal f (a) corresponding to the maximum value of the peak intensity on the a-axis based on the amplitude, phase, and frequency information.
- a step of generating an extrapolated signal by extracting a signal in a region smaller than a1 from the extended signal and a signal in a region larger than a2 in the extended signal, where a1 ⁇ a2. Synthesizing the extrapolated signal and the signals f (a) from a1 to a2 to generate a synthesized signal.
- the synthesized signal is further subjected to Fourier transform.
- the present invention includes repeating the signal generation method a plurality of times.
- optical spectrum interference signal and tomographic information signal in Fourier domain optical coherence tomography can be employed in the method (FT-IR).
- FD-OCT Fourier Domain-Optical Coherence Tomography
- FT-IR Fourier transform infrared spectroscopy Interference signal and spectral information signal
- optical intensity signal and analysis signal thereof optical intensity signal and analysis signal thereof
- audio signal and audio analysis signal, etc. can be employed in the method (FT-IR).
- the acquired optical spectrum interference signal I (k) is expressed by the following equation (1).
- k is the wave number
- x is the optical path length difference between the reference optical path and the measurement optical path
- S (k) is the spectral shape
- the tomographic information signal obtained by Fourier transforming the optical spectrum interference signal represented by the above formula (1) is
- the shape of the tomographic information signal reflects the Fourier transform shape of the spectrum shape.
- a tomographic information signal independent of the spectrum shape can be acquired by the following signal generation method.
- the present invention by taking an extrapolated signal in a wide band, it becomes possible to multiply the window function in a wide band, narrowing the width of the tomographic information signal of the peak intensity extrapolated from the extrapolated signal, Increase the resolution.
- the resolution ⁇ z of the tomographic information signal is expressed by the following equation (3).
- ⁇ c is the center wavelength of the Gaussian window
- ⁇ is the full width at half maximum of the Gaussian window.
- 301 indicates a spectrum obtained by the method of the present invention, and 301 indicates an extrapolation for a region other than the region from the wave number ks to ke indicating the spectral band of the light source, that is, a region less than ks and a region exceeding ke. It is obtained by adding a signal.
- 302 represents a spectrum corresponding to the prior art to which the present invention is not applied, and is obtained by performing a Fourier transform by multiplying a Gaussian window such that the end of the spectrum band of the light source output is 1 / e 2. is there.
- the full width at half maximum ⁇ k of the spectrum according to the present invention is about 2.4 times wider than 302, so that the resolution is improved about 2.4 times.
- the window function can be multiplied in the wide band, and the side noise around the tomographic information signal having the peak intensity extrapolated from the extrapolated signal can be reduced.
- FIG. 6 shows a signal spectrum corresponding to the prior art to which the present invention is not applied (a function in which the Gaussian function becomes zero at half value), and (B) shows a Fourier transform of the signal spectrum. The resulting signal is shown.
- FIG. 6C shows a signal spectrum to which the present invention is applied
- FIG. 6D shows a signal obtained by Fourier transforming the signal spectrum.
- the extrapolated signal can be extrapolated until the value of the bottom of the window function becomes sufficiently small, and the signal noise level obtained by Fourier transform as shown in FIG. (Side noise) can be reduced by 10 digits or more depending on the peak intensity value.
- the noise can be reduced.
- the noise level is improved by about one digit, and when the value of the Gaussian window function is extended to 0.001, the noise level is about 12 digits. Improved.
- the signal-to-noise ratio (SNR) of the image is desired to be 96 dB or higher, and the extension width is determined so that the noise level exceeds 96 dB.
- FIG. 1 (A) shows an optical interference spectrum signal obtained using a light source having a spectrum band in the range of wave number ks to ke. This will be explained by a specific example. Light output from a light source having a wide emission wave number or a light source whose emission wave number changes in a wide band in time is divided into at least two, and one of the lights is irradiated to the measurement object. The optical spectrum interference signal is shown by causing the scattered light or reflected light from the measurement object obtained by this to interfere with the other of the light.
- FIG. 1B is not an essential step in the present invention, but is performed as necessary, but a signal obtained by multiplying the optical spectrum interference signal of FIG. 1A by a window function for shaping an envelope. Is shown.
- FIG. 1 (C) shows an information signal related to a tomographic object to be measured, which is obtained by Fourier transforming the optical spectrum interference signal of FIG. 1 (B).
- FIG. 1 (C) shows an information signal related to a tomographic object to be measured, which is obtained by Fourier transforming the optical spectrum interference signal of FIG. 1 (B).
- FIG. 1 (C) shows an information signal related to a tomographic object to be measured, which is obtained by Fourier transforming the optical spectrum interference signal of FIG. 1 (B).
- FIG. 1 (C) shows an information signal related to a tomographic object to be measured, which is obtained by Fourier transforming the optical spectrum interference signal of FIG. 1 (B).
- the maximum peak intensity value 101 is detected from the information signal.
- the amplitude, phase, and frequency ⁇ of the optical spectrum interference signal corresponding to the maximum peak intensity value 101 are calculated, and the optical spectrum interference signal corresponding to the maximum peak intensity value 101 is calculated as the amplitude as shown in FIG. And a signal extended on the wave number space axis based on the phase and frequency information.
- the extended signal FIG. 2E
- FIG. 2F the extended signal
- the wave number ks and the wave number ke in the wave number band emitted by the light source are set to ks ⁇ ke, and the signal in the region larger than the wave number ke among the extended signal and the signal in the region smaller than ks.
- the signal is extracted to generate an extrapolated signal (FIG. 2G).
- a synthesized signal is generated by synthesizing the extracted signal obtained by extracting signals from wave number ks to wave number ke from the optical spectrum interference signal and the extrapolated signal (FIG. 2 (H)).
- the resulting synthesized signal is further Fourier transformed (FIG. 2 (I)).
- the information signal shown in FIG. 2 (I) has a narrower spectral width than that of the information signal shown in FIG. 1 (C), and noise is remarkably suppressed.
- the present invention includes a signal generation method having the following steps.
- a step of generating an extrapolated signal by extracting a signal in a region smaller than ks from the extended signal and a signal in a region larger than ke among the extended signals, where wave number ks ⁇ ke.
- the tomographic information signal is converted into a single peak by shaping the spectrum shape. And tomographic image quality can be improved.
- an optical spectrum interference signal is acquired (201).
- the envelope of the interference signal matches the spectral intensity of the light source output. Therefore, it is also possible to flatten the envelope of the interference signal and shape it into a rectangular wave shape by using the spectral intensity of the light source output separately obtained. It is also possible to multiply the rectangular wave shaped interference signal by a window function. Multiplying the window function can be done as needed.
- an information signal relating to the tomographic object to be measured is obtained by Fourier transforming the interference signal (202).
- the peak having the n-th peak intensity is detected from the information signal regarding the acquired fault (203).
- an extrapolated signal is generated corresponding to the maximum peak intensity of the obtained information signal (for example, an information signal related to a fault), a synthesized signal is generated, and a Fourier transform is performed up to a predetermined number N. Includes repeating.
- the detected amplitude value of the n-th peak intensity corresponds to the amplitude, phase, and frequency. Is calculated (205).
- the n-th peak intensity represents a vibration signal of an interference signal that vibrates at a frequency corresponding to the difference in optical path length between the reference optical path and the measurement optical path.
- a vibration waveform extended to the extrapolated wave number band is generated (206).
- a step of multiplying the extended signal by the window function may be added. From the signal extended on the wave number space axis (or from the signal multiplied by the window function), the signal in the region smaller than the wave number ks and the signal in the region where the wave number is larger than ke are extracted to generate an extrapolated signal. (207).
- the extracted signal obtained by extracting signals from wave number ks to wave number ke from the optical spectrum interference signal is synthesized with the extrapolated signal to generate a synthesized signal (208).
- the above composite signal is Fourier transformed to obtain a tomographic information signal (209).
- the processing up to this point increases the resolution of the detected peak intensity signal and suppresses side noise components around the peak intensity.
- step 203 only signal components distinguished from noise are extracted. Therefore, side noise components generated depending on the light source spectrum band can be suppressed by steps 201 to 209 and 210, and a tomographic information signal with high resolution and high sensitivity can be obtained.
- FIG. 3A shows an optical spectrum interference signal
- FIG. 3B shows an information signal obtained by Fourier transforming the optical spectrum interference signal.
- 351 is the maximum value of the peak intensity of the information signal (step 203 in FIG. 4).
- the n-th peak value is the first.
- FIG. 3C shows a synthesized signal (extrapolated signal) using a signal (extrapolated signal) extended on the wave number space axis based on the amplitude, phase and frequency of the optical spectrum interference signal corresponding to the maximum value 351 of the peak intensity. 208 steps).
- FIG. 3 (D) is an information signal obtained by Fourier transform of the combined signal, and 352 indicates the second peak intensity of the information signal obtained by first Fourier transform.
- 352 corresponds to the second peak value in the (n + 1) th peak value.
- 353 in FIG. 3 (F), 354 in FIG. 3 (H), and 355 in FIG. 3 (J) respectively indicate the third, fourth, and fifth peak intensities of the information signal obtained by the first Fourier transform. Show.
- an information signal in which noise components are remarkably suppressed is obtained by repeating the Fourier transform cycle of the present invention.
- N can also be regarded as the number of information signal peaks related to the fault. If signal generation is performed from this viewpoint, it is possible to increase the resolution and sensitivity of the entire information signal related to the tomography with a minimum number of repetitions.
- the N may be the number of peaks up to a peak that is equal to or lower than a predetermined intensity when the detected peak intensity is larger. This makes it possible to increase the resolution and sensitivity of a signal related to a tomogram of a certain value or more without repeating the calculation indefinitely.
- This method can take the following steps, for example.
- (a) The acquired interference signal is recorded so as to have an equal wave number interval, or data conversion is performed so as to have an equal wave number interval by signal processing.
- (b) The envelope of the interference signal matches the spectral intensity of the light source output. Therefore, using the spectral intensity of the light source output separately obtained, the envelope of the interference signal is flattened and shaped into a rectangular wave shape.
- An information signal related to the tomographic object to be measured is acquired by performing Fourier transform on the interference signal multiplied by the window function.
- a maximum peak value is detected from the acquired information signal regarding the fault, and a value obtained by multiplying the maximum peak value by a predetermined ratio m times is defined as an mth threshold value, and a peak intensity is equal to or greater than the mth threshold value. Detect all peaks.
- Each peak intensity equal to or greater than the mth threshold value represents an oscillation signal of an interference signal that vibrates at a frequency corresponding to the optical path length difference between the reference optical path and the measurement optical path. Therefore, the amplitude, phase, and frequency corresponding to each detected peak intensity value are calculated.
- a vibration waveform for each extrapolated wave number band is generated and added to obtain an extrapolated signal.
- M is determined in advance in the signal generation step (l). Alternatively, it may be repeated until the number of detected peaks reaches a predetermined number, instead of determining the Mth. This makes it possible to limit the number of calculations. Or you may repeat until a threshold value becomes below a predetermined value. This makes it possible to increase the resolution and sensitivity of the entire signal relating to a fault greater than a certain value without repeating the calculation indefinitely.
- the predetermined ratio for lowering the threshold value may be a ratio between the maximum value of the signal intensity obtained by acquiring the envelope waveform of the interference signal and Fourier transforming the envelope waveform and the second peak value. As a result, the predetermined ratio for lowering the threshold value can be maximized and the threshold value can be lowered efficiently.
- the signal generation method of the present invention is obtained by multiplying the signal f (a) by a window function in which the bottom of the waveform indicated by the function is 0 and extending the extended signal to a where the window function is 0. Includes form.
- the range in which the signal f ′ (a) is extended has a predetermined ratio of the maximum peak value of the signal obtained by Fourier transform of the window function in the extended range and the second peak value next to the peak value.
- the form which is more than the range used as a value is included.
- the signal f (a) is obtained by dividing the signal f (a) by the envelope information of the signal f (a) to shape the shape of the signal f (a) into a flat top shape. Can be multiplied by.
- side noise components proportional to the magnitude of the peak intensity value are sequentially suppressed by repeatedly generating and extrapolating an extrapolated signal in order from an information signal having a large peak intensity value. Since the side noise component to be suppressed is attached to the detected peak, the information signal buried in the side noise is not erased. Therefore, by suppressing the side noise components in order from the information signal having the largest peak intensity value, it is possible to correctly extract only the information signal and suppress the noise component over the entire information signal.
- the threshold value By reducing the threshold value as much as possible, many correct information signal peaks can be detected at one time. However, if the threshold value is too small, side noise is also detected as an information signal. On the other hand, the magnitude of the side noise is determined by a waveform obtained by Fourier transforming the envelope shape of the original signal. Therefore, the ratio of the maximum peak intensity of the waveform obtained by Fourier transform of the envelope shape and the peak value having the second magnitude is calculated and multiplied by the maximum peak intensity of the actually obtained information signal. By using as a threshold value, it is possible to correctly detect the peak intensity related to the information signal while avoiding side noise components.
- the extrapolated signal can theoretically widen the spectrum band infinitely.
- the spectrum band is infinite, the number of data becomes too large, and a process such as Fourier transform takes time.
- the shape of the information signal reflects the waveform shape obtained by Fourier transforming the envelope shape of the original signal. Therefore, in general, the shape of the information signal is shaped by multiplying the original signal by a window function such as a cos function, a Hanning function, or a Gaussian function that is symmetric and whose skirt approaches 0. Therefore, by setting the range of the extrapolated signal to a range in which the tail of the window function is sufficiently close to 0, the number of data can be limited, and the information signal can be increased in resolution and noise at high speed.
- a window function such as a cos function, a Hanning function, or a Gaussian function that is symmetric and whose skirt approaches 0. Therefore, by setting the range of the extrapolated signal to a range in which the tail of the window function is sufficiently close to 0, the number of data can be limited,
- the envelope shape of the original signal is not flat top, even if the original signal is directly multiplied by the window function, the influence of the envelope shape of the original signal remains and the information signal is also affected. Therefore, by dividing the optical spectrum interference signal by the envelope shape information and flattening the signal, it is possible to suppress the influence of the envelope shape of the original signal and perform shaping by the window function.
- the present invention includes an information signal processing device that generates an information signal F (b) by performing Fourier transform on the signal f (a).
- This information signal processing apparatus is configured to include the following.
- a processing unit that detects the maximum peak intensity of the information signal F (b).
- a processing unit that calculates the amplitude, phase, and frequency of the signal f (a) corresponding to the maximum peak intensity of the information signal F (b).
- a processing unit that generates a signal extended on the a-axis based on the amplitude, phase, and frequency information.
- a processing unit that generates an extrapolated signal by extracting a signal in a region smaller than a1 from the extended signal and a signal in a region larger than a2 out of the extended signal, where a1 ⁇ a2.
- a processing unit that combines the extrapolated signal and the signals f (a) from a1 to a2 to generate a combined signal.
- the present invention also includes an optical coherence tomography apparatus.
- the optical coherence tomography apparatus of the present invention includes a light source unit, a sample measurement unit that irradiates the sample with light from the light source unit, and transmits reflected light from the sample, and transmits light from the light source unit as reference light.
- a reference unit that causes the reflected light and the reference light to interfere with each other, a light detection unit that detects interference light from the interference unit, and light detected by the light detection unit.
- the image processing unit includes the signal generation device of the present invention.
- the optical coherence tomography apparatus shown in the present embodiment is a wavelength swept light source type optical coherence tomography apparatus (SweptwSource-Optical Coherence Tomography: SS-OCT) among FD-OCTs, and a schematic diagram is shown in FIG. In FIG. 7, the apparatus is configured using optical fibers, but may be configured in space without using optical fibers.
- SS-OCT wavelength swept light source type optical coherence tomography apparatus
- the imaging device (SS-OCT device) shown in FIG. 7 uses a light source 501 (light source unit) whose emission wavelength changes from 800 nm to 880 nm during 10 ⁇ s.
- a light source 501 for example, a light source that temporally changes the wavelength cut out from the broadband light source is used.
- the light output from the light source 501 is divided into two by an optical fiber coupler 510, and one of the divided lights is irradiated as a measurement light to the measurement object 515, and scattered light from the measurement object or The reflected light is again coupled to the optical fiber coupler 510 (interference unit) and propagated to the optical fiber coupler 518.
- the measurement light passes through a measurement light scanning optical system (specimen measurement unit) constituted by two galvanometer mirrors 512 and 513 orthogonal to each other. The measurement light is scanned over the measurement object by the measurement light scanning optical system.
- the delay line is a spatial optical system capable of changing the distance between the outgoing optical fiber end 516 and the incoming optical fiber end 517.
- the optical spectrum interference signal related to the tomography is obtained by the differential photodetector 519 (light detection unit) that detects the interference light by causing the scattered light or reflected light to interfere with the reference light.
- the differential photodetector 519 detects light with a response speed of 205 MHz. Therefore, optical spectrum interference data regarding 2050 faults is output during one wavelength sweep of the SS-OCT device.
- the output optical spectrum interference data relating to the tomographic image is taken into the computer 520 constituting the image processing unit for obtaining the tomographic image via the AD board.
- the light output from the wavelength swept light source is partially guided by the optical fiber coupler 504 to the Mach-Cender interferometer. Since the differential optical detection data of the optical spectrum interference signal from the Mach-Cender interferometer intersects with 0 at equal wave number intervals, it passes through another channel of the AD board as a wave number clock signal indicating a change in wave number over time. It is taken into the computer 520. Using the wave number clock signal, the optical spectrum interference data relating to the tomography is converted into data inside the computer so as to become interference signals of equal wave number intervals. At this time, data conversion is performed so that the number of data becomes 2048 points. Therefore, in order to obtain 2048 points or more where the wave number clock signal crosses 0, the optical path length difference between the two arms of the Mach-Cender interferometer is set to 8 mm or more.
- the light partially divided by the optical fiber coupler 502 is detected by the photodetector 503, and the intensity change of the light source is monitored.
- the intensity change information is also taken into the computer 520 by the AD board, used as an envelope of the optical spectrum interference signal, and used for flat-topping. By flattening the signal, it is possible to suppress the influence of the envelope shape of the original signal and perform shaping by the window function.
- the data acquisition by the AD board is synchronized using the trigger signal 521 output from the wavelength swept light source.
- the interference signal converted into data at the equal wave interval is subjected to fast Fourier transform, and a tomographic information signal on the optical axis irradiated with the measurement light is obtained.
- the Fourier-transformed signal includes the folding signal, half of the Fourier-transformed signal is extracted and used as the tomographic information signal.
- a reference optical path i.e., an optical path from an optical fiber coupler 510 that first branches light in the reference optical path to an optical fiber coupler 518 where interference occurs through a delay line
- a measurement optical path i.e., a measurement optical path.
- the frequency of the optical spectrum interference signal increases as the optical path length difference from the first optical fiber coupler 510 that splits the light to the optical fiber coupler 518 that is irradiated to the measurement object and causes interference is increased. Due to the structure of the apparatus, the optical path length difference is about 8 mm compared to 0 mm, and the amplitude of the optical spectrum interference signal is halved. Therefore, before acquiring the tomographic image, the delay line (distance between the emitted optical fiber end 516 and the incident optical fiber end 517) is set so that the tomographic information signal on the uppermost surface of the measurement object is near 0. Adjust.
- the 512 galvanometer mirrors constituting the measurement light scanning optical system are changed in angle in 1024 directions, optical spectrum interference signals are obtained in each direction, data conversion is performed at equal wave intervals, and fast Fourier transform is performed. A series of data processing is performed. Based on the intensity of the tomographic information signal in each direction of the galvano mirror 512, a tomographic image is obtained by adding shades in gray scale and arranging the tomographic information in 1024 directions.
- 505 and 508 are fiber couplers in a Mach-Zehnder interferometer for generating a wave number clock signal
- 506 and 507 are collimator lenses
- 509 is a differential photodetector
- 502 and 504 are fiber couplers.
- the signal generation method according to the present invention is applied to improve the image quality of the tomographic image.
- the signal generation method according to the present invention is applied for each direction of the acquired tomographic information signal.
- the maximum value of the peak intensity of the tomographic information signal is detected.
- the amplitude, phase and frequency of the optical spectrum interference signal corresponding to the maximum value of the peak intensity are calculated.
- the amplitude information is calculated from the magnitude of the detected peak value.
- the phase information is calculated from the real part and the imaginary part of the data obtained by Fourier transforming the optical spectrum interference signal.
- the frequency is calculated from the position of the detected peak. Based on the calculated amplitude, phase and frequency information, a sine wave signal extended on the wave number space axis is generated.
- the shape of the optical spectrum output from the light source is maximum at a wave number of 7496982.469 [m -1 ] (corresponding to a wavelength of 838.1 [nm]), and the wave number of 7199883 [m -1 ] (corresponding to a wavelength of 880 [nm]) It was a Gaussian function shape with an optical spectrum intensity at half value at a wave number 7785882 [m -1 ] (corresponding to a wavelength of 800 [nm]).
- wave number 8460395.981 [m -1 ] (corresponding to wavelength 742.6585376 [nm]) from which Gaussian function value falls to 1 / e 2 with respect to the maximum value, wave number 6533568.956 [m -1 ] (corresponding to wavelength 961.6773542 [nm]) ) Until the sine wave signal is extended.
- the extracted signal is added to the optical spectrum interference signal and further subjected to fast Fourier transform.
- Fig. 8 shows a schematic diagram of the spectral domain optical coherence tomography (Spectral Domain-Optical Coherence Tomography: SD-OCT) of this example.
- the apparatus is configured using optical fibers, but may be configured in space without using optical fibers.
- This SD-OCT apparatus uses a broadband light source 601 having a light emission wavelength from 800 [nm] to 880 [nm].
- a broadband LD Super-Luminescent Diode: SLD
- SC light super-continuum light
- the light output from the light source is divided into two by the optical fiber coupler 602, one of the divided lights is irradiated as the measurement light to the measurement object 607, and the scattered light or reflected light from the measurement object is irradiated.
- the optical fiber coupler 602 is coupled again.
- the measurement light passes through the measurement light scanning optical system constituted by two galvanometer mirrors 604 and 605 which are orthogonal to each other.
- the measurement light is scanned over the measurement object by the measurement light scanning optical system.
- the other of the divided lights is used as a reference light, and a delay line that can change the distance from the optical fiber end 608 to the mirror 609 is propagated.
- the reference light is applied to the mirror 609 of the delay line, and the reflected reference light is also coupled to the optical fiber coupler 602 again.
- the scattered light or reflected light and the reference light are caused to interfere with each other, and a spectroscope 610 acquires an optical spectrum interference signal.
- the detected optical spectrum interference signal is taken into the computer 611 constituting the image processing unit.
- the light detected by the adjacent pixels is not an equal wave number interval but an equal wavelength interval. Therefore, using the data prepared in advance, an internal interference signal with an equal wave number interval is used in the computer 611. Data is converted.
- the interference signal that has been data-converted at equal wave intervals is subjected to fast Fourier transform to obtain a tomographic information signal on the optical axis irradiated with the measurement light.
- the Fourier-transformed signal includes the folding signal, half of the Fourier-transformed signal is extracted and used as the tomographic information signal.
- the difference becomes longer, the frequency of the optical spectrum interference signal becomes higher.
- the amplitude of the optical spectrum interference signal becomes 1/2 when the optical path length difference is about 8 [mm] compared to 0 [mm]. Therefore, before acquiring the tomographic image, the delay line (distance from the optical fiber end 608 to the mirror 609) is adjusted so that the uppermost tomographic information signal of the measurement target comes near 0.
- One 604 of the galvano mirrors constituting the measurement light scanning optical system is changed in angle to 1024 directions, optical spectrum interference signals are obtained in each direction, data conversion is performed so as to have equal wave number intervals, and fast Fourier transform is performed. A series of data processing is performed. Based on the intensity of the tomographic information signal in each direction of the galvanometer mirror 604, a tomographic image is obtained by adding shades in gray scale and arranging the tomographic information in 1024 directions.
- the signal generation method according to the present invention is applied to improve the image quality of the tomographic image.
- the signal generation method according to the present invention is applied for each direction of the acquired tomographic information signal.
- the maximum value of the peak intensity of the tomographic information signal is detected.
- the amplitude, phase and frequency of the optical spectrum interference signal corresponding to the maximum value of the peak intensity are calculated.
- the amplitude information is calculated from the magnitude of the detected peak value.
- the phase information is calculated from the real part and the imaginary part of the data obtained by Fourier transforming the optical spectrum interference signal.
- the frequency is calculated from the position of the detected peak. Based on the calculated amplitude, phase and frequency information, a sine wave signal extended on the wave number space axis is generated.
- the shape of the optical spectrum output from the light source is maximum at a wave number of 7496982.469 [m -1 ] (corresponding to a wavelength of 838.1 [nm]), and the wave number of 7199883 [m -1 ] (corresponding to a wavelength of 880 [nm]) It was a Gaussian function shape with an optical spectrum intensity at half value at a wave number 7785882 [m -1 ] (corresponding to a wavelength of 800 [nm]).
- wave number 8460395.981 [m -1 ] (corresponding to wavelength 742.6585376 [nm]) from which Gaussian function value falls to 1 / e 2 with respect to the maximum value, wave number 6533568.956 [m -1 ] (corresponding to wavelength 961.6773542 [nm]) ) Until the sine wave signal is extended.
- the extracted signal is added to the optical spectrum interference signal and further subjected to fast Fourier transform.
- the signal generation method in all directions of the galvanometer mirrors 604 and 605, the image quality of the three-dimensional volume tomographic data can be improved.
- signal processing described in Example 1 is performed to improve image quality. Subsequently, for the second peak from the highest peak intensity value, extraction of the peak and calculation of the amplitude, phase, and frequency for the extracted peak are performed.
- the wave number 6533568.956 [m -1 ] (wavelength 961.6773542 [nm]) Sine wave signal is generated based on the calculated amplitude, phase and frequency.
- a signal in a region smaller than a wave number 7139893 [m -1 ] (corresponding to a wavelength of 880 [nm]) and a wave number 7785882 [m -1 ] (corresponding to a wavelength of 800 [nm]) of the extrapolated signal A signal in a larger area is extracted, and the extracted extrapolated signal and the optical spectrum interference signal are combined.
- the tomographic image quality can be further improved by the above signal processing as compared with the first embodiment.
- AD board used this time is because it is 12bit, the detected peak value may be as a signal generation method that repeats until the 10 -3.6 ( ⁇ 2- 12). In this, it becomes possible to reduce the noise level to 10 -3.6 ( ⁇ 2- 12).
- a peak is not detected one by one, but a threshold is provided, and the amplitude, phase, and frequency are calculated for each peak for peaks that are equal to or higher than the threshold.
- wave number 8460395.981 [m -1 ] (corresponding to wavelength 742.6585376 [nm]) from which the value of the Gaussian function falls to 1 / e 2 with respect to the maximum value
- wave number 6533568.956 [m -1 ] (wavelength 961.6773542 [nm])
- a sine wave signal is generated based on the calculated amplitude, phase and frequency.
- a signal in a region smaller than a wave number 7139893 [m -1 ] (corresponding to a wavelength of 880 [nm]) and a wave number 7785882 [m -1 ] (corresponding to a wavelength of 800 [nm]) of the extrapolated signal A signal in a larger area is extracted, and the extracted extrapolated signal and the optical spectrum interference signal are combined.
- the synthesized signal is fast Fourier transformed to obtain a tomographic information signal.
- the threshold value is obtained by acquiring the optical spectrum intensity waveform output from the light source, and detecting the maximum value and the second peak value of the signal intensity obtained by Fourier transforming the optical spectrum intensity waveform, A ratio between the peak value and the maximum peak value may be calculated and determined as a value obtained by multiplying the maximum peak value of the tomographic information signal obtained by Fourier transforming the optical spectrum interference signal. .
- the tomographic information signal is acquired in the same manner as in the first or second embodiment.
- a value obtained by multiplying the maximum peak value P 0 of the tomographic information signal obtained by Fourier transforming the optical spectrum interference signal by the reduction width 0.2 is set as a threshold value T m
- the reduction width 0.2 was obtained by referring to the ratio 0.2 between the first peak and the second peak of the sinc function obtained by Fourier transform of the rectangular wave.
- a signal generation method that repeats until the detected peak value reaches 10 ⁇ 3.6 may be used.
- the signal generation method may be repeated until the T m threshold value becomes 10 ⁇ 3.6 . In this, it becomes possible to reduce the noise level to 10 -3.6 ( ⁇ 2- 12).
- the reduction width is not 0.2, obtain the optical spectrum intensity waveform output from the light source, detect the maximum value and the second peak value of the signal intensity obtained by Fourier transform the optical spectrum intensity waveform, The ratio between the second peak value and the maximum peak value may be reduced.
- the tomographic information signal is acquired in the same manner as in the first or second embodiment.
- the window function may be a Gaussian function, a cos function or a Hanning function.
- the extrapolated signal is extended to a wave number where the skirt of the window function is zero.
- the number of data can be limited, and it becomes possible to increase the resolution and noise of the information signal at high speed.
- the window function is a Gaussian function
- the extended finite wave number range is set to a range in which the ratio between the maximum value of the waveform obtained by Fourier transform of the window function within the extended finite wave number range and the second peak value is 10 ⁇ 3.6 .
- the number of data can be limited, and it becomes possible to increase the resolution and noise of the information signal at high speed.
- the optical spectrum interference signal may be shaped by flattening it by dividing it by spectral intensity information obtained separately, and then multiplied by a window function.
- the extrapolated signal is multiplied by the same window function.
- the peak of the tomographic information signal after Fourier transform can be made unimodal and the tomographic image quality can be improved.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Signal Processing (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Health & Medical Sciences (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
Description
信号f(a)をフーリエ変換することで情報信号F(b)を生成する信号生成方法において、
前記情報信号F(b)のピーク強度の最大値を検出する工程と、
前記ピーク強度の最大値に対応する前記信号f(a)の振幅と位相と周波数とを算出する工程と、
前記ピーク強度の最大値に対応する前記信号f(a)を前記振幅及び位相及び周波数情報に基づきa軸上に延長した信号を生成する工程と、
a1 < a2として、前記延長した信号からa1よりも小さい領域の信号と、前記延長した信号のうちa2よりも大きい領域の信号とを抽出して外挿信号を生成する工程と、
前記外挿信号と前記a1からa2までの信号f(a)とを合成して合成信号を生成する工程と、
前記合成信号をさらにフーリエ変換する工程と、
を有することを特徴とする。
前記情報信号F(b)のピーク強度の最大値を検出する工程。
前記ピーク強度の最大値に対応する前記信号f(a)の振幅と位相と周波数とを算出する工程。
前記ピーク強度の最大値に対応する前記信号f(a)を前記振幅及び位相及び周波数情報に基づきa軸上に延長した信号f’(a)を生成する工程。
a1 < a2として、前記延長した信号からa1よりも小さい領域の信号と、前記延長した信号のうちa2よりも大きい領域の信号とを抽出して外挿信号を生成する工程。
前記外挿信号と前記a1からa2までの信号f(a)とを合成して合成信号を生成する工程。
そして前記合成信号をさらにフーリエ変換する工程である。
となり、光路長差に応じた周波数を有する正弦波のフーリエ変換信号とスペクトル形状をフーリエ変換した形状の畳みこみ演算となる。そのため、スペクトル帯域(例えば、光源の波長帯域に依存)が広ければ広いほど奥行き方向に高い分解能(層構造を分解表示できる能力)を有した断層情報信号が得られる。
(a) 取得された干渉信号は等波数間隔となるように記録する、あるいは信号処理により等波数間隔となるようにデータ変換する。
(b) 前記干渉信号のエンベロープは光源出力のスペクトル強度と一致している。そこで別に取得された光源出力のスペクトル強度を用いて、前記干渉信号のエンベロープをフラットトップ化し矩形波形状に整形する。
(c) 前記矩形波整形された干渉信号に窓関数を乗じる。
(d) 前記窓関数を乗じられた干渉信号をフーリエ変換することで前記測定対象の断層に関する情報信号を取得する。
(e) 前記取得された断層に関する情報信号から最大ピーク値を検出し、該最大ピーク値に所定の割合をm回乗じた値をm番目の閾値とし、ピーク強度が該m番目の閾値以上のピークを全て検出する。
(f) 前記m番目の閾値以上の各々のピーク強度は、参照光路と測定光路の光路長差に応じた周波数で振動する干渉信号の振動信号を表している。そこで、前記検出された各々のピーク強度の値が対応する振幅、位相、周波数を算出する。
(g) 算出された各々の振幅、位相、周波数情報に基づき、各々の外挿する波数帯域分の振動波形を生成し、足し合わせて外挿信号とする。
(h) 前記外挿信号に前記窓関数を乗じる。
(i) 前記窓関数を乗じられた外挿信号から波数ksよりも小さい領域の信号と前記外挿信号のうち波数がkeよりも大きい領域の信号とを抽出する。
(j) 前記抽出された外挿信号を、前記窓関数を乗じられた干渉信号に外挿した合成信号を生成する。
(k) 前記合成信号をフーリエ変換し断層情報信号を取得する。ここまでの処理により、検出されたピーク強度の信号は高分解能化され、該ピーク強度周辺のサイドノイズ成分は抑制される。
(l) 上記(e)~(k)の操作を、閾値を前記所定の割合だけ下げながらm=1番目、2番目、3番目・・・M番目まで繰り返す。
202 干渉信号をフーリエ変換して情報信号生成
203 情報信号のn番目のピーク強度検出
205 ピーク強度の最大値に対応する干渉信号の振幅、位相及び周波数を算出
206 振幅、位相及び周波数情報に基づき波数空間軸上に延長した信号生成
207 外挿信号の生成
208 合成信号の生成
209 合成信号をフーリエ変換
Claims (14)
- 信号f(a)をフーリエ変換することで情報信号F(b)を生成する情報信号生成方法において、
前記情報信号F(b)のピーク強度の最大値を検出する工程と、
前記情報信号F(b)のピーク強度の最大値に対応する前記信号f(a)の振幅と位相と周波数とを算出する工程と、
前記振幅、位相及び周波数情報に基づきa軸上に延長した信号f’(a)を生成する工程と、
a1 < a2として、前記延長した信号f’(a)からa1よりも小さい領域の信号と、前記延長した信号f’(a)のうちa2よりも大きい領域の信号とを抽出して外挿信号を生成する工程と、
前記外挿信号と前記a1からa2までの信号f(a)とを合成して合成信号を生成する工程と、
前記合成信号をさらにフーリエ変換する工程と、
を有することを特徴とする信号生成方法。 - 前記情報信号生成方法により得られた情報信号F(b)について、強度の第n番目(ここで、nは1以上の整数)のピーク値に次ぐ強度の第(n+1)番目のピーク値を検出する工程と、
前記第(n+1)番目のピーク値に対応する前記信号f(a)の振幅と位相と周波数とを算出する工程と、
前記振幅、位相及び周波数情報に基づきa軸上に延長した信号を生成する工程と、
a1 < a2として、前記延長した信号からa1よりも小さい領域の信号と、前記延長した信号のうちa2よりも大きい領域の信号とを抽出して外挿信号を生成する工程と、
前記外挿信号と前記a1からa2までの信号f(a)とを合成して合成信号を生成する工程と、
前記合成信号をさらにフーリエ変換する工程と、
をn=1からNまで繰り返す工程と、
を有することを特徴とする請求項1に記載の信号生成方法。 - 前記Nを前記情報信号F(b)のピークの数とすることを特徴とする請求項2に記載の信号生成方法。
- 前記Nは、大きい強度を示すピークから数えて所定の強度以下となるピークまでの数であることを特徴とする請求項2に記載の信号生成方法。
- 前記情報信号F(b)の強度に対して閾値を決める工程と、
前記情報信号F(b)で閾値を超える複数のピーク値を検出する工程と、
前記複数のピーク値のそれぞれに対応する前記信号f(a)の振幅と位相と周波数とを算出する工程と、
前記振幅、位相及び周波数情報に基づきa軸上に延長した信号を生成する工程と、
前記延長した信号からa1よりも小さい領域の信号と、前記延長した信号のうちa2よりも大きい領域の信号とを抽出して外挿信号を生成する工程と、
前記外挿信号と前記a1からa2までの信号f(a)とを合成して合成信号を生成する工程と、
前記合成信号をさらにフーリエ変換する工程と、
を有することを特徴とする請求項1に記載の信号生成方法。 - 前記閾値を所定の割合で下げる度に、
前記情報信号F(b)で閾値を超える複数のピーク値を検出する工程と、
前記複数のピーク値のそれぞれに対応する前記信号f(a)の振幅と位相と周波数とを算出する工程と、
前記振幅、位相及び周波数情報に基づきa軸上に延長した信号を生成する工程と、
前記延長した信号からa1よりも小さい領域の信号と、前記延長した信号のうちa2よりも大きい領域の信号とを抽出して外挿信号を生成する工程と、
前記外挿信号と前記a1からa2までの信号f(a)とを合成して合成信号を生成する工程と、
前記合成信号をさらにフーリエ変換する工程と、
を行うことを特徴とする請求項5に記載の信号生成方法。 - 前記閾値を決める工程は、f(a)のエンベロープ波形を取得し、該エンベロープ波形をフーリエ変換することで得られる信号の最大のピーク値とこれに次ぐ2番目のピーク値を検出する工程と、
前記2番目のピーク値と最大のピーク値の比を算出する工程とを有し、
前記比を前記情報信号F(b)の最大のピーク値に乗じて得られる値を、前記閾値とすることを特徴とする請求項5に記載の信号生成方法。 - 前記信号f(a)に、関数が示す波形の裾部が0となる窓関数を乗じ、前記延長した信号は前記窓関数が0となるaまで延長されたものであることを特徴とする請求項1に記載の信号生成方法。
- 前記信号f’(a)を延長する範囲は、前記延長した範囲の前記窓関数をフーリエ変換した信号の強度の最大のピーク値とこれに次ぐ2番目のピーク値との比が所定の値となる範囲以上であることを特徴とする請求項8に記載の信号生成方法。
- 前記信号f(a)は、該信号f(a)を該信号f(a)のエンベロープ情報で割ることで該信号f(a)の形状をフラットトップ形状に整形し、これに窓関数を乗じたものであることを特徴とする請求項1に記載の信号生成方法。
- 前記信号f(a)を光スペクトル干渉信号とし、前記情報信号F(b)を断層に関する情報信号として、
前記断層に関する情報信号のピーク強度の最大値を検出する工程と、
前記断層に関する情報信号のピーク強度の最大値に対応する前記光スペクトル干渉信号の振幅と位相と周波数とを算出する工程と、
前記振幅、位相及び周波数情報に基づき、前記a軸を波数軸として該波数軸上に延長した信号を生成する工程と、
波数ks < keとして、前記延長した信号からksよりも小さい領域の信号と、前記延長した信号のうちkeよりも大きい領域の信号とを抽出して外挿信号を生成する工程と、
前記外挿信号と前記ksからkeまでの光スペクトル干渉信号とを合成して合成信号を生成する工程と、
前記合成信号をさらにフーリエ変換する工程と、
を有することを特徴とする請求項1乃至10に記載の信号生成方法。 - 信号f(a)をフーリエ変換することで情報信号F(b)を生成する情報信号処理装置において、
前記情報信号F(b)のピーク強度の最大値を検出する処理部と、
前記情報信号F(b)のピーク強度の最大値に対応する前記信号f(a)の振幅と位相と周波数とを算出する処理部と、
前記振幅及び位相及び周波数情報に基づきa軸上に延長した信号を生成する処理部と、
a1 < a2として、前記延長した信号からa1よりも小さい領域の信号と、前記延長した信号のうちa2よりも大きい領域の信号とを抽出して外挿信号を生成する処理部と、
前記外挿信号と前記a1からa2までの信号f(a)とを合成して合成信号を生成する処理部と、
前記合成信号をさらにフーリエ変換する処理部と、
を有することを特徴とする信号生成装置。 - 前記信号f(a)を光スペクトル干渉信号、前記情報信号F(b)を断層に関する情報信号、
前記a軸を波数軸、とそれぞれすると共に、前記a1 < a2を、波数ks < keとした
ことを特徴とする請求項12に記載の信号生成装置。 - 光源部と、前記光源部からの光を検体に照射し、検体からの反射光を伝達させる検体測定部と、前記光源部からの光を参照光として伝達させる参照部と、前記反射光と前記参照光とを干渉させる干渉部と、前記干渉部からの干渉光を検出する光検出部と、前記光検出部で検出された光に基づいて前記検体の断層像を得る画像処理部と、を有する光干渉断層撮像装置であって、前記画像処理部は、請求項13に記載の断層に関する情報信号を生成する信号生成装置を備えることを特徴とする光干渉断層撮像装置。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013551079A JP5885758B2 (ja) | 2011-12-27 | 2011-12-27 | 情報信号生成方法 |
EP11878685.4A EP2799838B1 (en) | 2011-12-27 | 2011-12-27 | Information signal generating method |
CN201180076033.9A CN104011497B (zh) | 2011-12-27 | 2011-12-27 | 用于产生信息信号的方法 |
PCT/JP2011/080223 WO2013098942A1 (ja) | 2011-12-27 | 2011-12-27 | 情報信号生成方法 |
US13/721,403 US9600444B2 (en) | 2011-12-27 | 2012-12-20 | Method for generating information signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2011/080223 WO2013098942A1 (ja) | 2011-12-27 | 2011-12-27 | 情報信号生成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013098942A1 true WO2013098942A1 (ja) | 2013-07-04 |
Family
ID=48655395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/080223 WO2013098942A1 (ja) | 2011-12-27 | 2011-12-27 | 情報信号生成方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9600444B2 (ja) |
EP (1) | EP2799838B1 (ja) |
JP (1) | JP5885758B2 (ja) |
CN (1) | CN104011497B (ja) |
WO (1) | WO2013098942A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016188795A (ja) * | 2015-03-30 | 2016-11-04 | 株式会社Screenホールディングス | 信号処理方法および画像処理装置 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3109646A1 (de) * | 2015-06-23 | 2016-12-28 | Siemens Aktiengesellschaft | Verfahren zur analyse eines signals sowie vorrichtung zur durchführung des verfahrens |
DE102018200753A1 (de) * | 2018-01-18 | 2019-07-18 | Robert Bosch Gmbh | Verfahren und Vorrichtung zum Korrigieren eines Radarsignals und Radarvorrichtung |
US11357399B2 (en) * | 2018-08-02 | 2022-06-14 | Nidek Co., Ltd. | OCT apparatus |
CN114518162B (zh) * | 2022-01-24 | 2023-08-04 | 中国人民解放军海军工程大学 | 一种光纤水听器干涉信号强度补偿方法和系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63302307A (ja) * | 1987-06-02 | 1988-12-09 | Hitachi Ltd | 光学的膜厚測定方法 |
US5686922A (en) | 1995-09-29 | 1997-11-11 | Environmental Research Institute Of Michigan | Super spatially variant apodization (Super - SVA) |
JPH10510355A (ja) * | 1994-12-09 | 1998-10-06 | フォス エレクトリック アクティーゼルスカブ | 情報取得方法 |
JP2001008919A (ja) * | 1999-07-02 | 2001-01-16 | Hitachi Medical Corp | 磁気共鳴撮影装置及びケミカルシフト情報処理方法 |
JP2008216381A (ja) | 2007-02-28 | 2008-09-18 | Casio Comput Co Ltd | 音声分析合成装置、及びプログラム |
JP2008298767A (ja) * | 2007-05-02 | 2008-12-11 | Canon Inc | 光干渉断層計を用いた画像形成方法、及び光干渉断層装置 |
JP2009536740A (ja) * | 2006-05-10 | 2009-10-15 | ザ ジェネラル ホスピタル コーポレイション | サンプルの周波数領域画像形成を提供するためのプロセス、構成およびシステム |
JP2010223670A (ja) * | 2009-03-23 | 2010-10-07 | Toyama Univ | 光断層画像表示システム |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61111427A (ja) * | 1985-10-04 | 1986-05-29 | Hitachi Ltd | フ−リエ変換スペクトロメ−タの干渉信号処理方法 |
US5612784A (en) * | 1996-01-04 | 1997-03-18 | Bio-Rad Laboratories, Inc. | Digital signal processing technique for a FT-IR spectrometer using multiple modulations |
US6252668B1 (en) * | 1999-11-19 | 2001-06-26 | Zygo Corporation | Systems and methods for quantifying nonlinearities in interferometry systems |
EP1870029A1 (en) * | 2006-06-23 | 2007-12-26 | OPTOPOL Technology Spolka z o.o. | Apparatus and method for frequency domain optical coherence tomography |
JP5194963B2 (ja) * | 2008-04-03 | 2013-05-08 | 株式会社ニコン | 波形解析装置、波形解析プログラム、干渉計装置、パターン投影形状測定装置、及び波形解析方法 |
GB2472059B (en) * | 2009-07-23 | 2012-09-19 | Univ Loughborough | Apparatus for the absolute measurement of two dimensional optical path distributions using interferometry |
WO2011020092A1 (en) * | 2009-08-14 | 2011-02-17 | Pyrotek, Inc. | Waste heat system |
TWI417519B (zh) * | 2009-12-10 | 2013-12-01 | Ind Tech Res Inst | 干涉相位差量測方法及其系統 |
-
2011
- 2011-12-27 EP EP11878685.4A patent/EP2799838B1/en not_active Not-in-force
- 2011-12-27 WO PCT/JP2011/080223 patent/WO2013098942A1/ja active Application Filing
- 2011-12-27 CN CN201180076033.9A patent/CN104011497B/zh not_active Expired - Fee Related
- 2011-12-27 JP JP2013551079A patent/JP5885758B2/ja not_active Expired - Fee Related
-
2012
- 2012-12-20 US US13/721,403 patent/US9600444B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63302307A (ja) * | 1987-06-02 | 1988-12-09 | Hitachi Ltd | 光学的膜厚測定方法 |
JPH10510355A (ja) * | 1994-12-09 | 1998-10-06 | フォス エレクトリック アクティーゼルスカブ | 情報取得方法 |
US5686922A (en) | 1995-09-29 | 1997-11-11 | Environmental Research Institute Of Michigan | Super spatially variant apodization (Super - SVA) |
JP2001008919A (ja) * | 1999-07-02 | 2001-01-16 | Hitachi Medical Corp | 磁気共鳴撮影装置及びケミカルシフト情報処理方法 |
JP2009536740A (ja) * | 2006-05-10 | 2009-10-15 | ザ ジェネラル ホスピタル コーポレイション | サンプルの周波数領域画像形成を提供するためのプロセス、構成およびシステム |
JP2008216381A (ja) | 2007-02-28 | 2008-09-18 | Casio Comput Co Ltd | 音声分析合成装置、及びプログラム |
JP2008298767A (ja) * | 2007-05-02 | 2008-12-11 | Canon Inc | 光干渉断層計を用いた画像形成方法、及び光干渉断層装置 |
JP2010223670A (ja) * | 2009-03-23 | 2010-10-07 | Toyama Univ | 光断層画像表示システム |
Non-Patent Citations (1)
Title |
---|
See also references of EP2799838A4 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016188795A (ja) * | 2015-03-30 | 2016-11-04 | 株式会社Screenホールディングス | 信号処理方法および画像処理装置 |
Also Published As
Publication number | Publication date |
---|---|
JPWO2013098942A1 (ja) | 2015-04-30 |
EP2799838A1 (en) | 2014-11-05 |
EP2799838B1 (en) | 2016-07-27 |
CN104011497B (zh) | 2017-10-03 |
JP5885758B2 (ja) | 2016-03-15 |
US20130166239A1 (en) | 2013-06-27 |
CN104011497A (zh) | 2014-08-27 |
EP2799838A4 (en) | 2015-09-02 |
US9600444B2 (en) | 2017-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5680826B2 (ja) | 1以上のスペクトルを符号化する内視鏡技術によるデータ生成システム | |
US8457440B1 (en) | Method and system for background subtraction in medical optical coherence tomography system | |
US8504141B2 (en) | Optical tomographic image generating apparatus and optical tomographic image generating method | |
JP4389032B2 (ja) | 光コヒーレンストモグラフィーの画像処理装置 | |
JP5371315B2 (ja) | 光干渉断層撮像方法および光干渉断層撮像装置 | |
JP5885758B2 (ja) | 情報信号生成方法 | |
CN104568842B (zh) | 结合fdoct测量生物组织群延迟色散的装置及方法 | |
JP2006047264A (ja) | オプティカル・コヒーレント・トモグラフィー装置及びこれに用いる可変波長光発生装置並びに可変波長発光光源 | |
US10677581B2 (en) | Multiple reference OCT system | |
JP2012502262A (ja) | 半透明物質の画像化の方法および装置 | |
Cao et al. | Non-harmonic analysis applied to optical coherence tomography imaging | |
JP2008151697A (ja) | 光断層画像化装置 | |
US11118894B2 (en) | Systems, methods, and media for multiple reference arm spectral domain optical coherence tomography | |
CN107407601A (zh) | 用于补偿由光谱仪系统产生的干涉图的时间周期扰动的光谱仪系统和方法 | |
JP6214020B2 (ja) | 光断層イメージング法、その装置およびプログラム | |
JP2010223670A (ja) | 光断層画像表示システム | |
US20150173607A1 (en) | Method for acquiring optical tomographic image | |
JP5664564B2 (ja) | 光断層画像取得方法 | |
Rao et al. | Signal processing of spectral domain optical coherence tomography | |
JP5746741B2 (ja) | 画像生成装置、画像生成システム及び画像生成方法 | |
Liu et al. | Adaptive filtering of optical coherent tomography fringe data with ensemble empirical mode decomposition | |
JP2013057549A (ja) | 光断層画像取得方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11878685 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2013551079 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
REEP | Request for entry into the european phase |
Ref document number: 2011878685 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011878685 Country of ref document: EP |