A cold atom absorption imaging real-time despeckling method, device, equipment and storage medium
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- SOUTH CHINA NORMAL UNIV
- Filing Date
- 2026-03-30
- Publication Date
- 2026-07-14
Smart Images
Figure CN122391009A_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the fields of digital image processing and quantum precision measurement technology, and more specifically, to a method, apparatus, device, and storage medium for real-time destriating cold atom absorption imaging. Background Technology
[0002] Ultracold atomic gases, as highly pure and parameter-tunable quantum many-body systems, are ideal platforms for quantum simulation, quantum information processing, and precision measurement. Accurate measurement and analysis of physical quantities such as the spatial distribution and density fluctuations of ultracold atoms are crucial prerequisites for related research, and absorption imaging is currently the most mainstream detection technique in experiments.
[0003] Standard absorption imaging schemes typically require acquiring three images: the first is the absorption image with atoms present, recording the absorption signal after the resonant laser passes through the atomic cloud; the second is a reference image without atoms, usually acquired after the atoms have moved out of the imaging area; and the third is a dark-field image used to subtract camera dark noise. By processing these three images, the optical density of atoms can be calculated, thus obtaining information about the spatial distribution of atoms. However, in actual experiments, residual reflections and diffraction effects between optical elements, as well as environmental vibrations, can introduce complex interference fringe noise into the images. In particular, when mechanical vibrations or thermal drift cause relative displacement of the fringe patterns in the reference and absorption images, the backgrounds in the two images cannot be perfectly matched, ultimately introducing significant residual fringe structures into the optical density image, severely affecting the extraction of weak atomic signals.
[0004] Traditional post-processing methods (such as frequency domain filtering and statistical averaging) can alleviate noise problems to some extent, but they may lose some atomic structure information during denoising, especially when processing weak signals, making it difficult to balance denoising effectiveness with information preservation. While the fixed basis projection method can construct a basis set using a small number of reference images and approximate the background fringes of the target image through linear combination, this method relies solely on translation operations for basis set expansion, making it ill-suited for complex deformations. It also ignores the problem of fringe spatial frequency shifts caused by changes in object distance or optical system instability, resulting in an incomplete basis set and residual fringe noise. Furthermore, traditional methods typically require offline image processing after the experiment, making it impossible to monitor the atomic cloud state in real time during the experiment, thus limiting experimental efficiency and the timeliness of data acquisition.
[0005] Therefore, developing an imaging method that can both efficiently suppress dynamic stripe noise and have real-time processing capabilities is a technical challenge that urgently needs to be solved in the field of ultracold atom experiments. Summary of the Invention
[0006] In view of this, the present invention proposes a real-time destriating method, device, equipment and storage medium for cold atom absorption imaging. By constructing a multi-scale extended basis set generation strategy and an optimized real-time solver, the denoising process and experimental data acquisition are carried out simultaneously while ensuring the accuracy of noise suppression, thereby significantly improving imaging quality and experimental efficiency.
[0007] To achieve the above objectives, in one embodiment, the present invention provides a real-time destriating method for cold atom absorption imaging, comprising: S10: In each experimental cycle of the cold atom experiment, the original absorption image, reference image and dark field image corresponding to the cold atom experiment are acquired sequentially. S20, the constructed dynamic base set queue stores several frames of reference images acquired and forms an original multi-scale extended base set. The original multi-scale extended base set is masked to obtain an optimized multi-scale extended base set, and the acquired original absorption image is masked to obtain an optimized absorption image. S30, based on the obtained optimized multi-scale extended basis set and the optimized absorption image, a corrected reference image with the same size as the original absorption image is reconstructed; S40, based on the original absorption image, the dark field image and the obtained corrected reference image, calculate the corrected optical density and output the corresponding optical density image.
[0008] Furthermore, S20 includes, S201, sequentially perform validity checks on several frames of reference images acquired. When the current reference image is detected as valid, add it to the dynamic base set queue for storage. Specifically, when the reference image satisfies the following conditions: signal-to-noise ratio greater than a first threshold, image variance greater than a second threshold, and difference from the previous frame reference image less than a third threshold, the current reference image is determined to be valid. S202, perform multi-scale scaling and integer translation operations on each frame of reference image in the dynamic base set queue to generate the original multi-scale extended base set; S203, generate a mask matrix based on the preset position information of cold atomic clusters, and apply the mask matrix to the original multi-scale extended basis set and the acquired original absorption image respectively to obtain the optimized multi-scale extended basis set and the optimized absorption image.
[0009] Furthermore, in S202, the scaling factor set is used. The reference image is scaled using different scaling factors to obtain scaled images of different scales; for each scaled image of different scales, integer translation is performed along the "x" axis and "y" axis within a preset translation range to obtain multiple translated reference images, forming the original multi-scale extended base set.
[0010] Furthermore, in S203, based on the position information corresponding to the cold atom cluster, a rectangular region formed at a preset distance centered on the position information is set to 0, and the background region is set to 1, thereby generating the mask matrix.
[0011] Furthermore, S30 includes, S301, using the optimized multi-scale extended basis set and the optimized absorption image, calculate the autocorrelation matrix of the optimized multi-scale extended basis set and the cross-correlation vector between the optimized multi-scale extended basis set and the optimized absorption image; wherein, the autocorrelation matrix of the optimized multi-scale extended basis set is expressed as follows: in, Let each reference image in the optimized multi-scale extended basis set be unfolded into a one-dimensional vector matrix. To optimize the total number of multi-scale extended base sets, This represents the total number of valid pixels within the mask matrix; The cross-correlation vector between the optimized multi-scale basis set and the optimized absorption image is represented as follows: in, The optimized absorption image is expanded into a one-dimensional vector matrix; S302, a system of linear equations is constructed based on the autocorrelation matrix and the cross-correlation vector, and a corresponding set of linear combination coefficient vectors is obtained by solving the system; wherein, the constructed system of linear equations is represented as follows: S303, using the generated set of linear combination coefficient vectors, the original multi-scale extended basis set is linearly combined to reconstruct a corrected reference image with the same size as the original absorption image; the formula is as follows: in, For the obtained corrected reference image, This is the corresponding linear combination coefficient vector. It is the original multi-scale extended base set.
[0012] Furthermore, in S40, the formula for calculating the corrected optical density is as follows: in, This is the original absorption image. This is a dark field image. To correct the reference image, It is a very small regularization term.
[0013] Furthermore, in S20, a dynamic base set queue is constructed using a circular queue based on the first-in-first-out principle.
[0014] In one embodiment, the present invention also provides a real-time destriating device for cold atom absorption imaging, comprising: The image acquisition module is used to sequentially acquire the corresponding absorption image, reference image, and dark field image in each experimental cycle of the cold atom experiment. The dynamic base set management module is used to store several frames of reference images acquired by the constructed dynamic base set queue and form an original multi-scale extended base set. The original multi-scale extended base set is masked to obtain an optimized multi-scale extended base set, and the original absorption image is masked to obtain an optimized absorption image. The image reconstruction module is used to reconstruct a corrected reference image with the same size as the original absorption image based on the obtained optimized multi-scale extended base set and the optimized absorption image; The calculation and processing module is used to calculate the corrected optical density based on the original absorption image, the dark field image and the obtained corrected reference image, and output the corresponding optical density image.
[0015] In one embodiment, the present invention also provides a computer device including a memory and a processor, the memory storing a computer program, the processor executing the computer program to implement the steps of any of the methods described in the above embodiments.
[0016] In one embodiment, the present invention also provides a computer-readable storage medium having a computer program stored thereon, the computer program being executed by a processor as steps of any of the methods described in the above embodiments.
[0017] Compared with the prior art, the beneficial effects of the present invention include at least the following: 1. High-precision noise suppression: The multi-scale extended basis set constructed by multi-scale scaling and translation operations can more completely characterize the complex dynamic fringe noise caused by frequency drift and spatial displacement, significantly improving the accuracy of background reconstruction and effectively reducing residual fringes.
[0018] 2. Real-time processing capability: The introduction of a dynamic basis set update mechanism, an optimized stable biconjugate gradient solver, and a multi-threaded / graphics processor parallel acceleration architecture reduces the algorithm latency to the millisecond level, enabling real-time online denoising synchronized with experimental data acquisition and greatly improving experimental efficiency.
[0019] 3. High robustness and ease of use: The dynamic base set queue ensures the algorithm's adaptability to slow drift in the experimental environment. The interactive graphical user interface allows experimenters to adjust key parameters such as the mask in real time and observe the effects immediately, reducing the difficulty of operation and improving the user experience.
[0020] 4. Preservation of original signals: By eliminating atomic signal interference in the mask region and solving for noise patterns in the background region without atomic signals, this method effectively denoises while preserving the original atomic density information to the greatest extent, providing a high-quality data foundation for subsequent quantitative physical analysis. Attached Figure Description
[0021] The accompanying drawings, which form part of this invention, are used to provide a further understanding of the invention. The illustrative embodiments of the invention and their descriptions are used to explain the invention and do not constitute an improper limitation of the invention.
[0022] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0023] Furthermore, the accompanying drawings are not drawn to a 1:1 scale, and the relative dimensions of the various components are shown in the drawings only as examples and not necessarily to actual scale.
[0024] Figure 1 This is a schematic flowchart of a real-time stripe removal method for cold atom absorption imaging provided in an embodiment of the present invention; Figure 2 This is a schematic diagram illustrating the effect of applying a mask to a reference image in an embodiment of the present invention; Figure 3 This is a schematic diagram illustrating the effect of applying a mask to the absorption image in an embodiment of the present invention; Figure 4 This is a schematic diagram of the structure of a real-time destriating device for cold atom absorption imaging provided in an embodiment of the present invention; Figure 5 This is a comparison chart showing the difference between the optical density image before and after correction and the absorption image. Figure 6 An internal structural diagram of a computer device provided in an embodiment of the present invention. Detailed Implementation
[0025] To make the above-mentioned objects, features, and advantages of the present invention more apparent and understandable, specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. Many specific details are set forth in the following description to provide a thorough understanding of the present invention. However, the present invention can be practiced in many other ways different from those described herein, and those skilled in the art can make similar modifications without departing from the spirit of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below.
[0026] In one embodiment, please refer to Figure 1 This invention provides a real-time destriating method for cold atom absorption imaging. This embodiment illustrates the application of this method to a terminal. It is understood that this method can also be applied to a server, and further to a system including both a terminal and a server, and implemented through interaction between the terminal and the server. In this embodiment, the method includes the following steps: S10: In each experimental cycle of the cold atom experiment, the original absorption image, reference image and dark field image corresponding to the cold atom experiment are acquired sequentially. Specifically, in this embodiment, a cold atom experiment is first conducted in an ultra-high vacuum system environment, with the background pressure maintained at [insert pressure here]. Ultracold atomic gas, on the order of Pa, was prepared using existing magneto-optical trap technology and evaporative cooling technology, and then loaded into a scientific observation cavity. This embodiment uses ytterbium atoms (… 173 The invention uses the Yb) ultracold Fermi gas absorption imaging system as an example, but this method is also applicable to other alkali metal or alkaline earth metal atomic systems, and is not limited thereto. It is also understood that in this system, the probe light resonates with the cold atoms, and the wavelength of the probe light used for imaging is preferably 399 nm.
[0027] Within one experimental cycle, three images are acquired sequentially according to the standard absorption imaging procedure: absorption image (Contains atoms), reference image (No atoms) and dark field images (Without illumination), each image is used for subsequent calculations. Understandably, this involves absorbing the image. Image of atoms in their presence, used to record the absorption signal of a probe laser passing through an atomic cloud of cold atoms; reference image. The reference image is taken when there are no atoms, usually after the atoms have been removed from the shooting area; the dark field image is a dark field image taken under no-light conditions with all light sources turned off, used to subtract dark noise from the camera.
[0028] S20, the constructed dynamic base set queue stores several frames of reference images acquired and forms an original multi-scale extended base set. The original multi-scale extended base set is masked to obtain an optimized multi-scale extended base set, and the acquired original absorption image is masked to obtain an optimized absorption image. Specifically, when acquiring reference images showing the presence of atomic clusters, a dynamic base set queue is first constructed to store several frames of acquired reference images. Preferably, a circular queue based on the first-in-first-out principle can be used as the dynamic base set queue, and the maximum capacity of the dynamic base set queue is preset. In this embodiment, the maximum capacity .
[0029] Furthermore, step S20 specifically includes, S201, the validity of several frames of reference images acquired are sequentially checked. When the current reference image is detected as valid, it is added to the dynamic base set queue for storage; wherein, when the reference image meets the signal-to-noise ratio... Image variance and the degree of difference from the previous reference image At that time, it is determined that the current reference image detection is valid.
[0030] Understandably, when a current reference image is detected as valid, it is considered high-quality and can be added to the constructed dynamic basis set queue. If the queue is full, the oldest image is removed, thus ensuring that the constructed dynamic basis set queue always reflects the current experimental state. , , In this embodiment, a pre-set threshold is preferred. , , .
[0031] S202, perform multi-scale scaling and integer translation operations on each frame of reference image in the dynamic base set queue to generate the original multi-scale extended base set; Specifically, after storing several frames of valid reference images into the constructed dynamic basis set queue, each frame of reference image is further subjected to multi-scale scaling and integer translation operations to construct a corresponding multi-scale extended basis set. Multi-scale scaling is performed on each frame of reference image to simulate the spatial frequency changes of interference fringes caused by thermal drift or defocusing of the optical system. In this embodiment, before performing multi-scale scaling, a set of scaling factors is first defined. The reference image is scaled at different scales using different scaling factors from the scaling factor set. For example, in this embodiment, it can be taken as follows: That is, for each frame of the original reference image The image is scaled three times to obtain three scaled images at different scales; further, for each scaled image, within a preset translation range... Integer translation along the x-axis and y-axis is performed: Ishifted(a,b) = Iscaled(a+j,b+k), where j,k∈{-d,...,d}, resulting in multiple translated reference images. In this embodiment, the preset translation range is typically set to d=1 pixels. Therefore, it is understandable that after obtaining three scaled images at different scales, each scaled image is further translated by integers within the translation range [-1,1], i.e., by integer translation along the x-axis and y-axis within -1, 0, and 1 pixels, generating 9 translated images. It is understandable that for N=6 original reference images, a total of 6*3*9=162 images are generated, thus forming a multi-scale extended base set covering different fringe frequencies and displacements. Therefore, through multi-scale scaling and translation operations, the finite N original reference images can be expanded into an original multi-scale extended base set X containing multiple noise modes, and the scale of the original multi-scale extended base set is expanded to... .
[0032] S203, generate a mask matrix according to the preset position information of cold atom clusters, and apply the mask matrix to the original multi-scale extended basis set and the acquired original absorption image respectively to obtain the optimized multi-scale extended basis set and the optimized absorption image; Specifically, the experiment pre-sets the positional information corresponding to the cold atom clusters. After generating the multi-scale extended basis set, further calculations are performed based on the position information corresponding to the cold atom clusters. Generate mask matrix Defining atomic region masks in the image can eliminate interference from atomic signals on background noise estimation. Understandably, cold atomic clusters are typically located at the center of the image, such as... Figure 2 As shown, the effect of applying the mask matrix to the multi-scale extended basis set is presented; for example, Figure 3 As shown, the effect of applying the mask matrix to the absorption image to obtain the optimized absorption image is presented. Assume the position information corresponding to the cold atom clusters is... Size information (radius) (pixels), in Within the rectangular region (i.e., the region of interest containing the atomic cluster), 0s are set, and the background region is set to 1, generating a binary mask matrix. Furthermore, the generated binary mask matrix is applied to the original multi-scale extended basis set through Hadamard product. and the acquired absorption images The optimized multi-scale extended basis set after masking is obtained. and optimized absorption image after masking. Understandably, the optimized multi-scale extended basis set and optimized absorption image obtained after masking process shield the atomic signals, retaining only background pixels for subsequent analysis.
[0033] S30, based on the obtained optimized multi-scale extended basis set and the optimized absorption image, a corrected reference image with the same size as the original absorption image is reconstructed; Specifically, step S30 includes, S301, using the optimized multi-scale extended basis set and the optimized absorption image, calculate the autocorrelation matrix of the optimized multi-scale extended basis set and the cross-correlation vector between the optimized multi-scale extended basis set and the optimized absorption image; Specifically, the obtained optimized multi-scale extended base set Each reference image in the dataset is unfolded into a one-dimensional vector, forming the following vector matrix: in, To optimize the total number of multi-scale extended base sets, This represents the total number of valid pixels within the mask matrix.
[0034] The obtained optimized absorption image It can also be expanded into a one-dimensional vector. The autocorrelation matrix of the optimized multi-scale extended basis set is further calculated as follows: The optimized multi-scale basis set and cross-correlation vector are as follows: .
[0035] In this embodiment, to accelerate computation, PyCUDA is used to call the graphics processor to perform parallel matrix multiplication, which significantly reduces the time consumption.
[0036] S302, a system of linear equations is constructed based on the autocorrelation matrix and the cross-correlation vector, and a corresponding set of linear combination coefficient vectors is obtained by solving the system; wherein, the constructed system of linear equations is represented as follows: Specifically, after obtaining the optimized multi-scale extended basis set and the optimized absorption image, the autocorrelation matrix of the optimized multi-scale extended basis set is calculated. and the cross-correlation vector between the optimized multi-scale extended basis set and the optimized absorption image. Further based on the obtained autocorrelation matrix and cross-correlation vector Construct a system of linear equations Furthermore, the system of linear equations constituted... Solving this problem yields the corresponding linear combination coefficient vector. Understandably, when the size of the optimized multi-scale extended base set is 162, In this embodiment, the stable biconjugate gradient method is used to solve the linear equation system. This method has the advantages of fast convergence speed and high numerical stability, making it suitable for real-time processing. A tolerance of 10 can be set during the solution process. -8 The maximum number of iterations is 100. In other embodiments, the corresponding tolerance and maximum number of iterations are also set according to the specific circumstances.
[0037] S303, using the generated set of linear combination coefficient vectors, the original multi-scale extended basis set is linearly combined to reconstruct a corrected reference image with the same size as the original absorption image; Specifically, it preserves the original, unmasked, original multiscale extended basis set. Without changing, use the generated set of linear combination coefficient vectors A corrected reference image of the same size as the original absorption image is obtained by linearly combining the original multiscale extended basis set. The formula is as follows: Understandably, since the original multi-scale extended base set is used during combination, the reconstructed image is a corrected reference image that highly matches the background noise of the current absorption image. The corrected reference image has a stripe pattern that highly matches the background noise of the current absorption image across the entire image.
[0038] S40, based on the original absorption image, the dark field image and the obtained corrected reference image, calculate the corrected optical density and output the corresponding optical density image.
[0039] Specifically, based on the obtained corrected reference image of the same size as the original absorption image. Further combined with the acquired dark field images The corrected optical density is calculated and the corresponding optical density image is output. The calculation formula is as follows: in, This is the original absorption image. This is a dark field image. To correct the reference image, This is a very small regularization term used to avoid taking the logarithm of zero or negative values.
[0040] Understandably, the corrected optical density is calculated. Then, a corrected high signal-to-noise ratio atomic density distribution image is output, which will be subsequently... The data is displayed in real-time on the graphical user interface for monitoring by experimenters. In this embodiment, .
[0041] This invention provides a real-time destriating method for cold atom absorption imaging. By constructing a multi-scale extended basis set, it can more completely characterize the complex dynamic fringe noise caused by frequency drift and spatial displacement during cold atom absorption imaging, significantly improving the accuracy of background reconstruction. While ensuring noise suppression accuracy, it enables the denoising process to be carried out simultaneously with experimental data acquisition, thereby significantly improving imaging quality and experimental efficiency. By eliminating atomic signal interference in the mask region and solving for noise patterns in the background region without atomic signals, it effectively denoises while maximizing the preservation of the original atomic density information, providing a high-quality data foundation for subsequent quantitative physical analysis.
[0042] It should be understood that, although Figure 1 The steps in the flowchart are shown sequentially as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified herein, there is no strict order in which these steps are executed, and they can be performed in other orders. Figure 1 At least some of the steps in the process may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but may be executed at different times. The execution order of these steps or stages is not necessarily sequential, but may be executed in turn or alternately with other steps or at least some of the steps or stages in other steps.
[0043] In another embodiment, such as Figure 4 As shown, the present invention also provides a real-time destriating device for cold atom absorption imaging, comprising: The image acquisition module is used to sequentially acquire the corresponding absorption image, reference image, and dark field image in each experimental cycle of the cold atom experiment. The dynamic base set management module is used to store several frames of reference images acquired by the constructed dynamic base set queue and form an original multi-scale extended base set. The original multi-scale extended base set is masked to obtain an optimized multi-scale extended base set, and the original absorption image is masked to obtain an optimized absorption image. The image reconstruction module is used to reconstruct a corrected reference image with the same size as the original absorption image based on the obtained optimized multi-scale extended base set and the optimized absorption image; The calculation and processing module is used to calculate the corrected optical density based on the original absorption image, the dark field image and the obtained corrected reference image, and output the corresponding optical density image.
[0044] Specific limitations regarding the real-time destriating apparatus for cold atom absorption imaging provided by this invention can be found in the above-described limitations of the real-time destriating method for cold atom absorption imaging, and will not be repeated here. Each module in the aforementioned real-time destriating apparatus for cold atom absorption imaging can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the memory of a computer device as software, so that the processor can call and execute the operations corresponding to each module.
[0045] System Implementation The absorption imaging system of this invention is developed based on the PyQt5 framework, with the destriating device integrated as an independent module within the absorption imaging system software. The system employs a multi-threaded parallel architecture, with image acquisition, processing, reconstruction, and display running on separate threads to ensure a smooth interface. It utilizes a graphics processing unit (GPU) for accelerated computation, achieving real-time processing with millisecond-level latency. The graphical user interface supports interactive adjustment of core parameters such as mask parameters, scaling factors, and translation range, and displays the correction effect in real time. Users can adjust the mask position and size in real time via selection boxes on the interface; parameters such as scaling factors and translation range also support dynamic configuration, with immediate feedback on the adjustment effect.
[0046] Experimental results Applying the real-time destriating method provided in this invention to the 173Yb ultracold Fermi gas absorption imaging system, the horizontal stripe contrast in the corrected optical density image is significantly reduced compared to the original optical density image. Quantitative analysis shows that the background noise variance decreases from... Down to This represents a reduction of approximately 16.7%; the noise power in the Fourier spectrum is reduced by approximately 25%. For example... Figure 5 As shown in the figure, a comparison of the differences between the optical density image before and after correction and the absorption image is presented. It can be seen from the figure that the difference between the corrected optical density image and the absorption image is significantly reduced. The difference after correction is much smaller than the original difference, proving the effectiveness of the algorithm.
[0047] In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 6As shown, the computer device includes a processor, memory, communication interface, display screen, and input devices connected via a system bus. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, carrier networks, NFC (Near Field Communication), or other technologies. When executed by the processor, the computer program implements a real-time stripe removal method for cold atom absorption imaging. The display screen can be a liquid crystal display (LCD) or an e-ink display. The input devices can be a touch layer covering the display screen, buttons, a trackball, or a touchpad mounted on the computer device casing, or an external keyboard, touchpad, or mouse.
[0048] Those skilled in the art will understand that Figure 6 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.
[0049] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of the real-time destriating method for cold atom absorption imaging provided in the above embodiment.
[0050] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the steps of the real-time destriating method for cold atom absorption imaging provided in the above embodiment.
[0051] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the methods described above. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, or optical storage, etc. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.
[0052] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.
[0053] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the invention patent. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this patent application should be determined by the appended claims.
Claims
1. A real-time stripe removal method for cold atom absorption imaging, characterized in that, include: S10: In each experimental cycle of the cold atom experiment, the original absorption image, reference image and dark field image corresponding to the cold atom experiment are acquired sequentially. S20, the constructed dynamic base set queue stores several frames of reference images acquired and forms an original multi-scale extended base set. The original multi-scale extended base set is masked to obtain an optimized multi-scale extended base set, and the acquired original absorption image is masked to obtain an optimized absorption image. S30, based on the obtained optimized multi-scale extended basis set and the optimized absorption image, a corrected reference image with the same size as the original absorption image is reconstructed; S40, based on the original absorption image, the dark field image and the obtained corrected reference image, calculate the corrected optical density and output the corresponding optical density image.
2. The real-time destriating method for cold atom absorption imaging according to claim 1, characterized in that, S20 includes, S201, construct a dynamic base set queue, and sequentially perform validity checks on several frames of reference images acquired. When the current reference image is detected as valid, add it to the dynamic base set queue for storage. Specifically, when the reference image satisfies the following conditions: signal-to-noise ratio greater than a first threshold, image variance greater than a second threshold, and difference from the previous frame reference image less than a third threshold, the current reference image is determined to be valid. S202, perform multi-scale scaling and integer translation operations on each frame of reference image in the dynamic base set queue to generate the original multi-scale extended base set; S203, generate a mask matrix based on the preset position information of cold atomic clusters, and apply the mask matrix to the original multi-scale extended basis set and the acquired original absorption image respectively to obtain the optimized multi-scale extended basis set and the optimized absorption image.
3. The real-time destriating method for cold atom absorption imaging according to claim 2, characterized in that, In S202, the scaling factor set is used. The reference image is scaled using different scaling factors to obtain scaled images of different scales; for each scaled image of different scales, integer translation is performed along the "x" axis and "y" axis within a preset translation range to obtain multiple translated reference images, forming the original multi-scale extended base set.
4. The real-time destriating method for cold atom absorption imaging according to claim 2, characterized in that, In S203, based on the position information corresponding to the cold atom cluster, a rectangular area formed at a preset distance centered on the position information is set to 0, and the background area is set to 1, thereby generating the mask matrix.
5. The real-time destriating method for cold atom absorption imaging according to claim 1, characterized in that, S30 includes, S301, using the optimized multi-scale extended basis set and the optimized absorption image, calculate the autocorrelation matrix of the optimized multi-scale extended basis set and the cross-correlation vector between the optimized multi-scale extended basis set and the optimized absorption image; wherein, the autocorrelation matrix of the optimized multi-scale extended basis set is expressed as follows: in, Let each reference image in the optimized multi-scale extended basis set be unfolded into a one-dimensional vector matrix. To optimize the total number of multi-scale extended base sets, This represents the total number of valid pixels within the mask matrix; The cross-correlation vector between the optimized multi-scale basis set and the optimized absorption image is represented as follows: in, The optimized absorption image is expanded into a one-dimensional vector matrix; S302, a system of linear equations is constructed based on the autocorrelation matrix and the cross-correlation vector, and a corresponding set of linear combination coefficient vectors is obtained by solving the system; wherein, the constructed system of linear equations is represented as follows: S303, using the generated set of linear combination coefficient vectors, the original multi-scale extended basis set is linearly combined to reconstruct a corrected reference image with the same size as the original absorption image; the formula is as follows: in, For the obtained corrected reference image, This is the corresponding linear combination coefficient vector. It is the original multi-scale extended base set.
6. The real-time destriating method for cold atom absorption imaging according to claim 1, characterized in that, In S40, the formula for calculating the corrected optical density is as follows: in, This is the original absorption image. This is a dark field image. To correct the reference image, It is a very small regularization term.
7. The real-time destriating method for cold atom absorption imaging according to claim 1, characterized in that, In S20, the dynamic base set queue is constructed using a circular queue based on the first-in-first-out principle.
8. A real-time stripe removal device for cold atom absorption imaging, characterized in that, include: The image acquisition module is used to sequentially acquire the corresponding absorption image, reference image, and dark field image in each experimental cycle of the cold atom experiment. The dynamic base set management module is used to store several frames of reference images acquired by the constructed dynamic base set queue and form an original multi-scale extended base set. The original multi-scale extended base set is masked to obtain an optimized multi-scale extended base set, and the original absorption image is masked to obtain an optimized absorption image. The image reconstruction module is used to reconstruct a corrected reference image with the same size as the original absorption image based on the obtained optimized multi-scale extended base set and the optimized absorption image; The calculation and processing module is used to calculate the corrected optical density based on the original absorption image, the dark field image and the obtained corrected reference image, and output the corresponding optical density image.
9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 7.
10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.