CN108717685A - A kind of method and system of enhancing image resolution ratio - Google Patents
A kind of method and system of enhancing image resolution ratio Download PDFInfo
- Publication number
- CN108717685A CN108717685A CN201810456228.XA CN201810456228A CN108717685A CN 108717685 A CN108717685 A CN 108717685A CN 201810456228 A CN201810456228 A CN 201810456228A CN 108717685 A CN108717685 A CN 108717685A
- Authority
- CN
- China
- Prior art keywords
- image
- hot spot
- original image
- point spread
- spread function
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000002708 enhancing effect Effects 0.000 title claims abstract description 19
- 238000012545 processing Methods 0.000 claims abstract description 33
- 238000000605 extraction Methods 0.000 claims abstract description 25
- 238000006243 chemical reaction Methods 0.000 claims abstract description 19
- 230000000007 visual effect Effects 0.000 claims abstract description 8
- 230000003287 optical effect Effects 0.000 claims description 40
- 230000010287 polarization Effects 0.000 claims description 22
- 238000009826 distribution Methods 0.000 claims description 13
- 241000931526 Acer campestre Species 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- 230000008676 import Effects 0.000 claims description 6
- 239000004744 fabric Substances 0.000 claims 1
- 230000005764 inhibitory process Effects 0.000 abstract description 5
- 239000000284 extract Substances 0.000 abstract description 4
- 230000001737 promoting effect Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 5
- 239000007850 fluorescent dye Substances 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 101100117236 Drosophila melanogaster speck gene Proteins 0.000 description 4
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 4
- 229910052782 aluminium Inorganic materials 0.000 description 4
- 239000004411 aluminium Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- ZYGHJZDHTFUPRJ-UHFFFAOYSA-N coumarin Chemical compound C1=CC=C2OC(=O)C=CC2=C1 ZYGHJZDHTFUPRJ-UHFFFAOYSA-N 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 210000004027 cell Anatomy 0.000 description 2
- 210000000170 cell membrane Anatomy 0.000 description 2
- 229960000956 coumarin Drugs 0.000 description 2
- 235000001671 coumarin Nutrition 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004043 dyeing Methods 0.000 description 2
- 238000000386 microscopy Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000002203 pretreatment Methods 0.000 description 2
- 238000007334 copolymerization reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000002073 fluorescence micrograph Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Microscoopes, Condenser (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
The present invention provides a kind of method and system of enhancing image resolution ratio, wherein the method includes：Pending original image is obtained, and in the case where keeping the original image visual field size constant, the original image is extended to the pixel of specified quantity；Based on selected wavelet function, two-dimensional discrete wavelet conversion carried out to the image after expansion, and removes picture noise and extracts during the two-dimensional discrete wavelet conversion and specify image；Deconvolution processing is carried out to the specified image of extraction using preset algorithm, to obtain the image of resolution ratio enhancing.Technical solution provided by the present application can realize the inhibition to noise information while promoting image resolution ratio.
Description
Technical field
The present invention relates to technical field of image processing, more particularly to a kind of method and system of enhancing image resolution ratio.
Background technology
Deconvolution techniques be it is a kind of it is widely applied, without to target carry out in advance identification and signature analysis, have compared with
The image processing techniques of high universality can significantly improve the resolution ratio of imaging technique, be more than optical system used to obtain
The high-definition picture of system resolution ratio itself provides more sample informations for researcher.
Currently, the difference of difference or realization method based on algorithm, there are many kinds of types for deconvolution techniques, including most pass through
The Lucy-Richardson methods of allusion quotation and more flexible blind methods.The realization of Lucy-Richardson methods needs
Based on the accurate expection of the point spread function (Point Spread Function, PSF) to optical system, to realize to mesh
The deconvolution processing of logo image.The latter is not needing the accurate expection to the PSF of optical system for the former, only needs
Substantially to a range, blind algorithms that can go to approach actual PSF by built-in iterative algorithm to realize to target figure
The analysis of deconvoluting of picture.
However, deconvolution techniques used at present, including Lucy-Richardson and blind methods, exist excessive
The negative factor of processing.It is, deconvolution techniques while removing noise, improving resolution ratio, can additionally introduce a large amount of
Information introduces a large amount of undulatory patterns while removing random noise and shotnoise (shot noise).This is to image
The identification of middle effective information judges and the promotion of picture quality has apparent negative effect.
Invention content
The application's is designed to provide a kind of method and system of enhancing image resolution ratio, can promote image resolution
While rate, the inhibition to noise information is realized.
To achieve the above object, the application provides a kind of method of enhancing image resolution ratio, the method includes：
Pending original image is obtained, and in the case where keeping the original image visual field size constant, it will be described
Original image extends to the pixel of specified quantity；
Based on selected wavelet function, two-dimensional discrete wavelet conversion is carried out to the image after expansion, and small from two-dimensional discrete
Noise is removed during wave conversion and is extracted specifies image；
Deconvolution processing is carried out to the specified image of extraction using preset algorithm, to obtain the figure of resolution ratio enhancing
Picture.
Further, obtaining pending original image includes：
It imports the text data of the original image or imports the original image for specifying picture format.
Further, the pixel for the original image being extended to specified quantity includes：
According to the total quantity of the current pixel of the original image, determine with the total quantity of the current pixel most
The n times exponential value of close 2, described immediate 2 n times exponential value are more than or equal to the sum of the current pixel
Amount, and the N is integer；
The quantity for the pixel that the original image includes is extended to described immediate 2 n times exponential value.
Further, specified image is extracted from the result of two-dimensional discrete wavelet conversion includes：
Signal from the result of two-dimensional discrete wavelet conversion corresponding to scale of the extraction with peak signal energy carries out
Reduction, and the signal corresponding to smallest dimension signal and scale with least energy is eliminated, it is small from two-dimensional discrete to obtain
Noise and the specified image of extraction are removed during wave conversion.
Further, carrying out deconvolution processing to the specified image of extraction includes：
The point spread function of the fluorescence hot spot of calculating optical system, and based on the point spread function to the finger of extraction
Determine image and carries out deconvolution processing.
Further, the point spread function of calculating optical system includes：
If the unknown parameters of the fluorescence hot spot of the optical system, by the number for inputting microscope types, microcobjective
It is worth aperture and medium refraction index, and theoretical based on microscopical ideal hot spot, it is assumed that incident light is collimated light beam and no polarization, knot
Gaussian Profile is closed, the point spread function distribution of fluorescence hot spot is calculated；
Further, the point spread function of calculating optical system includes：
If the parameter of the optical system is it is known that by the wavelength of laser, polarization state, the numerical aperture of microcobjective and medium
Refractive index inputs hot spot simulation program, using numerical simulation, calculates different refractivity, numerical aperture, optical maser wavelength, polarization state
Under, the point spread function of laser is distributed, and calculates the point spread function distribution of fluorescence hot spot in turn.
To achieve the above object, the application also provides a kind of system of enhancing image resolution ratio, the system comprises：
Image expansion unit for obtaining pending original image, and is keeping the original image visual field size not
In the case of change, the original image is extended to the pixel of specified quantity；
Wavelet transform unit, for based on selected wavelet function, 2-d discrete wavelet change to be carried out to the image after expansion
It changes, and removes noise during the two-dimensional discrete wavelet conversion and extract and specify image；
Deconvolution processing unit, for carrying out deconvolution processing to the specified image of extraction using preset algorithm, with
Obtain the image of resolution ratio enhancing.
Further, the deconvolution processing unit includes：
Point spread function calculating module, the point spread function of the fluorescence hot spot for calculating optical system, and based on described
Point spread function carries out deconvolution processing to the specified image of extraction.
Further, the point spread function calculating module includes：
Unknown parameters processing module, it is aobvious by inputting if the unknown parameters of the fluorescence hot spot for the optical system
Micro mirror type, the numerical aperture of microcobjective and medium refraction index, and it is theoretical based on microscopical ideal hot spot, it is assumed that incident light
For collimated light beam and no polarization the point spread function distribution of fluorescence hot spot is calculated in conjunction with Gaussian Profile；
Parameter known treatment module, if the parameter for the optical system is it is known that by the wavelength of laser, polarization state, show
The numerical aperture and medium refraction index of speck mirror input hot spot simulation program, using numerical simulation, calculate different refractivity, numerical value
Under aperture, optical maser wavelength, polarization state, the point spread function of laser is distributed, and calculates the point spread function of fluorescence hot spot in turn
Distribution.
Therefore in the present invention, by the way that two-dimentional (2D) discrete wavelet analytical technology is combined with the methods of deconvoluting,
Pre-treatment is carried out to image with the analysis of 2D discrete wavelets, and then analysis of deconvoluting is carried out to image using the methods of deconvoluting again,
Finally it is aided with non-linear sharpening, to realize the very big inhibition to noise and be obviously improved to image resolution ratio.With existing skill
Art is compared, and the present invention at least has following advantageous effect：
1, noise reduction is extremely notable, and signal-to-noise ratio is high；
2, it can be obviously improved to graininess, point plaque-like, the resolution ratio of circle and Similar Round Object；
3, with strong points：It is analyzed by 2D discrete wavelets, sample is effectively extracted, in conjunction with deconvolution algorithm, improved
The clarity and resolution ratio of sample；
4, applied widely：Suitable for different wave length, the derived techniques of different fluorescent dye, different lens systems
Post processing of image, including inverted fluorescence microscope, Laser Scanning Confocal Microscope, stimulated radiation loss microscope etc..
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
It obtains it is clear that understand through the implementation of the invention.The purpose of the present invention and other advantages can be by the explanations write
Specifically noted structure is realized and is obtained in book, claims and attached drawing.
Below by drawings and examples, technical scheme of the present invention will be described in further detail.
Description of the drawings
Attached drawing is used to provide further understanding of the present invention, and a part for constitution instruction, the reality with the present invention
It applies example to be used to explain the present invention together, not be construed as limiting the invention.In the accompanying drawings：
Fig. 1 is the method flow diagram for enhancing image resolution ratio in the embodiment of the present invention；
Fig. 2 is the method flow diagram that PSF is calculated in the embodiment of the present invention；
Fig. 3 (a) is the original image in the present invention the first application example；
Fig. 3 (b) is treated image in the present invention the first application example；
Fig. 4 (a) is the original image in the present invention the second application example；
Fig. 4 (b) is treated image in the present invention the second application example；
Fig. 5 (a) is the original image in third application example of the present invention；
Fig. 5 (b) is treated image in third application example of the present invention.
Specific implementation mode
Hereinafter, preferred embodiments of the present invention will be described with reference to the accompanying drawings, it should be understood that preferred reality described herein
Apply example only for the purpose of illustrating and explaining the present invention and is not intended to limit the present invention.
Referring to Fig. 1, the application provides a kind of method of enhancing image resolution ratio, the method includes：
S1：Pending original image is obtained, and in the case where keeping the original image visual field size constant, by institute
State the pixel that original image extends to specified quantity；
S2：Based on selected wavelet function, two-dimensional discrete wavelet conversion carried out to the image after expansion, and from two-dimensional discrete
Noise is removed during wavelet transformation and is extracted specifies image；
S3：Deconvolution processing is carried out to the specified image of extraction using preset algorithm, to obtain resolution ratio enhancing
Image.
In the present embodiment, obtaining pending original image includes：
It imports the text data of the original image or imports the original image for specifying picture format.
In the present embodiment, the pixel for the original image being extended to specified quantity includes：
According to the total quantity of the current pixel of the original image, determine with the total quantity of the current pixel most
The n times exponential value of close 2, described immediate 2 n times exponential value are more than or equal to the sum of the current pixel
Amount, and the N is integer；
The quantity for the pixel that the original image includes is extended to described immediate 2 n times exponential value.
In the present embodiment, specified image is extracted from the result of two-dimensional discrete wavelet conversion includes：
Signal from the result of two-dimensional discrete wavelet conversion corresponding to scale of the extraction with peak signal energy carries out
Reduction, and the signal corresponding to smallest dimension signal and scale with least energy is eliminated, it is small from two-dimensional discrete to obtain
Noise and the specified image of extraction are removed during wave conversion.
In the present embodiment, carrying out deconvolution processing to the specified image of extraction includes：
The point spread function of the fluorescence hot spot of calculating optical system, and based on the point spread function to the finger of extraction
Determine image and carries out deconvolution processing.
In the present embodiment, the point spread function of the fluorescence hot spot of calculating optical system includes：
If the unknown parameters of the fluorescence hot spot of the optical system, by the number for inputting microscope types, microcobjective
It is worth aperture and medium refraction index, and theoretical based on microscopical ideal hot spot, it is assumed that incident light is collimated light beam and no polarization, knot
Gaussian Profile is closed, the point spread function distribution of fluorescence hot spot is calculated；
In the present embodiment, the point spread function of the fluorescence hot spot of calculating optical system includes：
If the parameter of the optical system is it is known that by the wavelength of laser, polarization state, the numerical aperture of microcobjective and medium
Refractive index inputs hot spot simulation program, using numerical simulation, calculates different refractivity, numerical aperture, optical maser wavelength, polarization state
Under, the point spread function of laser is distributed, and calculates the point spread function distribution of fluorescence hot spot in turn.
In a practical application scene, the technical solution of the application can turn to software with product, this software is for promotion
The resolution ratio of picture system and develop, developing instrument be Matlab softwares.This software includes mainly image expansion, 2D discrete wavelets
4 parts such as analysis, the image procossing based on deconvolution method and non-linear sharpening are constituted.
First, original image is read in from data file.Data file can be Document type data file, can also be
Or the image file of .GIF formats .GIFf.
Then, original image is expanded.Since the n times power that discrete wavelet is generally basede on 2 carries out, thus need former
Beginning image spreading to immediate 2 n times power is counted.Such as original image multiplies 100 pixels for 100.It is inserted by non-linear Gauss
Value is extended, and in the case where not influencing original image visual field size, obtains 128 images for multiplying 128 pixels.
In turn, by selecting different wavelet functions (such as db3), discrete wavelet processing is carried out to the image after expansion,
The signal of extraction energy out to out is restored, and eliminates the corresponding noise of energy smallest dimension, discrete small to obtain 2D
The image of wave extraction.
After the image for obtaining the extraction of 2D discrete wavelets, further deconvolution processing is carried out to image.It is built-in in this software
Four kinds of image algorithms, respectively Lucy-Richardson algorithms, Blind algorithms, regularized filter and Wiener
Filter algorithms.
Deconvolution method is dependent on the estimation or prediction to optical system fluorescence hot spot PSF.This program provides two kinds of moulds
Block estimates the PSF of the fluorescence hot spot of optical system, and by taking copolymerization is burnt and inverted fluorescence microscope as an example, specific flow is as schemed
Shown in 2：
(1) if the unknown parameters of exciting light or the non-optical profession of operator, as shown in the right side of fig 2, it is only necessary to input
Microscope types, the numerical aperture of microcobjective and medium refraction index, so it is theoretical based on microscopical ideal hot spot, it is assumed that enter
It is that collimated light beam (infinity focus) and no polarization (or circular polarization state) calculate ideal fluorescence light in conjunction with Gaussian Profile to penetrate light
Spot PSF distributions.
(2) if the parameter of exciting light it is known that and operator be familiar with basic optical principle, as shown on the left of Fig. 2, by light
Wavelength, polarization state, the numerical aperture of microcobjective and medium refraction index input hot spot simulation program, using numerical simulation, calculate
Under different refractivity, numerical aperture, optical maser wavelength, polarization state, the PSF of laser is distributed, and then calculates the PSF of fluorescence hot spot.
Microscopy Microscopy is lost for stimulated radiation, PSF function can also be by calculating exciting light and the exciting light that disappears
PSF so that estimate to obtain in conjunction with a variety of inhibition function modes.PSF is updated in software, is calculated based on four kinds of images above-mentioned
Method promotes the resolution ratio of image and reduces influence of noise, extracts true image information.Finally, it is further to promote image pair
Than degree, we are handled image using the non-linear sharpening algorithm of index, the super-resolution image that ultimately generates that treated.
Below from multiple technique effects for illustrating the application using example：
Application example 1：
After being illustrated in fig. 3 shown below to apply the image processing method, the handling result of confocal images.In Fig. 3 (b)
For the anodised aluminium porous structure of coated 102 fluorescent dyes of Coumarin in surface, the fluorescence speck in figure is (such as arrow institute
Show) by the fluorescent dye that is full of in the hole of the anodised aluminium of unprocessed mistake, also represent the structure in anodised aluminium hole.
These fluorescence hot spots are influenced by noise of detector and optical diffraction, it is difficult to effectively identify its exact position, size and shape.And
After the software of the application, as shown in Fig. 3 (b), while retaining useful signal, making an uproar for image is not only greatly inhibited
Sound, while greatly improving the resolution ratio to target.For example, in Fig. 3 (a), part shown in arrow is very fuzzy, and is linked to be
It is a piece of, it is difficult to distinguish.By the processing of the software, as shown in the arrow in Fig. 3 (b), the image to join together is by clearly area
It separates, and picture noise is minimum, the resolution ratio of original image is about 300nm, and the resolution ratio of treated image is more than
50nm。
Application example 2：
For the anodised aluminium porous structure of coated 102 fluorescent dyes of Coumarin in same surface, we utilize
Stimulated radiation loss microscope has also carried out scanning analysis to it.Fig. 4 (a) is that microscopical original scan figure is lost in stimulated radiation
Picture, although resolution ratio is promoted relative to Laser Scanning Confocal Microscope, noise is equally serious, it is difficult to effectively identify its accurate position
It sets, size and shape.And by the software of the application, as shown in Fig. 4 (b), while retaining useful signal, not only greatly
The noise of image is inhibited, while greatly improving the resolution ratio to target.The structure that cannot be distinguished in many Fig. 3 (b),
It can be distinguished in Fig. 4 (b), the resolution ratio of treated image is more than 40nm.
Application example 3：
Fig. 5 show the comparison before and after the processing of the inverted fluorescence microscope image of the cell of GFP dyeing.Fig. 5 (a) is dyeing
The original image of cell afterwards, it can be seen that the ambient noise of cell membrane produces seriously observation fluorescence speck (shown in arrow)
Negative effect so that image is very fuzzy.And in Fig. 5 (b), not only the bulk background of cell membrane fractions is successfully removed,
Meanwhile as shown by arrows, the spatial resolution of sample is greatly promoted.
The application also provides a kind of system of enhancing image resolution ratio, the system comprises：
Image expansion unit for obtaining pending original image, and is keeping the original image visual field size not
In the case of change, the original image is extended to the pixel of specified quantity；
Wavelet transform unit, for based on selected wavelet function, 2-d discrete wavelet change to be carried out to the image after expansion
It changes, and removes noise during the two-dimensional discrete wavelet conversion and extract and specify image；
Deconvolution processing unit, for carrying out deconvolution processing to the specified image of extraction using preset algorithm, with
Obtain the image of resolution ratio enhancing.
In the present embodiment, the deconvolution processing unit includes：
Point spread function calculating module, the point spread function of the fluorescence hot spot for calculating optical system, and based on described
Point spread function carries out deconvolution processing to the specified image of extraction.
In the present embodiment, the point spread function calculating module includes：
Unknown parameters processing module, it is aobvious by inputting if the unknown parameters of the fluorescence hot spot for the optical system
Micro mirror type, the numerical aperture of microcobjective and medium refraction index, and it is theoretical based on microscopical ideal hot spot, it is assumed that incident light
For collimated light beam and no polarization the point spread function distribution of fluorescence hot spot is calculated in conjunction with Gaussian Profile；
Parameter known treatment module, if the parameter for the optical system is it is known that by the wavelength of laser, polarization state, show
The numerical aperture and medium refraction index of speck mirror input hot spot simulation program, using numerical simulation, calculate different refractivity, numerical value
Under aperture, optical maser wavelength, polarization state, the point spread function of laser is distributed, and calculates the point spread function of fluorescence hot spot in turn
Distribution.
Therefore in the present invention, by the way that two-dimentional (2D) discrete wavelet analytical technology is combined with the methods of deconvoluting,
Pre-treatment is carried out to image with the analysis of 2D discrete wavelets, and then analysis of deconvoluting is carried out to image using the methods of deconvoluting again,
Finally it is aided with non-linear sharpening, to realize the very big inhibition to noise and be obviously improved to image resolution ratio.With existing skill
Art is compared, and the present invention at least has following advantageous effect：
1, noise reduction is extremely notable, and signal-to-noise ratio is high；
2, it can be obviously improved to graininess, point plaque-like, the resolution ratio of circle and Similar Round Object；
3, with strong points：It is analyzed by 2D discrete wavelets, sample is effectively extracted, in conjunction with deconvolution algorithm, improved
The clarity and resolution ratio of sample；
4, applied widely：Suitable for different wave length, the derived techniques of different fluorescent dye, different lens systems
Post processing of image, including inverted fluorescence microscope, Laser Scanning Confocal Microscope, stimulated radiation loss microscope etc..
Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art
God and range.In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technologies
Within, then the present invention is also intended to include these modifications and variations.
Claims (10)
1. a kind of method of enhancing image resolution ratio, which is characterized in that the method includes：
Pending original image is obtained, and in the case where keeping the original image visual field size constant, it will be described original
Image extends to the pixel of specified quantity；
Based on selected wavelet function, two-dimensional discrete wavelet conversion is carried out to the image after expansion, and become from 2-d discrete wavelet
Noise is removed during changing and is extracted specifies image；
Deconvolution processing is carried out to the specified image of extraction using preset algorithm, to obtain the image of resolution ratio enhancing.
2. according to the method described in claim 1, it is characterized in that, obtaining pending original image and including：
It imports the text data of the original image or imports the original image for specifying picture format.
3. according to the method described in claim 1, it is characterized in that, the original image to be extended to the pixel of specified quantity
Including：
According to the total quantity of the current pixel of the original image, determination is closest with the total quantity of the current pixel
2 n times exponential value, described immediate 2 n times exponential value is more than or equal to the total quantity of the current pixel,
And the N is integer；
The quantity for the pixel that the original image includes is extended to described immediate 2 n times exponential value.
4. according to the method described in claim 1, it is characterized in that, extracting specified figure from the result of two-dimensional discrete wavelet conversion
As including：
Signal from the result of two-dimensional discrete wavelet conversion corresponding to scale of the extraction with peak signal energy is restored,
And the signal corresponding to smallest dimension signal and scale with least energy is eliminated, to obtain from two-dimensional discrete wavelet conversion
During remove the specified image of noise and extraction.
5. according to the method described in claim 1, it is characterized in that, the specified image to extraction carries out deconvolution processing packet
It includes：
The point spread function of the fluorescence hot spot of calculating optical system, and based on the point spread function to the specified figure of extraction
As carrying out deconvolution processing.
6. according to the method described in claim 5, it is characterized in that, the point spread function packet of the fluorescence hot spot of calculating optical system
It includes：
If the unknown parameters of the fluorescence hot spot of the optical system, by the numerical aperture for inputting microscope types, microcobjective
Diameter and medium refraction index, and it is theoretical based on microscopical ideal hot spot, it is assumed that and incident light is collimated light beam and no polarization, in conjunction with height
This distribution calculates the point spread function distribution of fluorescence hot spot.
7. according to the method described in claim 5, it is characterized in that, the point spread function packet of the fluorescence hot spot of calculating optical system
It includes：
If the parameter of the optical system is it is known that by the wavelength of laser, polarization state, the numerical aperture of microcobjective and medium refraction
Rate inputs hot spot simulation program, using numerical simulation, calculates under different refractivity, numerical aperture, optical maser wavelength, polarization state,
The point spread function of laser is distributed, and calculates the point spread function distribution of fluorescence hot spot in turn.
8. a kind of system of enhancing image resolution ratio, which is characterized in that the system comprises：
Image expansion unit for obtaining pending original image, and is keeping the original image visual field size constant
In the case of, the original image is extended to the pixel of specified quantity；
Wavelet transform unit, for based on selected wavelet function, two-dimensional discrete wavelet conversion to be carried out to the image after expansion, and
Noise is removed during the two-dimensional discrete wavelet conversion and is extracted specifies image；
Deconvolution processing unit, for carrying out deconvolution processing to the specified image of extraction using preset algorithm, to obtain
The image of resolution ratio enhancing.
9. system according to claim 8, which is characterized in that the deconvolution processing unit includes：
Point spread function calculating module, the point spread function of the fluorescence hot spot for calculating optical system, and expanded based on the point
The specified image for dissipating function pair extraction carries out deconvolution processing.
10. system according to claim 9, which is characterized in that the point spread function calculating module includes：
Unknown parameters processing module, if the unknown parameters of the fluorescence hot spot for the optical system, by inputting microscope
Type, the numerical aperture of microcobjective and medium refraction index, and it is theoretical based on microscopical ideal hot spot, it is assumed that subject to incident light
Collimated optical beam and no polarization calculate the point spread function distribution of fluorescence hot spot in conjunction with Gaussian Profile；
Parameter known treatment module, if the parameter for the optical system is it is known that by the wavelength of laser, polarization state, micro- object
The numerical aperture and medium refraction index of mirror input hot spot simulation program, using numerical simulation, calculate different refractivity, numerical aperture
Under diameter, optical maser wavelength, polarization state, the point spread function of laser is distributed, and calculates the point spread function point of fluorescence hot spot in turn
Cloth.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810456228.XA CN108717685B (en) | 2018-05-14 | 2018-05-14 | Method and system for enhancing image resolution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810456228.XA CN108717685B (en) | 2018-05-14 | 2018-05-14 | Method and system for enhancing image resolution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108717685A true CN108717685A (en) | 2018-10-30 |
| CN108717685B CN108717685B (en) | 2021-07-06 |
Family
ID=63899828
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810456228.XA Active CN108717685B (en) | 2018-05-14 | 2018-05-14 | Method and system for enhancing image resolution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108717685B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110390270A (en) * | 2019-06-27 | 2019-10-29 | 北京集创北方科技股份有限公司 | Optical finger print recognition methods, device, electronic equipment and storage medium |
| CN112114045A (en) * | 2020-08-03 | 2020-12-22 | 中国科学院大学 | Resolution enhancement method for ultrasonic microscope |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004111609A2 (en) * | 2003-06-12 | 2004-12-23 | Predicant Biosciences, Inc. | Methods for accurate component intensity extraction from separations-mass spectrometry data |
| CN102708543A (en) * | 2012-04-19 | 2012-10-03 | 北京工商大学 | Fluorescent microscopic image restoring method based on blind deconvolution and sparse representation and device thereof |
| CN106770095A (en) * | 2016-11-30 | 2017-05-31 | 浙江大学 | A kind of super-resolution micro imaging method and device based on the modulation of non-linear hot spot |
-
2018
- 2018-05-14 CN CN201810456228.XA patent/CN108717685B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004111609A2 (en) * | 2003-06-12 | 2004-12-23 | Predicant Biosciences, Inc. | Methods for accurate component intensity extraction from separations-mass spectrometry data |
| CN102708543A (en) * | 2012-04-19 | 2012-10-03 | 北京工商大学 | Fluorescent microscopic image restoring method based on blind deconvolution and sparse representation and device thereof |
| CN106770095A (en) * | 2016-11-30 | 2017-05-31 | 浙江大学 | A kind of super-resolution micro imaging method and device based on the modulation of non-linear hot spot |
Non-Patent Citations (7)
| Title |
|---|
| KAY Y. LIU等: "A source wavelet deconvolution approach to improve the spatial resolution for radar-based breast imaging system", 《IEEE XPLORE》 * |
| M"HAMED BOULAKROUNE等: "Multiscale deconvolution using wavelet transform for improving the depth resolution in secondary ion mass spectrometry analysis", 《IEEE XPLORE》 * |
| RAFAEL C. GONZALEZ: "《数字图像处理》", 30 June 2011 * |
| 刘金琳: "超声图像复原方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 程鹏: "荧光显微图像的解卷积方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
| 肖昀: "点扫描成像高分辨率光学显微镜中点扩散函数的研究", 《中国博士学位论文全文数据库基础科学辑》 * |
| 魏通达: "共聚焦激光扫描光学显微成像关键技术研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110390270A (en) * | 2019-06-27 | 2019-10-29 | 北京集创北方科技股份有限公司 | Optical finger print recognition methods, device, electronic equipment and storage medium |
| CN112114045A (en) * | 2020-08-03 | 2020-12-22 | 中国科学院大学 | Resolution enhancement method for ultrasonic microscope |
| CN112114045B (en) * | 2020-08-03 | 2023-09-19 | 中国科学院大学 | Resolution enhancement method for ultrasonic microscope |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108717685B (en) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wallace et al. | A workingperson’s guide to deconvolution in light microscopy | |
| Sarder et al. | Deconvolution methods for 3-D fluorescence microscopy images | |
| Tai et al. | Single image defocus map estimation using local contrast prior | |
| US7224824B1 (en) | Computerized adaptive imaging | |
| AU2006303039B2 (en) | Method of processing an image | |
| Cannell et al. | Image enhancement by deconvolution | |
| JP2015192238A (en) | Image data generation device and image data generation method | |
| Holmes et al. | Blind deconvolution | |
| CN108717685A (en) | A kind of method and system of enhancing image resolution ratio | |
| JP6362062B2 (en) | Image generating apparatus and image generating method | |
| Rudnaya et al. | A derivative-based fast autofocus method in electron microscopy | |
| Kam | Microscopic differential interference contrast image processing by line integration (LID) and deconvolution | |
| Prigent et al. | SPITFIR (e): a supermaneuverable algorithm for fast denoising and deconvolution of 3D fluorescence microscopy images and videos | |
| Conchello et al. | Extended depth-of-focus microscopy via constrained deconvolution | |
| Pankajakshan | Blind deconvolution for confocal laser scanning microscopy | |
| TW201935074A (en) | Method and apparatus for extending depth of field during fluorescence microscopy imaging | |
| Shajkofci et al. | DeepFocus: a few-shot microscope slide auto-focus using a sample invariant CNN-based sharpness function | |
| JP2015191362A (en) | Image data generation apparatus and image data generation method | |
| Parton et al. | Lifting the fog: Image restoration by deconvolution | |
| Dogan et al. | A novel extended depth of field process based on nonsubsampled shearlet transform by estimating optimal range in microscopic systems | |
| Beckers et al. | Real-time, extended depth DIC microscopy | |
| Balasubramanian et al. | A ringing metric to evaluate the quality of images restored using iterative deconvolution algorithms | |
| Cao et al. | Dark-based Optical Sectioning assists Background Removal in Fluorescence Microscopy | |
| KR102578312B1 (en) | Hologram-based object three-dimensional mesurement device and method | |
| Mutahira et al. | Unorthodox approach toward microscopic shape from image focus using optical microscopy |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |