Method for eliminating electromagnetic interference noise in detection image and electronic equipment

Abstract

The invention discloses a method for eliminating electromagnetic interference noise in a detection image and electronic equipment. The method comprises the following steps: acquiring a detection signal of a detector to an astronomical environment; extracting a detection value of each pixel point in the overscan region from the detection signal; performing modeling analysis on the detection sequence according to the distribution attribute of the detection sequence in the time domain; performing mathematical modeling on the electromagnetic interference noise to construct an electromagnetic interference model of the electromagnetic interference noise in a time domain; determining values of related parameters in a sine function and a Gaussian function in the electromagnetic interference model according to the detection value of each pixel point in the overscan region; and based on an electromagnetic interference model, performing de-noising processing on the detection value of each pixel point in the photosensitive area to generate a detection image of an astronomical environment. According to the embodiment of the invention, the noise filtering effect can be improved, and the imaging quality of the detection image is improved.

Description

Technical Field

[0001] The present application relates to the field of space detection technology, and in particular to a method and electronic equipment for eliminating electromagnetic interference noise in a detection image. Background Art

[0002] When using a detector to detect the surrounding environment, the signal returned by the detector is mostly coupled with electromagnetic interference introduced by the detector itself or surrounding equipment. Specifically, it appears as roughly striped noise on the image, forming interference stripes. Affected by these interference stripes, the detector pixel reading will change, and the size and period of the impact are determined by the nature of the interference source.

[0003] In order to eliminate this electromagnetic interference, the following technical means can be used: using electromagnetic shielding technology to reduce electromagnetic interference, using grounding technology to eliminate electromagnetic interference, using wiring technology to improve electromagnetic interference, using filtering technology to reduce electromagnetic interference, etc. Among them, the usual treatment method is to perform electromagnetic shielding or circuit grounding on the detector in hardware, which can improve the problem of electromagnetic interference to a certain extent.

[0004] However, in some special cases, the above operations are no longer applicable. The main reasons are: First, in practice, it is difficult to achieve complete electromagnetic shielding in hardware, especially when the detector needs to operate simultaneously with complex circuits or high-power motor equipment, the electromagnetic shielding of the detector deteriorates. This is particularly obvious in astronomical observations. When the detector (such as CCD detector or CMOS detector) is working, the motor used to drive the telescope is always in operation, and the fan for cooling the detector is always in operation, which can lead to the deterioration of the electromagnetic shielding of the detector; Second, in special environments, the cost and maintenance investment to achieve effective electromagnetic shielding or grounding operations may be very high. For example, when using detectors for astronomical observations in harsh environments such as the North and South Poles, or when using space imaging satellites equipped with detectors for space exploration, the hardware electromagnetic shielding of the detectors is difficult, and the electromagnetic shielding is prone to failure during long-term working operations. If it fails, it will have a considerable impact on the imaging process of astronomical observations or space exploration.

[0005] In addition, unlike the above-mentioned means that rely on hardware, there is currently a method for improving detector image quality that can be implemented through software. The general principle is: the detector takes multiple photos of the target to obtain multiple images (including interference fringes), and performs "median superposition" processing on the multiple images. Finally, a detection image (excluding interference fringes) can be obtained, which can achieve the purpose of removing noise in the detection image.

[0006] However, this method of median superposition based on multiple images has the following disadvantages: First, this method requires multiple shots of the same target. Obviously, this method cannot be implemented with only one shot. If the number of shots is too small, the number of images is too small to achieve high-precision processing, which will form a residual effect. In the astronomical detection process, the detector takes a long time to shoot, and multiple shots will also waste noise time costs. Secondly, the median stacking method essentially requires that the same pixel is not covered by interference signals for most of the time, which requires strict control of the interference signal appearance pattern to finally obtain an image without interference fringes. In one case, if the position and intensity of the interference fringes on each of the multiple images taken are fixed, then the median stacking method will fail.

[0007] Due to the above-mentioned defects of the median superposition method, when using this method for environmental detection, high requirements are placed on the stability control of the detector, high requirements are placed on the computing power of the processor, and residual effects are prone to occur. The detection results have low reliability and poor practicality. Summary of the invention

[0008] In view of this, the embodiments of the present application provide a method and electronic device for eliminating electromagnetic interference noise in a detection image, which can improve the noise filtering effect and enhance the imaging quality of the detection image.

[0009] An embodiment of the present application provides a method for eliminating electromagnetic interference noise in a detection image, comprising: obtaining a detection signal of a detector for an astronomical environment, wherein the detection signal comprises a detection value of each pixel point in an overscan area and a detection value of each pixel point in a photosensitive area; extracting the detection value of each pixel point in the overscan area from the detection signal, and generating a detection sequence of electromagnetic interference noise according to a time domain output order of the detection value; modeling and analyzing the detection sequence according to the distribution properties of the detection sequence in the time domain to determine the two parts of the sine distribution and Gaussian distribution of the electromagnetic interference noise in the detection environment in the time domain; mathematically modeling the electromagnetic interference noise to construct an electromagnetic interference model of the electromagnetic interference noise in the time domain, the electromagnetic interference model comprising a sine function and a Gaussian function; determining the values ​​of relevant parameters in the sine function and the Gaussian function in the electromagnetic interference model according to the detection value of each pixel point in the overscan area; and denoising the detection value of each pixel point in the photosensitive area based on the electromagnetic interference model to generate a detection image of the astronomical environment.

[0010] Optionally, according to the method of an embodiment of the present application, the detector includes at least one of the following: a detector installed in an astronomical telescope, a detector carried in a space imaging satellite, a detector used in a harsh environment, a CCD detector, and a CMOS detector.

[0011] Optionally, according to the method of the embodiment of the present application, before mathematically modeling the electromagnetic interference noise, the method further includes: decomposing the detection sequence to obtain one or more pulse signal chains, where different pulse signal chains have different frequencies and amplitudes.

[0012] Optionally, according to the method of an embodiment of the present application, the method also includes: obtaining interference characteristic information corresponding to one or more hardware devices in the detection environment in which the detector is located; matching the one or more pulse signal chains obtained by decomposition processing with the interference characteristic information corresponding to one or more hardware devices, and determining the successfully matched hardware device as the interference source of the detector.

[0013] Optionally, according to the method of the embodiment of the present application, the electromagnetic interference model of the electromagnetic interference noise in the time domain adopts the following form:

[0014] Where t is the output order of the detection value, I(t) is the noise interference value corresponding to the t-th detection value, I0 is the interference value of the background noise in the detection environment, the subscript i represents the i-th pulse signal chain, N is a positive integer, and A i is the amplitude of the ith pulse signal chain, ω i and α i are the frequency and phase of the sinusoidal part in the i-th pulse signal chain, σ i and β i are the width and displacement of the Gaussian part in the i-th pulse signal chain, respectively.

[0015] Optionally, according to the method of an embodiment of the present application, the detection sequence is modeled and analyzed according to the distribution properties of the detection sequence in the time domain, including: decomposing the detection sequence to obtain one or more pulse signal chains, different pulse signal chains have different frequencies and amplitudes; modeling and analyzing each pulse signal chain obtained by the decomposition, wherein part or all of the pulse signals are respectively included in a sinusoidal part and a Gaussian part to construct a fitting function for each pulse signal chain; based on the fitting function of each pulse signal chain, a mathematical model of the detection sequence is constructed; and it is determined that the mathematical model of the detection sequence matches the distribution properties of the detection sequence.

[0016] Optionally, according to the method of the embodiment of the present application, the fitting function of a single pulse signal chain takes the following form:

[0017] Where, the subscript i represents the i-th pulse signal chain, t is the output order of the detection value, Ii(t) is the interference value of the i-th electromagnetic interference on the t-th output, and A i is the amplitude of the ith pulse signal chain, ω i and α iare the frequency and phase of the sinusoidal part in the i-th pulse signal chain, σ i and β i are the width and displacement of the Gaussian part in the i-th pulse signal chain, respectively.

[0018] An embodiment of the present application provides an electronic device, which includes a processor and a memory storing computer program instructions; when the processor executes the computer program instructions, the steps of the above method are implemented.

[0019] An embodiment of the present application provides a computer-readable storage medium, on which computer program instructions are stored. When the computer program instructions are executed by a processor, the steps of the method described above are implemented.

[0020] An embodiment of the present application provides a computer program product, wherein the computer program product includes computer program instructions, and when the computer program instructions are executed by a processor, the steps of the method described above are implemented.

[0021] According to an embodiment of the present application, the detection value of each pixel point in the overscan area of ​​the detector is obtained, and then a detection sequence of electromagnetic interference noise is generated according to the time domain output order of the detection value. After that, the detection sequence is modeled and analyzed according to the distribution properties of the detection sequence in the time domain. After determining that the electromagnetic interference noise in the detection environment in the time domain includes two parts, namely, the sinusoidal distribution and the Gaussian distribution, an electromagnetic interference model including a sinusoidal function and a Gaussian function is constructed, and the detection value of each pixel point in the overscan area is combined to determine the values ​​of relevant parameters in the sinusoidal function and the Gaussian function in the electromagnetic interference model, thereby obtaining an electromagnetic interference model including specific parameters. Based on the electromagnetic interference model, the actual fluctuation of the electromagnetic interference noise can be accurately represented, so that in the process of denoising, the interference value of the electromagnetic interference corresponding to each pixel point in the photosensitive area can be accurately determined. After denoising, a clear detection image of the astronomical environment can be generated. In the embodiment of the present application, for the denoising task of the detection image, there is no need to take multiple shots of the detector or to add additional hardware shielding equipment, which effectively reduces the denoising cost. Moreover, since there are sufficient detection values ​​in the overscan area, the electromagnetic interference model has sufficient accuracy even for a single shooting task, thereby ensuring the filtering effect of the electromagnetic interference noise and improving the imaging quality of the detection image, providing a clearer and more accurate image basis for the subsequent analysis of astronomical data. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the technical solution of the embodiments of the present application, the following briefly introduces the drawings in the embodiments of the present application.

[0023] Figure 1It is a flowchart of a method for eliminating electromagnetic interference noise in a detection image according to an embodiment of the present application.

[0024] Figure 2 A schematic diagram of interference noise provided in an embodiment of the present application.

[0025] Figure 3 It is a schematic diagram of a denoising effect of a detection image provided in an embodiment of the present application.

[0026] Figure 4 Schematic diagram of a pulse signal chain provided in an embodiment of the present application.

[0027] Figure 5 It is a schematic diagram of an electronic device used to implement the method for eliminating electromagnetic interference noise in a detection image according to an embodiment of the present application. DETAILED DESCRIPTION

[0028] The principles and spirit of the present application will be described below with reference to several exemplary embodiments. It should be understood that the purpose of providing these embodiments is to make the principles and spirit of the present application clearer and more thorough, so that those skilled in the art can better understand and implement the principles and spirit of the present application. The exemplary embodiments provided herein are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments herein, all other embodiments obtained by ordinary technicians of the art without creative work are within the scope of protection of this application.

[0029] It should be noted that the acquisition, storage, use, and processing of data in the embodiments of the present application are in compliance with the relevant provisions of national laws and regulations.

[0030] In this document, terms such as first, second, third, etc. are only used to distinguish one entity (or operation) from another entity (or operation), but not to require or imply any order or relationship between these entities (or operations).

[0031] The embodiments of the present application relate to terminal devices and / or servers. It is known to those skilled in the art that the embodiments of the present application can be implemented as a system, device, equipment, method, computer-readable storage medium or computer program product. Therefore, the present disclosure can be specifically implemented in at least one of the following forms: complete hardware, complete software, or a combination of hardware and software.

[0032] As people explore more and more space, there are more and more application environments for collecting image data through detectors. Detectors can capture light in the environment and convert it into electrical signals. Then, the detection image can be generated by combining the electrical signals. Detectors can usually be used, such as charge-coupled device (CCD) detectors, complementary metal-oxide-semiconductor (CMOS) detectors, etc., which are not listed here one by one.

[0033] In the related art, the working environment of the detector is susceptible to electromagnetic interference from the surrounding circuits, which causes striped noise to form in the detection image generated based on the electrical signal output by the detector. Since these surrounding circuits include circuits shared with the detector or other necessary electronic equipment, it is difficult to completely cut off the interference source to avoid noise. Although in the related art, improvement methods have been proposed from the hardware or software to reduce the striped noise in the detection image, the current improvement methods all have the problem of high denoising cost and poor denoising effect.

[0034] Based on this, the embodiment of the present application provides a method, device and electronic device for eliminating electromagnetic interference noise in the detection image, by obtaining the detection value of each pixel point in the overscan area of ​​the detector and the time domain output order of the detection value, constructing a suitable electromagnetic interference model, and determining accurate model parameters, so as to accurately represent the actual fluctuation of electromagnetic interference noise through the electromagnetic interference model. Then, in the process of denoising, the interference value of the electromagnetic interference corresponding to each pixel point in the photosensitive area can be accurately determined. After denoising, a clear detection image of the astronomical environment can be generated, providing clearer and more accurate image data for the subsequent analysis of astronomical data.

[0035] The following is an introduction to the method for eliminating electromagnetic interference noise in a detection image provided by an embodiment of the present application in conjunction with the accompanying drawings. Figure 1 A flowchart of a method for eliminating electromagnetic interference noise in a detection image provided in some embodiments of the present application. Figure 1 As shown, the method includes the following steps 101 to 106.

[0036] Step 101, obtaining a detection signal of the detector for the astronomical environment, wherein the detection signal includes a detection value of each pixel point in the overscan area and a detection value of each pixel point in the photosensitive area.

[0037] Step 102: extract the detection value of each pixel in the overscan area from the detection signal, and generate a detection sequence of electromagnetic interference noise according to the time domain output order of the detection value.

[0038] Step 103 , modeling and analyzing the detection sequence according to the distribution property of the detection sequence in the time domain, so as to determine the two parts of the sine distribution and the Gaussian distribution contained in the electromagnetic interference noise in the detection environment in the time domain.

[0039] Step 104 , mathematically modeling the electromagnetic interference noise to construct an electromagnetic interference model of the electromagnetic interference noise in the time domain, wherein the electromagnetic interference model includes a sine function and a Gaussian function.

[0040] Step 105: Determine the values ​​of relevant parameters in the sine function and the Gaussian function in the electromagnetic interference model according to the detection value of each pixel point in the overscan area.

[0041] Step 106, based on the electromagnetic interference model, denoising is performed on the detection value of each pixel point in the photosensitive area to generate a detection image of the astronomical environment.

[0042] The above steps are described in detail below in conjunction with specific embodiments.

[0043] Specifically, in step 101, the detection signal of the astronomical environment includes the detection value of each pixel point in the overscan area and the detection value of each pixel point in the photosensitive area.

[0044] Exemplarily, the detector includes at least one of the following: a detector installed in an astronomical telescope, a detector carried in a space imaging satellite, a detector used in a harsh environment, a CCD detector, or a CMOS detector. Among them, in a harsh environment such as an ice surface environment at high altitude or high altitude, specifically, a telescope array established at an Antarctic astronomical observatory.

[0045] The detector includes an overscan area and a photosensitive area. Optionally, the overscan area in the detector can be located around the photosensitive area, for example, at the beginning or end of each row of pixels, and the overscan area can also be located at other photosensitive positions, which is not specifically limited here.

[0046] When the detector performs a detection task, the pixels in the photosensitive area can receive ambient light, and the light collected by each pixel is converted into an electrical signal. After the output value corresponding to the electrical signal is output, the detection value of the pixel in the photosensitive area is obtained, and a detection image can be generated based on the detection value. Among them, the detection image may include the imaging of astronomical data, which includes but is not limited to celestial bodies such as stars, galaxies, or other objects that can be photographed by the detector.

