CN115994865A - Restoration method, acquisition restoration device and monitoring system for flutter blurred image - Google Patents

Restoration method, acquisition restoration device and monitoring system for flutter blurred image Download PDF

Info

Publication number
CN115994865A
CN115994865A CN202211315365.4A CN202211315365A CN115994865A CN 115994865 A CN115994865 A CN 115994865A CN 202211315365 A CN202211315365 A CN 202211315365A CN 115994865 A CN115994865 A CN 115994865A
Authority
CN
China
Prior art keywords
image
module
flutter
length
blurred
Prior art date
Legal status (The legal status 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 status listed.)
Pending
Application number
CN202211315365.4A
Other languages
Chinese (zh)
Inventor
李响
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Keshi Dalian Embedded Technology Development Co ltd
Original Assignee
Keshi Dalian Embedded Technology Development Co ltd
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
Application filed by Keshi Dalian Embedded Technology Development Co ltd filed Critical Keshi Dalian Embedded Technology Development Co ltd
Priority to CN202211315365.4A priority Critical patent/CN115994865A/en
Publication of CN115994865A publication Critical patent/CN115994865A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

The invention discloses a method for recovering a flutter blurred image, a device for collecting and recovering the flutter blurred image and a monitoring system, which comprise the following steps based on structural similarity and entropy: the method comprises the steps of collecting a flutter fuzzy image, constructing a fuzzy kernel through presetting a binary coding sequence and fuzzy lengths, calculating structural similarity indexes of the flutter fuzzy image and all decoded clear images, determining an image information entropy value in a certain range near the image with the maximum structural similarity, calculating an information entropy value of a reconstructed image, determining the minimum information entropy value as the fuzzy length of the most orderly reconstructed decoded image, and adopting an inverse filtering method to realize the clear restoration of the image. The invention starts from the aim of completely protecting the image information of the target, and simultaneously considers the method for removing the flutter fuzzy target from two angles of imaging and image post-processing. Modulating the exposure process by utilizing a pre-coding mode, protecting the original high-frequency information in the acquired image, and decoding the acquired image in the subsequent image restoration process to obtain a clear image.

Description

