Block rapid measurement system and method for phase contrast optical coherence elastography

Abstract

The invention discloses a body block rapid measurement system and method oriented to phase contrast optical coherence elastography. The method comprises the following steps: acquiring interference signals before and after deformation of a body block; performing field-of-view matching on the two industrial cameras; according to the matched same view field area, data clipping is carried out on the collected interference signals, and analysis signals of interference data before and after deformation are constructed; carrying out compressed sensing and reconstruction on the constructed analysis signal; fourier transform is carried out on the reconstructed signals, argument main values of the reconstructed signals are calculated, wrapped phase fields before and after deformation are obtained, and finally differential wrapped phases are calculated. According to the method, the analytic signal is constructed, the characteristic that the analytic signal is a unilateral spectrum after Fourier transform is utilized, and the advantage that the reconstruction precision of the compressed sensing for the signal with better sparsity is higher is combined, so that the sampling data volume is reduced, the deformation measurement speed of the body block is greatly improved, and meanwhile, the reconstruction precision of the signal can be improved when the height is under-sampled; and the obtained differential wrapped phase is clearer.

Description

Technical Field

[0001] The present invention relates to the technical field of optical imaging measurement, and in particular to a rapid volume measurement system and method for phase-contrast optical coherence elastic imaging. Background Art

[0002] Phase-Contrast Optical Coherence Elastography (PC-OCE) utilizes the high sensitivity of tomographic interferometry phase changes to achieve full-field, high-precision measurement of the internal deformation distribution of materials, making it the most promising method for characterizing internal mechanical properties.

[0003] The PC-OCE measurement steps generally include:

[0004] 1) Collect interference spectrum data before and after deformation of the block;

[0005] 2) Perform Fourier transform on the interference spectrum data before and after deformation along the wavenumber axis to demodulate the respective wrapped phases before and after deformation;

[0006] 3) Calculate the numerical difference of the two wrapped phases and use the wrapping operator to wrap the difference of the wrapped phases into [-π,π];

[0007] 4) The differential wrapped phase in step 3) is phase-unwrapped to obtain the deformation field distribution of the block.

[0008] In the PC-OCT system, the system depth range and resolution are:

[0009]

[0010] where δ z is the depth resolution, limited by the light source bandwidth Δλ and the light source center wavelength λ c ; Δ z It is limited by the sampling rate N of the interferometry system along the spectral axis, that is, the number of frames captured by the camera.

[0011] To meet the Nyquist sampling criteria and maintain measurement depth without sacrificing resolution, sufficient interferometric spectral data must be collected along the wavenumber axis. When measuring a block using PC-OCE, a large amount of data, both before and after deformation, must be collected in a short period of time. This makes it difficult for the camera to meet the required frame rate while maintaining measurement depth without sacrificing imaging resolution, thus limiting the speed of PC-OCE volume measurement.

[0012] At present, in order to ensure that resolution is not sacrificed and the measurement depth requirements are met while complying with the spectral Nyquist sampling theorem, existing technical solutions can be mainly divided into hardware solutions and software solutions. First, hardware solution: using a high-frame rate camera as the acquisition device, by increasing the number of shots per second, and subsequently processing the large amount of collected data using CPU thread parallel processing, combined with GPU acceleration to achieve full-process parallelization and increase data processing speed; Second, software solution: traditional compressed sensing method, randomly sampling real-number interference signals, first reducing the collection of a part of the data, and then restoring the original range of the signal through the compressed sensing algorithm, and finally calculating the differential wrapped phase after Fourier transform, to achieve a reduction in the amount of data that needs to be sampled.

[0013] While hardware solutions can ensure sufficient measurement range, they significantly increase the need for high-frame-rate cameras and computing hardware, significantly increasing system costs. Furthermore, the data volume remains enormous, making system integration difficult. Regarding software solutions, while traditional compressed sensing can reduce some data volume, it still cannot achieve a significant reduction. When undersampling is high, traditional compressed sensing methods cannot achieve good recovery results, resulting in blurred differential wrapped phases.

[0014] Therefore, the existing technology still has significant limitations. In order to improve the speed of data processing and reconstruction accuracy, it is necessary to explore a more effective phase-contrast optical coherence elastography rapid volume measurement system and method. Summary of the Invention

[0015] The purpose of the present invention is to overcome the deficiencies of the prior art and to provide a system and method for rapid volume measurement for phase-contrast optical coherence elastography.

[0016] To achieve the above objectives, the technical solutions provided by the present invention are:

[0017] A rapid volume measurement system for phase-contrast optical coherence elastography, comprising a tunable semiconductor wavelength scanning laser, a first 50:50 beam splitter prism, a second 50:50 beam splitter prism, a polarization beam splitter prism, a first 45° polarizer, a second 45° polarizer, a third 45° polarizer, a first industrial camera, a second industrial camera, a collimating lens, a quarter-wave plate, a first reflector, a second reflector, and a computer;

[0018] Wherein, the first 45° polarizer is provided between the collimating lens and the first 50:50 beam splitting prism;

[0019] The tunable semiconductor wavelength scanning laser is arranged on the side of the collimating lens away from the first 45° polarizer;

[0020] The first reflector and the second reflector are respectively placed in two directions of light splitting by the first 50:50 beam splitting prism;

[0021] The second 50:50 beam splitter prism is arranged on the light path after being reflected by the first reflector and on the light path after being reflected by the second reflector;

[0022] The 1 / 4 wave plate is arranged between the first reflector and the second 50:50 beam splitter prism;

[0023] The second 45° polarizer is provided between the first reflector and the quarter wave plate;

[0024] The third 45° polarizer is provided between the second reflector and the second 50:50 beam splitter prism;

[0025] The measured block is placed between the second reflector and the third 45° polarizer;

[0026] The polarization beam splitter prism is arranged in the light combining direction of the second 50:50 beam splitter prism;

[0027] The first industrial camera and the second industrial camera are respectively arranged in two directions of light splitting by the polarization beam splitter prism, and the lenses of the two cameras are facing the polarization beam splitter prism;

[0028] The computer is connected to the first industrial camera and the second industrial camera, receives images captured by the first industrial camera and the second industrial camera, and analyzes the images.

[0029] To achieve the above objectives, the present invention further provides a method for rapid volume measurement for phase-contrast optical coherence elastography, which is implemented using the above-mentioned rapid volume measurement system for phase-contrast optical coherence elastography, comprising:

[0030] The interference signals before and after the deformation of the block are collected by a first industrial camera and a second industrial camera in the block rapid measurement system for phase contrast optical coherence elastic imaging;

[0031] The computer performs the following operations:

[0032] Matching the fields of view of the first industrial camera and the second industrial camera;

[0033] According to the matching identical field of view areas, the interference signals collected by the first industrial camera and the second industrial camera are trimmed to obtain the interferometric data of the block before and after deformation in the same field of view, and the analytical signals of the interferometric data before and after deformation are constructed.

[0034] Perform compressed sensing on the constructed analytical signal and reconstruct it;

[0035] The Fourier transform is performed on the compressed sensed and reconstructed signals, and their principal values ​​of the argument are calculated to obtain the wrapped phase field before and after deformation, and finally the differential wrapped phase is calculated.

[0036] Furthermore, the process of collecting interference signals before and after the block deformation includes:

[0037] The light emitted by the tunable semiconductor wavelength scanning laser is connected to the optical path through an optical fiber. The light emitted is collimated into parallel light by a collimating lens, and then passes through the first 45° polarizer to become linearly polarized light and enters the first 50:50 splitter prism, and is equally divided into light rays emitting in two directions; the light emitting in the positive direction of the x-axis is reflected by the first reflector and then passes through the second 45° polarizer and the 1 / 4 wave plate to become circularly polarized light; the light along the negative direction of the y-axis is reflected by the second reflector, first passes through the block, and then passes through the third 45° polarizer to ensure that the light emitting to the second 50:50 splitter prism is linearly polarized; then the linearly polarized light and the circularly polarized light from the block are recombined at the second 50:50 splitter prism, and then separated into two orthogonal components by the polarization splitter prism, which are received by the first industrial camera and the second industrial camera, and finally the collected interference signal is transmitted to the computer 14.

[0038] Furthermore, matching the fields of view of the first industrial camera and the second industrial camera includes:

[0039] Using a shooting resolution board, an area in the image captured by the first industrial camera is selected as a template image, and the image captured by the second industrial camera is used as a reference image. The degree of matching between the template image and the reference image coverage area is calculated using a field of view matching algorithm. The template window is then translated from left to right and from top to bottom in the reference image, traversing the entire reference image pixel by pixel. Finally, the sub-image with the highest matching degree is found in the entire reference image as the matching result, achieving the field of view matching effect.

[0040] The formula of the field of view matching algorithm is as follows:

[0041]

[0042] Among them, NCC(u,v) is the normalized cross-correlation coefficient, T(x,y) is the pixel value of the template image at the coordinate (x,y), and I(x+u,y+v) is the pixel value of the target image at the coordinate (x+u,y+v). is the pixel mean of the template image, is the pixel mean of the target image corresponding to the window at the displacement (u, v).

[0043] Furthermore, constructing analytical signals of the interference data before and after deformation includes:

[0044] The interference signals collected by the first industrial camera and the second industrial camera are respectively expressed as:

[0045]

[0046] Where M is the number of surfaces involved in the interference; I R , I m is the DC component, Λ m is the optical path difference between the two beams of light that pass through the block and directly enter the first industrial camera and the second industrial camera, Δk is the laser wave number scanning range; k0 is the starting wave number of the light source, is the frequency generated when the two beams of light pass through the block and directly enter the first industrial camera and the second industrial camera interfere with each other, is the initial phase when the light directly incident on the first industrial camera and the second industrial camera interferes with the light passing through the block;

[0047] By eliminating the DC component I R , I m After that, two interference signals I x , I y , and then take the two signals as the real and imaginary parts of the analytical signal:

[0048] I(t)=I x (t)+I y (t)·j (6)

[0049] According to the matching identical field of view area, the interference signals collected by the first industrial camera and the second industrial camera are trimmed to obtain the real and imaginary parts of the analytical signals, thereby constructing the analytical signals of the interference data before and after deformation.

[0050] Furthermore, the constructed analytical signal is subjected to compressed sensing and reconstructed, including:

[0051] The constructed analytical signal is compressed and sensed. The premise of compressed sensing is to represent the signal sparsely. R N A set of vector bases in the space, forming a basis matrix Then for any interference data x∈R N Are expressed as:

[0052]

[0053] Where α is the sparse representation of x in the γ domain. For volume PC-OCE data, the Fourier basis is used as the sparse representation;

[0054] Then, a measurement matrix that is unrelated to the Fourier basis is used to project the block PC-OCE data to obtain the observation value y. Finally, x is reconstructed by solving the optimization problem. The expression is:

[0055] y=φx=φγα (8)

[0056] Where φγ is the sensing matrix; φ is the measurement matrix, which is also a random matrix;

[0057] By solving the NP-hard problem, the original signal x is reconstructed.

[0058] During PC-OCE volume measurement, the sheer volume of data leads to slow data processing and significant hardware strain. This significantly impacts the processing difficulty of PC-OCE volume measurements. Traditional compressed sensing methods can be used to reduce the amount of sampled data. However, in practice, this doesn't significantly reduce the number of samples collected from real-valued signals, and the differential wrapper phase calculated from the reconstructed signal becomes blurred.

[0059] Compared with the existing technology, the principles and advantages of this technical solution are as follows:

[0060] This technical solution directly constructs the analytical signal and utilizes its characteristic of a single-sided spectrum after Fourier transform. At the same time, it combines the advantage of compressed sensing that the better the reconstruction accuracy of the signal, the greater the sparsity. This effectively reduces the amount of sampled data and greatly improves the deformation measurement speed of the PC-OCE block. At the same time, it can also improve the reconstruction accuracy of the signal in the case of high undersampling, and the obtained differential wrapped phase is clearer. BRIEF DESCRIPTION OF THE DRAWINGS

[0061] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the services required for use in the embodiments or the prior art descriptions will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0062] Figure 1 Schematic diagram of a rapid volume measurement system for phase-contrast optical coherence elastography according to the present invention;

[0063] Figure 2 This is a flow chart of the method for rapid volume measurement for phase-contrast optical coherence elastography according to the present invention;

[0064] Figure 3 A comparison chart of the PC-OCE block data sections.

[0065] Reference numerals:

[0066] 1-tunable semiconductor wavelength scanning laser; 2-first 50:50 beam splitter prism; 3-second 50:50 beam splitter prism; 4-polarization beam splitter prism; 5-first 45° polarizer; 6-second 45° polarizer; 7-third 45° polarizer; 8-first industrial camera; 9-second industrial camera; 10-collimating lens; 11-1 / 4 wave plate; 12-first reflector; 13-second reflector; 14-computer; 15-block. DETAILED DESCRIPTION

[0067] The present invention will be further described below in conjunction with specific embodiments:

[0068] like Figure 1 As shown, the volume rapid measurement system for phase contrast optical coherence elastography described in this embodiment includes a tunable semiconductor wavelength scanning laser 1, a first 50:50 beam splitter prism 2, a second 50:50 beam splitter prism 3, a polarization beam splitter prism 4, a first 45° polarizer 5, a second 45° polarizer 6, a third 45° polarizer 7, a first industrial camera 8, a second industrial camera 9, a collimating lens 10, a 1 / 4 wave plate 11, a first reflector 12, a second reflector 13 and a computer 14.

[0069] Among them, the first 45° polarizer 5 is arranged between the collimating lens 10 and the first 50:50 beam splitting prism 2; the tunable semiconductor wavelength scanning laser 1 is arranged on the side of the collimating lens 10 away from the first 45° polarizer 5; the first reflector 12 and the second reflector 13 are respectively placed in the two directions of the light splitting of the first 50:50 beam splitting prism 2; the second 50:50 beam splitting prism 3 is arranged on the light path after the first reflector 12 is reflected and on the light path after the second reflector 13 is reflected; the 1 / 4 wave plate 11 is arranged between the first reflector 12 and the second 50:50 beam splitting prism 3; the second 45° polarizer 6 is arranged between the first reflector 12 and The first industrial camera 8 and the second industrial camera 9 are respectively arranged in the two directions of the light splitting of the polarization beam splitting prism 4 and the lenses of both cameras are facing the polarization beam splitting prism 4; the computer 14 is connected to the first industrial camera 8 and the second industrial camera 9 to receive the images captured by the first industrial camera 8 and the second industrial camera 9 and analyze the images.

[0070] like Figure 2 As shown, the working principle of this embodiment is as follows:

[0071] S1. Collect interference signals of the block 15 before and after deformation by the first industrial camera 8 and the second industrial camera 9. The process is as follows:

[0072] The light emitted by the tunable semiconductor wavelength scanning laser 1 is connected to the optical path through an optical fiber. The light emitted is collimated into parallel light by a collimating lens 10, and then passes through the first 45° polarizer 5 to become linearly polarized light and enters the first 50:50 beam splitter 2, and is equally divided into light rays emitting in two directions; the light emitting in the positive direction of the x-axis is reflected by the first reflector 12 and then passes through the second 45° polarizer 6 and the 1 / 4 wave plate 11 to become circularly polarized light; the light along the negative direction of the y-axis is reflected by the second reflector 13, first passes through the block 15, and then passes through the third 45° polarizer 7 to ensure that the light emitting to the second 50:50 beam splitter prism 3 is linearly polarized; then the linearly polarized light and the circularly polarized light from the block 15 are recombined at the second 50:50 beam splitter prism 3, and then separated into two orthogonal components by the polarization beam splitter 4, which are received by the first industrial camera 8 and the second industrial camera 9, and finally the collected interference signal is transmitted to the computer 14.

[0073] S2. The computer 14 performs field of view matching on the first industrial camera 8 and the second industrial camera 9. The matching process includes:

[0074] Using a shooting resolution board, an area is selected in the image captured by the first industrial camera 8 as a template image, and the image captured by the second industrial camera 9 is used as a reference image. The degree of matching between the template image and the reference image coverage area is calculated using a field of view matching algorithm. The template window is then translated from left to right and from top to bottom in the reference image, traversing the entire reference image pixel by pixel. Finally, the sub-image with the highest matching degree is found in the entire reference image as the matching result, thereby achieving the effect of field of view matching.

[0075] The formula of the field of view matching algorithm is as follows:

[0076]

[0077] Among them, NCC(u,v) is the normalized cross-correlation coefficient, T(x,y) is the pixel value of the template image at the coordinate (x,y), and I(x+u,y+v) is the pixel value of the target image at the coordinate (x+u,y+v). is the pixel mean of the template image, is the pixel mean of the target image corresponding to the window at the displacement (u, v).

[0078] S3. Based on the matched identical field of view, data clipping is performed on the interference signals collected by the first industrial camera 8 and the second industrial camera 9 to obtain interference data of the block 15 before and after deformation in the identical field of view, and analytical signals of the interference data before and after deformation are constructed.

[0079] This step specifically includes:

[0080] The interference signals collected by the first industrial camera 8 and the second industrial camera 9 are respectively expressed as:

[0081]

[0082] Where M is the number of surfaces involved in the interference; I R , I m is the DC component, Λ m is the optical path difference between the two beams of light that pass through the block 15 and directly enter the first industrial camera 8 and the second industrial camera 9, Δk is the laser wave number scanning range; k0 is the starting wave number of the light source, is the frequency generated when the two beams of light that pass through the block 15 and directly enter the first industrial camera 8 and the second industrial camera 9 interfere with each other, is the initial phase when the light directly incident on the first industrial camera 8 and the second industrial camera 9 interferes with the light passing through the block 15;

[0083] By eliminating the DC component I R , I m After that, two interference signals I x , I y , and then take the two signals as the real and imaginary parts of the analytical signal:

[0084] I(t)=I x (t)+I y (t)·j (6)

[0085] According to the matched identical field of view area, the interference signals collected by the first industrial camera 8 and the second industrial camera 9 are data clipped to obtain the real and imaginary parts of the analytical signals, thereby constructing analytical signals of the interference data before and after deformation.

[0086] S4, performing compressed sensing on the constructed analytical signal and reconstructing it;

[0087] This step specifically includes:

[0088] The constructed analytical signal is compressed and sensed. The premise of compressed sensing is to represent the signal sparsely. It is R N A set of vector bases in the space, forming a basis matrix Then for any interference data x∈R N Are expressed as:

[0089]

[0090] Where α is the sparse representation of x in the γ domain. For the volume 15PC-OCE data, the Fourier basis is used as the sparse representation;

[0091] Then, the measurement matrix that is unrelated to the Fourier basis is used to project the 15PC-OCE data of the block to obtain the observation value y. Finally, x is reconstructed by solving the optimization problem. The expression is:

[0092] y=φx=φγα (8)

[0093] Where A=φγ is the sensing matrix; φ is the measurement matrix, which is also a random matrix;

[0094] By solving the NP-hard problem, the original signal x is reconstructed.

[0095] S5. Perform Fourier transform on the compressed sensed and reconstructed signals, calculate their principal values ​​of the argument angles to obtain the wrapped phase field before and after deformation, and finally calculate the differential wrapped phase.

[0096] Figure 3 The figure shows one of the cross-sections of the block 15 data. From left to right, the compressed sensing sampling rate is gradually reduced, namely 2, 3, and 4 times undersampling. It can be seen that the differential wrapped phase of the traditional compressed sensing method will become more and more blurred. After using the present invention, it will be significantly better than the traditional compressed sensing method when four times undersampling (that is, only using 1 / 4 of the original data), greatly improving the data processing speed.

[0097] The embodiments described above are only preferred embodiments of the present invention and are not intended to limit the scope of implementation of the present invention. Therefore, any changes made based on the shape and principle of the present invention should be included in the scope of protection of the present invention.

Claims

1. A rapid volume measurement system for phase-contrast optical coherence elastography, characterized by: The device comprises a tunable semiconductor wavelength scanning laser, a first 50:50 beam splitter prism, a second 50:50 beam splitter prism, a polarization beam splitter prism, a first 45° polarizer, a second 45° polarizer, a third 45° polarizer, a first industrial camera, a second industrial camera, a collimating lens, a quarter wave plate, a first reflector, a second reflector, and a computer; Wherein, the first 45° polarizer is provided between the collimating lens and the first 50:50 beam splitting prism; The tunable semiconductor wavelength scanning laser is arranged on the side of the collimating lens away from the first 45° polarizer; The first reflector and the second reflector are respectively placed in two directions of light splitting by the first 50:50 beam splitting prism; The second 50:50 beam splitter prism is arranged on the light path after being reflected by the first reflector and on the light path after being reflected by the second reflector; The 1 / 4 wave plate is arranged between the first reflector and the second 50:50 beam splitter prism; The second 45° polarizer is provided between the first reflector and the quarter wave plate; The third 45° polarizer is provided between the second reflector and the second 50:50 beam splitter prism; The measured block is placed between the second reflector and the third 45° polarizer; The polarization beam splitter prism is arranged in the light combining direction of the second 50:50 beam splitter prism; The first industrial camera and the second industrial camera are respectively arranged in two directions of light splitting by the polarization beam splitter prism, and the lenses of the two cameras are facing the polarization beam splitter prism; The computer is connected to the first industrial camera and the second industrial camera, receives images captured by the first industrial camera and the second industrial camera, and analyzes the images.

2. A rapid volume measurement method for phase-contrast optical coherence elastography, characterized in that: The rapid volume measurement system for phase contrast optical coherence elastic imaging according to claim 1 is implemented, comprising: The interference signals before and after the deformation of the block are collected by a first industrial camera and a second industrial camera in the block rapid measurement system for phase contrast optical coherence elastic imaging; The computer performs the following operations: Matching the fields of view of the first industrial camera and the second industrial camera; According to the matching identical field of view areas, the interference signals collected by the first industrial camera and the second industrial camera are trimmed to obtain the interferometric data of the block before and after deformation in the same field of view, and the analytical signals of the interferometric data before and after deformation are constructed. Perform compressed sensing on the constructed analytical signal and reconstruct it; The Fourier transform is performed on the compressed sensed and reconstructed signals, and their principal values ​​of the argument are calculated to obtain the wrapped phase field before and after deformation, and finally the differential wrapped phase is calculated.

3. The method for rapid volume measurement based on phase contrast optical coherence elastography according to claim 2, characterized in that: The process of collecting interference signals before and after the deformation of the block includes: The light emitted by the tunable semiconductor wavelength scanning laser is connected to the optical path through an optical fiber. The emitted light is collimated into parallel light by a collimating lens, then passes through the first 45° polarizer to become linearly polarized light and enters the first 50:50 beam splitter prism, and is equally divided into light rays emitting in two directions; the light emitting in the positive direction of the x-axis is reflected by the first reflector and then passes through the second 45° polarizer and 1 / 4 wave plate to become circularly polarized light; the light along the negative direction of the y-axis is reflected by the second reflector, first passes through the block, and then passes through the third 45° polarizer to ensure that the light emitting to the second 50:50 beam splitter prism is linearly polarized; then the linearly polarized light and the circularly polarized light from the block are recombined at the second 50:50 beam splitter prism, and then separated into two orthogonal components by the polarization beam splitter prism, which are received by the first industrial camera and the second industrial camera, and finally the collected interference signal is transmitted to the computer.

4. The method for rapid volume measurement based on phase contrast optical coherence elastography according to claim 2, characterized in that: Match the field of view of the first industrial camera and the second industrial camera, including: Using a shooting resolution board, an area in the image captured by the first industrial camera is selected as a template image, and the image captured by the second industrial camera is used as a reference image. The degree of matching between the template image and the reference image coverage area is calculated using a field of view matching algorithm. The template window is then translated from left to right and from top to bottom in the reference image, traversing the entire reference image pixel by pixel. Finally, the sub-image with the highest matching degree is found in the entire reference image as the matching result, achieving the field of view matching effect. The formula of the field of view matching algorithm is as follows: Among them, NCC(u,v) is the normalized cross-correlation coefficient, T(x,y) is the pixel value of the template image at the coordinate (x,y), and I(x+u,y+v) is the pixel value of the target image at the coordinate (x+u,y+v). is the pixel mean of the template image, is the pixel mean of the target image corresponding to the window at the displacement (u, v).

5. The method for rapid volume measurement based on phase contrast optical coherence elastography according to claim 2, characterized in that: Construct analytical signals of interferometric data before and after deformation, including: The interference signals collected by the first industrial camera and the second industrial camera are respectively expressed as: Where M is the number of surfaces involved in the interference; I R , I m is the DC component, Λ m is the optical path difference between the two beams of light that pass through the block and directly enter the first industrial camera and the second industrial camera, Δk is the laser wave number scanning range; k0 is the starting wave number of the light source, is the frequency generated when the two beams of light pass through the block and directly enter the first industrial camera and the second industrial camera interfere with each other, is the initial phase when the light directly incident on the first industrial camera and the second industrial camera interferes with the light passing through the block; By eliminating the DC component I R , I m After that, two interference signals I x , I y , and then take the two signals as the real and imaginary parts of the analytical signal: I(t)=I x (t)+I y (t)·j (6) According to the matching identical field of view area, the interference signals collected by the first industrial camera and the second industrial camera are trimmed to obtain the real and imaginary parts of the analytical signals, thereby constructing the analytical signals of the interference data before and after deformation.

6. The method for rapid volume measurement based on phase contrast optical coherence elastography according to claim 2, characterized in that: Compressed sensing and reconstruction of the constructed analytical signal, including: The constructed analytical signal is compressed and sensed. The premise of compressed sensing is to represent the signal sparsely. It is R N A set of vector bases in the space, forming a basis matrix Then for any interference data x∈R N Are expressed as: Where α is the sparse representation of x in the γ domain. For volume PC-OCE data, the Fourier basis is used as the sparse representation; Then, a measurement matrix that is unrelated to the Fourier basis is used to project the block PC-OCE data to obtain the observation value y. Finally, x is reconstructed by solving the optimization problem. The expression is: y=φx=φγα (8) Where φγ is the sensing matrix; φ is the measurement matrix, which is also a random matrix; By solving the NP-hard problem, the original signal x is reconstructed.

Images