[0047] The overscan area in the detector is an area that is not illuminated. The detection value of each pixel in the overscan area is affected by the noise in the environment where the detector is located. In other words, the detection value of each pixel in the overscan area is based on the value generated by the noise in the environment where the detector is located, such as the background noise and the electromagnetic interference noise received by the readout channel during output.

[0048] It is understandable that the detection value of each pixel in the photosensitive area is also affected by the noise in the environment where the detector is located. Therefore, the detection image generated directly based on the detection value of each pixel in the photosensitive area will have noise. For example, the electromagnetic interference source in the detection environment is, for example, a periodically working electronic device such as a motor or an AC circuit. These electronic devices are prone to electromagnetic leakage. Therefore, when the detector outputs a signal through the readout channel, the output signal corresponding to each pixel in the photosensitive area and the overscan area will be affected.

[0049] As a specific example, the CCD detector used in the telescope is used to capture images. In the readout channel of the CCD detector, the interference pattern of the straight stripes and the inclined stripes can be clearly observed, where, for example, Figure 2 Stripe noise shown. The brightness of pixels affected by interfering stripes will increase or decrease by about 20-50 ADU.

[0050] In some embodiments, the readout period of the readout channel of the detector and the interference period of the electromagnetic interference source will affect the display position of the interference noise on the image. Therefore, the striped noise formed on the detected image may be inclined striped noise, vertical striped noise, or discontinuous striped noise. The display forms of the noise are not listed here one by one.

[0051] Next, step 102 is involved, in which the detection value of each pixel point in the overscan area is extracted from the detection signal, and a detection sequence of electromagnetic interference noise is generated according to the time domain output order of the detection value.

[0052] Specifically, each detector may correspond to one or more readout channels, and each readout channel sequentially outputs the corresponding electrical signal of each pixel based on the arrangement order of the pixels. For example, the readout channel of the detector sequentially outputs along the row direction of the pixels, or the signal readout channel sequentially outputs along the column direction of the pixels.

[0053] In some optional embodiments, the output mode of sequentially outputting along the row direction of pixels is consistent with the horizontal scanning of the pixel arrangement, which can improve the readout efficiency and facilitate the overall improvement of the image processing speed. The output mode of sequentially outputting along the column direction of pixels can be suitable for vertical priority data acquisition scenarios, for example, some detectors can be astronomical spectrometers or remote sensing equipment of specific modes. In the embodiments of the present application, the desired output mode can be selected according to the actual application scenario.

[0054] In some embodiments, the detector can output the detection value of each pixel in the overscan area and the detection value of each pixel in the photosensitive area through a readout channel, so that each detection value corresponds to a time domain output order, and the detection values ​​of different pixels correspond to different time domain output orders. Thus, the two-dimensional array is converted into a one-dimensional array.

[0055] Optionally, if the detection signal outputs detection values ​​through multiple readout channels, a one-dimensional array is generated for each readout channel, thereby obtaining a one-dimensional array corresponding to the multiple readout channels. Subsequently, any one of the one-dimensional arrays can be selected for analysis, and an electromagnetic interference model can be constructed, and relevant parameters of the model can be determined.

[0056] In some embodiments, a one-dimensional array may include the detection value of each pixel in the overscan area and the detection value of each pixel in the photosensitive area. By combining the position of each pixel in the overscan area and the position of each pixel in the photosensitive area, the detection value corresponding to each pixel in the overscan area is extracted from the array, thereby obtaining a detection sequence of electromagnetic interference noise.

[0057] The detection sequence of electromagnetic interference noise arranged in the order of time domain output includes the interference information in the environment where the detector is located. Therefore, the electromagnetic interference model can be accurately constructed by analyzing the detection sequence of electromagnetic interference noise. Among them, the distribution properties of the detection sequence in the time domain include periodic distribution such as the frequency of pulse occurrence and the attenuation of pulse amplitude.

[0058] Specifically, it involves step 103 and step 104, and the detection sequence is modeled and analyzed according to the distribution properties of the detection sequence in the time domain to determine the two parts of the sine distribution and Gaussian distribution contained in the electromagnetic interference noise in the detection environment in the time domain. Then, the electromagnetic interference noise is mathematically modeled to construct an electromagnetic interference model including a sine function and a Gaussian function, that is, an electromagnetic interference model of the electromagnetic interference noise in the time domain. In the electromagnetic interference model, the sinusoidal part and the Gaussian part are in a product operation relationship.

[0059] As a specific example, the operational relationship between the sinusoidal part and the Gaussian part in the electromagnetic interference model can be expressed as: , the sine function part is: ; The Gaussian function part is: .

[0060] Based on the electromagnetic interference model constructed in the embodiment of the present application, the sine function in the electromagnetic interference model can be used to characterize the periodic partial features in the interference information, for example, the periodic interference noise caused by the electronic equipment working in the detection environment. The Gaussian function in the electromagnetic interference model can control the sine wave corresponding to the sine function to fluctuate within the envelope range corresponding to the Gaussian function term, so that the amplitude variation range of the sine wave can be controlled within the envelope corresponding to the Gaussian function term. After determining the specific parameters, the fluctuation of the electromagnetic interference noise can be accurately represented.

[0061] After the electromagnetic interference model is determined, the next step involves step 105, in which the values ​​of relevant parameters in the sine function and the Gaussian function in the electromagnetic interference model are determined according to the detection values ​​of each pixel point in the overscan area.

[0062] For example, the parameters related to the sine function include pulse amplitude, frequency, phase, etc., where the amplitude can represent the interference intensity of the interference source and the frequency can represent the interference interval. The parameters related to the Gaussian function include the width and displacement of the Gaussian part, where the width of the Gaussian part can also be called the extension degree of the pulse.

[0063] When calculating the parameter value, the time domain output order of the detection value in the detection sequence of the electromagnetic interference noise and the size of the detection value can be fitted to determine the value of the relevant parameter in the sine function and the Gaussian function. Among them, the fitting processing method includes but is not limited to the least squares method and other optimization algorithms. Through the fitting processing, the value of the relevant parameter can be adjusted to the preferred value, so that it can accurately represent the fluctuation of the electromagnetic interference noise in one dimension, which is conducive to the subsequent accurate removal of the noise in the image.

[0064] In some embodiments, the detection value of each pixel in the overscan area is a value generated based on the noise in the environment where the detector is located. The noise in the environment where the detector is located also includes background noise. Since the background noise is relatively stable and smooth, the detection value corresponding to the background noise is usually a constant I0. Optionally, the background noise is determined by, for example, selecting a detection value with a relatively stable fluctuation range, and determining the interference value of the background noise by calculating the mean or standard deviation. In the embodiments of the present application, the method for obtaining the interference value of the background noise is not specifically limited.

[0065] Optionally, before determining the values ​​of relevant parameters in the sine function and the Gaussian function, the detection value corresponding to the background noise in the detection sequence of the electromagnetic interference noise can be removed first, and then the parameter values ​​of the relevant parameters in the sine function and the Gaussian function can be determined based on the detection sequence without the background noise.

[0066] Exemplarily, the detection sequence of electromagnetic interference noise without removing the background noise is, for example, {I0, I0, X1, X2, X3, I0, X4, X5, I0, I0, I0}, and the detection sequence after removing the background noise is, for example, {0, 0, X1-I0, X2-I0, X3-I0, 0, X4-I0, X5-I0, 0, 0, 0}.

[0067] The electromagnetic interference model includes at least a sine function part and a Gaussian function part. After determining the values ​​of the relevant parameters in the sine function and the Gaussian function in the electromagnetic interference model, next, step 106 is involved, and the detection value of each pixel point in the photosensitive area can be denoised to remove the interference stripes in the image. The detection value after denoising can be used to generate a clear detection image.

[0068] Specifically, the time domain output order corresponding to each pixel point in the photosensitive area is input into the electromagnetic interference model, and the electromagnetic interference value corresponding to each pixel point in each photosensitive area is determined by the electromagnetic interference model. According to the electromagnetic interference value corresponding to each pixel point in the photosensitive area, the detection value corresponding to each pixel point in the photosensitive area can be denoised to remove the interference fringes in the detection image. For example, Figure 3 is a schematic diagram of a denoising effect of a set of detection images provided in an embodiment of the present application, and a photosensitive area image Figure 3 (a), after denoising, the detection image is obtained as follows Figure 3 (b) shown.

[0069] According to the embodiment of the present application, for the denoising task of the detection image, there is no need for the detector to take multiple shots or to add additional hardware shielding equipment, which effectively reduces the denoising cost. Moreover, since there are sufficient detection values ​​in the overscan area, even for a single shooting task, the electromagnetic interference model has sufficient accuracy to ensure the filtering effect of electromagnetic interference noise, which is beneficial to improving the imaging quality of the detection image and provides a clearer and more accurate image basis for the subsequent analysis of astronomical data.

[0070] Compared with the traditional denoising method, it is necessary to take multiple images and then use median superposition to remove interference fringes. This denoising method also requires that the position and intensity of the interference fringes in each image are fixed. Therefore, once the interference source in the detection environment changes, for example, the interference frequency changes, the position of the same interference signal in each image will also change. This makes the denoising effect of removing interference fringes by median superposition poor. The denoising method provided in the embodiment of the present application can not only be applied to each detection image separately, but also even if the position of the same interference signal in the image changes, it can accurately locate the interference value corresponding to each pixel, thereby achieving accurate denoising of the original detection image.

[0071] In some optional embodiments, the electromagnetic interference model may further include a constant part for representing a corresponding detection value of the background noise in the detection environment.

[0072] As a specific example, the electromagnetic interference model of electromagnetic interference noise in the time domain can be in the following form: (1) Where t is the output order of the detection value, I(t) is the noise interference value corresponding to the t-th detection value, I0 is the interference value of the background noise in the detection environment, the subscript i represents the i-th pulse signal chain, N is a positive integer, and A i is the amplitude of the ith pulse signal chain, ω i and α i are the frequency and phase of the sinusoidal part in the i-th pulse signal chain, σ i and β i are the width and displacement of the Gaussian part in the i-th pulse signal chain, respectively.

[0073] Based on the three parts of constant, sinusoidal distribution and Gaussian distribution, it can not only remove the interference noise caused by equipment such as electromagnetic interference sources, but also remove the background noise. It can further eliminate the internal random noise and system bias, achieve more comprehensive removal of noise in the detection image, and improve the imaging quality of the detection image.

[0074] In some embodiments, the source of electromagnetic interference in the detection environment is a periodically working electronic device such as a motor or an AC circuit. Therefore, the detection environment may include multiple electronic devices that interfere with the detector noise. Before mathematically modeling the electromagnetic interference noise, the detection sequence can be decomposed to obtain one or more pulse signal chains, and different pulse signal chains have different frequencies and amplitudes. A fitting function can be established for each pulse signal chain separately, and then a mathematical model of the detection sequence is constructed based on the fitting functions of each pulse signal chain.

[0075] Exemplarily, taking the detection environment including two interference sources as an example, the detection sequence is {I0, X, I0, Y, I0, X, I0, Y, I0}, where X corresponds to an interference source, including multiple detection values, Y corresponds to another interference source, including multiple detection values, and I0 is the noise value corresponding to the background noise. By decomposing the detection sequence, two pulse signal chains can be obtained, namely {I0, X, I0, I0, I0, X, I0, I0, I0} and {I0, I0, I0, Y, I0, I0, I0, Y, I0}. The two pulse signal chains correspond to a frequency and amplitude respectively.

[0076] In some embodiments, due to the large number of electronic devices in the detection environment, it is not easy to locate which device is causing the interference. Optionally, the interference characteristic information corresponding to one or more hardware devices in the detection environment where the detector is located can be obtained; the one or more pulse signal chains obtained by the decomposition process are matched with the interference characteristic information corresponding to one or more hardware devices, and the successfully matched hardware devices are determined as the interference source of the detector.

[0077] Specifically, the hardware device is an electronic device that works periodically, such as a motor or an AC circuit. The interference characteristic information of the hardware device, for example, the power supply frequency of the AC motor. For each pulse signal chain obtained by the decomposition process, the periodic information corresponding to the pulse signal chain is extracted, such as pulse amplitude, frequency, phase, etc. Based on the matching of the periodic information corresponding to the pulse signal chain with the interference characteristic information of the hardware device, it is possible to identify which specific devices cause interference to the readout channel of the detector. For example, if the power supply frequency of the AC motor is consistent with the pulse signal chain and the AC motor, it can be considered that the AC motor matches the pulse signal chain and the AC motor causes interference to the readout channel of the detector.

[0078] In one example, the matching results can also help troubleshoot potential hardware failures. For example, each pulse signal chain matches a hardware device in the detection environment, but there are still hardware devices in the detection environment that have not been successfully matched. Based on this, it is possible to check whether the unmatched hardware device has a fault to ensure the smooth progress of the entire detection task.

[0079] As a specific example, after fitting the detection sequence, it is determined that the main frequency of the electromagnetic interference is 15.89kHz. Based on the method provided in the embodiment of the present application, the frequencies of each interference source can be accurately found. Through further analysis, it is determined that the electromagnetic interference noise is actually caused by the superposition effect of four interference sources. These interference sources not only have different amplitudes, but the frequency difference between different interference sources is even only 0.1Hz. Since the ceramic resonator is at this level, the crystal oscillator has higher accuracy, and the RC oscillator is far inferior to this, this tiny frequency difference indicates that the interference source may use a ceramic resonator. Furthermore, interference positioning can be performed based on the ceramic resonator, for example, to find out the equipment associated with the ceramic resonator. In this example, it is speculated that the interference source comes from the four power supply systems of the right ascension and declination drive motors. After further inspection of the power supply system by the engineer, it is determined that the motor drive system has a ground wire break fault. As a result, the fault can be repaired in time to ensure the smooth progress of the detection.

[0080] In some embodiments of the present application, in order to further improve the accuracy of the values ​​of relevant parameters in the sine function and the Gaussian function in the electromagnetic interference model, in the step of constructing the electromagnetic interference model, the mathematical model corresponding to each pulse signal chain can be established respectively, that is, the fitting function of the pulse signal chain. Specifically, the detection sequence is modeled and analyzed according to the distribution properties of the detection sequence in the time domain, including: decomposing the detection sequence to obtain one or more pulse signal chains, and different pulse signal chains have different frequencies and amplitudes; modeling and analyzing each pulse signal chain obtained by decomposition, wherein part or all of the pulse signals are respectively included in the sine part and the Gaussian part to construct the fitting function of each pulse signal chain; based on the fitting function of each pulse signal chain, the mathematical model of the detection sequence is constructed; and it is determined that the mathematical model of the detection sequence matches the distribution property of the detection sequence.

[0081] Exemplarily, the fitting function of a single pulse signal chain takes the following form: (2) Where, the subscript i represents the i-th pulse signal chain, t is the output order of the detection value, Ii(t) is the interference value of the i-th electromagnetic interference on the background noise of the t-th output, and A i is the amplitude of the ith pulse signal chain, ω i and α i are the frequency and phase of the sinusoidal part in the i-th pulse signal chain, σ i and β i are the width and displacement of the Gaussian part in the i-th pulse signal chain, respectively.

[0082] For example, Figure 4 is a schematic diagram of a pulse signal chain provided in an embodiment of the present application, Figure 4 Three different pulse signal chains are shown, where the three pulse signal chains have different amplitudes, intervals and frequencies. For each pulse signal chain, by establishing a separate fitting function to represent its fluctuation, it is possible to more flexibly respond to various interference sources in the detection environment. Afterwards, a summation operation relationship between each fitting sub-function is established to determine the electromagnetic interference model, so that the interference situation of the entire detection environment can be more accurately represented, which is conducive to improving the accuracy of the subsequent determination of the corresponding interference value of each pixel point in the photosensitive area.

[0083] In some embodiments, in the step of denoising the detection values ​​of each pixel point in the photosensitive area based on the electromagnetic interference model and generating a detection image of the astronomical environment, the interference value corresponding to each pixel point in the photosensitive area and the detection value of each pixel point in the photosensitive area can be converted into pixel values ​​respectively, and then the pixel value corresponding to the detection value and the pixel value corresponding to the interference value are subtracted to obtain the target pixel value, and the detection image based on the target pixel value is, that is, the detection image after denoising.

[0084] Optionally, during the denoising process, the detection value of each pixel point can be directly subtracted from the interference value corresponding to each pixel point to obtain the target detection value, and then the target detection value is converted into a target pixel value, and then a detection image is constructed based on the target pixel value, that is, the detection image after denoising.

[0085] The detection image obtained by processing according to the embodiment of the present application can accurately remove the interference noise in the photosensitive area image, while also accurately retaining the image information of the astronomical target.

[0086] Corresponding to the method embodiment of the present application, the embodiment of the present application also provides a device for eliminating electromagnetic interference noise in a detection image, including a data acquisition module and a data processing module. The data acquisition module is used to acquire the detection signal of the detector to the astronomical environment, and the data processing module is used to process the detection signal according to the method for eliminating electromagnetic interference noise in the detection image provided in the embodiment of the present application to generate a clear detection image.

[0087] It can be understood that the device for eliminating electromagnetic interference noise in the detection image of the embodiment of the present application can correspond to the execution entity of the method for eliminating electromagnetic interference noise in the detection image provided by the embodiment of the present application. The specific details of the operation and / or function of each module / unit of the device for eliminating electromagnetic interference noise in the detection image can be found in the description of the corresponding parts of the method for eliminating electromagnetic interference noise in the detection image provided by the above-mentioned embodiment of the present application. For the sake of brevity, they will not be repeated here.

[0088] The electronic device in the embodiment of the present application may be a user terminal device, a server, other computing devices, or a cloud server. Figure 5 A schematic diagram of the hardware structure of an electronic device according to an embodiment of the present application is shown. The electronic device may include a processor 501 and a memory 502 storing computer program instructions. When the processor 501 executes the computer program instructions, the process or function of any of the above-mentioned embodiments of the method is implemented.

[0089] Specifically, the processor 501 may include a central processing unit (CPU), or an application-specific integrated circuit (ASIC), or may be configured to implement one or more integrated circuits of the embodiments of the present application. The memory 502 may include a large-capacity memory for data or instructions. For example, the memory 502 may be at least one of the following: a hard disk drive (HDD), a read-only memory (ROM), a random access memory (RAM), a floppy disk drive, a flash memory, an optical disk, a magneto-optical disk, a tape, a universal serial bus (USB) drive, or other physical / tangible memory storage devices. For another example, the memory 502 may include a removable or non-removable (or fixed) medium. For another example, the memory 502 may be inside or outside the integrated gateway disaster recovery device. The memory 502 may be a non-volatile solid-state memory. In other words, the memory 502 generally includes a tangible (non-transitory) computer-readable storage medium (such as a memory device) encoded with computer-executable instructions, and when the software is executed (such as by one or more processors), the operations described in the method of the embodiments of the present application may be performed. The processor 501 implements the process or function of any method in the above embodiments by reading and executing the computer program instructions stored in the memory 502 .

[0090] In one example, Figure 5The electronic device shown may also include a communication interface 503 and a bus 510. Among them, the processor 501, the memory 502, and the communication interface 503 are connected through the bus 510 and complete the communication between each other. The communication interface 503 is mainly used to realize the communication between the modules, devices, units and / or devices in the embodiment of the present application. The bus 510 includes hardware, software or both, and can couple the components of the online data traffic billing device to each other. For example, the bus may include at least one of the following: an accelerated graphics port (AGP) or other graphics bus, an enhanced industrial standard architecture (EISA) bus, a front-end bus (FSB), a hypertransport (HT) interconnect, an industrial standard architecture (ISA) bus, an infinite bandwidth interconnect, a low pin count (LPC) bus, a memory bus, a micro channel architecture (MCA) bus, a peripheral component interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a serial advanced technology attachment (SATA) bus, a video electronics standard association local (VLB) bus or other suitable buses. The bus 510 may include one or more buses. Although the embodiments of the present application describe or illustrate a specific bus, the embodiments of the present application may consider any suitable bus or interconnection method.

[0091] In combination with the method in the above embodiments, an embodiment of the present application also provides a computer-readable storage medium, on which computer program instructions are stored. When the computer program instructions are executed by a processor, the process or function of any method in the above embodiments is implemented.

[0092] In addition, an embodiment of the present application further provides a computer program product, on which computer program instructions are stored. When the computer program instructions are executed by a processor, the process or function of any one of the methods in the above embodiments is implemented.

[0093] The above exemplarily describes the flowcharts and / or block diagrams of the methods, devices, systems and computer program products of the embodiments of the present application, and describes the relevant various aspects. It should be understood that each box or combination thereof in the flowchart and / or block diagram can be implemented by computer program instructions, or by dedicated hardware that performs a specified function or action, or by a combination of dedicated hardware and computer instructions. For example, these computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer or other programmable data processing device to form a machine that enables these instructions executed by such a processor to enable the implementation of the functions / actions specified in each box or combination thereof in the flowchart and / or block diagram. Such a processor can be a general-purpose processor, a special-purpose processor, a special application processor or a field programmable logic circuit.

[0094] The functional blocks shown in the structural block diagram of the embodiment of the present application can be implemented as hardware, software, firmware or a combination thereof. When implemented in hardware, it can be, for example, an electronic circuit, an application-specific integrated circuit (ASIC), appropriate firmware, a plug-in, a function card, etc.; when implemented in software, it is a program or code segment used to perform the required task. The program or code segment can be stored in a memory, or transmitted on a transmission medium or a communication link via a data signal carried in a carrier. The code segment can be downloaded via a computer network such as the Internet, an intranet, etc.

[0095] It should be noted that the present application is not limited to the specific configurations and processes described above or shown in the figures. The above is only a specific implementation mode of the present application. Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working process of the described system, device, module or unit can refer to the corresponding process in the method embodiment without further description. It should be understood that the scope of protection of the present application is not limited to this. Any technician familiar with the technical field can think of various equivalent modifications or substitutions within the technical scope disclosed in the present application, and these modifications or substitutions should be included in the scope of protection of the present application.

Claims

1. A method for eliminating electromagnetic interference noise in a detection image, characterized in that: include: Acquire a detection signal of the detector to the astronomical environment, wherein the detection signal includes a detection value of each pixel point in the overscan area and a detection value of each pixel point in the photosensitive area; Extracting the detection value of each pixel point in the overscan area from the detection signal, and generating a detection sequence of electromagnetic interference noise according to the time domain output order of the detection value; Modeling and analyzing the detection sequence according to the distribution properties of the detection sequence in the time domain to determine the two parts of the sine distribution and the Gaussian distribution contained in the electromagnetic interference noise in the detection environment in the time domain; Mathematically modeling the electromagnetic interference noise to construct an electromagnetic interference model of the electromagnetic interference noise in the time domain, wherein the electromagnetic interference model includes a sine function and a Gaussian function; Determining values ​​of relevant parameters in the sine function and the Gaussian function in the electromagnetic interference model according to the detection values ​​of each pixel point in the overscan area; Based on the electromagnetic interference model, the detection value of each pixel point in the photosensitive area is denoised to generate a detection image of the astronomical environment.

2. The method according to claim 1, characterized in that The detector includes at least one of the following: a detector installed in an astronomical telescope, a detector carried in a space imaging satellite, a detector used in a harsh environment, a CCD detector, and a CMOS detector.

3. The method according to claim 1, characterized in that Before mathematically modeling the electromagnetic interference noise, the method further includes: The detection sequence is decomposed to obtain one or more pulse signal chains, and different pulse signal chains have different frequencies and amplitudes.

4. The method according to claim 3, characterized in that The method further comprises: Obtain interference characteristic information corresponding to one or more hardware devices in the detection environment where the detector is located; The one or more pulse signal chains obtained by the decomposition process are matched with the interference characteristic information corresponding to the one or more hardware devices, and the hardware devices with successful matching are determined as interference sources of the detector.

5. The method according to claim 3, characterized in that: The electromagnetic interference model of the electromagnetic interference noise in the time domain is in the following form: Where t is the output order of the detection value, I(t) is the noise interference value corresponding to the t-th detection value, I0 is the interference value of the background noise in the detection environment, the subscript i represents the i-th pulse signal chain, N is a positive integer, and A i is the amplitude of the ith pulse signal chain, ω i and α i are the frequency and phase of the sinusoidal part in the i-th pulse signal chain, σ i and β i are the width and displacement of the Gaussian part in the i-th pulse signal chain, respectively.

6. The method according to claim 1, characterized in that The modeling and analyzing the detection sequence according to the distribution property of the detection sequence in the time domain includes: Decomposing the detection sequence to obtain one or more pulse signal chains, wherein different pulse signal chains have different frequencies and amplitudes; Modeling and analyzing each pulse signal chain obtained by decomposition, wherein part or all of the pulse signals are respectively made to include a sinusoidal part and a Gaussian part, so as to construct a fitting function of each pulse signal chain; Constructing a mathematical model of the detection sequence based on the fitting functions of each pulse signal chain; A mathematical model of the detection sequence is determined to match a distribution property of the detection sequence.

7. The method according to claim 6, characterized in that The fit function for a single pulse signal chain takes the following form: Where, the subscript i represents the i-th pulse signal chain, t is the output order of the detection value, Ii(t) is the interference value of the i-th electromagnetic interference on the t-th output, and A i is the amplitude of the ith pulse signal chain, ω i and α i are the frequency and phase of the sinusoidal part in the i-th pulse signal chain, σ i and β i are the width and displacement of the Gaussian part in the i-th pulse signal chain, respectively.

8. An electronic device, characterized in that: The electronic device comprises: a processor and a memory storing computer program instructions; when the electronic device executes the computer program instructions, the method according to any one of claims 1 to 7 is implemented.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer program instructions, and when the computer program instructions are executed by a processor, the method according to any one of claims 1 to 7 is implemented.

10. A computer program product, characterized in that, It comprises computer program instructions, which, when executed by a processor, implement the method according to any one of claims 1 to 7.

Images