WO2022184167A1 - Imaging method and apparatus, device, and storage medium - Google Patents

Imaging method and apparatus, device, and storage medium Download PDF

Info

Publication number
WO2022184167A1
WO2022184167A1 PCT/CN2022/079318 CN2022079318W WO2022184167A1 WO 2022184167 A1 WO2022184167 A1 WO 2022184167A1 CN 2022079318 W CN2022079318 W CN 2022079318W WO 2022184167 A1 WO2022184167 A1 WO 2022184167A1
Authority
WO
WIPO (PCT)
Prior art keywords
light intensity
intensity value
pixel point
image
time
Prior art date
Application number
PCT/CN2022/079318
Other languages
French (fr)
Chinese (zh)
Inventor
熊瑞勤
赵菁
黄铁军
Original Assignee
脉冲视觉(北京)科技有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from CN202110240363.2A external-priority patent/CN113034634B/en
Priority claimed from CN202110253200.8A external-priority patent/CN113067979A/en
Application filed by 脉冲视觉(北京)科技有限公司 filed Critical 脉冲视觉(北京)科技有限公司
Publication of WO2022184167A1 publication Critical patent/WO2022184167A1/en

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/60Control of cameras or camera modules

Definitions

  • the present application relates to the technical field of image processing, and in particular, to an imaging method, apparatus, device, and storage medium.
  • Pulse cameras represent visual information in the form of pulse arrays and can continuously record changes in light intensity without the concept of frame rate and exposure time, breaking through the limitations of traditional cameras. The texture details in the scene can be reconstructed.
  • the image reconstruction algorithm of the pulse camera includes two types: the image reconstruction algorithm based on the pulse interval and the image reconstruction algorithm based on the time window averaging.
  • the image reconstruction algorithm based on the pulse interval utilizes the characteristic that the pulse interval decreases with the increase of the light intensity, and reconstructs the light intensity in a short period of time with the help of the two pulses before and after.
  • this algorithm can describe the outline of high-speed moving objects, the reconstructed signal is usually not stable enough, and the pixel value has obvious fluctuations in the time direction;
  • the image reconstruction algorithm based on the time window average uses the pulse firing frequency to increase with the increase of light intensity. This feature is used to estimate the average light intensity within a time window.
  • the algorithm improves the signal-to-noise ratio of the reconstructed image to a certain extent, when there is motion of the object, the average in the time window will cause motion blur in the reconstructed image.
  • Embodiments of the present application provide an imaging method, apparatus, device, and storage medium, which can improve the quality of reconstructed images and obtain clear images in high-speed motion scenes.
  • a brief summary is given below. This summary is not intended to be an extensive review, nor is it intended to identify key/critical elements or delineate the scope of protection of these embodiments. Its sole purpose is to present some concepts in a simplified form as a prelude to the detailed description that follows.
  • the embodiments of the present application provide an imaging method, comprising:
  • filtering in the time dimension is performed on each pixel point of the image frame in the image sequence that meets the specified condition to update the light intensity value of each pixel point, and obtain the updated light intensity value;
  • An updated image is generated according to the updated light intensity value of each pixel.
  • the specified condition includes time k, which is within a preset time period.
  • filtering in the time dimension is performed on each pixel of the image frame in the image sequence that satisfies the specified condition to update the light intensity value of each pixel to obtain the updated light intensity values, including:
  • each pixel point is filtered in the time dimension to obtain the updated light intensity value.
  • filtering in the time dimension is performed on each pixel point to obtain the updated light intensity value, including:
  • the following formula is used to filter each pixel point in the time dimension to obtain the updated light intensity value
  • I k (x, y) represents the updated light intensity value at the (x, y) pixel point in the image frame at time k
  • C represents the regularization parameter
  • r represents the radius of the filtering time window.
  • the regularization parameter C is determined by the following formula:
  • is the setting parameter.
  • the setting parameters are determined by the following formula:
  • the imaging method further includes:
  • the corresponding pixel is an integer pixel, directly take the light intensity value at the corresponding pixel;
  • the light intensity value of the corresponding pixel is obtained by interpolation.
  • filtering in the time dimension is performed on each pixel of the image frame in the image sequence that satisfies the specified condition to update the light intensity value of each pixel to obtain the updated light intensity values, including:
  • the autoregressive model of the pixels in the image sequence in the time direction is adaptively learned to determine the model parameters of the autoregressive model
  • each pixel point is filtered in the time dimension to obtain the updated light intensity value.
  • the formula of the autoregressive model is:
  • I k represents the theoretical light intensity value at time k
  • I t represents the theoretical light intensity value at time t
  • T k represents the set of adjacent moments at time k
  • P k ⁇ t (x, y) is the image frame at time k
  • ⁇ t is the model parameter
  • ⁇ k (x, y) is the model error on the (x, y) pixel point.
  • adaptive learning is performed on the autoregressive model of the pixel points in the image sequence in the time direction to determine the model parameters of the autoregressive model, including:
  • the autoregressive model is adaptively learned by the least square method to solve the objective function to determine the model parameters.
  • the objective function is:
  • W represents the set of pixels in the spatial range centered on (x, y); Represents the estimated light intensity value at time k; represents the estimated light intensity value at time t.
  • filtering in the time dimension is performed on each pixel point to obtain the updated light intensity value, including:
  • the following formula is used to filter each pixel point in the time dimension to obtain the updated light intensity value
  • I′ k represents the updated light intensity value at time k.
  • the imaging method further includes: obtaining the relative motion between the image sequences according to any one or more of optical flow method, pixel point matching, pixel point motion alignment, and pixel point relative position offset estimation relation.
  • the imaging method before obtaining the reconstructed image sequence according to the pulse array within a preset time period, the imaging method further includes:
  • the pixel-based pulse array of the object to be imaged is obtained.
  • a reconstructed image sequence is obtained according to the pulse array within a preset time period, including:
  • the image sequence is reconstructed using the pulse interval algorithm.
  • an imaging device including:
  • a reconstruction module configured to obtain a reconstructed image sequence according to the pulse array within a preset time period
  • the update module is used to filter each pixel point of the image frame that meets the specified condition in the image sequence in the time dimension according to the relative motion relationship between the image sequences to update the light intensity value of each pixel point, and obtain the updated light intensity value ;
  • the imaging module is configured to generate an updated image according to the updated light intensity value of each pixel point.
  • embodiments of the present application provide an electronic device, including a processor and a memory storing program instructions, where the processor is configured to execute the imaging method provided by the foregoing embodiments when executing the program instructions.
  • the embodiments of the present application provide an imaging device, including the electronic device provided by the above embodiments.
  • the embodiments of the present application provide a computer-readable medium on which computer-readable instructions are stored, and the computer-readable instructions can be executed by a processor to implement the imaging method provided by the foregoing embodiments.
  • a reconstructed image sequence is obtained according to a pulse array within a preset time period; according to the relative motion relationship between the image sequences, each pixel of an image frame that meets a specified condition in the image sequence is analyzed. Filtering in the time dimension is performed to update the light intensity value of each pixel point to obtain the updated light intensity value; an updated image is generated according to the updated light intensity value of each pixel point.
  • the method proposes to filter the image pixels along the moving trajectory of the object in the temporal dimension, which can be applied to high-speed moving scenes to generate images with high signal-to-noise ratio and no blur, and achieve better visual effects. The problem of clear imaging of high-speed moving objects.
  • FIG. 1 is a schematic flowchart of an imaging method according to an exemplary embodiment of the present application.
  • FIG. 2 is a schematic diagram of a pulse array according to an exemplary embodiment of the present application.
  • FIG. 3 is a schematic diagram of motion alignment according to an exemplary embodiment of the present application.
  • FIG. 4 is a schematic diagram of calculating a light intensity value according to an interpolation method according to an exemplary embodiment of the present application.
  • FIG. 5 is a schematic flowchart of an imaging method according to another exemplary embodiment of the present application.
  • FIG. 6 is a schematic flowchart of an imaging method according to another exemplary embodiment of the present application.
  • FIG. 7 is a schematic flowchart of an imaging method according to another exemplary embodiment of the present application.
  • FIG. 8 is a schematic flowchart of an imaging method according to another exemplary embodiment of the present application.
  • FIG. 9 is a schematic diagram of an imaging picture taken by different imaging methods according to an exemplary embodiment of the present application.
  • FIG. 10 is a schematic diagram of an imaging picture taken by different imaging methods according to another exemplary embodiment of the present application.
  • FIG. 11 is a schematic structural diagram of an imaging device according to an exemplary embodiment of the present application.
  • FIG. 12 is a schematic structural diagram of an electronic device according to an exemplary embodiment of the present application.
  • FIG. 13 is a schematic diagram of a computer storage medium according to an exemplary embodiment of the present application.
  • the imaging algorithm of the pulse camera based on the pulse signal can describe the contour of the high-speed moving object
  • the reconstructed signal is usually not stable and perfect, and the pixel value in the time direction has obvious fluctuations, or the reconstructed image has motion blur. question. Therefore, a new imaging method is urgently needed to solve the above problem of unclear imaging of pulse cameras in high-speed motion scenes.
  • pulse arrays have correlations in the temporal direction, and that correlations behave differently in different spatial locations of the image. For example, in a smooth motion area with relatively continuous motion, the temporal correlation of the image content along the motion trajectory is relatively strong; while for an area with relatively complex local motion, the image content at a certain moment may be different from the image content at the previous and subsequent times. There is a corresponding relationship, and even occlusion or disappearance may occur. At this time, the correlation of the time direction of the image content along the motion trajectory is relatively weak.
  • an embodiment of the present application proposes an imaging method.
  • the imaging method includes the following contents.
  • S110 Obtain a reconstructed image sequence according to the pulse array within a preset time period.
  • the imaging method further includes: acquiring a pixel-based pulse array of the object to be imaged according to the pulse signal within a preset time period.
  • the embodiments of the present application may use a bionic pulse camera for imaging, and each pixel in the bionic pulse camera includes a photoreceptor, an integrator, and a threshold comparator.
  • the photoreceptor is responsible for converting the optical signal into an electrical signal
  • the integrator accumulates the converted photo-generated charges
  • the threshold comparator repeatedly compares the accumulated charges at an ultra-high frequency. When the accumulated charges reach a preset threshold , the pixel is pulsed and the integrator is cleared.
  • Figure 2 is a schematic diagram of a pulse array.
  • the pulse camera represents visual information in the form of "H*W*T” pulse array
  • H*W is the spatial resolution of the pulse camera
  • T is the sampling times.
  • the pulse array consists of two symbols, "0" and "1", where “1" (solid dot) indicates that there is a pulse firing at this spatiotemporal position; "0” (open dot) means that there is no pulse firing at this spatiotemporal position.
  • the bionic pulse camera can continuously record changes in light intensity without the concept of frame rate and exposure time, breaking through the limitations of traditional cameras.
  • a basic image sequence is reconstructed according to the acquired pulse array within a preset time period.
  • a pulse interval algorithm can be used to reconstruct the image sequence, or other algorithms can be used to estimate the basic reconstructed image sequence according to the pulse array.
  • the image sequence is reconstructed using a pulse interval algorithm, and the pulse interval is the time interval between two adjacent pulses at the same pixel position, which is equivalent to the integration time used by the integrator to accumulate the next pulse.
  • the integrator of the pixel in the strong illumination area can accumulate the photogenerated charge to the pulse emission threshold in a shorter time, and the pulse interval is small; the integrator of the pixel in the weak illumination area takes a longer time to accumulate the charge to Threshold of pulse release, pulse interval is larger.
  • the basic image sequence is reconstructed
  • the basic image sequence can be reconstructed according to the pulse array using the pulse interval algorithm.
  • S120 According to the relative motion relationship between the image sequences, perform filtering in the time dimension on each pixel of the image frame in the image sequence that meets the specified condition to update the light intensity value of each pixel to obtain an updated light intensity value.
  • the specified condition may be a certain moment, or an image frame with a certain feature.
  • the specified condition includes time k, which is within a preset time period.
  • the motion alignment of the image sequences can be performed to obtain the corresponding relationship between pixels in different image frames in the image sequence.
  • the camera position is fixed, with the movement of the object, a certain point in the object has different pixel positions in the image frames at different times. Therefore, for a certain pixel in the image frame at time k, corresponding pixels in the image frame at other times can be obtained.
  • the motion trajectory of each pixel in different image frames can be obtained.
  • the light intensity value of the pixel point at different times can be determined. For example, for each pixel in the image frame at time k, filtering in the time dimension can be performed on the pixel based on the corresponding light intensity values of the pixel in the image frames at different times to update the light intensity of the pixel. Intensity value, get updated light intensity value. An updated image can be generated according to the updated light intensity value of each pixel in the image frame at time k.
  • the imaging method provided by the embodiment of the present application utilizes the temporal direction correlation of the pulse array to improve the quality of the reconstructed image.
  • filtering in the time dimension is performed on each pixel point of the image frame in the image sequence that meets the specified condition to update the light intensity value of each pixel point, so as to obtain the updated light Intensity value, and then generate an updated image according to the updated light intensity value of each pixel, so that the quality of the reconstructed image can be improved, a clear image in a high-speed motion scene can be obtained, and a better visual effect can be achieved.
  • step S120 includes: determining the corresponding pixel points of each pixel point in the image frame at time t of the image sequence according to the relative motion relationship; Dimensional filtering to get updated light intensity values.
  • the corresponding pixel of the pixel in the image frame at other times may be determined based on the relative motion relationship.
  • the light intensity value of the corresponding pixel point can be considered as the light intensity value corresponding to the pixel point in the image frame at time k at other time points in the image frame.
  • the imaging method further includes: obtaining the relative relationship between the image sequences according to any one or more of optical flow method, pixel point matching, pixel point motion alignment, and pixel point relative position offset estimation. sports relationship.
  • the relative motion relationship between image sequences is obtained according to the optical flow method.
  • the current frame is calculated using an optical flow algorithm such as Horn–Schunck
  • Horn–Schunck The optical flow field between the r frame images before and after it is used to determine the relative motion relationship between the images.
  • the optical flow field can be understood as a motion field, and the apparent motion of the image brightness pattern is the optical flow.
  • the optical flow method is mostly based on the principle of light intensity consistency, that is, it is assumed that the same object has the same light intensity value.
  • the (x, y) point of the first image corresponds to the (x', y') point of the second image
  • UV_t ⁇ k(x, y) is the optical flow at the (x, y) point (describes the two relative motion at that point between frames).
  • the optical flow method can be used to solve the optical flow (relative motion) between two frames based on the assumption of light intensity consistency, and further, based on the optical flow (relative motion), the correspondence between images between different frames can be determined, That is, motion alignment, get each pixel (x, y) in the current frame in other the corresponding position on the frame
  • Figure 3 is a schematic diagram of motion alignment. As shown in Figure 3, the positions of free-falling triangles on different images are different. Based on the optical flow information, we can determine the correspondence between the images and find the same object point in the scene at different time points. The position on the object determines the trajectory of the object.
  • the imaging method further includes: Obtain the light intensity value of the pixel point of the current frame at the corresponding pixel point on other frame images (obtain the light intensity value of the corresponding pixel point).
  • acquiring the light intensity value of the corresponding pixel point includes: if the corresponding pixel point is an integer pixel point, directly obtaining the light intensity value at the corresponding pixel point; if the corresponding pixel point is a sub-pixel point, using an interpolation method Get the light intensity value of the corresponding pixel point.
  • the light intensity value at the corresponding pixel is directly taken; if the pixel of the current frame corresponds to the pixel in other frame images is a sub-pixel point, then the light intensity value at the corresponding pixel point is obtained by interpolation.
  • Figure 4 is a schematic diagram of calculating the light intensity value using the interpolation method. As shown in Figure 4, assuming that after the motion is aligned, the corresponding point of the current point on other frame images is not at an integer position, then it is necessary to use the surrounding whole point. light intensity value to estimate the light intensity value at that point. Taking bilinear interpolation as an example, P is a non-integer point, and the light intensity value of this point is unknown, then the brightness information of the surrounding integer points Q 11 , Q 12 , Q 21 , and Q 22 can be used to interpolate the light intensity of point P. light intensity value.
  • the specific methods are as follows:
  • f(.) represents the light intensity value of a certain point.
  • filtering each pixel point in the time dimension to obtain the updated light intensity value includes: according to the light intensity value of the corresponding pixel point, through the following formula, for each pixel point The pixel points are filtered in the time dimension to obtain the updated light intensity value,
  • I k (x, y) represents the updated light intensity value at the (x, y) pixel point in the image frame at time k
  • C represents the regularization parameter
  • r represents the radius of the filtering time window.
  • Weights The selection of is gradually reduced with the increase of the time interval between two frames, which can be calculated by the following expression or other similar expressions.
  • the following formula can be used to filter each pixel point in the time dimension to obtain the updated light intensity value
  • T represents the total number of frames in the image sequence, 1 ⁇ k ⁇ T. It can be determined according to the method provided in the foregoing embodiment, or determined by other methods according to the actual situation.
  • the motion trajectory of the object can be obtained, and the signal is filtered along the time dimension along the motion trajectory of the object, and the image frame at time k can be updated (x, y) The light intensity value on the pixel point.
  • motion trajectory filtering in the time dimension can be performed on the basic reconstructed image sequence to improve the stability of the signal, thereby improving the quality of the reconstructed image.
  • the motion trajectory filtering along the time dimension realizes the utilization of the time correlation of the pulse array, and the motion alignment ensures the accuracy of the time direction correlation used for the moving object.
  • FIG. 5 is an imaging method provided by another embodiment of the present application, specifically an imaging method based on a bionic pulse camera.
  • FIG. 5 is an example of the embodiment of FIG. 1 . As shown in FIG. 5 , the method specifically includes the following contents.
  • S510 Obtain a reconstructed image sequence according to the pulse array within a preset time period.
  • S530 Perform motion alignment on the image sequence according to the relative motion relationship to obtain the corresponding relationship between pixels in different image frames.
  • S540 According to the correspondence between the pixels of different image frames, perform filtering on the pixels of the image sequence in the time dimension, and update the light intensity values of the pixels of the image to be reconstructed to generate an updated image.
  • one frame of images or multiple frames of images in the image sequence can be reconstructed.
  • the method includes the following contents.
  • a bionic pulse camera is used for imaging, and a pixel-based pulse array of the object to be imaged is obtained according to the pulse signal within a preset time period. Then, based on the pulse array, a pulse interval algorithm is used to obtain a reconstructed image sequence (basic reconstruction).
  • the optical flow estimation algorithm such as Horn-Schunck is used to calculate the optical flow field between the current frame and its previous and previous frame images, so as to determine the relative motion relationship between the images.
  • the optical flow method can be used to solve the optical flow (relative motion) between two frames based on the assumption of light intensity consistency, and further, based on the optical flow (relative motion), the correspondence between images between different frames can be determined, i.e. motion alignment. Get each pixel (x, y) in the current frame in other the corresponding position on the frame
  • the method further includes: acquiring the light intensity values of the pixels of the current frame at the corresponding pixels of other frame images.
  • the motion trajectory of the object can be obtained, and the signal is filtered in the time dimension along the motion trajectory of the object, and the light intensity value of the pixel is updated to generate an updated image.
  • the image sequence reconstruction is performed on the scene at different times according to the characteristics of the pulse array; then, on the basis of the reconstructed image sequence, the optical flow estimation algorithm is used to determine the motion trajectory of the object along the time direction, along the respective The motion track performs motion alignment on the image sequence; finally, the reconstructed image sequence is subjected to time dimension motion track filtering to improve the stability of the signal, thereby improving the quality of the reconstructed image, obtaining clear images in high-speed motion scenes, and achieving better visual effects. , which solves the problem that high-speed moving objects cannot be clearly imaged in the prior art.
  • the imaging method based on the bionic pulse camera provided in the embodiment of the present application can be applied to a high-speed motion scene, and can generate an image with a high signal-to-noise ratio and no blur.
  • step S120 includes: performing adaptive learning on the autoregressive model of the pixels in the image sequence in the time direction according to the relative motion relationship to determine the model parameters of the autoregressive model; The corresponding pixel points in the image frame at time t of the sequence; according to the model parameters and the light intensity value of the corresponding pixel point, each pixel point is filtered in the time dimension to obtain the updated light intensity value.
  • the autoregressive model can be set in advance as required.
  • the formula of the autoregressive model is:
  • I k represents the theoretical light intensity value (image theoretical value) at time k
  • I t represents the theoretical light intensity value (image theoretical value) at time t
  • T k represents the set of adjacent moments at time k
  • P k ⁇ t ( x, y) is the corresponding pixel position of the pixel (x, y) in the image frame at time k in the image frame at time t
  • ⁇ t is the model parameter
  • ⁇ k (x, y) is the pixel at (x, y) ) model error on pixels, which may correspond to some image details or noise.
  • adaptive learning is performed on the autoregressive model of the pixel points in the image sequence in the time direction, so as to determine the model parameters of the autoregressive model, including: according to the relative motion relationship, through the least two The method of multiplying to solve the objective function performs adaptive learning on the autoregressive model to determine the model parameters.
  • the pixels on the image sequence can be adaptively selected, the autoregressive model can be adaptively learned, and the model parameters of the autoregressive model can be determined.
  • the objective function is solved by the least squares method to obtain the model parameters of the autoregressive model.
  • the objective function can be:
  • W represents the set of pixels in the spatial range centered on (x, y); Represents the estimated light intensity value at time k (image estimated value); Indicates the estimated light intensity value (image estimated value) at time t.
  • the temporal correlation of the image content along the motion trajectory is different. Therefore, it is necessary to adjust the parameters of the autoregressive model adaptively according to the content of the local region of the image sequence, so as to determine the suitable local Filter weights for regions.
  • the filter weight is related to the content of the pulse array, and the filter weight at different spatial positions of an image may be different.
  • the present application assumes that pixels with similar spatial positions have relatively consistent temporal correlations, so the model parameters can be adaptively selected according to the data at the local spatial positions, and can be obtained by optimizing the above objective function.
  • each pixel point in the time dimension to obtain the updated light intensity value, including: according to the model parameter and the light intensity value of the corresponding pixel point, Through the following formula, each pixel is filtered in the time dimension to obtain the updated light intensity value,
  • I′ k represents the updated light intensity value at time k, or is called image adaptive imaging value.
  • this formula is used to filter the pixel points in the time dimension, and the obtained result is used as the updated light intensity value on the (x, y) pixel point in the image frame at time k, which can be finally obtained.
  • Adaptive imaging of the object is obtained.
  • an autoregressive model of each pixel in the time direction may be established according to the motion trajectory of each pixel of the image sequence.
  • the motion trajectory of each pixel of the image sequence may be determined according to the relative motion relationship between the image sequences, and reference may be made to the description in the foregoing embodiments for the determination of the motion trajectory. Then the pixels on the image sequence are adaptively selected, and the autoregressive model is adaptively learned to determine the model parameters.
  • FIG. 7 is an imaging method provided by another embodiment of the present application, specifically an adaptive imaging method based on a pulse signal.
  • FIG. 7 is an example of the embodiment of FIG. 1 . As shown in FIG. 7 , the method specifically includes the following contents.
  • S710 Obtain a reconstructed image sequence according to the pulse array within a preset time period.
  • S720 Determine the motion trajectory of each pixel of the image sequence according to the relative motion relationship between the image sequences.
  • S730 Establish an autoregressive model of each pixel in the time direction according to the motion trajectory of each pixel in the image sequence.
  • an appropriate autoregressive model can be established according to the motion trajectory of each pixel point.
  • S740 Adaptively select pixels on the image sequence, perform adaptive learning on the autoregressive model, and determine model parameters of the autoregressive model.
  • S750 Perform filtering in the time dimension on the pixels of the image sequence according to the model parameters, and update the light intensity values of the pixels of the image to be reconstructed to generate an updated image.
  • one frame of images or multiple frames of images in the image sequence can be reconstructed.
  • the method includes the following contents.
  • a basic reconstructed image sequence (reconstructed basic image) is estimated from the pulse array based on the pulse interval information.
  • the basic reconstructed image is generated by the pulse array, and a pulse spacing algorithm or other reconstruction method can be used.
  • the basic image sequence is reconstructed: in, represents the reconstructed image sequence, represents each reconstructed image.
  • the motion trajectory of each pixel point of the image is determined (establishing the motion trajectory) according to the basic reconstructed image sequence by means of methods such as optical flow.
  • an autoregressive model in the time dimension is established along the direction of the motion trajectory (establishing an autoregressive model), and the model parameters are adaptively adjusted according to the temporal correlation of neighboring pixels (adaptive learning model parameters).
  • the motion trajectory filtering of the time dimension (time dimension filtering) is performed on the image pixels to reconstruct a higher quality image.
  • the present application establishes an autoregressive model between image pixels at different times along the motion trajectory of the object, and adaptively adjusts the model parameters according to the content of the image sequence, so as to accurately utilize the temporal direction correlation of the image signal, thereby improving the reconstructed image. the quality of.
  • the adaptive imaging method based on pulse signals of the present application obtains a pixel-based pulse array of an object to be imaged according to the pulse signals within a preset time period; obtains a reconstructed image sequence within a preset time period according to the pulse array;
  • the relative motion relationship between the image sequences determines the motion trajectory of each pixel in the reconstructed image sequence; according to the motion trajectory of each pixel, the autoregressive model of each pixel in the time direction is established; the pixels on the reconstructed image sequence are selected , the autoregressive model is adaptively learned, and the model parameters of the autoregressive model are determined; according to the model parameters, the pixel points of the reconstructed image sequence are filtered in the time dimension, and the light intensity value of the pixel points is updated to generate an updated image.
  • the present application adaptively utilizes the temporal direction correlation of the pulse array to improve the quality of the reconstructed image.
  • an autoregressive model along the time direction of the motion trajectory is established; then, the model parameters are adaptively adjusted according to the correlation structure of the local spatial position of the pulse array , in order to improve the accuracy of the autoregressive model; finally, the basic reconstructed image is filtered in the time dimension by the established local adaptive autoregressive model to improve the stability of the signal, thereby improving the quality of the reconstructed image.
  • the locally adaptive autoregressive model ensures the accuracy of the used time-direction correlation, effectively reduces the influence of outliers, and ensures the effectiveness and robustness of the filtering algorithm.
  • the imaging method provided by the present application such as an imaging method based on a biomimetic pulse camera, or an adaptive imaging method based on a pulse signal, will be further described through specific application scenarios.
  • FIG. 9 is a schematic diagram of a car traveling at a speed of 100 km/h shot by different imaging methods according to an exemplary embodiment.
  • (a) is an initial pulse matrix image
  • (b) is an image reconstructed using the pulse interval method
  • (c) is an image reconstructed using the imaging method provided by the embodiment of the present application. It can be seen from the pictures that, according to the image reconstruction algorithm of the embodiment of the present application, a clear image of a vehicle traveling at a high speed can be obtained, and a better visual effect can be achieved.
  • FIG. 10 is a schematic diagram of a free-falling ragdoll photographed by different imaging methods according to an exemplary embodiment.
  • (a) is an initial pulse matrix image
  • (b) is an image reconstructed using the pulse interval method
  • (c) is an image reconstructed using the imaging method provided by the embodiments of the present application. It can be seen from the images that the embodiments of the present application can clearly reconstruct the motion details in the scene, and the reconstructed images not only avoid motion blur but also have a high signal-to-noise ratio.
  • An embodiment of the present application further provides an imaging apparatus, which is used for executing the imaging method of the foregoing embodiment.
  • the apparatus includes: a reconstruction module 1110 , an update module 1120 and an imaging module 1130 .
  • the reconstruction module 1110 is used for obtaining the reconstructed image sequence according to the pulse array in the preset time period; the updating module 1120 is used for, according to the relative motion relationship between the image sequences, for the image frames that meet the specified conditions in the image sequence Each pixel of the pixel is filtered in the time dimension to update the light intensity value of each pixel to obtain an updated light intensity value; the imaging module 1130 is configured to generate an updated image according to the updated light intensity value of each pixel.
  • the specified condition includes time k, which is within a preset time period.
  • the updating module 1120 is configured to: determine the corresponding pixel points of each pixel point in the image frame at time t of the image sequence according to the relative motion relationship; Filtering in the time dimension to get updated light intensity values.
  • the update module 1120 is configured to: perform filtering on each pixel in the time dimension according to the light intensity value of the corresponding pixel point by the following formula to obtain the updated light intensity value,
  • I k (x, y) represents the updated light intensity value at the (x, y) pixel point in the image frame at time k
  • C represents the regularization parameter
  • r represents the radius of the filtering time window.
  • the regularization parameter C is determined by the following formula:
  • is the setting parameter.
  • the setting parameters are determined by the following formula:
  • the imaging device further includes: a first acquiring module 1140, configured to acquire the light intensity value of the corresponding pixel point.
  • the first obtaining module 1140 is configured to: if the corresponding pixel is an integer pixel, directly obtain the light intensity value at the corresponding pixel; if the corresponding pixel is a sub-pixel, obtain the corresponding pixel by interpolation The light intensity value of the pixel.
  • the updating module 1120 is configured to: perform adaptive learning on the autoregressive model of the pixel points in the image sequence in the time direction according to the relative motion relationship, so as to determine the model parameters of the autoregressive model; The corresponding pixel points in the image frame at time t of the image sequence; according to the model parameters and the light intensity value of the corresponding pixel point, each pixel point is filtered in the time dimension to obtain the updated light intensity value.
  • the formula of the autoregressive model is:
  • I k represents the theoretical light intensity value at time k
  • I t represents the theoretical light intensity value at time t
  • T k represents the set of adjacent moments at time k
  • P k ⁇ t (x, y) is the image frame at time k
  • ⁇ t is the model parameter
  • ⁇ k (x, y) is the model error on the (x, y) pixel point.
  • the update module 1120 is configured to: perform adaptive learning on the autoregressive model by solving the objective function by the least squares method according to the relative motion relationship, so as to determine the model parameters.
  • the objective function is:
  • W represents the set of pixels in the spatial range centered on (x, y); Represents the estimated light intensity value at time k; represents the estimated light intensity value at time t.
  • the update module 1120 is configured to: filter each pixel in the time dimension to obtain the updated light intensity value by using the following formula according to the model parameters and the light intensity value of the corresponding pixel point,
  • I′ k represents the updated light intensity value at time k.
  • the imaging device further includes: a second acquisition module 1150, configured to: according to any one or more of optical flow method, pixel point matching, pixel point motion alignment and pixel point relative position offset estimation, Obtain the relative motion relationship between image sequences.
  • a second acquisition module 1150 configured to: according to any one or more of optical flow method, pixel point matching, pixel point motion alignment and pixel point relative position offset estimation, Obtain the relative motion relationship between image sequences.
  • the imaging device further includes: a third acquisition module 1160, configured to: acquire a pixel-based pulse array of the object to be imaged according to the pulse signal within a preset time period.
  • the reconstruction module 1110 is configured to: reconstruct an image sequence by using a pulse interval algorithm according to the pulse array.
  • the imaging apparatus provided in the above-mentioned embodiments executes the imaging method
  • only the division of the above-mentioned functional modules is used as an example.
  • the internal structure of the device is divided into different functional modules to complete all or part of the functions described above.
  • the imaging apparatus and the imaging method embodiments provided by the above embodiments belong to the same concept, and the specific implementation process and effects thereof are described in the description of the method embodiment section, which will not be repeated here.
  • the embodiment of the present application further provides an electronic device corresponding to the imaging method provided by the foregoing embodiments, so as to execute the foregoing imaging method.
  • FIG. 12 shows a schematic diagram of an electronic device provided by some embodiments of the present application.
  • the electronic device includes: a processor 1210 and a memory 1220 .
  • the memory 1220 stores a computer program that can be executed on the processor 1210.
  • the processor 1210 runs the computer program, the imaging method provided by any of the foregoing embodiments of the present application is executed.
  • the electronic device also includes a bus 1230 and a communication interface 1240 .
  • the processor 1210 , the communication interface 1240 and the memory 1220 are connected through the bus 1230 .
  • the memory 1220 may include a high-speed random access memory (RAM: Random Access Memory), and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.
  • RAM Random Access Memory
  • non-volatile memory such as at least one disk memory.
  • the communication connection between the network element of the system and at least one other network element is implemented through at least one communication interface 1240 (which may be wired or wireless), which may use the Internet, a wide area network, a local network, a metropolitan area network, and the like.
  • the bus 1230 may be an ISA bus, a PCI bus, an EISA bus, or the like.
  • the bus can be divided into address bus, data bus, control bus and so on.
  • the memory 1220 is used to store a program, and the processor 1210 executes the program after receiving the execution instruction.
  • the imaging method based on the bionic pulse camera disclosed in any of the foregoing embodiments of the present application may be applied to the processor 1210, or Implemented by processor 1210 .
  • the processor 1210 may be an integrated circuit chip with signal processing capability.
  • each step of the above-mentioned method can be completed by an integrated logic circuit of hardware in the processor 1210 or an instruction in the form of software.
  • the above-mentioned processor 900 can be a general-purpose processor, including a central processing unit (Central Processing Unit, referred to as CPU), a network processor (Network Processor, referred to as NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), off-the-shelf programmable gate array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware component.
  • CPU Central Processing Unit
  • NP Network Processor
  • DSP digital signal processor
  • ASIC application-specific integrated circuit
  • FPGA off-the-shelf programmable gate array
  • the processor 1210 may implement or execute the methods, steps, and logical block diagrams disclosed in the embodiments of this application.
  • a general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
  • the steps of the methods disclosed in the embodiments of the present application may be directly executed by a hardware decoding processor, or executed by a combination of hardware and software modules in the decoding processor.
  • the software modules may be located in random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, registers and other storage media mature in the art.
  • the storage medium is located in the memory 1220, and the processor 1210 reads the information in the memory 1220, and completes the steps of the above method in combination with its hardware.
  • the electronic device provided by the embodiment of the present application and the imaging method provided by the embodiment of the present application are based on the same inventive concept, and have the same beneficial effects as the method adopted, operated or realized.
  • Embodiments of the present application further provide an imaging device, including the electronic device provided by the foregoing embodiments.
  • the imaging device may be a camera, a monitor, or the like.
  • Embodiments of the present application also provide a computer-readable storage medium corresponding to the imaging method provided by the foregoing embodiments.
  • the computer-readable storage medium shown is an optical disc 1300, on which a computer program ( That is, a program product), when the computer program is executed by the processor, the imaging method provided by any of the foregoing embodiments will be executed.
  • examples of computer-readable storage media may also include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory Memory (RAM), Read-Only Memory (ROM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Flash Memory or other optical and magnetic storage media will not be repeated here.
  • PRAM phase-change memory
  • SRAM static random access memory
  • DRAM dynamic random access memory
  • RAM random access memory
  • ROM Read-Only Memory
  • EEPROM Electrically Erasable Programmable Read-Only Memory
  • Flash Memory or other optical and magnetic storage media will not be repeated here.
  • the computer-readable storage medium provided by the above embodiments of the present application and the imaging methods provided by the embodiments of the present application are based on the same inventive concept, and have the same beneficial effects as the methods used, run or implemented by the application programs stored thereon.

Abstract

The present application discloses an imaging method and apparatus, device, and storage medium. The imaging method comprises: obtaining reconstructed image sequences according to a pulse array within a preset time period; according to relative motion relationships between the image sequences, performing filtering in the time dimension on each pixel point of image frames in the image sequences that meets a specified condition to update a light intensity value of each pixel point to obtain updated light intensity values; and generating updated images according to the updated light intensity value of each pixel point. The imaging method provided by the embodiments of the present application solves the problem of difficulty of traditional cameras in imaging a high-speed moving object clearly, and enables the generation of a clear image having a high signal-to-noise ratio.

Description

成像方法、装置、设备及存储介质Imaging method, apparatus, device and storage medium 技术领域technical field
本申请涉及图像处理技术领域,特别涉及一种成像方法、装置、设备及存储介质。The present application relates to the technical field of image processing, and in particular, to an imaging method, apparatus, device, and storage medium.
发明背景Background of the Invention
传统数字相机通常以固定帧率进行拍摄成像,每一帧图像按照以下方式生成:在一定的曝光时间窗内,图像传感器的每个像素对入射光进行光电转换和累计,曝光结束后经过模数转换得到该像素的光照总量。通常曝光时间的长度由光照强度决定,当光线较弱时往往需要增加曝光时间来抑制暗电流噪声的影响。这种方式无法对高速运动物体进行有效成像,往往导致高速运动物体的成像模糊。近年来,生物视网膜中央凹独特的神经元连接结构和神经节细胞的积分发放模型为视觉采样提供了新的思路。通过对视网膜中央凹的模拟和抽象,一种包含光感受器、积分器和阈值比较器的脉冲相机被提出。脉冲相机以脉冲阵列的形式表示视觉信息,能够持续记录光强的变化,不存在帧率和曝光时间的概念,突破了传统相机的局限性,既可以实现对高速运动过程的捕捉和记录,又可以重构出场景中纹理细节。Traditional digital cameras usually shoot and image at a fixed frame rate, and each frame of image is generated in the following way: within a certain exposure time window, each pixel of the image sensor performs photoelectric conversion and accumulation of incident light, and after exposure, the analog-to-digital Convert to get the total amount of light for this pixel. Usually the length of the exposure time is determined by the light intensity. When the light is weak, it is often necessary to increase the exposure time to suppress the influence of dark current noise. This method cannot effectively image high-speed moving objects, which often leads to blurred imaging of high-speed moving objects. In recent years, the unique neuronal connection structure of the biological fovea and the integral firing model of ganglion cells have provided new ideas for visual sampling. Through the simulation and abstraction of the fovea, a pulsed camera including a photoreceptor, an integrator, and a threshold comparator is proposed. Pulse cameras represent visual information in the form of pulse arrays and can continuously record changes in light intensity without the concept of frame rate and exposure time, breaking through the limitations of traditional cameras. The texture details in the scene can be reconstructed.
现有技术中,脉冲相机的图像重建算法包括两种:基于脉冲间隔的图像重建算法和基于时间窗平均的图像重建算法。其中,基于脉冲间隔的图像重建算法利用脉冲间隔随着光强的增加而减小这一特性,借助前后两个脉冲重构出一小段时间内的光强。该算法虽然能刻画出高速运动物体的轮廓,但重建信号通常不够稳定,在时间方向看像素值有明显的波动;基于时间窗平均的图像重建算法利用脉冲发放频率随着光强的增加而增加这一特性来估计时间窗内的平均光强。该算法虽然在一定程度上提高了重建图像的信噪比,但是当物体存在运动时,时间窗内的平均会导致重建图像存在运动模糊。In the prior art, the image reconstruction algorithm of the pulse camera includes two types: the image reconstruction algorithm based on the pulse interval and the image reconstruction algorithm based on the time window averaging. Among them, the image reconstruction algorithm based on the pulse interval utilizes the characteristic that the pulse interval decreases with the increase of the light intensity, and reconstructs the light intensity in a short period of time with the help of the two pulses before and after. Although this algorithm can describe the outline of high-speed moving objects, the reconstructed signal is usually not stable enough, and the pixel value has obvious fluctuations in the time direction; the image reconstruction algorithm based on the time window average uses the pulse firing frequency to increase with the increase of light intensity. This feature is used to estimate the average light intensity within a time window. Although the algorithm improves the signal-to-noise ratio of the reconstructed image to a certain extent, when there is motion of the object, the average in the time window will cause motion blur in the reconstructed image.
因此,现有技术中利用脉冲相机进行图像重建的方法,在高速运动场景下无法得到清晰的图像。Therefore, in the prior art method for image reconstruction using a pulse camera, a clear image cannot be obtained in a high-speed motion scene.
发明内容SUMMARY OF THE INVENTION
本申请实施例提供了一种成像方法、装置、设备及存储介质,可以提升重建图像的质量,得到高速运动场景下清晰的图像。为了对披露的实施例的一些方面有一个基本的理解,下面给出了简单的概括。该概括部分不是泛泛评述,也不是要确定关键/重要组成元素或描绘这些实施例的保护范围。其唯一目的是用简单的形式呈现一些概念,以此作为后面的详细说明的序言。Embodiments of the present application provide an imaging method, apparatus, device, and storage medium, which can improve the quality of reconstructed images and obtain clear images in high-speed motion scenes. In order to provide a basic understanding of some aspects of the disclosed embodiments, a brief summary is given below. This summary is not intended to be an extensive review, nor is it intended to identify key/critical elements or delineate the scope of protection of these embodiments. Its sole purpose is to present some concepts in a simplified form as a prelude to the detailed description that follows.
第一方面,本申请实施例提供了一种成像方法,包括:In the first aspect, the embodiments of the present application provide an imaging method, comprising:
根据预设时间段内的脉冲阵列,得到重构的图像序列;Obtain a reconstructed image sequence according to the pulse array within a preset time period;
根据图像序列之间的相对运动关系,对图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新各像素点的光强值,得到更新光强值;According to the relative motion relationship between the image sequences, filtering in the time dimension is performed on each pixel point of the image frame in the image sequence that meets the specified condition to update the light intensity value of each pixel point, and obtain the updated light intensity value;
根据各像素点的更新光强值生成更新图像。An updated image is generated according to the updated light intensity value of each pixel.
在一个实施例中,指定条件包括k时刻,k时刻位于预设时间段内。In one embodiment, the specified condition includes time k, which is within a preset time period.
在一个实施例中,根据图像序列之间的相对运动关系,对图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新各像素点的光强值,得到更新光强值,包括:In one embodiment, according to the relative motion relationship between the image sequences, filtering in the time dimension is performed on each pixel of the image frame in the image sequence that satisfies the specified condition to update the light intensity value of each pixel to obtain the updated light intensity values, including:
根据相对运动关系,确定各像素点在图像序列的t时刻的图像帧中的对应像素点;According to the relative motion relationship, determine the corresponding pixel points of each pixel point in the image frame at time t of the image sequence;
根据对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值。According to the light intensity value of the corresponding pixel point, each pixel point is filtered in the time dimension to obtain the updated light intensity value.
在一个实施例中,根据对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值,包括:In one embodiment, according to the light intensity value of the corresponding pixel point, filtering in the time dimension is performed on each pixel point to obtain the updated light intensity value, including:
根据对应像素点的光强值,通过如下公式,对各像素点进行时间维度上的滤波以得到更新光强值,According to the light intensity value of the corresponding pixel point, the following formula is used to filter each pixel point in the time dimension to obtain the updated light intensity value,
Figure PCTCN2022079318-appb-000001
Figure PCTCN2022079318-appb-000001
其中,I k(x,y)表示k时刻的图像帧内(x,y)像素点处的更新光强值,C表示正则化参数,
Figure PCTCN2022079318-appb-000002
表示t时刻的图像帧在重建更新图像时的滤波权重,
Figure PCTCN2022079318-appb-000003
表示k时刻的图像帧中的像素点在t时刻的图像帧中的对应像素点处的光强值,r表示滤波时间窗口半径。
Among them, I k (x, y) represents the updated light intensity value at the (x, y) pixel point in the image frame at time k, C represents the regularization parameter,
Figure PCTCN2022079318-appb-000002
Represents the filter weight of the image frame at time t when reconstructing the updated image,
Figure PCTCN2022079318-appb-000003
represents the light intensity value of the pixel in the image frame at time k at the corresponding pixel point in the image frame at time t, and r represents the radius of the filtering time window.
在一个实施例中,正则化参数C通过如下公式确定:In one embodiment, the regularization parameter C is determined by the following formula:
Figure PCTCN2022079318-appb-000004
Figure PCTCN2022079318-appb-000004
在一个实施例中,
Figure PCTCN2022079318-appb-000005
通过如下公式确定:
In one embodiment,
Figure PCTCN2022079318-appb-000005
Determined by the following formula:
Figure PCTCN2022079318-appb-000006
Figure PCTCN2022079318-appb-000006
其中,σ为设定参数。Among them, σ is the setting parameter.
在一个实施例中,设定参数通过如下公式确定:In one embodiment, the setting parameters are determined by the following formula:
σ=(2*r+1)/3。σ=(2*r+1)/3.
在一个实施例中,该成像方法还包括:In one embodiment, the imaging method further includes:
若对应像素点为整像素点,则直接取对应像素点处的光强值;If the corresponding pixel is an integer pixel, directly take the light intensity value at the corresponding pixel;
若对应像素点为亚像素点,则通过插值法得到对应像素点的光强值。If the corresponding pixel is a sub-pixel, the light intensity value of the corresponding pixel is obtained by interpolation.
在一个实施例中,根据图像序列之间的相对运动关系,对图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新各像素点的光强值,得到更新光强值,包括:In one embodiment, according to the relative motion relationship between the image sequences, filtering in the time dimension is performed on each pixel of the image frame in the image sequence that satisfies the specified condition to update the light intensity value of each pixel to obtain the updated light intensity values, including:
根据相对运动关系,对图像序列中像素点在时间方向上的自回归模型进行自适应学习,以确定自回归模型的模型参数;According to the relative motion relationship, the autoregressive model of the pixels in the image sequence in the time direction is adaptively learned to determine the model parameters of the autoregressive model;
确定各像素点在图像序列的t时刻的图像帧中的对应像素点;Determine the corresponding pixel points of each pixel point in the image frame at time t of the image sequence;
根据模型参数以及对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值。According to the model parameters and the light intensity value of the corresponding pixel point, each pixel point is filtered in the time dimension to obtain the updated light intensity value.
在一个实施例中,自回归模型的公式为:In one embodiment, the formula of the autoregressive model is:
Figure PCTCN2022079318-appb-000007
Figure PCTCN2022079318-appb-000007
其中,I k表示k时刻的理论光强值;I t表示t时刻的理论光强值;T k表示k时刻相邻时刻的集合;P k→t(x,y)为k时刻的图像帧内的像素点(x,y)在t时刻的图像帧内的对应像素位置;α t为模型参数;Δ k(x,y)为在(x,y)像素点上的模型误差。 Among them, I k represents the theoretical light intensity value at time k; I t represents the theoretical light intensity value at time t; T k represents the set of adjacent moments at time k; P k→t (x, y) is the image frame at time k The corresponding pixel position of the pixel point (x, y) in the image frame at time t; α t is the model parameter; Δ k (x, y) is the model error on the (x, y) pixel point.
在一个实施例中,根据相对运动关系,对图像序列中像素点在时间方向上的自回归模型进行自适应学习,以确定自回归模型的模型参数,包括:In one embodiment, according to the relative motion relationship, adaptive learning is performed on the autoregressive model of the pixel points in the image sequence in the time direction to determine the model parameters of the autoregressive model, including:
根据相对运动关系,通过最小二乘法求解目标函数的方式,对自回归模型进行自适应学习,以确定模型参数。According to the relative motion relationship, the autoregressive model is adaptively learned by the least square method to solve the objective function to determine the model parameters.
在一个实施例中,目标函数为:In one embodiment, the objective function is:
Figure PCTCN2022079318-appb-000008
Figure PCTCN2022079318-appb-000008
其中,W表示以(x,y)为中心的空间范围内像素点集合;
Figure PCTCN2022079318-appb-000009
表示k时刻的估计光强值;
Figure PCTCN2022079318-appb-000010
表示t时刻的估计光强值。
Among them, W represents the set of pixels in the spatial range centered on (x, y);
Figure PCTCN2022079318-appb-000009
Represents the estimated light intensity value at time k;
Figure PCTCN2022079318-appb-000010
represents the estimated light intensity value at time t.
在一个实施例中,根据模型参数以及对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值,包括:In one embodiment, according to the model parameters and the light intensity value of the corresponding pixel point, filtering in the time dimension is performed on each pixel point to obtain the updated light intensity value, including:
根据模型参数以及对应像素点的光强值,通过如下公式,对各像素点进行时间维度上的滤波以得到更新光强值,According to the model parameters and the light intensity value of the corresponding pixel point, the following formula is used to filter each pixel point in the time dimension to obtain the updated light intensity value,
Figure PCTCN2022079318-appb-000011
Figure PCTCN2022079318-appb-000011
其中,I′ k表示k时刻的更新光强值。 Among them, I′ k represents the updated light intensity value at time k.
在一个实施例中,该成像方法还包括:根据光流法、像素点匹配、像素点运动对齐和像素点相对位置偏移估计中的任一项或多项,得到图像序列之间的相对运动关系。In one embodiment, the imaging method further includes: obtaining the relative motion between the image sequences according to any one or more of optical flow method, pixel point matching, pixel point motion alignment, and pixel point relative position offset estimation relation.
在一个实施例中,在根据预设时间段内的脉冲阵列,得到重构的图像序列之前,该成像方法还包括:In one embodiment, before obtaining the reconstructed image sequence according to the pulse array within a preset time period, the imaging method further includes:
根据预设时间段内的脉冲信号,获取待成像物体基于像素点的脉冲阵列。According to the pulse signal within a preset time period, the pixel-based pulse array of the object to be imaged is obtained.
在一个实施例中,根据预设时间段内的脉冲阵列,得到重构的图像序列,包括:In one embodiment, a reconstructed image sequence is obtained according to the pulse array within a preset time period, including:
根据脉冲阵列,采用脉冲间隔算法重构出图像序列。According to the pulse array, the image sequence is reconstructed using the pulse interval algorithm.
第二方面,本申请实施例提供了一种成像装置,包括:In a second aspect, an embodiment of the present application provides an imaging device, including:
重构模块,用于根据预设时间段内的脉冲阵列,得到重构的图像序列;a reconstruction module, configured to obtain a reconstructed image sequence according to the pulse array within a preset time period;
更新模块,用于根据图像序列之间的相对运动关系,对图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新各像素点的光强值,得到更新光强值;The update module is used to filter each pixel point of the image frame that meets the specified condition in the image sequence in the time dimension according to the relative motion relationship between the image sequences to update the light intensity value of each pixel point, and obtain the updated light intensity value ;
成像模块,用于根据各像素点的更新光强值生成更新图像。The imaging module is configured to generate an updated image according to the updated light intensity value of each pixel point.
第三方面,本申请实施例提供了一种电子设备,包括处理器和存储有程序指令的存储器,处理器被配置为在执行程序指令时,执行上述实施例提供的成像方法。In a third aspect, embodiments of the present application provide an electronic device, including a processor and a memory storing program instructions, where the processor is configured to execute the imaging method provided by the foregoing embodiments when executing the program instructions.
第四方面,本申请实施例提供了一种成像设备,包括上述实施例提供的电子设备。In a fourth aspect, the embodiments of the present application provide an imaging device, including the electronic device provided by the above embodiments.
第五方面,本申请实施例提供了一种计算机可读介质,其上存储有计算机可读指令,计算机可读指令可被处理器执行以实现上述实施例提供的成像方法。In a fifth aspect, the embodiments of the present application provide a computer-readable medium on which computer-readable instructions are stored, and the computer-readable instructions can be executed by a processor to implement the imaging method provided by the foregoing embodiments.
本申请实施例提供的成像方法,根据预设时间段内的脉冲阵列,得到重构的图像序列;根据图像序列之间的相对运动关系,对图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新各像素点的光强值,得到更新光强值;根据各像素点的更新光强值生成更新图像。该方法提出了沿物体运动轨迹对图像像素进行时间维度滤波,能够应用于高速运动的场景,产生信噪比高、不模糊的图像,实现较好的视觉效果,解决了现有技术中无法对高速运动的物体进行清晰成像的问题。In the imaging method provided by the embodiment of the present application, a reconstructed image sequence is obtained according to a pulse array within a preset time period; according to the relative motion relationship between the image sequences, each pixel of an image frame that meets a specified condition in the image sequence is analyzed. Filtering in the time dimension is performed to update the light intensity value of each pixel point to obtain the updated light intensity value; an updated image is generated according to the updated light intensity value of each pixel point. The method proposes to filter the image pixels along the moving trajectory of the object in the temporal dimension, which can be applied to high-speed moving scenes to generate images with high signal-to-noise ratio and no blur, and achieve better visual effects. The problem of clear imaging of high-speed moving objects.
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本申请。It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not limiting of the present application.
附图简要说明Brief Description of Drawings
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本申请的实施例,并与说明书一起用于解释本申请的原理。The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the application and together with the description serve to explain the principles of the application.
图1是根据本申请一示例性实施例示出的成像方法的流程示意图。FIG. 1 is a schematic flowchart of an imaging method according to an exemplary embodiment of the present application.
图2是根据本申请一示例性实施例示出的脉冲阵列的示意图。FIG. 2 is a schematic diagram of a pulse array according to an exemplary embodiment of the present application.
图3是根据本申请一示例性实施例示出的运动对齐的示意图。FIG. 3 is a schematic diagram of motion alignment according to an exemplary embodiment of the present application.
图4是根据本申请一示例性实施例示出的根据插值法计算光强值的示意图。FIG. 4 is a schematic diagram of calculating a light intensity value according to an interpolation method according to an exemplary embodiment of the present application.
图5是根据本申请另一示例性实施例示出的成像方法的流程示意图。FIG. 5 is a schematic flowchart of an imaging method according to another exemplary embodiment of the present application.
图6是根据本申请另一示例性实施例示出的成像方法的流程示意图。FIG. 6 is a schematic flowchart of an imaging method according to another exemplary embodiment of the present application.
图7是根据本申请另一示例性实施例示出的成像方法的流程示意图。FIG. 7 is a schematic flowchart of an imaging method according to another exemplary embodiment of the present application.
图8是根据本申请另一示例性实施例示出的成像方法的流程示意图。FIG. 8 is a schematic flowchart of an imaging method according to another exemplary embodiment of the present application.
图9是根据本申请一示例性实施例示出的不同成像方法拍摄的成像图片的示意图。FIG. 9 is a schematic diagram of an imaging picture taken by different imaging methods according to an exemplary embodiment of the present application.
图10是根据本申请另一示例性实施例示出的不同成像方法拍摄的成像图片的示意图。FIG. 10 is a schematic diagram of an imaging picture taken by different imaging methods according to another exemplary embodiment of the present application.
图11是根据本申请一示例性实施例示出的成像装置的结构示意图。FIG. 11 is a schematic structural diagram of an imaging device according to an exemplary embodiment of the present application.
图12是根据本申请一示例性实施例示出的电子设备的结构示意图。FIG. 12 is a schematic structural diagram of an electronic device according to an exemplary embodiment of the present application.
图13是根据本申请一示例性实施例示出的计算机存储介质的示意图。FIG. 13 is a schematic diagram of a computer storage medium according to an exemplary embodiment of the present application.
实施本发明的方式MODES OF IMPLEMENTING THE INVENTION
为了能够更加详尽地了解本申请实施例的特点与技术内容,下面结合附图对本申请实施例的实现进行详细阐述,所附附图仅供参考说明之用,并非用来限定本申请实施例。在以下的技术描述中,为方便解释起见,通过多个细节以提供对所披露实施例的充分理解。然而,在没有这些细节的情况下,一个或一个以上实施例仍然可以实施。在其它情况下,为简化附图,熟知的结构和装置可以简化展示。显然,本申请所描述的实施例仅是本申请的一部分实施例,而不是所有实施例的穷举。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。In order to have a more detailed understanding of the features and technical contents of the embodiments of the present application, the implementation of the embodiments of the present application will be described in detail below with reference to the accompanying drawings. In the following technical description, for the convenience of explanation, numerous details are provided to provide a thorough understanding of the disclosed embodiments. However, one or more embodiments may be practiced without these details. In other instances, well-known structures and devices may be shown simplified in order to simplify the drawings. Obviously, the embodiments described in the present application are only a part of the embodiments of the present application, rather than an exhaustive list of all the embodiments. It should be noted that the embodiments in the present application and the features of the embodiments may be combined with each other in the case of no conflict.
在实现本申请的过程中,申请人发现传统相机无法对高速运动物体进行有效成像,往往导致高速运动物体的成像模糊。因此,受启发于生物视网膜的视觉采样机制,采集脉冲阵列信号的新型摄像机进入人们视野,包括基于光照强度变化发放脉冲信号的传感器,如动态视觉传感器(Dynamic Vision Sensor,DVS)、基于异步时间的图像传感器(Asynchronous Time-based image Sensor,ATIS)、动态有源像素视觉传感器(Dynamic and Active Pixel Vision Sensor,DAVIS)等,基于光照强度累积强度发放信号的传感器,如光强累积传感器等。这种摄像机的传感器采集一定时间,一定区域内光信号的信息,具有高动态范围,高时间分辨率等优点。In the process of realizing the present application, the applicant found that traditional cameras cannot effectively image high-speed moving objects, which often leads to blurred images of high-speed moving objects. Therefore, inspired by the visual sampling mechanism of the biological retina, new cameras that collect pulse array signals have entered people's field of vision, including sensors that emit pulse signals based on changes in light intensity, such as Dynamic Vision Sensor (DVS), asynchronous time-based sensors. Image sensor (Asynchronous Time-based image Sensor, ATIS), dynamic active pixel vision sensor (Dynamic and Active Pixel Vision Sensor, DAVIS), etc., sensors that emit signals based on the cumulative intensity of light intensity, such as light intensity cumulative sensors, etc. The sensor of this camera collects the information of the light signal in a certain area for a certain time, and has the advantages of high dynamic range and high time resolution.
而基于脉冲信号的脉冲相机的成像算法,虽然能刻画出高速运动物体的轮廓,但重建信号通常不够稳定和完善,在时间方向看像素值有明显的波动,或者会出现重建图像存在运动模糊的问题。因此,目前亟需一种新的成像方法来解决以上在高速运动场景下脉冲相机的成像不清晰问题。While the imaging algorithm of the pulse camera based on the pulse signal can describe the contour of the high-speed moving object, the reconstructed signal is usually not stable and perfect, and the pixel value in the time direction has obvious fluctuations, or the reconstructed image has motion blur. question. Therefore, a new imaging method is urgently needed to solve the above problem of unclear imaging of pulse cameras in high-speed motion scenes.
申请人还发现,脉冲阵列具有时间方向上的相关性,而这种相关性在图像的不同空间位置上的表现也不尽相同。比如,在运动比较连续的平滑运动区域,图像内容沿着运动轨迹的时间方向相关性比较强;而对于具有比较复杂的局部运动的区域,某一时刻的图像内容可能与前后时刻的图像内容不具有对应关系,甚至可能出现遮挡或消失,此时,图像内容沿着运动轨迹的时间方向相关性比较弱。Applicants have also found that pulse arrays have correlations in the temporal direction, and that correlations behave differently in different spatial locations of the image. For example, in a smooth motion area with relatively continuous motion, the temporal correlation of the image content along the motion trajectory is relatively strong; while for an area with relatively complex local motion, the image content at a certain moment may be different from the image content at the previous and subsequent times. There is a corresponding relationship, and even occlusion or disappearance may occur. At this time, the correlation of the time direction of the image content along the motion trajectory is relatively weak.
为了合理地利用脉冲阵列时间方向上的相关性以提升重建图像的质量,基于以上,本申请实施例提出了一种成像方法。In order to reasonably utilize the correlation in the time direction of the pulse array to improve the quality of the reconstructed image, based on the above, an embodiment of the present application proposes an imaging method.
如图1所示,该成像方法包括如下内容。As shown in FIG. 1, the imaging method includes the following contents.
S110:根据预设时间段内的脉冲阵列,得到重构的图像序列。S110: Obtain a reconstructed image sequence according to the pulse array within a preset time period.
在一种可能的实现方式中,在执行步骤S110之前,该成像方法还包括:根据预设时间段内的脉冲信号,获取待成像物体基于像素点的脉冲阵列。In a possible implementation manner, before step S110 is performed, the imaging method further includes: acquiring a pixel-based pulse array of the object to be imaged according to the pulse signal within a preset time period.
具体地,为了打破传统相机的局限,本申请实施例可采用仿生式脉冲相机进行拍摄成像,仿生式脉冲相机中每个像素均包含光感受器、积分器和阈值比较器。其中,光感受器负责将光信号转换为电信号,积分器对转换后的光生电荷进行累积,阈值比较器以超高的频率反复地对累积的电荷进行比较,当累积后的电荷达到预设阈值时,该像素处进行脉冲发放,并清空积分器。图2是一种脉冲阵列的示意图,如图2所示,脉冲相机以“H*W*T”脉冲阵列的形式表示视觉信息,“H*W”为脉冲相机的空间分辨率,“T”为采样次数。我们称作单个像素输出的信号序列为“脉冲序列”,脉冲阵列在某个时刻的截面为“脉冲矩阵”。脉冲阵列由“0”和“1”两种符号组成,其中“1”(实心点)表示在该时空位置处有脉冲发放;“0”(空心点)表示该时空位置无脉冲发放。通过该仿生式脉冲相机,能够持续记录光强的变化,不存在帧率和曝光时间的概念,突破了传统相机的局限性。Specifically, in order to break the limitations of traditional cameras, the embodiments of the present application may use a bionic pulse camera for imaging, and each pixel in the bionic pulse camera includes a photoreceptor, an integrator, and a threshold comparator. Among them, the photoreceptor is responsible for converting the optical signal into an electrical signal, the integrator accumulates the converted photo-generated charges, and the threshold comparator repeatedly compares the accumulated charges at an ultra-high frequency. When the accumulated charges reach a preset threshold , the pixel is pulsed and the integrator is cleared. Figure 2 is a schematic diagram of a pulse array. As shown in Figure 2, the pulse camera represents visual information in the form of "H*W*T" pulse array, "H*W" is the spatial resolution of the pulse camera, and "T" is the sampling times. We call the signal sequence output by a single pixel as a "pulse sequence", and the cross-section of the pulse array at a certain moment is a "pulse matrix". The pulse array consists of two symbols, "0" and "1", where "1" (solid dot) indicates that there is a pulse firing at this spatiotemporal position; "0" (open dot) means that there is no pulse firing at this spatiotemporal position. The bionic pulse camera can continuously record changes in light intensity without the concept of frame rate and exposure time, breaking through the limitations of traditional cameras.
进一步地,获取到待成像物体基于像素点的脉冲阵列之后,根据获取到的预设时间段内的脉冲阵列,重构出基本的图像序列。在一种可能的实现方式中,根据脉冲阵列,可以采用脉冲间隔算法重构出图像序列,也可以采用其他算法根据脉冲阵列估计出基本的重建图像序列。Further, after acquiring the pixel-based pulse array of the object to be imaged, a basic image sequence is reconstructed according to the acquired pulse array within a preset time period. In a possible implementation manner, according to the pulse array, a pulse interval algorithm can be used to reconstruct the image sequence, or other algorithms can be used to estimate the basic reconstructed image sequence according to the pulse array.
在一个实施例中,采用脉冲间隔算法重构出图像序列,脉冲间隔为同一像素位置两个相邻脉冲所间隔的时间,等价于积分器累积得到后一个脉冲所用的积分时间。通常而言,强光照区 域像素的积分器可以在较短的时间内将光生电荷累积到脉冲发放阈值,脉冲间隔较小;弱光照区域像素的积分器则要花费更长的时间将电荷累积到脉冲发放阈值,脉冲间隔较大。根据脉冲间的平均光强与脉冲间隔的大小成反比这一特性,重构出基本图像序列
Figure PCTCN2022079318-appb-000012
In one embodiment, the image sequence is reconstructed using a pulse interval algorithm, and the pulse interval is the time interval between two adjacent pulses at the same pixel position, which is equivalent to the integration time used by the integrator to accumulate the next pulse. Generally speaking, the integrator of the pixel in the strong illumination area can accumulate the photogenerated charge to the pulse emission threshold in a shorter time, and the pulse interval is small; the integrator of the pixel in the weak illumination area takes a longer time to accumulate the charge to Threshold of pulse release, pulse interval is larger. According to the characteristic that the average light intensity between pulses is inversely proportional to the size of the pulse interval, the basic image sequence is reconstructed
Figure PCTCN2022079318-appb-000012
根据该步骤,可以根据脉冲阵列,采用脉冲间隔算法重构出基本的图像序列。According to this step, the basic image sequence can be reconstructed according to the pulse array using the pulse interval algorithm.
S120:根据图像序列之间的相对运动关系,对图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新各像素点的光强值,得到更新光强值。S120: According to the relative motion relationship between the image sequences, perform filtering in the time dimension on each pixel of the image frame in the image sequence that meets the specified condition to update the light intensity value of each pixel to obtain an updated light intensity value.
具体地,指定条件可以是某个时刻,或者是具有某特征的图像帧。例如,指定条件包括k时刻,k时刻位于预设时间段内。Specifically, the specified condition may be a certain moment, or an image frame with a certain feature. For example, the specified condition includes time k, which is within a preset time period.
根据图像序列之间的相对运动关系可以将图像序列进行运动对齐,得到图像序列中不同图像帧内各像素点间的对应关系。当相机位置固定时,随着物体的运动,物体中的某一点在不同时刻的图像帧中的像素点位置不同。因此,针对k时刻的图像帧内的某一像素点,可以获得其他时刻图像帧内的对应像素点。根据图像序列中不同图像帧内各像素点间的对应关系,可以获取不同图像帧内各像素点的运动轨迹。According to the relative motion relationship between the image sequences, the motion alignment of the image sequences can be performed to obtain the corresponding relationship between pixels in different image frames in the image sequence. When the camera position is fixed, with the movement of the object, a certain point in the object has different pixel positions in the image frames at different times. Therefore, for a certain pixel in the image frame at time k, corresponding pixels in the image frame at other times can be obtained. According to the correspondence between each pixel in different image frames in the image sequence, the motion trajectory of each pixel in different image frames can be obtained.
基于像素点的运动轨迹,可以确定像素点在不同时刻的光强值。例如,针对k时刻的图像帧内的每个像素点,可以基于该像素点在不同时刻的图像帧中对应的光强值对该像素点进行时间维度上的滤波,以更新该像素点的光强值,得到更新光强值。根据k时刻的图像帧内的各个像素点的更新光强值可以生成更新图像。Based on the motion trajectory of the pixel point, the light intensity value of the pixel point at different times can be determined. For example, for each pixel in the image frame at time k, filtering in the time dimension can be performed on the pixel based on the corresponding light intensity values of the pixel in the image frames at different times to update the light intensity of the pixel. Intensity value, get updated light intensity value. An updated image can be generated according to the updated light intensity value of each pixel in the image frame at time k.
S130:根据各像素点的更新光强值生成更新图像。S130: Generate an updated image according to the updated light intensity value of each pixel.
本申请实施例提供的成像方法,利用脉冲阵列的时间方向相关性,提升重建图像的质量。The imaging method provided by the embodiment of the present application utilizes the temporal direction correlation of the pulse array to improve the quality of the reconstructed image.
本申请实施例中,通过根据图像序列之间的相对运动关系,对图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新各像素点的光强值,得到更新光强值,进而根据各像素点的更新光强值生成更新图像,从而可以提升重建图像的质量,得到高速运动场景下清晰的图像,实现较好的视觉效果,解决了现有技术中无法对高速运动的物体进行清晰成像的问题。In the embodiment of the present application, according to the relative motion relationship between the image sequences, filtering in the time dimension is performed on each pixel point of the image frame in the image sequence that meets the specified condition to update the light intensity value of each pixel point, so as to obtain the updated light Intensity value, and then generate an updated image according to the updated light intensity value of each pixel, so that the quality of the reconstructed image can be improved, a clear image in a high-speed motion scene can be obtained, and a better visual effect can be achieved. The problem of clear imaging of moving objects.
根据本申请一实施例,步骤S120包括:根据相对运动关系,确定各像素点在图像序列的t时刻的图像帧中的对应像素点;根据对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值。According to an embodiment of the present application, step S120 includes: determining the corresponding pixel points of each pixel point in the image frame at time t of the image sequence according to the relative motion relationship; Dimensional filtering to get updated light intensity values.
具体地,针对k时刻的图像帧内的每个像素点,可以基于相对运动关系确定该像素点在其他时刻图像帧中的对应像素点。对应像素点的光强值可以认为是k时刻的图像帧内的该像素点在其他时刻图像帧中对应的光强值。Specifically, for each pixel in the image frame at time k, the corresponding pixel of the pixel in the image frame at other times may be determined based on the relative motion relationship. The light intensity value of the corresponding pixel point can be considered as the light intensity value corresponding to the pixel point in the image frame at time k at other time points in the image frame.
根据本申请一实施例,该成像方法还包括:根据光流法、像素点匹配、像素点运动对齐和像素点相对位置偏移估计中的任一项或多项,得到图像序列之间的相对运动关系。According to an embodiment of the present application, the imaging method further includes: obtaining the relative relationship between the image sequences according to any one or more of optical flow method, pixel point matching, pixel point motion alignment, and pixel point relative position offset estimation. sports relationship.
例如,根据光流法得到图像序列之间的相对运动关系。For example, the relative motion relationship between image sequences is obtained according to the optical flow method.
在一种可能的实现方式中,利用Horn–Schunck等光流算法计算当前帧
Figure PCTCN2022079318-appb-000013
与其前后r帧图像间的光流场,从而确定图像间的相对运动关系。
In a possible implementation, the current frame is calculated using an optical flow algorithm such as Horn–Schunck
Figure PCTCN2022079318-appb-000013
The optical flow field between the r frame images before and after it is used to determine the relative motion relationship between the images.
其中,计算光流场是为了找到图像序列间的对应关系。光流场可以理解为就是运动场,图像亮度模式的表观运动就是光流。Among them, the calculation of the optical flow field is to find the correspondence between the image sequences. The optical flow field can be understood as a motion field, and the apparent motion of the image brightness pattern is the optical flow.
在一个示例性场景中,假设所拍摄场景中有一个三角形正在自由下落,那么保持相机不动的话,它在不同时刻所拍摄的图像中的位置不同,利用光流法来找到同一物体在不同图像上的对应关系,光流法大多基于光强一致性原理,即假设相同物体光强值一样。例如,第一张图的(x,y)点对应于第二张图的(x’,y’)点,UV_t^k(x,y)就是(x,y)点的光流(描述两帧间该点处的相对运动)。In an exemplary scene, assuming that a triangle is in free fall in the captured scene, then if the camera is kept still, its position in the images captured at different times is different, and the optical flow method is used to find the same object in different images. The optical flow method is mostly based on the principle of light intensity consistency, that is, it is assumed that the same object has the same light intensity value. For example, the (x, y) point of the first image corresponds to the (x', y') point of the second image, and UV_t^k(x, y) is the optical flow at the (x, y) point (describes the two relative motion at that point between frames).
根据上述步骤,可以基于光强一致性假设,利用光流法求解两帧之间的光流(相对运动),进一步地,基于光流(相对运动),可以确定不同帧间图像的对应关系,即运动对齐,得到当前帧内每个像素点(x,y)在其他
Figure PCTCN2022079318-appb-000014
帧上的对应位置
Figure PCTCN2022079318-appb-000015
According to the above steps, the optical flow method can be used to solve the optical flow (relative motion) between two frames based on the assumption of light intensity consistency, and further, based on the optical flow (relative motion), the correspondence between images between different frames can be determined, That is, motion alignment, get each pixel (x, y) in the current frame in other
Figure PCTCN2022079318-appb-000014
the corresponding position on the frame
Figure PCTCN2022079318-appb-000015
图3是运动对齐的示意图,如图3所示,自由下落的三角形在不同图像上的位置不同,基于光流信息我们可以确定图像间的对应关系,找到场景中同一物体点在不同时刻点图像上的位置,确定物体的运动轨迹。Figure 3 is a schematic diagram of motion alignment. As shown in Figure 3, the positions of free-falling triangles on different images are different. Based on the optical flow information, we can determine the correspondence between the images and find the same object point in the scene at different time points. The position on the object determines the trajectory of the object.
在一个实施例中,得到不同帧各像素点间的对应关系之后(根据相对运动关系,确定各像素点在图像序列的t时刻的图像帧中的对应像素点之后),该成像方法还包括:获取当前帧的像素点在其他帧图像上的对应像素点处的光强值(获取对应像素点的光强值)。In one embodiment, after obtaining the corresponding relationship between the pixels in different frames (according to the relative motion relationship, after determining the corresponding pixels of each pixel in the image frame at time t of the image sequence), the imaging method further includes: Obtain the light intensity value of the pixel point of the current frame at the corresponding pixel point on other frame images (obtain the light intensity value of the corresponding pixel point).
在一个实施例中,获取对应像素点的光强值包括:若对应像素点为整像素点,则直接取对应像素点处的光强值;若对应像素点为亚像素点,则通过插值法得到对应像素点的光强值。In one embodiment, acquiring the light intensity value of the corresponding pixel point includes: if the corresponding pixel point is an integer pixel point, directly obtaining the light intensity value at the corresponding pixel point; if the corresponding pixel point is a sub-pixel point, using an interpolation method Get the light intensity value of the corresponding pixel point.
具体地,若当前帧的像素点在其他帧图像上的对应像素点为整像素点,则直接取对应像素点处的光强值;若当前帧的像素点在其他帧图像上的对应像素点为亚像素点,则通过插值法得到对应像素点处的光强值。Specifically, if the pixel of the current frame corresponds to the pixel in other frame images, the light intensity value at the corresponding pixel is directly taken; if the pixel of the current frame corresponds to the pixel in other frame images is a sub-pixel point, then the light intensity value at the corresponding pixel point is obtained by interpolation.
图4是一种利用插值法计算光强值的示意图,如图4所示,假设运动对齐后,当前点在其他帧图像上的对应点未处于整数位置,那么需利用其周围整点上的光强值来估计出该点处的光强值。以双线性插值举例,P为非整数点,该点的光强值未知,那么可以利用其周围整点Q 11,Q 12,Q 21,Q 22上的亮度信息,来插值出P点的光强值。具体做法如下: Figure 4 is a schematic diagram of calculating the light intensity value using the interpolation method. As shown in Figure 4, assuming that after the motion is aligned, the corresponding point of the current point on other frame images is not at an integer position, then it is necessary to use the surrounding whole point. light intensity value to estimate the light intensity value at that point. Taking bilinear interpolation as an example, P is a non-integer point, and the light intensity value of this point is unknown, then the brightness information of the surrounding integer points Q 11 , Q 12 , Q 21 , and Q 22 can be used to interpolate the light intensity of point P. light intensity value. The specific methods are as follows:
Figure PCTCN2022079318-appb-000016
where R 1=(x,y 1)
Figure PCTCN2022079318-appb-000016
where R 1 =(x,y 1 )
Figure PCTCN2022079318-appb-000017
where R 2=(x,y 2)
Figure PCTCN2022079318-appb-000017
where R 2 =(x,y 2 )
Figure PCTCN2022079318-appb-000018
Figure PCTCN2022079318-appb-000018
其中,f(.)表示某点的光强值。Among them, f(.) represents the light intensity value of a certain point.
根据本申请一实施例,根据对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值,包括:根据对应像素点的光强值,通过如下公式,对各像素点进行时间维度上的滤波以得到更新光强值,According to an embodiment of the present application, according to the light intensity value of the corresponding pixel point, filtering each pixel point in the time dimension to obtain the updated light intensity value includes: according to the light intensity value of the corresponding pixel point, through the following formula, for each pixel point The pixel points are filtered in the time dimension to obtain the updated light intensity value,
Figure PCTCN2022079318-appb-000019
Figure PCTCN2022079318-appb-000019
其中,I k(x,y)表示k时刻的图像帧内(x,y)像素点处的更新光强值,C表示正则化参数,
Figure PCTCN2022079318-appb-000020
表示t时刻的图像帧在重建更新图像(重建k时刻的图像帧)时的滤波权重,
Figure PCTCN2022079318-appb-000021
表示k时刻的图像帧中的像素点在t时刻的图像帧中的对应像素点处的光强值,r表示滤波时间窗口半径。
Among them, I k (x, y) represents the updated light intensity value at the (x, y) pixel point in the image frame at time k, C represents the regularization parameter,
Figure PCTCN2022079318-appb-000020
Represents the filter weight of the image frame at time t when reconstructing the updated image (reconstructing the image frame at time k),
Figure PCTCN2022079318-appb-000021
represents the light intensity value of the pixel in the image frame at time k at the corresponding pixel point in the image frame at time t, and r represents the radius of the filtering time window.
权重
Figure PCTCN2022079318-appb-000022
的选取随着两帧间时间间隔的增大而逐渐减小,可以采用下述表达式或其他类似表达式计算得出。
Weights
Figure PCTCN2022079318-appb-000022
The selection of is gradually reduced with the increase of the time interval between two frames, which can be calculated by the following expression or other similar expressions.
Figure PCTCN2022079318-appb-000023
Figure PCTCN2022079318-appb-000023
其中,σ可为设定参数,其取值由时间窗口大小决定,例如σ=(2*r+1)/3,r为滤波时间窗口半径。Wherein, σ can be a set parameter, and its value is determined by the size of the time window, for example, σ=(2*r+1)/3, and r is the radius of the filter time window.
根据本申请一实施例,
Figure PCTCN2022079318-appb-000024
According to an embodiment of the present application,
Figure PCTCN2022079318-appb-000024
可选地,根据对应像素点的光强值,可通过如下公式,对各像素点进行时间维度上的滤波以得到更新光强值,Optionally, according to the light intensity value of the corresponding pixel point, the following formula can be used to filter each pixel point in the time dimension to obtain the updated light intensity value,
Figure PCTCN2022079318-appb-000025
Figure PCTCN2022079318-appb-000025
其中,T表示图像序列的总帧数,
Figure PCTCN2022079318-appb-000026
1≤k≤T。
Figure PCTCN2022079318-appb-000027
可以根据上述实施例提供的方法确定,或根据实际情况通过其他方法确定。
where T represents the total number of frames in the image sequence,
Figure PCTCN2022079318-appb-000026
1≤k≤T.
Figure PCTCN2022079318-appb-000027
It can be determined according to the method provided in the foregoing embodiment, or determined by other methods according to the actual situation.
在本申请实施例中,根据不同帧各像素点间的对应关系,可得到物体的运动轨迹,沿着物体的运动轨迹对信号进行时间维度的滤波,可更新k时刻的图像帧内(x,y)像素点上的光强值。In the embodiment of the present application, according to the corresponding relationship between the pixels in different frames, the motion trajectory of the object can be obtained, and the signal is filtered along the time dimension along the motion trajectory of the object, and the image frame at time k can be updated (x, y) The light intensity value on the pixel point.
本实施例可以对基本重建图像序列进行时间维度的运动轨迹滤波以提升信号的稳定性,从而提升重建图像的质量。其中,沿着时间维度的运动轨迹滤波实现了对脉冲阵列的时间相关性的利用,运动对齐保证了针对运动物体而言其所利用的时间方向相关性的准确性。In this embodiment, motion trajectory filtering in the time dimension can be performed on the basic reconstructed image sequence to improve the stability of the signal, thereby improving the quality of the reconstructed image. Among them, the motion trajectory filtering along the time dimension realizes the utilization of the time correlation of the pulse array, and the motion alignment ensures the accuracy of the time direction correlation used for the moving object.
图5为本申请另一实施例提供的成像方法,具体为基于仿生式脉冲相机的成像方法,图5是图1实施例的例子,相同之处不再赘述,此处着重描述不同之处。如图5所示,该方法具体包括以下内容。FIG. 5 is an imaging method provided by another embodiment of the present application, specifically an imaging method based on a bionic pulse camera. FIG. 5 is an example of the embodiment of FIG. 1 . As shown in FIG. 5 , the method specifically includes the following contents.
S510:根据预设时间段内的脉冲阵列,得到重构的图像序列。S510: Obtain a reconstructed image sequence according to the pulse array within a preset time period.
S520:根据光流法得到图像序列之间的相对运动关系。S520: Obtain the relative motion relationship between the image sequences according to the optical flow method.
S530:根据相对运动关系将图像序列进行运动对齐,得到不同图像帧各像素点间的对应关系。S530: Perform motion alignment on the image sequence according to the relative motion relationship to obtain the corresponding relationship between pixels in different image frames.
S540:根据不同图像帧各像素点间的对应关系,对图像序列的像素点进行时间维度上的滤波,更新待重构图像的像素点的光强值以生成更新图像。S540: According to the correspondence between the pixels of different image frames, perform filtering on the pixels of the image sequence in the time dimension, and update the light intensity values of the pixels of the image to be reconstructed to generate an updated image.
具体地,根据本申请实施例提供的成像方法,可以重构图像序列中的一帧图像或多帧图像。Specifically, according to the imaging method provided by the embodiments of the present application, one frame of images or multiple frames of images in the image sequence can be reconstructed.
为了便于理解本申请实施例提供的基于仿生式脉冲相机的成像方法,下面结合附图6进行说明。如图6所示,该方法包括如下内容。In order to facilitate the understanding of the imaging method based on the bionic pulse camera provided by the embodiment of the present application, the following description is made with reference to FIG. 6 . As shown in Figure 6, the method includes the following contents.
首先,采用仿生式脉冲相机进行拍摄成像,根据预设时间段内的脉冲信号,得到待成像物体基于像素点的脉冲阵列。然后根据脉冲阵列,采用脉冲间隔算法得到重构的图像序列(基础重建)。First, a bionic pulse camera is used for imaging, and a pixel-based pulse array of the object to be imaged is obtained according to the pulse signal within a preset time period. Then, based on the pulse array, a pulse interval algorithm is used to obtain a reconstructed image sequence (basic reconstruction).
进一步地,利用Horn–Schunck等光流估计算法计算当前帧与其前后帧图像间的光流场,从而确定图像间的相对运动关系。Further, the optical flow estimation algorithm such as Horn-Schunck is used to calculate the optical flow field between the current frame and its previous and previous frame images, so as to determine the relative motion relationship between the images.
根据上述步骤,可以基于光强一致性假设,利用光流法求解两帧之间的光流(相对运动),进一步地,基于光流(相对运动),可以确定不同帧间图像的对应关系,即运动对齐。得到当前帧内每个像素点(x,y)在其他
Figure PCTCN2022079318-appb-000028
帧上的对应位置
Figure PCTCN2022079318-appb-000029
According to the above steps, the optical flow method can be used to solve the optical flow (relative motion) between two frames based on the assumption of light intensity consistency, and further, based on the optical flow (relative motion), the correspondence between images between different frames can be determined, i.e. motion alignment. Get each pixel (x, y) in the current frame in other
Figure PCTCN2022079318-appb-000028
the corresponding position on the frame
Figure PCTCN2022079318-appb-000029
得到不同帧各像素点间的对应关系之后,该方法还包括:获取当前帧的像素点在其他帧图像上的对应像素点处的光强值。After obtaining the correspondence between the pixels of different frames, the method further includes: acquiring the light intensity values of the pixels of the current frame at the corresponding pixels of other frame images.
最后,根据不同帧各像素点间的对应关系,可得到物体的运动轨迹,沿着物体的运动轨迹对信号进行时间维度的滤波,更新像素点的光强值以生成更新图像。Finally, according to the correspondence between the pixels in different frames, the motion trajectory of the object can be obtained, and the signal is filtered in the time dimension along the motion trajectory of the object, and the light intensity value of the pixel is updated to generate an updated image.
本申请实施例中,首先根据脉冲阵列的特性对不同时刻的场景进行图像序列重建;而后在重建图像序列的基础上采用光流估计算法确定物体沿着时间方向的运动轨迹,沿着物体各自的运动轨迹对图像序列进行运动对齐;最后对重建图像序列进行时间维度的运动轨迹滤波以提升信号的稳定性,从而提升重建图像的质量,得到高速运动场景下清晰的图像,实现较好的视觉效果,解决了现有技术中无法对高速运动的物体进行清晰成像的问题。本申请实施例提供的基于仿生式脉冲相机的成像方法,能够应用于高速运动的场景,可以产生信噪比高、不模糊的图像。In the embodiment of the present application, firstly, the image sequence reconstruction is performed on the scene at different times according to the characteristics of the pulse array; then, on the basis of the reconstructed image sequence, the optical flow estimation algorithm is used to determine the motion trajectory of the object along the time direction, along the respective The motion track performs motion alignment on the image sequence; finally, the reconstructed image sequence is subjected to time dimension motion track filtering to improve the stability of the signal, thereby improving the quality of the reconstructed image, obtaining clear images in high-speed motion scenes, and achieving better visual effects. , which solves the problem that high-speed moving objects cannot be clearly imaged in the prior art. The imaging method based on the bionic pulse camera provided in the embodiment of the present application can be applied to a high-speed motion scene, and can generate an image with a high signal-to-noise ratio and no blur.
根据本申请一实施例,步骤S120包括:根据相对运动关系,对图像序列中像素点在时间方向上的自回归模型进行自适应学习,以确定自回归模型的模型参数;确定各像素点在图像序列的t时刻的图像帧中的对应像素点;根据模型参数以及对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值。According to an embodiment of the present application, step S120 includes: performing adaptive learning on the autoregressive model of the pixels in the image sequence in the time direction according to the relative motion relationship to determine the model parameters of the autoregressive model; The corresponding pixel points in the image frame at time t of the sequence; according to the model parameters and the light intensity value of the corresponding pixel point, each pixel point is filtered in the time dimension to obtain the updated light intensity value.
具体地,自回归模型可以根据需要提前设置。Specifically, the autoregressive model can be set in advance as required.
在一实施例中,自回归模型的公式为:In one embodiment, the formula of the autoregressive model is:
Figure PCTCN2022079318-appb-000030
Figure PCTCN2022079318-appb-000030
其中,I k表示k时刻的理论光强值(图像理论值);I t表示t时刻的理论光强值(图像理论值);T k表示k时刻相邻时刻的集合;P k→t(x,y)为k时刻的图像帧内的像素点(x,y)在t时刻的图像帧内的对应像素位置;α t为模型参数;Δ k(x,y)为在(x,y)像素点上的模型误差,可能对应一些图像细节或噪声。 Among them, I k represents the theoretical light intensity value (image theoretical value) at time k; I t represents the theoretical light intensity value (image theoretical value) at time t; T k represents the set of adjacent moments at time k; P k→t ( x, y) is the corresponding pixel position of the pixel (x, y) in the image frame at time k in the image frame at time t; α t is the model parameter; Δ k (x, y) is the pixel at (x, y) ) model error on pixels, which may correspond to some image details or noise.
根据本申请一实施例,根据相对运动关系,对图像序列中像素点在时间方向上的自回归模型进行自适应学习,以确定自回归模型的模型参数,包括:根据相对运动关系,通过最小二乘法求解目标函数的方式,对自回归模型进行自适应学习,以确定模型参数。According to an embodiment of the present application, according to the relative motion relationship, adaptive learning is performed on the autoregressive model of the pixel points in the image sequence in the time direction, so as to determine the model parameters of the autoregressive model, including: according to the relative motion relationship, through the least two The method of multiplying to solve the objective function performs adaptive learning on the autoregressive model to determine the model parameters.
具体地,可自适应选择图像序列上的像素点,对自回归模型进行自适应学习,确定自回归模型的模型参数。例如,通过最小二乘法求解目标函数,得到自回归模型的模型参数,目标函数可以为:Specifically, the pixels on the image sequence can be adaptively selected, the autoregressive model can be adaptively learned, and the model parameters of the autoregressive model can be determined. For example, the objective function is solved by the least squares method to obtain the model parameters of the autoregressive model. The objective function can be:
Figure PCTCN2022079318-appb-000031
Figure PCTCN2022079318-appb-000031
其中,W表示以(x,y)为中心的空间范围内像素点集合;
Figure PCTCN2022079318-appb-000032
表示k时刻的估计光强值(图像估计值);
Figure PCTCN2022079318-appb-000033
表示t时刻的估计光强值(图像估计值)。
Among them, W represents the set of pixels in the spatial range centered on (x, y);
Figure PCTCN2022079318-appb-000032
Represents the estimated light intensity value at time k (image estimated value);
Figure PCTCN2022079318-appb-000033
Indicates the estimated light intensity value (image estimated value) at time t.
在这一步骤中,不同图像区域中,图像内容沿着运动轨迹的时间方向相关性不同,因此,需要根据图像序列局部区域的内容对自回归模型的参数进行自适应地调整,从而确定适合局部区域的滤波权重。在对脉冲阵列中的像素进行时间方向上的滤波时,滤波权重与脉冲阵列的内容相关,一帧图像不同空间位置上的滤波权重可能存在差异。In this step, in different image regions, the temporal correlation of the image content along the motion trajectory is different. Therefore, it is necessary to adjust the parameters of the autoregressive model adaptively according to the content of the local region of the image sequence, so as to determine the suitable local Filter weights for regions. When filtering the pixels in the pulse array in the temporal direction, the filter weight is related to the content of the pulse array, and the filter weight at different spatial positions of an image may be different.
因此,本申请假设空间位置相近的像素点具有相对一致的时间相关性,故模型参数可以根据局部空间位置上的数据自适应地选择,可以通过优化以上目标函数得出。Therefore, the present application assumes that pixels with similar spatial positions have relatively consistent temporal correlations, so the model parameters can be adaptively selected according to the data at the local spatial positions, and can be obtained by optimizing the above objective function.
根据本申请一实施例,根据模型参数以及对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值,包括:根据模型参数以及对应像素点的光强值,通过如下公式,对各像素点进行时间维度上的滤波以得到更新光强值,According to an embodiment of the present application, according to the model parameters and the light intensity value of the corresponding pixel point, filtering each pixel point in the time dimension to obtain the updated light intensity value, including: according to the model parameter and the light intensity value of the corresponding pixel point, Through the following formula, each pixel is filtered in the time dimension to obtain the updated light intensity value,
Figure PCTCN2022079318-appb-000034
Figure PCTCN2022079318-appb-000034
其中,I′ k表示k时刻的更新光强值,或者称为图像自适应成像值。具体地,在得到模型参数的基础上,利用该公式对像素点进行时间维度上的滤波,得到的结果作为k时刻的图像帧内(x,y)像素点上的更新光强值,最终可得到物体的自适应成像。 Among them, I′ k represents the updated light intensity value at time k, or is called image adaptive imaging value. Specifically, on the basis of obtaining the model parameters, this formula is used to filter the pixel points in the time dimension, and the obtained result is used as the updated light intensity value on the (x, y) pixel point in the image frame at time k, which can be finally obtained. Adaptive imaging of the object is obtained.
在其他实施例中,可以根据图像序列的各像素点的运动轨迹,建立各像素点在时间方向上的自回归模型。例如,可根据图像序列之间的相对运动关系,确定图像序列的各像素点的运动轨迹,运动轨迹的确定可参见上述实施例中的描述。然后自适应选择图像序列上的像素点,对自回归模型进行自适应学习,以确定模型参数。In other embodiments, an autoregressive model of each pixel in the time direction may be established according to the motion trajectory of each pixel of the image sequence. For example, the motion trajectory of each pixel of the image sequence may be determined according to the relative motion relationship between the image sequences, and reference may be made to the description in the foregoing embodiments for the determination of the motion trajectory. Then the pixels on the image sequence are adaptively selected, and the autoregressive model is adaptively learned to determine the model parameters.
图7为本申请另一实施例提供的成像方法,具体为基于脉冲信号的自适应成像方法,图7是图1实施例的例子,相同之处不再赘述,此处着重描述不同之处。如图7所示,该方法具体包括以下内容。FIG. 7 is an imaging method provided by another embodiment of the present application, specifically an adaptive imaging method based on a pulse signal. FIG. 7 is an example of the embodiment of FIG. 1 . As shown in FIG. 7 , the method specifically includes the following contents.
S710:根据预设时间段内的脉冲阵列,得到重构的图像序列。S710: Obtain a reconstructed image sequence according to the pulse array within a preset time period.
S720:根据图像序列之间的相对运动关系,确定图像序列各像素点的运动轨迹。S720: Determine the motion trajectory of each pixel of the image sequence according to the relative motion relationship between the image sequences.
S730:根据图像序列各像素点的运动轨迹,建立各像素点在时间方向上的自回归模型。S730: Establish an autoregressive model of each pixel in the time direction according to the motion trajectory of each pixel in the image sequence.
具体地,可以根据各像素点的运动轨迹,建立合适的自回归模型。Specifically, an appropriate autoregressive model can be established according to the motion trajectory of each pixel point.
S740:自适应选择图像序列上的像素点,对自回归模型进行自适应学习,确定自回归模型的模型参数。S740: Adaptively select pixels on the image sequence, perform adaptive learning on the autoregressive model, and determine model parameters of the autoregressive model.
S750:根据模型参数,对图像序列的像素点进行时间维度上的滤波,更新待重构图像的像素点的光强值以生成更新图像。S750: Perform filtering in the time dimension on the pixels of the image sequence according to the model parameters, and update the light intensity values of the pixels of the image to be reconstructed to generate an updated image.
具体地,根据本申请实施例提供的成像方法,可以重构图像序列中的一帧图像或多帧图像。Specifically, according to the imaging method provided by the embodiments of the present application, one frame of images or multiple frames of images in the image sequence can be reconstructed.
为了便于理解本申请实施例提供的基于脉冲信号的自适应成像方法,下面结合附图8进行说明。如图8所示,该方法包括如下内容。In order to facilitate the understanding of the adaptive imaging method based on the pulse signal provided by the embodiment of the present application, the following description is made with reference to FIG. 8 . As shown in Figure 8, the method includes the following contents.
首先,根据脉冲间隔信息从脉冲阵列中估计出基本的重建图像序列(重建基本图像)。由脉冲阵列生成基本的重建图像,可以采用脉冲间隔算法或其他重建方法。主要根据脉冲间的平均光流与脉冲间隔的大小成反比这一特性,重构出基本图像序列:
Figure PCTCN2022079318-appb-000035
其中,
Figure PCTCN2022079318-appb-000036
表示重建图像序列,
Figure PCTCN2022079318-appb-000037
表示每一张重建图像。然后,借助光流等方法根据基本重建图像序列确定图像各像素点的运动轨迹(确立运动轨迹)。
First, a basic reconstructed image sequence (reconstructed basic image) is estimated from the pulse array based on the pulse interval information. The basic reconstructed image is generated by the pulse array, and a pulse spacing algorithm or other reconstruction method can be used. According to the characteristic that the average optical flow between pulses is inversely proportional to the size of the pulse interval, the basic image sequence is reconstructed:
Figure PCTCN2022079318-appb-000035
in,
Figure PCTCN2022079318-appb-000036
represents the reconstructed image sequence,
Figure PCTCN2022079318-appb-000037
represents each reconstructed image. Then, the motion trajectory of each pixel point of the image is determined (establishing the motion trajectory) according to the basic reconstructed image sequence by means of methods such as optical flow.
接下来,沿着运动轨迹方向建立时间维度上的自回归模型(建立自回归模型),并根据邻域像素点的时间相关性自适应地调整模型参数(自适应学习模型参数)。Next, an autoregressive model in the time dimension is established along the direction of the motion trajectory (establishing an autoregressive model), and the model parameters are adaptively adjusted according to the temporal correlation of neighboring pixels (adaptive learning model parameters).
最后,根据得到的模型对图像像素进行时间维度的运动轨迹滤波(时间维度滤波),以重建出更高质量的图像。Finally, according to the obtained model, the motion trajectory filtering of the time dimension (time dimension filtering) is performed on the image pixels to reconstruct a higher quality image.
本申请沿着物体的运动轨迹,建立不同时刻图像像素点之间的自回归模型,并根据图像序列内容来自适应地调整模型参数,以准确地利用图像信号的时间方向相关性,从而提升重建图像的质量。The present application establishes an autoregressive model between image pixels at different times along the motion trajectory of the object, and adaptively adjusts the model parameters according to the content of the image sequence, so as to accurately utilize the temporal direction correlation of the image signal, thereby improving the reconstructed image. the quality of.
本申请的基于脉冲信号的自适应成像方法,根据预设时间段内的脉冲信号,得到待成像物体基于像素点的脉冲阵列;根据脉冲阵列,得到预设时间段内的重建图像序列;根据重建图像序列之间的相对运动关系,确定重建图像序列的各像素点的运动轨迹;根据各像素点的运动轨迹,建立各像素点在时间方向上的自回归模型;选择重建图像序列上的像素点,对自回归模型进行自适应学习,确定自回归模型的模型参数;根据模型参数,对重建图像序列的像素点进行时间维度上滤波,更新像素点的光强值以生成更新图像。The adaptive imaging method based on pulse signals of the present application obtains a pixel-based pulse array of an object to be imaged according to the pulse signals within a preset time period; obtains a reconstructed image sequence within a preset time period according to the pulse array; The relative motion relationship between the image sequences determines the motion trajectory of each pixel in the reconstructed image sequence; according to the motion trajectory of each pixel, the autoregressive model of each pixel in the time direction is established; the pixels on the reconstructed image sequence are selected , the autoregressive model is adaptively learned, and the model parameters of the autoregressive model are determined; according to the model parameters, the pixel points of the reconstructed image sequence are filtered in the time dimension, and the light intensity value of the pixel points is updated to generate an updated image.
本申请自适应地利用了脉冲阵列的时间方向相关性,提升了重建图像的质量。为了合理地利用脉冲阵列的时间方向相关性,首先,建立了沿着运动轨迹的时间方向上的自回归模型;然后,根据脉冲阵列局部空间位置上所具有的相关性结构自适应地调整模型参数,以提升自回归模型的准确性;最后,利用所建立的局部自适应的自回归模型对基本重建图像进行时间维度的滤波以提升信号的稳定性,从而提升重建图像的质量。The present application adaptively utilizes the temporal direction correlation of the pulse array to improve the quality of the reconstructed image. In order to reasonably utilize the temporal correlation of the pulse array, firstly, an autoregressive model along the time direction of the motion trajectory is established; then, the model parameters are adaptively adjusted according to the correlation structure of the local spatial position of the pulse array , in order to improve the accuracy of the autoregressive model; finally, the basic reconstructed image is filtered in the time dimension by the established local adaptive autoregressive model to improve the stability of the signal, thereby improving the quality of the reconstructed image.
其中,局部自适应的自回归模型保证了所利用的时间方向相关性的准确性,有效地降低了异常值的影响,保证了滤波算法的有效性和鲁棒性。Among them, the locally adaptive autoregressive model ensures the accuracy of the used time-direction correlation, effectively reduces the influence of outliers, and ensures the effectiveness and robustness of the filtering algorithm.
以下,通过具体场景应用来进一步说明本申请提供的成像方法,如基于仿生式脉冲相机的成像方法,或基于脉冲信号的自适应成像方法。Hereinafter, the imaging method provided by the present application, such as an imaging method based on a biomimetic pulse camera, or an adaptive imaging method based on a pulse signal, will be further described through specific application scenarios.
在一个示例性场景中,在脉冲相机所拍摄的真实脉冲阵列数据上进行实验,测试了本实施例提供的成像方法在高速运动场景下的重建性能。图9是根据一示例性实施例示出的不同成像方法拍摄的以100km/h速度行驶的汽车的示意图。如图9所示,其中(a)是初始脉冲矩阵图像,(b)是采用脉冲间隔法重建的图像,(c)是采用本申请实施例提供的成像方法重建的图像。由图片可知,根据本申请实施例的图像重建算法,可以得到高速行驶的汽车的清晰图像,实现较好的视觉效果。In an exemplary scene, an experiment is performed on real pulse array data captured by a pulse camera to test the reconstruction performance of the imaging method provided in this embodiment in a high-speed motion scene. FIG. 9 is a schematic diagram of a car traveling at a speed of 100 km/h shot by different imaging methods according to an exemplary embodiment. As shown in FIG. 9 , (a) is an initial pulse matrix image, (b) is an image reconstructed using the pulse interval method, and (c) is an image reconstructed using the imaging method provided by the embodiment of the present application. It can be seen from the pictures that, according to the image reconstruction algorithm of the embodiment of the present application, a clear image of a vehicle traveling at a high speed can be obtained, and a better visual effect can be achieved.
图10是根据一示例性实施例示出的不同成像方法所拍摄的自由下落的布娃娃的示意图。其中(a)是初始脉冲矩阵图像,(b)是采用脉冲间隔法重建的图像,(c)是采用本申请实施例提供的成像方法重建的图像。根据图像可知,本申请实施例能清晰地重建出场景中的运动细节,重建出的图像不仅避免了运动模糊而且具有高信噪比。FIG. 10 is a schematic diagram of a free-falling ragdoll photographed by different imaging methods according to an exemplary embodiment. (a) is an initial pulse matrix image, (b) is an image reconstructed using the pulse interval method, and (c) is an image reconstructed using the imaging method provided by the embodiments of the present application. It can be seen from the images that the embodiments of the present application can clearly reconstruct the motion details in the scene, and the reconstructed images not only avoid motion blur but also have a high signal-to-noise ratio.
本申请实施例还提供一种成像装置,该装置用于执行上述实施例的成像方法。对于本实施例的成像装置中未披露的细节,请参照上述实施例中的成像方法部分的描述。如图11所示,该装置包括:重构模块1110,更新模块1120以及成像模块1130。An embodiment of the present application further provides an imaging apparatus, which is used for executing the imaging method of the foregoing embodiment. For details not disclosed in the imaging device of this embodiment, please refer to the description of the imaging method in the above-mentioned embodiment. As shown in FIG. 11 , the apparatus includes: a reconstruction module 1110 , an update module 1120 and an imaging module 1130 .
重构模块1110,用于根据预设时间段内的脉冲阵列,得到重构的图像序列;更新模块1120,用于根据图像序列之间的相对运动关系,对图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新各像素点的光强值,得到更新光强值;成像模块1130,用于根据各像素点的更新光强值生成更新图像。The reconstruction module 1110 is used for obtaining the reconstructed image sequence according to the pulse array in the preset time period; the updating module 1120 is used for, according to the relative motion relationship between the image sequences, for the image frames that meet the specified conditions in the image sequence Each pixel of the pixel is filtered in the time dimension to update the light intensity value of each pixel to obtain an updated light intensity value; the imaging module 1130 is configured to generate an updated image according to the updated light intensity value of each pixel.
在一个实施例中,指定条件包括k时刻,k时刻位于预设时间段内。In one embodiment, the specified condition includes time k, which is within a preset time period.
在一个实施例中,更新模块1120用于:根据相对运动关系,确定各像素点在图像序列的t时刻的图像帧中的对应像素点;根据对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值。In one embodiment, the updating module 1120 is configured to: determine the corresponding pixel points of each pixel point in the image frame at time t of the image sequence according to the relative motion relationship; Filtering in the time dimension to get updated light intensity values.
在一个实施例中,更新模块1120用于:根据对应像素点的光强值,通过如下公式,对各像素点进行时间维度上的滤波以得到更新光强值,In one embodiment, the update module 1120 is configured to: perform filtering on each pixel in the time dimension according to the light intensity value of the corresponding pixel point by the following formula to obtain the updated light intensity value,
Figure PCTCN2022079318-appb-000038
Figure PCTCN2022079318-appb-000038
其中,I k(x,y)表示k时刻的图像帧内(x,y)像素点处的更新光强值,C表示正则化参数,
Figure PCTCN2022079318-appb-000039
表示t时刻的图像帧在重建更新图像时的滤波权重,
Figure PCTCN2022079318-appb-000040
表示k时刻的图像帧中的像素点在t时刻的图像帧中的对应像素点处的光强值,r表示滤波时间窗口半径。
Among them, I k (x, y) represents the updated light intensity value at the (x, y) pixel point in the image frame at time k, C represents the regularization parameter,
Figure PCTCN2022079318-appb-000039
Represents the filter weight of the image frame at time t when reconstructing the updated image,
Figure PCTCN2022079318-appb-000040
represents the light intensity value of the pixel in the image frame at time k at the corresponding pixel point in the image frame at time t, and r represents the radius of the filtering time window.
在一个实施例中,正则化参数C通过如下公式确定:In one embodiment, the regularization parameter C is determined by the following formula:
Figure PCTCN2022079318-appb-000041
Figure PCTCN2022079318-appb-000041
在一个实施例中,
Figure PCTCN2022079318-appb-000042
通过如下公式确定:
In one embodiment,
Figure PCTCN2022079318-appb-000042
Determined by the following formula:
Figure PCTCN2022079318-appb-000043
Figure PCTCN2022079318-appb-000043
其中,σ为设定参数。Among them, σ is the setting parameter.
在一个实施例中,设定参数通过如下公式确定:In one embodiment, the setting parameters are determined by the following formula:
σ=(2*r+1)/3。σ=(2*r+1)/3.
在一个实施例中,成像装置还包括:第一获取模块1140,用于获取对应像素点的光强值。In one embodiment, the imaging device further includes: a first acquiring module 1140, configured to acquire the light intensity value of the corresponding pixel point.
在一个实施例中,第一获取模块1140用于:若对应像素点为整像素点,则直接取对应像素点处的光强值;若对应像素点为亚像素点,则通过插值法得到对应像素点的光强值。In one embodiment, the first obtaining module 1140 is configured to: if the corresponding pixel is an integer pixel, directly obtain the light intensity value at the corresponding pixel; if the corresponding pixel is a sub-pixel, obtain the corresponding pixel by interpolation The light intensity value of the pixel.
在一个实施例中,更新模块1120用于:根据相对运动关系,对图像序列中像素点在时间方向上的自回归模型进行自适应学习,以确定自回归模型的模型参数;确定各像素点在图像序列的t时刻的图像帧中的对应像素点;根据模型参数以及对应像素点的光强值,对各像素点进行时间维度上的滤波以得到更新光强值。In one embodiment, the updating module 1120 is configured to: perform adaptive learning on the autoregressive model of the pixel points in the image sequence in the time direction according to the relative motion relationship, so as to determine the model parameters of the autoregressive model; The corresponding pixel points in the image frame at time t of the image sequence; according to the model parameters and the light intensity value of the corresponding pixel point, each pixel point is filtered in the time dimension to obtain the updated light intensity value.
在一个实施例中,自回归模型的公式为:In one embodiment, the formula of the autoregressive model is:
Figure PCTCN2022079318-appb-000044
Figure PCTCN2022079318-appb-000044
其中,I k表示k时刻的理论光强值;I t表示t时刻的理论光强值;T k表示k时刻相邻时刻的集合;P k→t(x,y)为k时刻的图像帧内的像素点(x,y)在t时刻的图像帧内的对应像素位置;α t为模型参数;Δ k(x,y)为在(x,y)像素点上的模型误差。 Among them, I k represents the theoretical light intensity value at time k; I t represents the theoretical light intensity value at time t; T k represents the set of adjacent moments at time k; P k→t (x, y) is the image frame at time k The corresponding pixel position of the pixel point (x, y) in the image frame at time t; α t is the model parameter; Δ k (x, y) is the model error on the (x, y) pixel point.
在一个实施例中,更新模块1120用于:根据相对运动关系,通过最小二乘法求解目标函数的方式,对自回归模型进行自适应学习,以确定模型参数。In one embodiment, the update module 1120 is configured to: perform adaptive learning on the autoregressive model by solving the objective function by the least squares method according to the relative motion relationship, so as to determine the model parameters.
在一个实施例中,目标函数为:In one embodiment, the objective function is:
Figure PCTCN2022079318-appb-000045
Figure PCTCN2022079318-appb-000045
其中,W表示以(x,y)为中心的空间范围内像素点集合;
Figure PCTCN2022079318-appb-000046
表示k时刻的估计光强值;
Figure PCTCN2022079318-appb-000047
表示t时刻的估计光强值。
Among them, W represents the set of pixels in the spatial range centered on (x, y);
Figure PCTCN2022079318-appb-000046
Represents the estimated light intensity value at time k;
Figure PCTCN2022079318-appb-000047
represents the estimated light intensity value at time t.
在一个实施例中,更新模块1120用于:根据模型参数以及对应像素点的光强值,通过如下公式,对各像素点进行时间维度上的滤波以得到更新光强值,In one embodiment, the update module 1120 is configured to: filter each pixel in the time dimension to obtain the updated light intensity value by using the following formula according to the model parameters and the light intensity value of the corresponding pixel point,
Figure PCTCN2022079318-appb-000048
Figure PCTCN2022079318-appb-000048
其中,I′ k表示k时刻的更新光强值。 Among them, I′ k represents the updated light intensity value at time k.
在一个实施例中,成像装置还包括:第二获取模块1150,用于:根据光流法、像素点匹配、像素点运动对齐和像素点相对位置偏移估计中的任一项或多项,得到图像序列之间的相对运动关系。In one embodiment, the imaging device further includes: a second acquisition module 1150, configured to: according to any one or more of optical flow method, pixel point matching, pixel point motion alignment and pixel point relative position offset estimation, Obtain the relative motion relationship between image sequences.
在一个实施例中,成像装置还包括:第三获取模块1160,用于:根据预设时间段内的脉冲信号,获取待成像物体基于像素点的脉冲阵列。In one embodiment, the imaging device further includes: a third acquisition module 1160, configured to: acquire a pixel-based pulse array of the object to be imaged according to the pulse signal within a preset time period.
在一个实施例中,重构模块1110用于:根据脉冲阵列,采用脉冲间隔算法重构出图像序列。In one embodiment, the reconstruction module 1110 is configured to: reconstruct an image sequence by using a pulse interval algorithm according to the pulse array.
需要说明的是,上述实施例提供的成像装置在执行成像方法时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的成像装置与成像方法实施例属于同一构思,其具体实现过程以及效果详见方法实施例部分的描述,这里不再赘述。It should be noted that, when the imaging apparatus provided in the above-mentioned embodiments executes the imaging method, only the division of the above-mentioned functional modules is used as an example. The internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the imaging apparatus and the imaging method embodiments provided by the above embodiments belong to the same concept, and the specific implementation process and effects thereof are described in the description of the method embodiment section, which will not be repeated here.
本申请实施例还提供一种与前述实施例所提供的成像方法对应的电子设备,以执行上述成像方法。The embodiment of the present application further provides an electronic device corresponding to the imaging method provided by the foregoing embodiments, so as to execute the foregoing imaging method.
请参考图12,其示出了本申请的一些实施例所提供的一种电子设备的示意图。如图12所示,电子设备包括:处理器1210以及存储器1220。存储器1220中存储有可在处理器1210上运行的计算机程序,处理器1210运行计算机程序时执行本申请前述任一实施例所提供的成像方法。进一步地,电子设备还包括总线1230和通信接口1240。处理器1210、通信接口1240和存储器1220通过总线1230连接。Please refer to FIG. 12 , which shows a schematic diagram of an electronic device provided by some embodiments of the present application. As shown in FIG. 12 , the electronic device includes: a processor 1210 and a memory 1220 . The memory 1220 stores a computer program that can be executed on the processor 1210. When the processor 1210 runs the computer program, the imaging method provided by any of the foregoing embodiments of the present application is executed. Further, the electronic device also includes a bus 1230 and a communication interface 1240 . The processor 1210 , the communication interface 1240 and the memory 1220 are connected through the bus 1230 .
其中,存储器1220可能包含高速随机存取存储器(RAM:Random Access Memory),也可能还包括非不稳定的存储器(non-volatile memory),例如至少一个磁盘存储器。通过至少一个通信接口1240(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接,可以使用互联网、广域网、本地网、城域网等。The memory 1220 may include a high-speed random access memory (RAM: Random Access Memory), and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory. The communication connection between the network element of the system and at least one other network element is implemented through at least one communication interface 1240 (which may be wired or wireless), which may use the Internet, a wide area network, a local network, a metropolitan area network, and the like.
总线1230可以是ISA总线、PCI总线或EISA总线等。总线可以分为地址总线、数据总线、控制总线等。其中,存储器1220用于存储程序,处理器1210在接收到执行指令后,执行程序,前述本申请实施例任一实施方式揭示的基于仿生式脉冲相机的成像方法可以应用于处理器1210中,或者由处理器1210实现。The bus 1230 may be an ISA bus, a PCI bus, an EISA bus, or the like. The bus can be divided into address bus, data bus, control bus and so on. The memory 1220 is used to store a program, and the processor 1210 executes the program after receiving the execution instruction. The imaging method based on the bionic pulse camera disclosed in any of the foregoing embodiments of the present application may be applied to the processor 1210, or Implemented by processor 1210 .
处理器1210可能是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器1210中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器900可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。处理器1210可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。本申请实施例所公开的方法的步骤可以直接由硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器1220,处理器1210读取存储器1220中的信息,结合其硬件完成上述方法的步骤。The processor 1210 may be an integrated circuit chip with signal processing capability. In the implementation process, each step of the above-mentioned method can be completed by an integrated logic circuit of hardware in the processor 1210 or an instruction in the form of software. The above-mentioned processor 900 can be a general-purpose processor, including a central processing unit (Central Processing Unit, referred to as CPU), a network processor (Network Processor, referred to as NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), off-the-shelf programmable gate array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware component. The processor 1210 may implement or execute the methods, steps, and logical block diagrams disclosed in the embodiments of this application. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the methods disclosed in the embodiments of the present application may be directly executed by a hardware decoding processor, or executed by a combination of hardware and software modules in the decoding processor. The software modules may be located in random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, registers and other storage media mature in the art. The storage medium is located in the memory 1220, and the processor 1210 reads the information in the memory 1220, and completes the steps of the above method in combination with its hardware.
本申请实施例提供的电子设备与本申请实施例提供的成像方法出于相同的发明构思,具有与其采用、运行或实现的方法相同的有益效果。The electronic device provided by the embodiment of the present application and the imaging method provided by the embodiment of the present application are based on the same inventive concept, and have the same beneficial effects as the method adopted, operated or realized.
本申请实施例还提供一种成像设备,包括上述实施例提供的电子设备。例如,成像设备可以是相机、监视器等。Embodiments of the present application further provide an imaging device, including the electronic device provided by the foregoing embodiments. For example, the imaging device may be a camera, a monitor, or the like.
本申请实施例还提供一种与前述实施例所提供的成像方法对应的计算机可读存储介质,请参考图13,其示出的计算机可读存储介质为光盘1300,其上存储有计算机程序(即程序产品),计算机程序在被处理器运行时,会执行前述任意实施例所提供的成像方法。Embodiments of the present application also provide a computer-readable storage medium corresponding to the imaging method provided by the foregoing embodiments. Please refer to FIG. 13 , the computer-readable storage medium shown is an optical disc 1300, on which a computer program ( That is, a program product), when the computer program is executed by the processor, the imaging method provided by any of the foregoing embodiments will be executed.
需要说明的是,计算机可读存储介质的例子还可以包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器 (RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他光学、磁性存储介质,在此不再一一赘述。It should be noted that examples of computer-readable storage media may also include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory Memory (RAM), Read-Only Memory (ROM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Flash Memory or other optical and magnetic storage media will not be repeated here.
本申请的上述实施例提供的计算机可读存储介质与本申请实施例提供的成像方法出于相同的发明构思,具有与其存储的应用程序所采用、运行或实现的方法相同的有益效果。The computer-readable storage medium provided by the above embodiments of the present application and the imaging methods provided by the embodiments of the present application are based on the same inventive concept, and have the same beneficial effects as the methods used, run or implemented by the application programs stored thereon.
以上所述,仅为本申请较佳的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。The above are only the preferred specific embodiments of the present application, but the protection scope of the present application is not limited to this. Substitutions should be covered within the protection scope of this application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims (20)

  1. 一种成像方法,其特征在于,包括:An imaging method, characterized in that, comprising:
    根据预设时间段内的脉冲阵列,得到重构的图像序列;Obtain a reconstructed image sequence according to the pulse array within a preset time period;
    根据所述图像序列之间的相对运动关系,对所述图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新所述各像素点的光强值,得到更新光强值;According to the relative motion relationship between the image sequences, filtering in the time dimension is performed on each pixel of the image frame in the image sequence that satisfies the specified condition to update the light intensity value of each pixel to obtain the updated light intensity value;
    根据所述各像素点的所述更新光强值生成更新图像。An updated image is generated according to the updated light intensity value of each pixel point.
  2. 根据权利要求1所述的方法,其特征在于,所述指定条件包括k时刻,所述k时刻位于所述预设时间段内。The method according to claim 1, wherein the specified condition includes time k, and the time k is within the preset time period.
  3. 根据权利要求2所述的方法,其特征在于,所述根据所述图像序列之间的相对运动关系,对所述图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新所述各像素点的光强值,得到更新光强值,包括:The method according to claim 2, wherein, according to the relative motion relationship between the image sequences, filtering in the time dimension is performed on each pixel of the image frame in the image sequence that meets the specified condition to Update the light intensity value of each pixel to obtain the updated light intensity value, including:
    根据所述相对运动关系,确定所述各像素点在所述图像序列的t时刻的图像帧中的对应像素点;According to the relative motion relationship, determine the corresponding pixel points of each pixel point in the image frame at time t of the image sequence;
    根据所述对应像素点的光强值,对所述各像素点进行时间维度上的滤波以得到所述更新光强值。According to the light intensity value of the corresponding pixel point, filtering in the time dimension is performed on each pixel point to obtain the updated light intensity value.
  4. 根据权利要求3所述的方法,其特征在于,所述根据所述对应像素点的光强值,对所述各像素点进行时间维度上的滤波以得到所述更新光强值,包括:The method according to claim 3, wherein, according to the light intensity value of the corresponding pixel point, filtering the pixel points in the time dimension to obtain the updated light intensity value, comprising:
    根据所述对应像素点的光强值,通过如下公式,对所述各像素点进行时间维度上的滤波以得到所述更新光强值,According to the light intensity value of the corresponding pixel point, the following formula is used to filter each pixel point in the time dimension to obtain the updated light intensity value,
    Figure PCTCN2022079318-appb-100001
    Figure PCTCN2022079318-appb-100001
    其中,I k(x,y)表示所述k时刻的图像帧内(x,y)像素点处的更新光强值,C表示正则化参数,
    Figure PCTCN2022079318-appb-100002
    表示所述t时刻的图像帧在重建所述更新图像时的滤波权重,
    Figure PCTCN2022079318-appb-100003
    表示所述k时刻的图像帧中的像素点在所述t时刻的图像帧中的对应像素点处的光强值,r表示滤波时间窗口半径。
    Among them, I k (x, y) represents the updated light intensity value at the (x, y) pixel point in the image frame at the k moment, C represents the regularization parameter,
    Figure PCTCN2022079318-appb-100002
    represents the filter weight of the image frame at time t when reconstructing the updated image,
    Figure PCTCN2022079318-appb-100003
    represents the light intensity value of the pixel point in the image frame at time k at the corresponding pixel point in the image frame at time t, and r represents the radius of the filtering time window.
  5. 根据权利要求4所述的方法,其特征在于,所述正则化参数C通过如下公式确定:The method according to claim 4, wherein the regularization parameter C is determined by the following formula:
    Figure PCTCN2022079318-appb-100004
    Figure PCTCN2022079318-appb-100004
  6. 根据权利要求4或5所述的方法,其特征在于,
    Figure PCTCN2022079318-appb-100005
    通过如下公式确定:
    The method according to claim 4 or 5, wherein,
    Figure PCTCN2022079318-appb-100005
    Determined by the following formula:
    Figure PCTCN2022079318-appb-100006
    Figure PCTCN2022079318-appb-100006
    其中,σ为设定参数。Among them, σ is the setting parameter.
  7. 根据权利要求6所述的方法,其特征在于,所述设定参数通过如下公式确定:The method according to claim 6, wherein the set parameter is determined by the following formula:
    σ=(2*r+1)/3。σ=(2*r+1)/3.
  8. 根据权利要求3至7中任一项所述的方法,其特征在于,还包括:The method according to any one of claims 3 to 7, further comprising:
    若所述对应像素点为整像素点,则直接取所述对应像素点处的光强值;If the corresponding pixel is an integer pixel, directly take the light intensity value at the corresponding pixel;
    若所述对应像素点为亚像素点,则通过插值法得到所述对应像素点的光强值。If the corresponding pixel point is a sub-pixel point, the light intensity value of the corresponding pixel point is obtained by an interpolation method.
  9. 根据权利要求2所述的方法,其特征在于,所述根据所述图像序列之间的相对运动关系,对所述图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新所述各像素点的光强值,得到更新光强值,包括:The method according to claim 2, wherein, according to the relative motion relationship between the image sequences, filtering in the time dimension is performed on each pixel of the image frame in the image sequence that meets the specified condition to Update the light intensity value of each pixel to obtain the updated light intensity value, including:
    根据所述相对运动关系,对所述图像序列中像素点在时间方向上的自回归模型进行自适应学习,以确定所述自回归模型的模型参数;According to the relative motion relationship, adaptive learning is performed on the autoregressive model of the pixel points in the image sequence in the time direction, so as to determine the model parameters of the autoregressive model;
    确定所述各像素点在所述图像序列的t时刻的图像帧中的对应像素点;Determine the corresponding pixel points of each pixel point in the image frame at time t of the image sequence;
    根据所述模型参数以及所述对应像素点的光强值,对所述各像素点进行时间维度上的滤波以得到所述更新光强值。According to the model parameter and the light intensity value of the corresponding pixel point, filtering in the time dimension is performed on each pixel point to obtain the updated light intensity value.
  10. 根据权利要求9所述的方法,其特征在于,所述自回归模型的公式为:The method according to claim 9, wherein the formula of the autoregressive model is:
    Figure PCTCN2022079318-appb-100007
    Figure PCTCN2022079318-appb-100007
    其中,I k表示所述k时刻的理论光强值;I t表示所述t时刻的理论光强值;T k表示所述k时刻相邻时刻的集合;P k→t(x,y)为所述k时刻的图像帧内的像素点(x,y)在所述t时刻的图像帧内的对应像素位置;α t为所述模型参数;Δ k(x,y)为在(x,y)像素点上的模型误差。 Among them, I k represents the theoretical light intensity value at time k; I t represents the theoretical light intensity value at time t; T k represents the set of adjacent moments at time k; P k→t (x, y) is the corresponding pixel position of the pixel point (x, y) in the image frame at the time k at the time t; α t is the model parameter; Δ k (x, y) is the pixel at (x, y) , y) the model error at the pixel point.
  11. 根据权利要求10所述的方法,其特征在于,所述根据所述相对运动关系,对所述图像序列中像素点在时间方向上的自回归模型进行自适应学习,以确定所述自回归模型的模型参数,包括:The method according to claim 10, wherein the autoregressive model of the pixels in the image sequence in the time direction is adaptively learned according to the relative motion relationship to determine the autoregressive model model parameters, including:
    根据所述相对运动关系,通过最小二乘法求解目标函数的方式,对所述自回归模型进行自适应学习,以确定所述模型参数。According to the relative motion relationship, the autoregressive model is adaptively learned by means of least squares method to solve the objective function, so as to determine the model parameters.
  12. 根据权利要求11所述的方法,其特征在于,所述目标函数为:The method according to claim 11, wherein the objective function is:
    Figure PCTCN2022079318-appb-100008
    Figure PCTCN2022079318-appb-100008
    其中,W表示以(x,y)为中心的空间范围内像素点集合;
    Figure PCTCN2022079318-appb-100009
    表示所述k时刻的估计光强值;
    Figure PCTCN2022079318-appb-100010
    表示所述t时刻的估计光强值。
    Among them, W represents the set of pixels in the spatial range centered on (x, y);
    Figure PCTCN2022079318-appb-100009
    represents the estimated light intensity value at time k;
    Figure PCTCN2022079318-appb-100010
    represents the estimated light intensity value at the time t.
  13. 根据权利要求10至12中任一项所述的方法,其特征在于,所述根据所述模型参数以及所述对应像素点的光强值,对所述各像素点进行时间维度上的滤波以得到所述更新光强值,包括:The method according to any one of claims 10 to 12, wherein, according to the model parameter and the light intensity value of the corresponding pixel point, filtering the pixel points in the time dimension to Get the updated light intensity value, including:
    根据所述模型参数以及所述对应像素点的光强值,通过如下公式,对所述各像素点进行时间维度上的滤波以得到所述更新光强值,According to the model parameters and the light intensity value of the corresponding pixel point, the following formula is used to filter each pixel point in the time dimension to obtain the updated light intensity value,
    Figure PCTCN2022079318-appb-100011
    Figure PCTCN2022079318-appb-100011
    其中,I′ k表示所述k时刻的更新光强值。 Wherein, I′ k represents the updated light intensity value at time k.
  14. 根据权利要求1至13中任一项所述的方法,其特征在于,还包括:The method according to any one of claims 1 to 13, further comprising:
    根据光流法、像素点匹配、像素点运动对齐和像素点相对位置偏移估计中的任一项或多项,得到所述图像序列之间的所述相对运动关系。The relative motion relationship between the image sequences is obtained according to any one or more of the optical flow method, pixel point matching, pixel point motion alignment, and pixel point relative position offset estimation.
  15. 根据权利要求1至14中任一项所述的方法,其特征在于,在所述根据预设时间段内的脉冲阵列,得到重构的图像序列之前,还包括:The method according to any one of claims 1 to 14, characterized in that before obtaining the reconstructed image sequence according to the pulse array within a preset time period, the method further comprises:
    根据所述预设时间段内的脉冲信号,获取待成像物体基于像素点的所述脉冲阵列。According to the pulse signal within the preset time period, the pixel-based pulse array of the object to be imaged is acquired.
  16. 根据权利要求1至15中任一项所述的方法,其特征在于,所述根据预设时间段内的脉冲阵列,得到重构的图像序列,包括:The method according to any one of claims 1 to 15, wherein the obtaining a reconstructed image sequence according to the pulse array within a preset time period comprises:
    根据所述脉冲阵列,采用脉冲间隔算法重构出所述图像序列。From the pulse array, the image sequence is reconstructed using a pulse spacing algorithm.
  17. 一种成像装置,其特征在于,包括:An imaging device, characterized in that, comprising:
    重构模块,用于根据预设时间段内的脉冲阵列,得到重构的图像序列;a reconstruction module, configured to obtain a reconstructed image sequence according to the pulse array within a preset time period;
    更新模块,用于根据所述图像序列之间的相对运动关系,对所述图像序列中满足指定条件的图像帧的各像素点进行时间维度上的滤波以更新所述各像素点的光强值,得到更新光强值;an update module, configured to perform filtering in the time dimension on each pixel point of the image frame in the image sequence that meets the specified condition according to the relative motion relationship between the image sequences to update the light intensity value of each pixel point , get the updated light intensity value;
    成像模块,用于根据所述各像素点的所述更新光强值生成更新图像。An imaging module, configured to generate an updated image according to the updated light intensity value of each pixel point.
  18. 一种电子设备,其特征在于,包括处理器和存储有程序指令的存储器,所述处理器被配置为在执行所述程序指令时,执行如权利要求1至16中任一项所述的成像方法。An electronic device, characterized by comprising a processor and a memory storing program instructions, the processor is configured to perform the imaging according to any one of claims 1 to 16 when executing the program instructions method.
  19. 一种成像设备,其特征在于,包括如权利要求18所述的电子设备。An imaging device, characterized by comprising the electronic device as claimed in claim 18 .
  20. 一种计算机可读介质,其特征在于,其上存储有计算机可读指令,所述计算机可读指令可被处理器执行以实现如权利要求1至16中任一项所述的成像方法。A computer-readable medium having computer-readable instructions stored thereon, the computer-readable instructions being executable by a processor to implement the imaging method according to any one of claims 1 to 16.
PCT/CN2022/079318 2021-03-04 2022-03-04 Imaging method and apparatus, device, and storage medium WO2022184167A1 (en)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN202110240363.2A CN113034634B (en) 2021-03-04 2021-03-04 Adaptive imaging method, system and computer medium based on pulse signal
CN202110240363.2 2021-03-04
CN202110253200.8 2021-03-04
CN202110253200.8A CN113067979A (en) 2021-03-04 2021-03-04 Imaging method, device, equipment and storage medium based on bionic pulse camera

Publications (1)

Publication Number Publication Date
WO2022184167A1 true WO2022184167A1 (en) 2022-09-09

Family

ID=83154756

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/079318 WO2022184167A1 (en) 2021-03-04 2022-03-04 Imaging method and apparatus, device, and storage medium

Country Status (1)

Country Link
WO (1) WO2022184167A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116012675A (en) * 2023-02-14 2023-04-25 荣耀终端有限公司 Model training method, image processing method and electronic equipment
CN116456202A (en) * 2023-04-25 2023-07-18 北京大学 Pulse camera and color imaging method and device thereof

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200128245A1 (en) * 2016-01-22 2020-04-23 Peking University Imaging method and device
US20200397350A1 (en) * 2019-06-20 2020-12-24 Ethicon Llc Super resolution and color motion artifact correction in a pulsed hyperspectral imaging system
WO2021005608A1 (en) * 2019-07-09 2021-01-14 Voyage81 Real-time image generation in moving scenes
CN113034634A (en) * 2021-03-04 2021-06-25 北京大学 Adaptive imaging method, system and computer medium based on pulse signal
CN113067979A (en) * 2021-03-04 2021-07-02 北京大学 Imaging method, device, equipment and storage medium based on bionic pulse camera

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200128245A1 (en) * 2016-01-22 2020-04-23 Peking University Imaging method and device
US20200397350A1 (en) * 2019-06-20 2020-12-24 Ethicon Llc Super resolution and color motion artifact correction in a pulsed hyperspectral imaging system
WO2021005608A1 (en) * 2019-07-09 2021-01-14 Voyage81 Real-time image generation in moving scenes
CN113034634A (en) * 2021-03-04 2021-06-25 北京大学 Adaptive imaging method, system and computer medium based on pulse signal
CN113067979A (en) * 2021-03-04 2021-07-02 北京大学 Imaging method, device, equipment and storage medium based on bionic pulse camera

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
ZHAO JING; XIONG RUIQIN; HUANG TIEJUN: "High-Speed Motion Scene Reconstruction for Spike Camera via Motion Aligned Filtering", 2019 IEEE INTERNATIONAL SYMPOSIUM ON CIRCUITS AND SYSTEMS (ISCAS), IEEE, 12 October 2020 (2020-10-12), pages 1 - 5, XP033933546, ISSN: 2158-1525, ISBN: 978-1-7281-3320-1, DOI: 10.1109/ISCAS45731.2020.9181055 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116012675A (en) * 2023-02-14 2023-04-25 荣耀终端有限公司 Model training method, image processing method and electronic equipment
CN116012675B (en) * 2023-02-14 2023-08-11 荣耀终端有限公司 Model training method, image processing method and electronic equipment
CN116456202A (en) * 2023-04-25 2023-07-18 北京大学 Pulse camera and color imaging method and device thereof
CN116456202B (en) * 2023-04-25 2023-12-15 北京大学 Pulse camera and color imaging method and device thereof

Similar Documents

Publication Publication Date Title
Chen et al. Reblur2deblur: Deblurring videos via self-supervised learning
WO2022184167A1 (en) Imaging method and apparatus, device, and storage medium
CN113067979A (en) Imaging method, device, equipment and storage medium based on bionic pulse camera
EP3046320B1 (en) Method for generating an hdr image of a scene based on a tradeoff between brightness distribution and motion
Xu et al. Motion deblurring with real events
CN102077572B (en) Method and apparatus for motion blur and ghosting prevention in imaging system
CN113034634B (en) Adaptive imaging method, system and computer medium based on pulse signal
CN102970464B (en) Image processing apparatus and image processing method
CN111709895A (en) Image blind deblurring method and system based on attention mechanism
CN114584703B (en) Imaging method, device, equipment and storage medium of bionic pulse camera
US11303793B2 (en) System and method for high-resolution, high-speed, and noise-robust imaging
CN112288642A (en) Ghost detection method, image fusion method and corresponding device
Zhao et al. Super resolve dynamic scene from continuous spike streams
CN112837245A (en) Dynamic scene deblurring method based on multi-mode fusion
Li et al. Real-world deep local motion deblurring
Zhou et al. DeLiEve-Net: Deblurring low-light images with light streaks and local events
Zhang et al. Self-supervised image restoration with blurry and noisy pairs
US8213496B2 (en) Image processing device, image processing method, and image processing program
Huang et al. Inverting the imaging process by learning an implicit camera model
CN116091337B (en) Image enhancement method and device based on event signal nerve coding mode
CN111798484A (en) Continuous dense optical flow estimation method and system based on event camera
Nayak et al. Morphology based iterative back-projection for super-resolution reconstruction of image
Zhang et al. Bracketing is All You Need: Unifying Image Restoration and Enhancement Tasks with Multi-Exposure Images
US11749004B2 (en) Event detector and method of generating textural image based on event count decay factor and net polarity
CN114630009A (en) Noise and tailing removing method for pulse sequence type image sensor

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22762629

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 19-01-2024)