Restoration method, acquisition restoration device and monitoring system for flutter blurred image
Technical Field
The invention relates to the technical field of image restoration, in particular to a restoration method, an acquisition restoration device and a monitoring system for a flutter blurred image.
Background
With the development of digital storage media, consumer electronic devices such as digital cameras and smart phones have become a main source for acquiring emerging image information, and image clarity is one of the primary factors for guaranteeing reliable and effective information acquisition. Today, where digital cameras are popular, image information is collected anytime and anywhere, and each collected image contains characteristics of a target. The image is optimal when the acquired image can correctly reflect the characteristics of the target area. When the acquired image is unclear or the acquired image can not reflect the specific characteristics of the target area at the time, the effective information is lost, and the image is degraded. If the difference of the video monitoring images between day and night illumination is large, aviation and remote sensing images are easily affected by the atmosphere, and underwater imaging is affected by multipath interference such as refraction and the like. Therefore, when the image is degraded, the effective information is missing or wrong, so that the image quality is seriously affected, and how to protect the original information of the image is a problem to be solved in front-end image acquisition and back-end computer image processing.
The coded exposure imaging technique is a computational imaging means proposed by Raskar et al in 2006, and the core idea is to control the opening and closing of the camera shutter by presetting a specific binary coding sequence during the exposure of the camera. Compared with the traditional camera exposure mode, the method is equivalent to setting a broadband filter on a time domain, so that as much as possible of medium-high frequency information is reserved in the image acquisition process from the perspective of frequency domain rate analysis, the zero point part on the frequency domain is eliminated, the reversibility of restoration is realized, and the disease state problem of the restoration of the blurred image is greatly improved.
Image degradation is a comprehensive characterization of image blur, distortion, noise, etc., and as a result, the loss of image information. The reason for image degradation is various, and the image is generally subjected to the processes of acquisition, imaging, transmission, storage and the like, and each part is likely to cause image degradation. Motion blur is a problem frequently encountered in the optical imaging process, and the reason for the generation is that the acquired image is blurred due to the relative displacement between the acquired object and the camera in the camera exposure process, so that the image resolution is reduced, and the imaging quality is greatly influenced. The motion blurred image restoration technology is to restore the blurred image into a clear image by modeling and mathematical solving an imaging physical process by utilizing the existing motion blurred image on the premise of not re-acquiring a target scenery, and has important application value in the fields of civil, military and the like.
The Chinese patent CN110097509B is used for acquiring a fuzzy image by utilizing a coding exposure mode, performing primary extraction of a target by adopting a background difference method, and then comprehensively coding exposure motion fuzzy superposition characteristics and motion priori information to realize accurate extraction of a motion fuzzy target region, and performing PSF accurate estimation and restoration reconstruction by combining a student-t restoration algorithm, wherein a restoration result can be obtained after 2-3 iterations. In the present invention, it is necessary to input a scene background image and a local moving object blurred image. The Chinese patent publication CN202011448781.2 utilizes a deep learning method to estimate the fuzzy core, solves the problem that the fuzzy core is difficult to estimate in the traditional method, and searches out the image block with the highest content of the high-frequency information of the image by using a genetic algorithm. But does not solve the matching problem between coding and motion blur length. Because of the important role of the accuracy of blur length estimation in image restoration reconstruction, the present invention utilizes a method of jointly estimating blur length based on image Structural Similarity (SSIM) and image information Entropy (Entropy) to restore the encoded exposure image.
Disclosure of Invention
In order to effectively solve the technical problems, the invention adopts the following technical scheme based on structural similarity and entropy: a method for restoring a flutter blurred image,
step one: collecting flutter blurred images;
step two: by presetting the code length to be
Figure 477909DEST_PATH_IMAGE001
Binary coding sequence->
Figure 129470DEST_PATH_IMAGE002
And blur Length->
Figure 846891DEST_PATH_IMAGE003
Construction of fuzzy core->
Figure 242100DEST_PATH_IMAGE004
Said->
Figure 728576DEST_PATH_IMAGE001
For binary coding sequences->
Figure 816618DEST_PATH_IMAGE002
Is defined by the formula->
Figure 83651DEST_PATH_IMAGE005
Calculating to obtain a reconstructed image, wherein->
Figure 954655DEST_PATH_IMAGE006
For flutter blurred image +.>
Figure 623534DEST_PATH_IMAGE007
To decode an image;
step three: calculating a flutter blurred image
Figure 820160DEST_PATH_IMAGE008
And all decoded pictures +.>
Figure 308910DEST_PATH_IMAGE007
Structural similarity index->
Figure 311501DEST_PATH_IMAGE009
Length of blur at the time of searching for maximum +.>
Figure 772570DEST_PATH_IMAGE010
The formula is: />
Figure 936835DEST_PATH_IMAGE011
In the formula->
Figure 110284DEST_PATH_IMAGE012
Search range for fuzzy length;
step four: determining an image information entropy value in a certain range near an image with the maximum structural similarity, wherein a search interval is
Figure 650987DEST_PATH_IMAGE013
Wherein->
Figure 966561DEST_PATH_IMAGE014
Is positive integer->
Figure 567307DEST_PATH_IMAGE015
,/>
Figure 30649DEST_PATH_IMAGE016
Is made of->
Figure 312726DEST_PATH_IMAGE017
The total blur length limit range is determined to be +.>
Figure 810704DEST_PATH_IMAGE016
Within the limits according to the formula ∈ ->
Figure 254454DEST_PATH_IMAGE018
And
Figure 470672DEST_PATH_IMAGE019
calculating spatial entropy of each decoded image>
Figure 353178DEST_PATH_IMAGE020
And spectral entropy after image DCT (discrete cosine transform)>
Figure 643345DEST_PATH_IMAGE021
Wherein->
Figure 320314DEST_PATH_IMAGE022
Wherein->
Figure 695931DEST_PATH_IMAGE023
Is a characteristic binary group->
Figure 382128DEST_PATH_IMAGE024
Frequency of occurrence, < >>
Figure 526801DEST_PATH_IMAGE025
Gray representing a pixelMetric value->
Figure 640251DEST_PATH_IMAGE026
Representing neighborhood gray scale mean,/->
Figure 565481DEST_PATH_IMAGE027
Is the scale of the image; if signal->
Figure 727472DEST_PATH_IMAGE028
DCT transformation into (1)
Figure 54548DEST_PATH_IMAGE029
,/>
Figure 276582DEST_PATH_IMAGE030
Is a generalized frequency domain variable, normalized DCT transform coefficient is +.>
Figure 423530DEST_PATH_IMAGE031
Step five: formula (VI)
Figure 451529DEST_PATH_IMAGE032
Calculating the information entropy value of the reconstructed image>
Figure 570795DEST_PATH_IMAGE033
Wherein->
Figure 26047DEST_PATH_IMAGE034
For the picture space entropy value->
Figure 863553DEST_PATH_IMAGE020
Weight of occupied->
Figure 429663DEST_PATH_IMAGE035
Spectral entropy values for the coefficients of the DCT-transform (discrete cosine transform) of an image>
Figure 465752DEST_PATH_IMAGE021
The weight is occupied;
step six: taking out
Figure 29589DEST_PATH_IMAGE033
Minimum value, determined as most ordered reconstructed decoded picture +.>
Figure 151129DEST_PATH_IMAGE036
Is +.>
Figure 458613DEST_PATH_IMAGE037
The formula is->
Figure 614788DEST_PATH_IMAGE038
;/>
Step seven: the image is restored clearly by adopting an inverse filtering method, and the formula is
Figure 411843DEST_PATH_IMAGE039
Wherein->
Figure 958362DEST_PATH_IMAGE036
For the most ordered reconstruction of the decoded pictures +.>
Figure 866275DEST_PATH_IMAGE008
For flutter blurred image +.>
Figure 814639DEST_PATH_IMAGE040
Is a fuzzy kernel.
Further, the step one is that the flutter blurred image is collected by setting exposure time
Figure 782595DEST_PATH_IMAGE041
Exposure time +.>
Figure 878727DEST_PATH_IMAGE041
Average divide into->
Figure 528014DEST_PATH_IMAGE001
Time slots, & gt>
Figure 658781DEST_PATH_IMAGE001
For the number of time slots (i.e. the code length), each time slot time is +.>
Figure 729463DEST_PATH_IMAGE042
Each time slot is +.>
Figure 47311DEST_PATH_IMAGE043
Collecting charges, outputting independently, and then superposing to form a flutter blurred image, wherein the total time required for collecting the image is
Figure 562606DEST_PATH_IMAGE044
Further, the first step sets an exposure time
Figure 485563DEST_PATH_IMAGE041
Exposure time +.>
Figure 795322DEST_PATH_IMAGE041
Average divide into->
Figure 538150DEST_PATH_IMAGE001
Time slots, each time slot time is +>
Figure 591556DEST_PATH_IMAGE042
,/>
Figure 369019DEST_PATH_IMAGE001
The transmission process is driven out only once for the number of time slots, i.e
Figure 849679DEST_PATH_IMAGE045
Further, the second step is to preset the length to be
Figure 407700DEST_PATH_IMAGE001
Binary coding sequence->
Figure 936901DEST_PATH_IMAGE002
And blur Length->
Figure 896767DEST_PATH_IMAGE046
Construction of fuzzy core->
Figure 486011DEST_PATH_IMAGE040
The method comprises changing image convolution operation into product operation, adding zero vector +.>
Figure 265748DEST_PATH_IMAGE047
An improved mathematical model of the relative displacement of the acquired image in the spatial domain is built, which corresponds to the blurred pixel length of the actual target projected into the image plane>
Figure 660958DEST_PATH_IMAGE048
Greater than the pre-encoded length->
Figure 413013DEST_PATH_IMAGE001
It is necessary to add/subtract the code sequence>
Figure 235475DEST_PATH_IMAGE002
Then a plurality of 0 s are added to make the length equal to the length of the blurred pixel, so as to form a length of +.>
Figure 440192DEST_PATH_IMAGE049
Is->
Figure 639092DEST_PATH_IMAGE008
In the formula->
Figure 245654DEST_PATH_IMAGE003
For the length of ambiguity after zero padding, +.>
Figure 239018DEST_PATH_IMAGE048
For blurring image +.>
Figure 931030DEST_PATH_IMAGE008
In the direction of movementA length; at this time, the blur kernel for reconstruction +.>
Figure 933621DEST_PATH_IMAGE040
According to the blurring length after zero padding +.>
Figure 457006DEST_PATH_IMAGE003
The Toeplitz matrix is formed by single direction shift.
On the other hand, the invention also provides a device for collecting and recovering the flutter blurred image, which uses a method for recovering the flutter blurred image to collect and recover the flutter blurred image, and comprises the following steps: the device comprises an image acquisition module, a core control module, a communication module, a storage module, a clock module, a transmission and display module; the image acquisition module consists of a CCD image sensor, a time sequence driving circuit, a signal conditioning and converting circuit and the like; the time sequence driving circuit is divided into a horizontal time sequence driving module and a vertical time sequence driving module, and the time sequence driving circuit provides proper driving level for the CCD image sensor; the signal conditioning and converting circuit comprises analog signal amplification filtering and analog-to-digital conversion; the pre-coding sends an instruction to the image acquisition module through the core control module to enable the image acquisition module to acquire a flutter blurred image of the motion; the image data and the prefabricated codes are stored in the storage module, and the whole circuit is ensured by the clock module and the communication module; the CCD image sensor includes a horizontal shift register and a vertical shift register.
Further, the core control module comprises a time sequence module, a decoding reconstruction module, an image reading module, a DDR control module, an FPGA gigabit network core, an exposure code and a reference clock.
Further, the driving level is a three-state level.
Further, the image acquisition module further comprises a laser generator and a laser receiver, and the CCD image sensor receives laser emitted by the laser generator through the laser receiver to perform automatic focusing so as to solve the problem of image blurring caused by inaccurate focusing.
Further, the image acquisition module further comprises an ultrasonic generator and an ultrasonic receiver, and the CCD image sensor receives ultrasonic waves sent by the ultrasonic generator through the ultrasonic receiver to perform automatic focusing so as to solve the problem of image blurring caused by focusing inaccuracy.
Further, the image acquisition module further comprises an infrared generator and an infrared receiver, and the CCD image sensor receives infrared light emitted by the infrared generator through the infrared receiver to perform automatic focusing so as to solve the problem of image blurring caused by inaccurate focusing.
Furthermore, the acquisition and restoration device for the flutter blurred image further comprises a three-dimensional model and a visualization module, wherein the image acquisition module is connected with the three-dimensional model and the visualization module, and the three-dimensional model and the visualization module are used for carrying out dynamic display in a three-dimensional visualization mode.
Further, the core control module further includes a target tracking module that designates a specific dither blur image as a target image when a plurality of dither blur images are acquired.
On the other hand, the invention also provides a monitoring system, which uses a plurality of the acquisition and restoration devices of the flutter blurred images, the monitoring system also comprises a remote data transmission module and a blockchain module, the blockchain module is divided into a monitoring node and a management node, the remote data transmission module is respectively installed on the plurality of the acquisition and restoration devices of the flutter blurred images to form a plurality of monitoring nodes, and the monitoring nodes reconstruct the decoding images generated by the nodes in the most order
Figure 558955DEST_PATH_IMAGE050
Timestamping, generating simultaneously a pair of public and private keys, said most ordered reconstructed decoded image to be timestamped +.>
Figure 800580DEST_PATH_IMAGE050
Uploading to the blockchain module, performing digital signature by using a private key, and transmitting a public key to the management node; the management node is a central control system provided with the remote data transmission module, and is connected withReceiving the public key of the monitoring node to verify the private key, and checking the most orderly reconstructed decoded image +.>
Figure 278966DEST_PATH_IMAGE050
The beneficial effects of the invention are as follows:
1. the invention discloses a coding exposure imaging restoration method for generating blurring by flutter of a target relative to a camera in a single motion direction.
2. The invention starts from the aim of completely protecting the target image information, and provides a method for removing the blurring of a moving target from two angles of imaging and image post-processing. Modulating the exposure process by utilizing a pre-coding mode, protecting the original high-frequency information in the acquired image, and decoding the acquired image in the subsequent image restoration process to obtain a clear image.
3. The method omits manual selection intervention or the use time of external measurement equipment, and breaks the image restriction of natural image rules. The complexity of equipment used in the prior art is obviously reduced by estimating through a combined image restoration algorithm of image structure similarity and entropy.
4. The image acquisition and restoration device disclosed by the invention can solve the restoration problem of motion blurred images, realize clear restored images of cameras with higher frame rates with lower cost, and utilize the non-reference image evaluation index test, and compared with general exposure, the average value of the improvement of the related indexes is about 2 times of the related index value of the general common exposure.
5. The monitoring system of the invention uses the remote data transmission and control module to uplink the image information to the block chain module, and uses the advantages of the block chain technology such as decentralization, encryption, non-falsification and the like, so that the monitoring system can be widely applied to the industries such as traffic monitoring and the like.
Drawings
In order to better express the technical scheme of the invention, the following description of the invention is given by way of the accompanying drawings:
fig. 1 is a schematic diagram of an image structure similarity structure according to a first embodiment;
FIG. 2 is a fuzzy image fuzzy kernel construction diagram of a pre-coding exposure of the first embodiment;
FIG. 3 is a diagram showing the relationship between one-dimensional relative motion and convolution calculation according to the first embodiment;
fig. 4 is a diagram of a process of acquisition and reconstruction decoding of a flutter blurred image according to the first embodiment;
fig. 5 is a block diagram of a flutter blurred image acquisition and restoration device according to the third embodiment;
FIG. 6 is a flow chart of a codec procedure of a three-dither blurred image according to an embodiment;
FIG. 7 is a graph showing the experimental results of the present invention;
reference numerals illustrate: 1. the device comprises an image acquisition module 11, a CCD image sensor 12, a time sequence driving circuit 13, a signal conditioning and converting circuit 2, a core control module 21, a time sequence module 22, a decoding reconstruction module 23, an image reading module 24, a DDR control module 25, an FPGA gigabit network core module 26, an exposure coding module 27, a reference clock 3, a communication module 4, a storage module 5 and a clock module.
Detailed Description
Preferred embodiments of the present invention are described below with reference to the accompanying drawings. It should be understood by those skilled in the art that these embodiments are merely for explaining the technical principles of the present invention, and are not intended to limit the scope of the present invention.
It should be noted that, in the description of the present invention, terms such as "upper," "lower," "left," "right," "inner," "outer," and the like indicate directions or positional relationships based on the directions or positional relationships shown in the drawings, which are merely for convenience of description, and do not indicate or imply that the apparatus or elements must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the present invention.
Furthermore, it should be noted that, in the description of the present invention, unless explicitly specified and limited otherwise, the terms "mounted," "connected," and "connected" are to be construed broadly, and may be either fixedly connected, detachably connected, or integrally connected, for example; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present invention can be understood by those skilled in the art according to the specific circumstances.
The working principle of the invention is as follows: and a Toeplitz matrix is utilized to form a fuzzy core, and the problem of image restoration of a target and a camera in a single motion direction is solved according to the construction form of the fuzzy core.
The coded exposure will be for a complete exposure time
Figure 656858DEST_PATH_IMAGE051
Is divided into->
Figure 257603DEST_PATH_IMAGE001
The number of the time slots is consistent with the preset code length, and whether each time slot is exposed is consistent with the code word of the corresponding bit. Taking a one-dimensional signal as an example, if the acquisition length is +.>
Figure 658629DEST_PATH_IMAGE048
Is>
Figure 3022DEST_PATH_IMAGE052
The signal and length are +>
Figure 173104DEST_PATH_IMAGE001
Is>
Figure 944751DEST_PATH_IMAGE053
The convolution operation (namely binary code sequence) can be expressed as
Figure 833072DEST_PATH_IMAGE054
/>
From the relative relation of signals, the above is relative to one-dimensional signals
Figure 715578DEST_PATH_IMAGE055
Movement in one direction causes data misalignmentAnd (5) superposition. The method can be understood as using Toeplitz matrix to divide the clear one-dimensional signal +.>
Figure 68062DEST_PATH_IMAGE055
And (2) fuzzy core->
Figure 948293DEST_PATH_IMAGE002
The matrix convolution relationship between the two is converted into a product relationship of the matrix. If->
Figure 386227DEST_PATH_IMAGE002
Any one of the data is +.>
Figure 4247DEST_PATH_IMAGE056
The method is continuous data superposition in the traditional sense; while->
Figure 148921DEST_PATH_IMAGE002
The data in (a) is->
Figure 262370DEST_PATH_IMAGE056
Or->
Figure 390863DEST_PATH_IMAGE057
When represented as discontinuous or intermittent data superposition, where
Figure 615171DEST_PATH_IMAGE058
From the viewpoint of shutter exposure, the code exposure differs from the normal exposure in that the shutter is not kept in a one-time open state, but is controlled to open and close by a preset code at a certain frequency. If you get
Figure 942248DEST_PATH_IMAGE056
Indicating the slot exposure; otherwise, if->
Figure 633123DEST_PATH_IMAGE057
Indicating that the time slot is not exposed.
As established in fig. 3 is a static binary code sequenceColumn of
Figure 45650DEST_PATH_IMAGE002
One-dimensional clear signal shifted from time-sharing>
Figure 11332DEST_PATH_IMAGE055
Is a mathematical model of the acquisition of (a). Taking the one-dimensional relative motion of the imaging system and the target object as an example, the encoding exposure imaging and image restoration process will be described. If blur image +.>
Figure 192914DEST_PATH_IMAGE008
Length in the direction of movement is +.>
Figure 648167DEST_PATH_IMAGE059
Is>
Figure 220093DEST_PATH_IMAGE055
Exposing by an imaging system, whether the imaging system exposes and has a length of +>
Figure 786204DEST_PATH_IMAGE001
Is a preset binary coding sequence +.>
Figure 25555DEST_PATH_IMAGE002
And keep the same. Will +_in FIG. 3>
Figure 651709DEST_PATH_IMAGE001
Is set to +.>
Figure 773248DEST_PATH_IMAGE060
. If the two signals move relatively, the acquisition length is equal to
Figure 346312DEST_PATH_IMAGE061
Signal of->
Figure 502487DEST_PATH_IMAGE008
The signal is a convolution model between the two. When->
Figure 971646DEST_PATH_IMAGE056
At +.>
Figure 580482DEST_PATH_IMAGE008
Middle is accumulated with time->
Figure 426078DEST_PATH_IMAGE055
Information in (a); while->
Figure 436759DEST_PATH_IMAGE057
When (I)>
Figure 607978DEST_PATH_IMAGE008
Without any newly added information. Thus, the acquisition signal +.>
Figure 704109DEST_PATH_IMAGE008
Whether or not to accumulate and +.>
Figure 87817DEST_PATH_IMAGE062
Whether it is 1. In the restoration process, coding by time-sharing dislocation ∈>
Figure 953005DEST_PATH_IMAGE002
Restoring the relatively stationary one-dimensional signal +.>
Figure 295125DEST_PATH_IMAGE055
. The acquisition encoding exposure and image reconstruction decoding process is shown in fig. 4.
The same exposure code and different blur lengths can form different blur kernels, and the quality of the restored reconstructed decoded image is also different. In order to find a clear restored image from a plurality of reconstructed decoded images, the structural similarity of the image quality evaluation function is taken as a comparison basis, so that high-quality reconstruction of the image is formed, and a system block diagram of the high-quality reconstruction is shown in figure 1.
In the context of figure 1 of the drawings,
Figure 878553DEST_PATH_IMAGE008
flutter blurred image indicative of degradation caused by flutter, < >>
Figure 393848DEST_PATH_IMAGE063
Representing the original sharp image. Similar to human visual sense, the structural similarity will be from luminance +.>
Figure 785646DEST_PATH_IMAGE064
Contrast->
Figure 95405DEST_PATH_IMAGE065
Structural->
Figure 103812DEST_PATH_IMAGE066
And independently judging the structural similarity of the two images at three angles. These three quantities can be expressed as
Figure 157219DEST_PATH_IMAGE067
、/>
Figure 200261DEST_PATH_IMAGE068
、/>
Figure 680921DEST_PATH_IMAGE069
There is a compound of the formula (I),
Figure 238941DEST_PATH_IMAGE070
,/>
Figure 520142DEST_PATH_IMAGE071
,/>
Figure 214428DEST_PATH_IMAGE072
Figure 131569DEST_PATH_IMAGE073
,/>
Figure 114568DEST_PATH_IMAGE074
,/>
Figure 509777DEST_PATH_IMAGE075
、/>
Figure 261833DEST_PATH_IMAGE076
for blurring image +.>
Figure 84295DEST_PATH_IMAGE008
And original sharp image->
Figure 351328DEST_PATH_IMAGE063
An average value of the pixels; />
Figure 222332DEST_PATH_IMAGE077
For blurring image +.>
Figure 891211DEST_PATH_IMAGE008
And original sharp image->
Figure 884575DEST_PATH_IMAGE063
Standard deviation of each; />
Figure 842167DEST_PATH_IMAGE078
For blurring image +.>
Figure 579179DEST_PATH_IMAGE008
And original sharp image->
Figure 774668DEST_PATH_IMAGE063
Is a cross-correlation function of (2); to avoid the situation that the molecular denominator is zero, three very small positive numbers are defined, e.g.>
Figure 204512DEST_PATH_IMAGE079
、/>
Figure 446137DEST_PATH_IMAGE080
、/>
Figure 924523DEST_PATH_IMAGE081
. In this embodiment, <' > a->
Figure 567994DEST_PATH_IMAGE082
For a range of pixel values of a grey image, e.g. 8 bits for a grey image, i.e +.>
Figure 840844DEST_PATH_IMAGE083
。/>
Figure 569765DEST_PATH_IMAGE084
、/>
Figure 648580DEST_PATH_IMAGE085
Is a default value.
The image is blurred by the flutter
Figure 818661DEST_PATH_IMAGE008
And original sharp image->
Figure 590308DEST_PATH_IMAGE063
The similarity function between can be expressed as:
Figure 540947DEST_PATH_IMAGE086
in this case, parameters
Figure DEST_PATH_IMAGE087
For adjusting->
Figure 626714DEST_PATH_IMAGE064
、/>
Figure 979198DEST_PATH_IMAGE065
、/>
Figure 859429DEST_PATH_IMAGE066
The ratio of the three parameters is as
Figure 172730DEST_PATH_IMAGE088
There is->
Figure 593347DEST_PATH_IMAGE089
The formula shows the original sharp image +.>
Figure 65917DEST_PATH_IMAGE063
And flutter blurred image +.>
Figure 851470DEST_PATH_IMAGE008
Structural similarity of (3). />
Figure 776701DEST_PATH_IMAGE090
The value range of (2) is +.>
Figure 204271DEST_PATH_IMAGE091
。/>
Figure 265768DEST_PATH_IMAGE090
The larger the value of (c), the more similar the two image qualities are.
However, in practical experiments, clear images cannot be obtained, and flutter blurred images
Figure 550119DEST_PATH_IMAGE008
Is the only data that can be obtained, and the corresponding reconstructed image can be obtained by different blur lengths. Because the acquired code exposure blurred image is derived from the target object, similarity evaluation can be performed by reconstructing the image and the code exposure blurred image. Therefore, the dither blurred image +.>
Figure 900329DEST_PATH_IMAGE006
And encoding exposure decoding image +.>
Figure 928328DEST_PATH_IMAGE007
The structural similarity index of (2) is defined as +.>
Figure 41734DEST_PATH_IMAGE092
. However, the dither blurred image +.>
Figure 231407DEST_PATH_IMAGE006
Is a degraded image, with which the decoded image is evaluated +.>
Figure 6596DEST_PATH_IMAGE007
Deviations in the structural similarity index are caused, and therefore, the order of the images needs to be calculated to finally determine the restored image.
The structural similarity method can only find the reconstructed image which is most similar to the coded exposure blurred image, but the image is not necessarily the most ordered image which accords with the natural statistics rule. The original clear image should be an orderly natural image, so that the restored reconstructed image is orderly used as an index for judging the image quality. To avoid deviations from structural similarity alone, information entropy is introduced. Entropy is a measure representing the amount of information. And determining the image range most similar to the blurred image by combining the structural similarity, and then searching the most ordered image information to determine the final restored image.
When a system is more ordered, the information entropy of the system is lower; conversely, the more chaotic the system, the higher the entropy of the system. The order degree of gray values among image pixels in the search interval is calculated by using the spatial entropy; meanwhile, in order to evaluate the flat characteristic of the reconstructed decoded signal, the signal is transformed into the frequency domain using discrete cosine transform and its spectral entropy is calculated. And determining the optimal fuzzy length by adjusting the weight estimation occupied by the spatial entropy and the spectral entropy, and reconstructing and restoring the clear image. The ordering degree formula of the spatial entropy representation image is as follows:
Figure 572707DEST_PATH_IMAGE093
wherein->
Figure 812058DEST_PATH_IMAGE094
Representing the maximum value in the gray scale range of the image; />
Figure 172632DEST_PATH_IMAGE095
Representing gray scale +.>
Figure 559751DEST_PATH_IMAGE096
Is a possibility of (1). The one-dimensional image entropy cannot reflect the spatial characteristics of the gray distribution of the image, and two-dimensional entropy of the image needs to be introduced, wherein the spatial distribution formula of the neighborhood gray mean value of the characteristic image is as follows:
Figure 867236DEST_PATH_IMAGE097
in the formula->
Figure 757832DEST_PATH_IMAGE098
Wherein->
Figure 554886DEST_PATH_IMAGE099
Is a characteristic binary group->
Figure 101405DEST_PATH_IMAGE100
Frequency of occurrence, < >>
Figure 274898DEST_PATH_IMAGE096
Representing the gray value of a pixel +.>
Figure 223262DEST_PATH_IMAGE101
Representing a neighborhood gray level average; />
Figure 191218DEST_PATH_IMAGE102
Is the scale of the image.
To obtain an accurate solution, spectral entropy is introduced to detect the flatness of the signal spectrum. When the signal has strong correlation in the spatial domain, the transformation into the frequency domain appears as a concentrated convergence of a specific region, where the spatial domain signal is transformed into the frequency domain using a Discrete Cosine Transform (DCT). Since the image pixel value is a real number, the DCT is also a real number operation, and the operation speed is faster than the complex number operation in the Fourier transform. If the signal is
Figure 287350DEST_PATH_IMAGE103
DCT transform to +.>
Figure 936637DEST_PATH_IMAGE104
,/>
Figure 801825DEST_PATH_IMAGE105
Is a generalized frequency domain variable, the normalized DCT coefficient is:/>
Figure 878365DEST_PATH_IMAGE106
Its spectral entropy can be expressed as +.>
Figure 461793DEST_PATH_IMAGE107
1-4, a method for restoring a flutter blurred image:
step one: the acquisition of the flutter blurred image is completed by using the CCD image sensor 11, and the exposure time is set
Figure 649192DEST_PATH_IMAGE108
Exposure time +.>
Figure 634466DEST_PATH_IMAGE108
Average divide into->
Figure 944224DEST_PATH_IMAGE109
Time slots, each time slot time is +>
Figure 952632DEST_PATH_IMAGE110
,/>
Figure 6038DEST_PATH_IMAGE109
For the number of time slots, each time slot +.>
Figure 845818DEST_PATH_IMAGE111
Collecting charges, outputting independently, and then superposing to form a flutter blurred image, wherein the total time required for collecting the image is
Figure 264161DEST_PATH_IMAGE112
Step two: by presetting the code length to be
Figure 556603DEST_PATH_IMAGE001
Binary coding sequence->
Figure 351383DEST_PATH_IMAGE002
And blur Length->
Figure 45670DEST_PATH_IMAGE046
Construction of fuzzy core->
Figure 697231DEST_PATH_IMAGE040
Said->
Figure 680230DEST_PATH_IMAGE109
For binary coding sequences->
Figure 75440DEST_PATH_IMAGE002
Is defined by the formula->
Figure 624233DEST_PATH_IMAGE113
Calculating to obtain a reconstructed image, wherein->
Figure 649958DEST_PATH_IMAGE008
For flutter blurred image +.>
Figure 916991DEST_PATH_IMAGE063
To decode a sharp image;
as shown in fig. 2, by presetting the length to be
Figure 53574DEST_PATH_IMAGE001
Binary coding sequence->
Figure 456874DEST_PATH_IMAGE002
And blur Length->
Figure 715817DEST_PATH_IMAGE046
Construction of fuzzy core->
Figure 401970DEST_PATH_IMAGE040
The method is that the image convolution operation is changed into a product operationFuzzy core +.>
Figure 138982DEST_PATH_IMAGE040
The shift to the right and down in the main diagonal forms a Toeplitz-like matrix, as shown in fig. 2 (a), where the blur length of the object projected onto the image plane is equal to the code length. It is in fact difficult to ensure that the blur length of the object projected onto the image plane is equal to the encoding length. If the length of the blurred pixels of the actual target projected into the image plane is large, a clear restored image cannot be obtained by using the model.
In order not to change the original one-dimensional signal
Figure 600050DEST_PATH_IMAGE052
Frequency information of (2) in the original one-dimensional signal +.>
Figure 29894DEST_PATH_IMAGE052
The tail is added with zero vector of specific length +.>
Figure 5940DEST_PATH_IMAGE047
An improved mathematical model of the relative displacement of the acquired image in the spatial domain is built, as shown in fig. 2 (b). This situation corresponds to a blurred pixel length of the actual target projected into the image plane being larger than the pre-coding length +.>
Figure 484326DEST_PATH_IMAGE001
After encoding, a plurality of 0 s are added to be equal to the length of the blurred pixel, so that the length is +.>
Figure 127797DEST_PATH_IMAGE114
Is->
Figure 400647DEST_PATH_IMAGE008
In the formula->
Figure 129568DEST_PATH_IMAGE046
For the length of ambiguity after zero padding, +.>
Figure 208383DEST_PATH_IMAGE048
For blurring image +.>
Figure 644043DEST_PATH_IMAGE008
Length in the direction of movement; at this time, the blur kernel for reconstruction +.>
Figure 150111DEST_PATH_IMAGE040
According to the blurring length after zero padding +.>
Figure 366329DEST_PATH_IMAGE046
The shift constitutes a Toeplitz matrix form.
One-dimensional signal
Figure 186517DEST_PATH_IMAGE002
And->
Figure 539001DEST_PATH_IMAGE052
The process of the relative shift superposition is equivalent to the convolution process of signals. When there is a flutter blurred image +.>
Figure 419232DEST_PATH_IMAGE008
Length in the direction of movement is +.>
Figure 857167DEST_PATH_IMAGE048
Is>
Figure 277784DEST_PATH_IMAGE052
When generating a displacement, a one-dimensional signal with a length m is +.>
Figure 688037DEST_PATH_IMAGE115
The convolution process will produce a length of +>
Figure 535907DEST_PATH_IMAGE116
Is->
Figure 398821DEST_PATH_IMAGE006
. If at this time->
Figure 623129DEST_PATH_IMAGE115
For binary coding sequence, one-dimensional signal is selectively superimposed according to coding rule>
Figure 950205DEST_PATH_IMAGE052
Is provided. Thus, the imaging process of the coded exposure is similar to the front-to-back accumulated smear at the time of ordinary exposure, and is the accumulation of the exposed image in a plurality of coded time slots.
Step three: calculating a flutter blurred image
Figure 172239DEST_PATH_IMAGE117
And all decoded sharp pictures +.>
Figure 584766DEST_PATH_IMAGE118
Is of the structural similarity index of (2)
Figure 347185DEST_PATH_IMAGE119
Length of blur at the time of searching for maximum +.>
Figure 466451DEST_PATH_IMAGE120
The formula is: />
Figure 921703DEST_PATH_IMAGE121
In the formula->
Figure 759209DEST_PATH_IMAGE122
Search range for fuzzy length;
step four: determining an image information entropy value in a certain range near an image with the maximum structural similarity, wherein a search interval is
Figure 325320DEST_PATH_IMAGE123
Wherein->
Figure 361409DEST_PATH_IMAGE124
Is positive integer->
Figure 925245DEST_PATH_IMAGE125
,/>
Figure 312364DEST_PATH_IMAGE126
Is made of->
Figure 354270DEST_PATH_IMAGE127
The total blur length limit range is determined to be +.>
Figure 510445DEST_PATH_IMAGE126
Within the limits according to the formula ∈ ->
Figure 307499DEST_PATH_IMAGE128
And
Figure 854018DEST_PATH_IMAGE129
calculating spatial entropy of each decoded image>
Figure 27511DEST_PATH_IMAGE130
And spectral entropy after image DCT (discrete cosine transform)>
Figure 38192DEST_PATH_IMAGE131
Wherein->
Figure 943831DEST_PATH_IMAGE132
Wherein->
Figure 39963DEST_PATH_IMAGE133
Is a characteristic binary group->
Figure 417812DEST_PATH_IMAGE134
Frequency of occurrence, < >>
Figure 548579DEST_PATH_IMAGE025
Representing the gray value of a pixel +.>
Figure 687436DEST_PATH_IMAGE026
Representing neighborhood gray scale mean,/->
Figure 208547DEST_PATH_IMAGE135
Is the scale of the image; if signal->
Figure 458263DEST_PATH_IMAGE136
DCT transformation into (1)
Figure 443536DEST_PATH_IMAGE137
,/>
Figure 690978DEST_PATH_IMAGE138
Is a generalized frequency domain variable, normalized DCT transform coefficient is +.>
Figure 761702DEST_PATH_IMAGE139
Step five: formula (VI)
Figure 752792DEST_PATH_IMAGE140
Calculating the information entropy value of the reconstructed image>
Figure 592572DEST_PATH_IMAGE141
Wherein->
Figure 73232DEST_PATH_IMAGE142
For the picture space entropy value->
Figure 303356DEST_PATH_IMAGE130
Weight of occupied->
Figure 160454DEST_PATH_IMAGE143
Spectral entropy values for the coefficients of the DCT-transform (discrete cosine transform) of an image>
Figure 854740DEST_PATH_IMAGE131
The weight is occupied;
step six: taking out
Figure 443985DEST_PATH_IMAGE141
Minimum value, determined as most ordered reconstructed decoded picture +.>
Figure 489301DEST_PATH_IMAGE144
Is +.>
Figure 822193DEST_PATH_IMAGE145
The formula is->
Figure 370986DEST_PATH_IMAGE146
;/>
Step seven: the image is restored clearly by adopting an inverse filtering method, and the formula is
Figure 459028DEST_PATH_IMAGE147
Wherein->
Figure 663745DEST_PATH_IMAGE144
For the most ordered reconstruction of the decoded pictures +.>
Figure 597065DEST_PATH_IMAGE117
For flutter blurred image +.>
Figure 203627DEST_PATH_IMAGE148
Is a fuzzy kernel.
The second embodiment differs from the first embodiment in that the exposure time is set in the first step
Figure 462570DEST_PATH_IMAGE041
Exposure time +.>
Figure 951320DEST_PATH_IMAGE041
Average divide into->
Figure 891595DEST_PATH_IMAGE149
Time slots, each time slot time is +>
Figure 414980DEST_PATH_IMAGE150
,/>
Figure 516928DEST_PATH_IMAGE149
The transmission process is driven out only once for the number of time slots, i.e. +.>
Figure 758554DEST_PATH_IMAGE151
The method uses only one charge to drive and transfer time
Figure 299256DEST_PATH_IMAGE152
And meanwhile, the occurrence probability of electronic noise is reduced. The total exposure time of the method is kept consistent with the common exposure time. I.e. under equal conditions the pixel displacement of the object and camera projection in the image plane coincides with a normal camera. The method is characterized in that a motion blur superposition process is arranged in front, the motion process is converted into a charge accumulation process conforming to a preset coding rule, and in the process, because the sum of time slots is the total exposure time, the image acquisition interval time is the smallest in the scheme, namely the minimum blur length is obtained under the same condition, so that the clear image can be restored and reconstructed conveniently.
Embodiment three, a device for acquiring and recovering a flutter blurred image, is shown in fig. 5-6. As shown in fig. 5, the device consists of an image acquisition module 1, a core control module 2, a communication module 3, a storage module 4, a clock module 5, a three-dimensional model and visualization module, and a transmission and display module (not shown in the figure); the image acquisition module 1 consists of a CCD image sensor 11, a time sequence driving circuit 12, a signal conditioning and converting circuit 13 and the like; the time sequence driving circuit 2 is divided into a horizontal time sequence driving module and a vertical time sequence driving module, and provides proper driving level for the CCD image sensor 11, wherein the driving level is a three-state level; the signal conditioning and converting circuit 13 includes analog signal amplification filtering and analog-to-digital conversion. The CCD image sensor 11 includes a horizontal shift register and a vertical shift register (not shown). The core control module 2 comprises a time sequence module 21, a decoding reconstruction module 22, an image reading module 23, a DDR control module 24, an FPGA gigabit network core module 25, an exposure coding module 26 and a reference clock 27; the image acquisition module 1 dynamically displays in a three-dimensional visualization mode through a three-dimensional model and a visualization module. The pre-coding sends an instruction to the image acquisition module 1 through the core control module 2 to enable the image acquisition module to acquire a flutter blurred image of the motion; the image data and the prefabricated codes are stored in the storage module 4, and the whole circuit is ensured by the clock module 5 and the communication module 3.
To solve the image blur caused by the out-of-focus, the image acquisition module 1 may use a generator and a receiver, through which the CCD image sensor 11 receives the signal from the generator for auto-focusing. The generator and receiver may use commercially available types of products, not limited to laser generators and receivers, ultrasonic generators and receivers, infrared generators and receivers.
The core control module 2 can be finished by XC6SLX45T-3FG484C in Spartan6 series manufactured by Xilinx company; the reference clock 27 is selected from
Figure 614831DEST_PATH_IMAGE153
A crystal oscillator; the driving time sequence required by the CCD image sensor 11 is generated by a time sequence module 21, a time sequence driving signal is generated by a core processing module (a field programmable logic device FPGA is used in the embodiment of the invention), and a special chip CXD3400N is selected because two levels and three levels are required to be ensured together for driving the CCD image sensor 11. The image reading module 23 reads the image data and saves the data in the DDR data memory 24 through the DDR management core of the FPGA. The decoding and reconstructing module 22 in the core control module 2 is reserved, and can reconstruct the clear restored image by intra-chip decoding by using a relevant restoring method after reading the flutter blurred image in the DDR data memory 24. The FPGA gigabit network core module 25 adopted for data communication with the upper computer is a high-speed bidirectional transmission channel, and reliably transmits image data and status commands. The physical layer transceiver chip adopts a physical layer gigabit Ethernet transceiver 88E1111; the CCD image sensor 11 adopts a high-sensitivity and low-noise ICX204AL black-and-white image sensor; the signal conditioning and converting circuit 13 uses a chip AD9949. Other modules are all commercially available universal devices.
As shown in fig. 6, the dither blurred image codec program flow chart. Preset binary coding sequence
Figure 215577DEST_PATH_IMAGE053
When a control signal for shutter change is supplied to the substrate SUB,is also used as a blur kernel in the decoding module>
Figure 678919DEST_PATH_IMAGE004
Is a structure of (a). Therefore, when the dither matrix image conforming to the preset coding +.>
Figure 960996DEST_PATH_IMAGE006
After being acquired, the decoding module uses the inverse filtering method +.>
Figure 458973DEST_PATH_IMAGE154
And obtaining a decoded clear image. Meanwhile, the decoding module is used as an optional module in the embedded coding camera, and the blurred image can be restored by being transmitted to an upper computer for processing and restoring or by downloading a complex reconstruction method.
The fourth embodiment is a monitoring system, which uses a plurality of vibration blurred image acquisition and restoration devices, the monitoring system further comprises a remote data transmission module and a blockchain module, the blockchain module is divided into a monitoring node and a management node, the plurality of vibration blurred image acquisition and restoration devices are respectively provided with the remote data transmission module to form a plurality of monitoring nodes, and the monitoring nodes reconstruct the most orderly decoded image generated by the monitoring nodes
Figure 902724DEST_PATH_IMAGE155
Timestamping, generating simultaneously a pair of public and private keys, said most ordered reconstructed decoded image to be timestamped +.>
Figure 118942DEST_PATH_IMAGE155
Uploading to the blockchain module, performing digital signature by using a private key, and transmitting a public key to the management node; the management node is a central control system provided with the remote data transmission module, receives the public key of the monitoring node to verify the private key, and checks the most orderly reconstructed decoded image
Figure 1447DEST_PATH_IMAGE155
Fig. 7 is a graph comparing experimental results of other methods using the prior art with experimental results of using the method and apparatus of the present invention.
The above examples illustrate only a few embodiments of the invention, which are described in detail and are not to be construed as limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention.

Claims (10)

1. A method for restoring a flutter blurred image is characterized by comprising the following steps of
Step one: collecting flutter blurred images;
step two: by presetting the code length to be
Figure 506860DEST_PATH_IMAGE001
Binary coding sequence->
Figure 27971DEST_PATH_IMAGE002
And blur Length->
Figure 277687DEST_PATH_IMAGE003
Construction of fuzzy core->
Figure 200643DEST_PATH_IMAGE004
The said
Figure 510402DEST_PATH_IMAGE001
For binary coding sequences->
Figure 581126DEST_PATH_IMAGE002
Is defined by the formula->
Figure 572216DEST_PATH_IMAGE005
Calculating to obtain a reconstructed image, wherein->
Figure 411996DEST_PATH_IMAGE006
For flutter blurred image +.>
Figure 830339DEST_PATH_IMAGE007
To decode an image;
step three: calculating a flutter blurred image
Figure 122780DEST_PATH_IMAGE006
And all decoded pictures +.>
Figure 979877DEST_PATH_IMAGE007
Structural similarity index->
Figure 611847DEST_PATH_IMAGE008
Length of blur at the time of searching for maximum +.>
Figure 263408DEST_PATH_IMAGE009
The formula is: />
Figure 246408DEST_PATH_IMAGE010
In the formula->
Figure 641617DEST_PATH_IMAGE011
Search range for fuzzy length;
step four: determining an image information entropy value in a certain range near an image with the maximum structural similarity, wherein a search interval is
Figure 128093DEST_PATH_IMAGE012
Wherein->
Figure 216135DEST_PATH_IMAGE013
Is positive integer->
Figure 483168DEST_PATH_IMAGE014
,/>
Figure 619751DEST_PATH_IMAGE015
Is made of->
Figure 23051DEST_PATH_IMAGE016
The total blur length limit range is determined to be +.>
Figure 219677DEST_PATH_IMAGE015
Within the limits according to the formula ∈ ->
Figure 974006DEST_PATH_IMAGE017
And
Figure 648701DEST_PATH_IMAGE018
calculating spatial entropy of each decoded image>
Figure 172087DEST_PATH_IMAGE019
And spectral entropy after image DCT (discrete cosine transform)>
Figure 601931DEST_PATH_IMAGE020
Wherein->
Figure 509801DEST_PATH_IMAGE021
Wherein->
Figure 50504DEST_PATH_IMAGE022
Is a characteristic binary group->
Figure 631658DEST_PATH_IMAGE023
Frequency of occurrence, < >>
Figure 966824DEST_PATH_IMAGE024
Representing the gray value of a pixel +.>
Figure 633429DEST_PATH_IMAGE025
Representing a neighborhoodGray scale mean value->
Figure 712243DEST_PATH_IMAGE026
Is the scale of the image; if signal->
Figure 210221DEST_PATH_IMAGE027
DCT transformation into (1)
Figure 653972DEST_PATH_IMAGE028
,/>
Figure 870189DEST_PATH_IMAGE029
Is a generalized frequency domain variable, normalized DCT transform coefficient is +.>
Figure 690378DEST_PATH_IMAGE030
Step five: formula (VI)
Figure 42862DEST_PATH_IMAGE031
Calculating the information entropy value of the reconstructed image>
Figure 985410DEST_PATH_IMAGE032
Wherein->
Figure 361028DEST_PATH_IMAGE033
For the picture space entropy value->
Figure 781645DEST_PATH_IMAGE019
Weight of occupied->
Figure 191897DEST_PATH_IMAGE034
Spectral entropy values for the coefficients of the DCT-transform (discrete cosine transform) of an image>
Figure 39768DEST_PATH_IMAGE020
The weight is occupied;
step six: taking out
Figure 964998DEST_PATH_IMAGE032
Minimum value, determined as most ordered reconstructed decoded picture +.>
Figure 126989DEST_PATH_IMAGE035
Is +.>
Figure 454066DEST_PATH_IMAGE036
The formula is->
Figure 676099DEST_PATH_IMAGE037
Step seven: the image is restored clearly by adopting an inverse filtering method, and the formula is
Figure 88626DEST_PATH_IMAGE038
Wherein->
Figure 851046DEST_PATH_IMAGE039
For the most ordered reconstruction of the decoded pictures +.>
Figure 970312DEST_PATH_IMAGE006
For flutter blurred image +.>
Figure 425564DEST_PATH_IMAGE004
Is a fuzzy kernel.
2. The method of claim 1, wherein the step one, the acquisition of the dither blurred image, sets the exposure time
Figure 263070DEST_PATH_IMAGE040
Exposure time +.>
Figure 829180DEST_PATH_IMAGE040
Average divide into->
Figure 802952DEST_PATH_IMAGE041
Time slots each of which is
Figure 429106DEST_PATH_IMAGE042
Each time slot is +.>
Figure 816225DEST_PATH_IMAGE043
After the collected charges are singly output, the collected charges are overlapped to form a flutter blurred image, and the total time required for collecting the images is +.>
Figure 858130DEST_PATH_IMAGE044
3. The method of recovering a dither blur image according to claim 1, wherein the first step sets an exposure time
Figure 14305DEST_PATH_IMAGE045
Exposure time +.>
Figure 14622DEST_PATH_IMAGE045
Average divide into->
Figure 623458DEST_PATH_IMAGE046
Time slots, each time slot time is +>
Figure 796950DEST_PATH_IMAGE047
The transmission of the driving out is done only once, i.e. +.>
Figure 479736DEST_PATH_IMAGE048
4. The method of claim 1, wherein the second step is performed by presetting the length to be
Figure 447692DEST_PATH_IMAGE001
Binary coding sequence->
Figure 481507DEST_PATH_IMAGE002
And blur Length->
Figure 927531DEST_PATH_IMAGE003
Construction of fuzzy core->
Figure 990122DEST_PATH_IMAGE004
The method comprises changing image convolution operation into product operation, adding zero vector +_of specific length at the tail of the original signal>
Figure 128980DEST_PATH_IMAGE049
Establishing an improved mathematical model of the relative displacement of the acquired image in the space domain, forming a mathematical model with a length of +.>
Figure 712408DEST_PATH_IMAGE050
Is->
Figure 899806DEST_PATH_IMAGE006
In the formula->
Figure 885080DEST_PATH_IMAGE003
For the length of ambiguity after zero padding, +.>
Figure 398101DEST_PATH_IMAGE051
For blurring image +.>
Figure 468825DEST_PATH_IMAGE006
Length in the direction of movement; at this time, the blur kernel for reconstruction +.>
Figure 522232DEST_PATH_IMAGE004
According to the blurring length after zero padding +.>
Figure 34116DEST_PATH_IMAGE003
The shift constitutes a Toeplitz matrix form.
5. A dither blurred image acquisition and restoration apparatus characterized in that the dither blurred image acquisition and restoration is realized using the dither blurred image restoration method of any one of claims 1 to 4, comprising: the device comprises an image acquisition module, a core control module, a communication module, a storage module, a clock module, a transmission and display module; the image acquisition module consists of a CCD image sensor, a time sequence driving circuit and a signal conditioning and converting circuit; the time sequence driving circuit is divided into a horizontal time sequence driving module and a vertical time sequence driving module, and provides proper driving level for the CCD image sensor; the signal conditioning and converting circuit comprises analog signal amplification filtering and analog-to-digital conversion; the CCD image sensor includes a horizontal shift register and a vertical shift register.
6. The device for acquiring and recovering the flutter blurred image according to claim 5, wherein the core control module comprises a time sequence module, a decoding reconstruction module, an image reading module, a DDR control module, an FPGA gigabit network core, an exposure code and a reference clock.
7. The apparatus of claim 5, wherein the driving level is a three-state level.
8. The device for acquiring and recovering the flutter blurred image of claim 5, further comprising a three-dimensional model and a visualization module, wherein the image acquisition module is connected with the three-dimensional model and the visualization module.
9. The apparatus of claim 5, wherein the core control module further comprises a target tracking module.
10. A monitoring system, characterized in that the flutter blurred image acquisition and restoration device as defined in any one of claims 5-9 is used, the monitoring system further comprises a remote data transmission module and a blockchain module, the blockchain module is divided into a monitoring node and a management node, the remote data transmission module is respectively installed on a plurality of the flutter blurred image acquisition and restoration devices to form a plurality of monitoring nodes, and the monitoring nodes reconstruct the most orderly decoded images generated by the nodes
Figure 514776DEST_PATH_IMAGE039
Timestamping, generating simultaneously a pair of public and private keys, said most ordered reconstructed decoded image to be timestamped +.>
Figure 744900DEST_PATH_IMAGE039
Uploading to the blockchain module, performing digital signature by using a private key, and transmitting a public key to the management node; the management node is a central control system provided with the remote data transmission module, receives the public key of the monitoring node to verify the private key, and views the most orderly reconstructed decoded image +.>
Figure 601997DEST_PATH_IMAGE039
。/>
CN202211315365.4A 2022-10-26 2022-10-26 Restoration method, acquisition restoration device and monitoring system for flutter blurred image Pending CN115994865A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211315365.4A CN115994865A (en) 2022-10-26 2022-10-26 Restoration method, acquisition restoration device and monitoring system for flutter blurred image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211315365.4A CN115994865A (en) 2022-10-26 2022-10-26 Restoration method, acquisition restoration device and monitoring system for flutter blurred image

