Single snapshot multi-frequency demodulation method

Abstract

A single snapshot multi-frequency demodulation method for a modulated image obtained by modulating and summing one or more original components at different frequencies in a time domain or spatial domain, especially for a modulated image including multiple frequency components. AC and DC component values of each pixel at each frequency are extracted sequentially, and then an original AC and DC component image corresponding to each frequency is obtained. The method can be used in the time or spatial domain, can decompose multiple frequency component images using single measurement, has the advantages of fast speed, higher demodulation precision and good de-noising effect, meets the requirements for acquiring multiple pieces of frequency information at a time and overcomes inevitable errors in multiple measurements. Further, multiple pieces of image information can also be transmitted once using the demodulation method, so that parallel real-time transmission of the information in the communication field is realized.

Claims

1. A single multi-frequency snapshot demodulation method, comprising: in a time domain or spatial domain, for a modulated image including one or more frequencies, firstly extracting AC and DC component values of each pixel at each frequency in sequence, and then obtaining an original AC and DC component image corresponding to each frequency for the modulated image, selecting a T.sub.1T.sub.2 matrix as a kernel, wherein T1 and T2 are respectively least common multiples of transverse and longitudinal components in each AC component cycle, then in the image, selecting kernel units having the kernel size from the image by using each pixel as the upper left corner of the matrix in a predetermined order, and performing integral summation on the kernel units and corresponding patterns respectively.

2. The method of claim 1, wherein the integral summation further comprises: multiplying each kernel unit in the image by cosinoidal and sinusoidal kernel patterns at the same frequency and in the same direction, and then performing integral summation to obtain cosinoidal and sinusoidal modulated patterns having the same sizes as the original image.

3. The method of claim 1, wherein for a modulated image having AC components at different frequencies f.sub.1, f.sub.2, . . . , the relation between the components is as shown in formula (1): $\begin{array}{cc}f\left(x,y\right)=\underset{i=1}{\overset{k}{.Math.}}{A}_i\cos \left(2f_i.Math.\left(m_i.Math.x+n_i.Math.y\right)+{}_i\right)+B& \left(1\right)\\ \sqrt{m_i^2+n_i^2}=1& \left(2\right)\end{array}$ wherein, k1 is the number of AC modulation components, A.sub.i, f.sub.i and .sub.i are respectively the amplitude, frequency and spatial initial phase of each AC component, m.sub.i[1,1] and n.sub.i[1,1] jointly determine the direction of AC modulation stripes and satisfy formula (2), B is a DC component, and x and y are spatial coordinates.

4. The method of claim 3, wherein the amplitude of the AC component of the modulated pattern is further demodulated using formula (3), and the DC component is demodulated using formula (4): $\begin{array}{cc}{A}_i=\frac{\sqrt{\begin{array}{c}{\left[\underset{}{}f\left(x,y\right).Math.\cos \left(2f_i.Math.\left(m_i.Math.x+n_i.Math.y\right)\right)\mathrm{dxdy}\right]}^2+\\ {\left[\underset{}{}f\left(x,y\right).Math.\sin \left(2f_i.Math.\left(m_i.Math.x+n_i.Math.y\right)\right)\mathrm{dxdy}\right]}^2\end{array}}}{\underset{}{}{\cos }^2\left(2f_i.Math.\left(m_i.Math.x+n_i.Math.y\right)\right)\mathrm{dxdy}}& \left(3\right)\\ B=\frac{1}{T_1{T}_2}\underset{}{}f\left(x,y\right)\mathrm{dxdy}& \left(4\right)\end{array}$ wherein, is a kernel unit corresponding to each pixel in the image, and T.sub.1 and T.sub.2 are respectively least common multiples of transverse and longitudinal components in each AC component cycle.

5. The method of claim 1, wherein when the AC components of the image at different frequencies or in different directions are extracted, the kernel patterns are also different, and have the same frequencies and directions as the extracted AC modulated patterns.

6. The method of claim 1, wherein the predetermined order is from left to right and from top to bottom, starting from the upper left corner of the image.

7. A spatial frequency domain imaging method, comprising projecting modulated light including one or more frequencies onto a sample, and collecting the light reflected from the sample, transmitted light or excited fluorescent light by a CCD camera, wherein the light collected by the CCD camera is decomposed into a DC portion and an AC portion; and the method is characterized in that the light collected by the CCD camera is demodulated using the single multi-frequency snapshot demodulation method of claim 1.

8. A signal transmission method, comprising the modulation and demodulation method for spatial frequency domain information of claim 1.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) To describe the technical solutions of the embodiments of the present disclosure more clearly, here is a brief introduction of the accompanying drawings necessary for the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only some exemplary embodiments of the present disclosure.

(2) FIG. 1 is a diagram for demodulating an image including two spatial frequencies using a single snapshot multi-frequency demodulation method.

(3) FIG. 2 is a contrast diagram of demodulation on a pattern including one frequency using the three-phase shifting standard method, the Hilbert transform method and the single snapshot multi-frequency demodulation method.

(4) FIG. 3 is a contrast diagram of demodulation on a pattern including two frequencies using the Hilbert transform method and the single snapshot demodulation method.

(5) FIG. 4 is a contrast diagram of demodulation on a pattern including three frequencies using a single snapshot demodulation method.

DETAILED DESCRIPTION

(6) A clear and complete description will be made to the technical solutions in the embodiments of the present disclosure below in combination with the drawings of the present disclosure. Obviously, the embodiments described are only part of the embodiments of the present disclosure, not all of them. All other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present disclosure without creative efforts shall fall within the protection scope of the present disclosure.

(7) The following introduces a single snapshot multi-frequency demodulation method, taking spatial frequency domain imaging as an example. Compared with the standard three-phase shifting method, the method has the advantages of fast speed, higher demodulation precision and good de-noising effect, meets the requirements for acquiring multiple pieces of frequency information at a time and overcomes inevitable errors in multiple measurements.

(8) The single multi-frequency snapshot demodulation method is to extract AC and DC component values of each pixel at each frequency for a modulated image obtained by modulating and summing one or more original components at different frequencies in a spatial domain, and then to obtain an original AC and DC component image corresponding to each frequency. As shown in FIG. 1, an image including two spatial frequencies is demodulated using the single snapshot multi-frequency demodulation method, and the image measured for a single time includes components of two different spatial frequencies (f.sub.1, f.sub.2). The single snapshot multi-frequency demodulation method can accurately demodulate spatial modulated images of different frequencies.

(9) The specific steps are as follows:

(10) For an image having AC components of two different frequencies, if the frequencies are respectively f.sub.1 and f.sub.2, the relation between the components is shown in formula (10). In order to extract the AC components of the different frequencies in the image, firstly a T.sub.1T.sub.2 matrix is selected as a kernel, wherein T.sub.1 and T.sub.2 are respectively least common multiples of transverse and longitudinal components in each AC component cycle; then, kernel units having the kernel size are selected from the image by using each pixel as the upper left corner of the matrix in sequence from left to right and from top to bottom, and integral summation is performed on the kernel units and corresponding kernel patterns respectively. Further, when the AC components of the image at different frequencies or in different directions are extracted, the kernel patterns are also different, but should be ensured to have the same frequencies and directions as the extracted components. Each kernel unit in the image is multiplied by cosinoidal and sinusoidal patterns at the same frequency and in the same direction, and then integral summation is performed to obtain cosinoidal and sinusoidal modulated patterns having the same sizes as the original patterns.

(11) Then, the amplitude of the AC component of the modulated pattern is demodulated using formulas (11) and (12), and the DC component is demodulated using formula (13).
f(x,y)=a cos(2f.sub.1.Math.x+.sub.1)+b cos(2f.sub.2.Math.(x+y)+.sub.2)+c(10)

(12) In formula (10), a and b are the amplitudes of AC components, f.sub.1 and f.sub.2 are frequencies of AC modulated patterns, .sub.1 and .sub.2 are spatial initial phases, c is a DC component, and x and y are spatial coordinates.

(13) $\begin{array}{cc}a=\frac{\sqrt{\begin{array}{c}{\left[\underset{}{}f\left(x,y\right)\cos \left(2f_1.Math.x\right)\mathrm{dxdy}\right]}^2+\\ {\left[\underset{}{}f\left(x,y\right)\sin \left(2f_1.Math.x\right)\mathrm{dxdy}\right]}^2\end{array}}}{\underset{}{}{\cos }^2\left(2f_1.Math.x\right)\mathrm{dxdy}}& \left(11\right)\\ b=\frac{\sqrt{\begin{array}{c}{\left[\underset{}{}f\left(x,y\right)\cos \left(2f_2.Math.\left(x+y\right)\right)\mathrm{dxdy}\right]}^2+\\ {\left[\underset{}{}f\left(x,y\right)\sin \left(2f_2.Math.\left(x+y\right)\right)\mathrm{dxdy}\right]}^2\end{array}}}{\underset{}{}{\cos }^2\left(2f_2.Math.\left(x+y\right)\right)\mathrm{dxdy}}& \left(12\right)\\ c=\frac{1}{T_1{T}_2}\underset{}{}f\left(x,y\right)\mathrm{dxdy}& \left(13\right)\end{array}$

(14) Wherein, is a kernel unit corresponding to each pixel in the image, and T.sub.1 and T.sub.2 are respectively least common multiples of transverse and longitudinal components in each AC component cycle.

(15) Further, for a modulated image having AC components at different frequencies f.sub.1, f.sub.2, . . . , the relation between the components is as shown in formula (14):

(16) $\begin{array}{cc}f\left(x,y\right)=\underset{i=1}{\overset{k}{.Math.}}{A}_i\cos \left(2f_i.Math.\left(m_i.Math.x+n_i.Math.y\right)+{}_i\right)+B& \left(14\right)\end{array}$

(17) Wherein, k1 is the number of AC modulation components, A.sub.i, f.sub.i and .sub.i are respectively the amplitude, frequency and spatial initial phase of each AC component, m.sub.i[1,1] and n.sub.i[1,1] jointly determine the direction of AC modulation stripes, B is a DC component, and x and y are spatial coordinates. In addition, the amplitude of the AC component of the modulated pattern is demodulated using formula (15), and the DC component is demodulated using formula (16):

(18) $\begin{array}{cc}{A}_i=\frac{\sqrt{\begin{array}{c}{\left[\underset{}{}f\left(x,y\right).Math.\cos \left(2f_i.Math.\left(m_i.Math.x+n_i.Math.y\right)\right)\mathrm{dxdy}\right]}^2+\\ {\left[\underset{}{}f\left(x,y\right).Math.\sin \left(2f_i.Math.\left(m_i.Math.x+n_i.Math.y\right)\right)\mathrm{dxdy}\right]}^2\end{array}}}{\underset{}{}{\cos }^2\left(2f_i.Math.\left(m_i.Math.x+n_i.Math.y\right)\right)\mathrm{dxdy}}& \left(15\right)\\ B=\frac{1}{T_1{T}_2}\underset{}{}f\left(x,y\right)\mathrm{dxdy}& \left(16\right)\end{array}$

(19) Wherein, is a kernel unit corresponding to each pixel in the image, and T.sub.1 and T.sub.2 are respectively least common multiples of transverse and longitudinal components in each AC component cycle.

(20) Contrast Experiment:

(21) FIG. 2 shows the results of demodulation on a pattern including one frequency using the three-phase shifting standard method, the Hilbert transform method and the single snapshot multi-frequency demodulation method respectively. The original image includes three components, which are respectively a DC component, an AC component and a noise signal. The advantages and disadvantages of the methods will be compared through several different demodulation methods. According to the results obtained, the three methods can all fulfill the purpose of demodulating AC signals from noise and DC signals. However, regarding the noise suppression effect, the single snapshot multi-frequency demodulation method is superior to the three-phase shifting standard method and the Hilbert transform method, and can well suppress noise interference.

(22) As shown in FIG. 3, a pattern including two frequencies is demodulated using the Hilbert transform method and the single snapshot demodulation method respectively: in FIG. 3(a), the original image is obtained by transversely modulating (a Chinese character) (f.sub.y=0.2), longitudinally modulating (a Chinese character) (f.sub.x=0.1), then superposing and , and adding a DC component and noise signals. Since it is known that the three-phase shifting standard method can only demodulate a single-frequency AC component, only the Hilbert transform method and the single snapshot multi-frequency demodulation method are compared here. FIG. 3(b) shows the result of demodulation on FIG. 3(a) using the Hilbert transform method, and it can be seen from the result that the demodulated AC components cannot be used for separating signals of different frequencies when the AC components have two or more frequencies. However, it can be seen from FIGS. 3(c) and 3(d) that the single snapshot multi-frequency demodulation method can well demodulate AC modulation signals of different frequencies.

(23) As shown in FIG. 4, a pattern including three frequencies is demodulated using the single snapshot demodulation method respectively: in FIG. 4(a), the original image is obtained by transversely modulating (a Chinese character) (f.sub.y=0.1), longitudinally modulating (f.sub.x=0.25), obliquely modulating 45 (f.sub.x=0.25, f.sub.y=0.25), then superposing , and , and adding a DC component and noise signals. Since it is known that the three-phase shifting standard method can only demodulate a single-frequency AC component and the Hilbert transform method cannot distinguish multiple frequency modulation signals, FIG. 4 only shows the results obtained by the single snapshot multi-frequency demodulation method. From FIG. 4(b), FIG. 4(c) and FIG. 4(d), it can be seen that the single snapshot multi-frequency demodulation method can well demodulate AC modulation signals of different frequencies.

(24) From the experimental results, it can be clearly obtained that the single snapshot multi-frequency demodulation method not only can well demodulate a single-frequency AC component from noise and suppress noise interference, but also can demodulate AC components of two or more frequencies. In this way, real-time imaging of the MTF function at different spatial frequencies can be realized, the optical parameters of tissues can be quickly resolved, and the problem of noise interference caused by multiple times of imaging is thus solved.

(25) According to the single snapshot multi-frequency demodulation method provided by the above embodiments of the present disclosure, for a modulated image obtained by modulating and summing one or more original components at different frequencies in a spatial domain, especially for a modulated image including multiple frequency components, AC and DC component values of each pixel at each frequency are extracted sequentially through the method, and then an original AC and DC component image corresponding to each frequency is obtained. The method has the advantages of fast speed, higher demodulation precision and good de-noising effect, meets the requirements for acquiring multiple pieces of frequency information at a time and overcomes inevitable errors in multiple measurements. At the same time, multiple pieces of image information can also be transmitted at a time using the demodulation method, so that parallel real-time transmission of the information in the communication field is realized.

(26) It could be understood by a person skilled in the art that the steps, measures or schemes of various operations, methods or processes discussed in the present disclosure can be alternated, changed, combined or omitted. Further, other steps, measures or schemes of various operations, methods or processes discussed in the present disclosure can also be alternated, changed, rearranged, decomposed, combined or omitted. Further, the steps, measures or schemes of various operations, methods or processes of the prior art, which are the same as those in the present disclosure, can also be alternated, changed, rearranged, decomposed, combined or omitted.

(27) Described above are only part of the embodiments of the present disclosure. It should be pointed out that a person of ordinary skill in the art may further make improvements and adjustments without departing from the principle of the present disclosure, and these improvements and adjustments fall within the protection scope of the present disclosure.