Publications (1)

Publication Number Publication Date
CN115994865A true CN115994865A (en) 2023-04-21

Family

ID=85990993

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211315365.4A Pending CN115994865A (en) 2022-10-26 2022-10-26 Restoration method, acquisition restoration device and monitoring system for flutter blurred image

Country Status (1)

Country Link
CN (1) CN115994865A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117812275A (en) * 2024-02-28 2024-04-02 哈尔滨学院 Image optimization communication method for volleyball auxiliary training

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117812275A (en) * 2024-02-28 2024-04-02 哈尔滨学院 Image optimization communication method for volleyball auxiliary training
CN117812275B (en) * 2024-02-28 2024-05-28 哈尔滨学院 Image optimization communication method for volleyball auxiliary training

Similar Documents

Publication Publication Date Title
Kuznietsov et al. Semi-supervised deep learning for monocular depth map prediction
US9225889B1 (en) Photographic image acquisition device and method
CN102726037B (en) Image processing apparatus, camera head and image processing method
CN103139469B (en) Utilize multiresolution process to generate the system and method for robust depth map
Shu et al. Imaging via three-dimensional compressive sampling (3DCS)
EP4055555A1 (en) Noise reconstruction for image denoising
CN110880163B (en) Low-light color imaging method based on deep learning
Niu et al. Low cost edge sensing for high quality demosaicking
CN115994865A (en) Restoration method, acquisition restoration device and monitoring system for flutter blurred image
KR20200084419A (en) Appratus for generating moire removing model, method for removing moire and imaging device for removing moire
KR20160004912A (en) Method and apparatus for image capturing and simultaneous depth extraction
CN104376547A (en) Motion blurred image restoration method
US20120148108A1 (en) Image processing apparatus and method therefor
CN116703752A (en) Image defogging method and device of near infrared fused transducer structure
CN117571128B (en) High-resolution polarized spectrum image imaging method and system
CN104539851B (en) High-speed imaging system and method based on pixel optimization coding exposure
US20090316994A1 (en) Method and filter for recovery of disparities in a video stream
CN111583345B (en) Method, device and equipment for acquiring camera parameters and storage medium
CN113658128A (en) Image blurring degree determining method, data set constructing method and deblurring method
CN116389912B (en) Method for reconstructing high-frame-rate high-dynamic-range video by fusing pulse camera with common camera
CN111210390A (en) Motion blur restoration method based on Golay sequence complementary code word set
US20230080120A1 (en) Monocular depth estimation device and depth estimation method
WO2022207110A1 (en) Noise reconstruction for image denoising
Reshetouski et al. Lensless mismatched aspect ratio imaging
CN113573076A (en) Method and apparatus for video encoding

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination