CN103345727A - Method for reconstructing binary optical image spectrum - Google Patents
Method for reconstructing binary optical image spectrum Download PDFInfo
- Publication number
- CN103345727A CN103345727A CN2013102376378A CN201310237637A CN103345727A CN 103345727 A CN103345727 A CN 103345727A CN 2013102376378 A CN2013102376378 A CN 2013102376378A CN 201310237637 A CN201310237637 A CN 201310237637A CN 103345727 A CN103345727 A CN 103345727A
- Authority
- CN
- China
- Prior art keywords
- image
- spectrum
- magnification
- standard deviation
- point spread
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 62
- 230000003287 optical effect Effects 0.000 title claims abstract description 55
- 238000012952 Resampling Methods 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 2
- 238000001914 filtration Methods 0.000 abstract 2
- 230000000694 effects Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 230000007812 deficiency Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
A method for reconstructing a binary optical image spectrum includes the following detailed steps that (1), original binary optical image spectrum data are read in; (2), according to parameters and the structure of an optical system, whether the optical system is an equal magnification system or not is judged, and if the optical system is an unequal magnification system, images of all wave bands with unequal magnifications are converted into equal magnification images; (3) standard errors of point spread functions of all the wave bands are calculated; (4), a noise level is analyzed, and if noise does not exist, a traditional inverse filtering method is used for restoring spatial dimension images; (5), for the images with the noise and the small standard errors of the point spread functions, a Jansson-VanCitter algorithm is used for restoring the spatial dimension images; for the images with the noise and the large stand errors of the point spread functions, a Norbert wiener inverse filtering method is adopted for restoring the spatial dimension images; (6), a linear deconvolution method is adopted for reconstructing spectrum dimension data; (7), the result of image spectrum reconstruction is obtained.
Description
(1) technical field
The present invention relates to a kind of method for reconstructing of binary optical image spectrum, belong to high-spectral data and rebuild and processing technology field, it is applicable to the treatment theory methods and applications technical research of binary optical system data.
(2) background technology
High spectrum image has been realized the breakthrough raising of spectral resolution, has realized " spectrum unification ".But it is too much to be directed to the high spectrum image wave band, thereby the energy that causes each wave band to obtain is less, has produced low this problem of signal to noise ratio (S/N ratio).Binary optical elements can need not carry out light splitting by slit in the system simultaneously at space peacekeeping spectrum dimensional imaging along the optical axis direction chromatic dispersion, has strengthened projectile energy greatly, thereby has improved signal to noise ratio (S/N ratio).But simultaneously, this light that causes each wave band also can become vague image in adjacent band except becoming the sharply defined image in self focal plane, and namely each width of cloth image of binary optical image all is the stack of the defocusing amount of the sharply defined image of this wave band and adjacent band, cannot directly use.The rebuilding spectrum technology of binary optical image can be removed the influence of defocusing amount, obtains the spectroscopic data of true available binary optical image.
To the rebuilding spectrum Study on Technology mainly from the research of point spread function and research two aspects of Image Restoration Algorithm are carried out.
Make J
0Be the zeroth order Bessel function of the first kind, λ is wavelength, ρ=||
v||=(ξ
2+ η
2)
1/2Be the polar coordinates of the face that becomes of lens back,
v=(ξ η) is the coordinate of putting on the back focal plane,
NA, M are respectively numerical aperture and the enlargement factors of lens, z
dBe lens to the distance of receiving plane,
Be the polar coordinates on the two dimensional sample plane, z
0Be the degree of depth of three-dimensional samples, P (ρ, z
0, T) be the pupil function of lens, then the model of the point spread function of actual light path is
But parameter is too much in this model, much is difficult to obtain, and comprises complex calculation such as integration, differentiate, can greatly strengthen operand, therefore is difficult in practice use.And Gauss model parameter similar with it is few, calculates easyly, can be used for replacing point spread function is simulated.Image restoration relates to disposal routes such as denoising, 2D signal deconvolution, and common method has nearest neighbor method, linear solution convolution algorithm, nonlinear iteration algorithm etc. at present.For the model of different known conditions, different imaging characteristicses, should select for use different algorithms to carry out rebuilding spectrum work, and these methods also have its different characteristics separately.It is very fast to receive the liftering algorithm speed as arest neighbors, liftering, dimension, but recovering quality is not high; And need the algorithm speed of iteration slower, but the recovery result is generally better.
At present, the image spectrum data re-establishing method based on binary optical mainly contains following deficiency: the optical system magnification of (1) general supposing the system is identical, does not consider not wait the situation of magnification system, makes method for reconstructing not be suitable for the magnification system that do not wait; (2) do not carry out the selection of method for reconstructing according to the different noise level of data, occur noise easily and amplify, flood problems such as image spectrum useful information; (3) do not carry out the selection of distinct methods according to image spectrum data standard difference difference, influence reconstructed image spectroscopic data quality thereby occur ringing effect easily.
(3) summary of the invention
The method for reconstructing that the purpose of this invention is to provide a kind of binary optical image spectrum, it overcome the existing rebuilding spectrum method scope of application narrower, be difficult to deficiency that the image spectrum data of different system are handled, can determine best rebuilding spectrum algorithm voluntarily according to systematic parameter, it is the method for reconstructing of the binary optical image spectrum that a kind of applicability is strong, reliability is high.
Technical solution of the present invention is: a kind of imaging system for different parameters, can determine automatically to obtain the method that best binary optical image spectrum is rebuild effect according to its characteristics.This method mainly is based on the standard deviation of system point spread function, noise level, the isoparametric judgement of magnification, method by image resampling will not wait that the magnification image is converted into etc. magnification image, utilize methods such as linear solution convolution, nonlinear iteration that image is carried out the reconstruction of space peacekeeping spectrum dimension, thereby obtain the image spectrum reconstructed results.
The method for reconstructing of a kind of binary optical image spectrum of the present invention, its concrete steps are as follows:
(1) reads in original binary optical image spectrum data;
(2) carry out whether waiting the magnification system to judge according to optical system parameter and structure, if do not wait the magnification optical system, each band image that does not wait magnification such as is converted at the magnification image;
(3) calculate the standard deviation of each wave band point spread function;
(4) noise level is analyzed, if noiseless uses traditional liftering method to carry out space dimension image restoration;
(5) for the image spectrum data that have noise level, carry out the selection of space dimension restored method according to the difference of point spread function standard deviation;
(6) carry out the reconstruction of spectrum dimension data according to the binary optical image-forming principle;
(7) obtain the image spectrum reconstructed results.
Wherein, whether what carry out in the step (2) is to wait the magnification system to judge, its objective is magnification such as will not wait that the magnification image is converted into, at first carrying out the identical imaging region of different-waveband selects, adopting method for resampling will not wait each band image of magnification to be converted to minimum magnification then is the image of benchmark, makes this algorithm applicable to waiting magnification and not waiting two kinds of optical systems of magnification.
Wherein, the standard deviation of each wave band point spread function of calculating described in the step (3) refers to that sharp two-dimentional Gauss model calculates point spread function:
Wherein, μ
1Be the mean value of x direction, σ
1Be the standard deviation of x direction, μ
2Be the mean value of y direction, σ
2Be the standard deviation of y direction, equate in the standard deviation of x direction and y direction, when the center is identical that this formula can be reduced to:
Wherein, μ is the mean value of this Gauss model, can adopt the actual image point coordinate figure; σ is the standard deviation of Gauss model, the maximal value substitution formula of point spread function matrix in the experiment can be tried to achieve.
Wherein, noise level is analyzed described in the step (4) is to adopt the signal noise ratio (snr) of image based on the image statistics parameter to calculate the analysis of realization noise level, and carries out the selection of image spectrum method for reconstructing according to noise level.
Wherein, step (5) described " for the image spectrum data that have noise level; carry out the selection of space dimension restored method according to the difference of point spread function standard deviation ", refer to for having noise and the less image of point spread function standard deviation, use Jansson-Van Citter method to carry out space dimension image restoration; For the bigger image of the standard deviation that has noise and point spread function, adopt dimension to receive the liftering method and carry out the image restoration of space dimension.
Wherein, " the carrying out the reconstruction of spectrum dimension data according to the binary optical image-forming principle " described in the step (6) refers to carry out according to the binary optical image-forming principle reconstruction of spectrum dimension data, and this spectrum dimension data method for reconstructing adopts the linear solution convolution method.
The present invention's advantage compared with prior art is: it is comparatively independent, single to have overcome traditional rebuilding spectrum algorithm, only be applicable to the shortcoming of the optical system with certain features, can carry out rebuilding spectrum to the resulting binary optical image of the system of difference spread function, noise level, system's magnification simultaneously, all can obtain effect preferably.It has following advantage: the method that (1) utilizes image resampling magnification such as is converted into each band image, and it is magnification is unified to minimum magnification, make method can be applied to wait magnification and the optical system that does not wait magnification simultaneously, and the gridiron pattern effect of effective removal of images; (2) by the judgement to noise level, the real image that has noise is adopted different image spectrum reconstruction algorithm with muting image, suppressed noise effect effectively, thereby realized that image spectrum is rebuild fast and accurately; (3) by the judgement to the point spread function standard deviation, can take diverse ways to carry out image spectrum for different system and rebuild, eliminated " ring " effect of rebuilding back image spectrum data effectively.
(4) description of drawings
Fig. 1 is FB(flow block) of the present invention
(5) embodiment
Be 450nm-900nm with the spectral band scope, the wave band number is 30 data instance, and Fig. 1 is seen in the method for reconstructing specific implementation process explanation of a kind of binary optical image spectrum that the present invention relates to, the method for reconstructing of a kind of binary optical image spectrum of the present invention, the specific implementation step is as follows:
(1) reads in original binary optical image spectrum data; Read in the binary optical image spectrum data of 30 wave bands in the 450nm-900nm wavelength band;
(2) carry out whether waiting the magnification system to judge according to optical system parameter and structure, if do not wait the magnification optical system, each band image that does not wait magnification such as is converted at the magnification image; Judge that by structure, the parameter of optical system whether optical system the magnification optical system such as is, inequality through judging each band image magnification, for not waiting magnification optical system, at first determine the scope of same area in each band image according to the magnification relation, be of a size of benchmark with same area in the minimum magnification image again, same area in each band image is carried out image resampling, obtaining waiting the image of magnification, is benchmark can prevent the gridiron pattern effect under the prerequisite that guarantees arithmetic speed appearance with minimum magnification;
(3) calculate the standard deviation of each wave band point spread function; The point spread function matrix of each wave band data of reading in of statistics with the maximum point of the maximal value substitution Gauss model of matrix, can be determined Gauss model, utilizes the Gauss model statistical parameter to calculate the standard deviation of system point spread function;
(4) noise level is analyzed, if noiseless uses traditional liftering method to carry out space dimension image restoration; Judge according to Analysis signal-to-noise ratio (SNR), comprised noise in the raw data of reading in, therefore, adopt the algorithm in the step (5) to carry out image restoration;
(5) for the image spectrum data that have noise level, carry out the selection of space dimension restored method according to the difference of point spread function standard deviation; If the point spread function standard deviation is less, then can use Jansson-Van Citter algorithm to carry out the image restoration of space dimension, reach the effect of image restoration by the method for nonlinear iteration, though speed is slower, but this method can not amplified noise, and recovery result's sharpness is higher; But the method has strict requirement to the standard deviation of point spread function, when the standard deviation of point spread function is big, can produce very serious ringing owing to repeatedly image being blured in iterative process; Standard deviation for point spread function is bigger, adopt dimension to receive the liftering method and carry out image restoration, the standard deviation of the raw data point spread function that reads in through judgement is 1.34, therefore adopt dimension to receive the liftering method and carry out space dimension image restoration, this method has added the factor of considering noise on the basis of traditional liftering method, noise can be exaggerated, and owing to there is not iteration can not produce ringing;
(6) carry out the reconstruction of spectrum dimension data according to the binary optical image-forming principle; Tie up at spectrum, it all is stacks of the defocusing amount of the sharply defined image of this wave band and adjacent band that the binary optical image does not have the image of each wave band, the response power that each wave band light can cause is to be determined at the point spread function of spectrum dimension by system, therefore, the spectrum dimension is carried out rebuilding spectrum can adopt the linear solution convolution method, the point spread function of spectrum dimension is the one dimension Gauss model;
(7) obtain the image spectrum reconstructed results; The result that will carry out the reconstruction of space peacekeeping spectrum dimension exports, thereby obtains the image spectrum data reconstructed results of 450-900nm wavelength band.
Claims (6)
1. the method for reconstructing of a binary optical image spectrum, it is characterized in that: it comprises following concrete steps:
(1) reads in original binary optical image spectrum data;
(2) carry out whether waiting the magnification system to judge according to optical system parameter and structure, if do not wait the magnification optical system, each band image that does not wait magnification such as is converted at the magnification image;
(3) calculate the standard deviation of each wave band point spread function;
(4) noise level is analyzed, if noiseless uses traditional liftering method to carry out space dimension image restoration;
(5) for the image spectrum data that have noise level, carry out the selection of space dimension restored method according to the difference of point spread function standard deviation;
(6) carry out the reconstruction of spectrum dimension data according to the binary optical image-forming principle;
(7) obtain the image spectrum reconstructed results.
2. the method for reconstructing of a kind of binary optical image spectrum according to claim 1, it is characterized in that: whether what carry out in its step (2) is to wait the magnification system to judge, its objective is magnification such as will not wait that the magnification image is converted into, at first carrying out the identical imaging region of different-waveband selects, adopting method for resampling will not wait each band image of magnification to be converted to minimum magnification then is the image of benchmark, makes this algorithm applicable to waiting magnification and not waiting two kinds of optical systems of magnification.
3. the method for reconstructing of a kind of binary optical image spectrum according to claim 1, it is characterized in that: the standard deviation of each wave band point spread function of calculating described in the step (3) refers to utilize Gauss model to calculate the standard deviation of point spread function:
Wherein, μ
1Be the mean value of x direction, σ
1Be the standard deviation of x direction, μ
2Be the mean value of y direction, σ
2Be the standard deviation of y direction, equate in the standard deviation of x direction and y direction, when the center is identical that this formula can be reduced to:
Wherein, μ is the mean value of this Gauss model, can adopt the actual image point coordinate figure; σ is the standard deviation of Gauss model, the maximal value substitution formula of point spread function matrix in the experiment can be tried to achieve.
4. the method for reconstructing of a kind of binary optical image spectrum according to claim 1, it is characterized in that: noise level is analyzed described in its step (4), be to adopt the signal noise ratio (snr) of image based on the image statistics parameter to calculate the analysis of realization noise level, and carry out the selection of image spectrum method for reconstructing according to noise level.
5. the method for reconstructing of a kind of binary optical image spectrum according to claim 1, it is characterized in that: its step (5) described " for the image spectrum data that have noise level; carry out the selection of space dimension restored method according to the difference of point spread function standard deviation ", refer to for having noise and the less image of point spread function standard deviation, use Jansson-Van Citter method to carry out space dimension image restoration; For the bigger image of the standard deviation that has noise and point spread function, adopt dimension to receive the liftering method and carry out the image restoration of space dimension.
6. the method for reconstructing of a kind of binary optical image spectrum according to claim 1, it is characterized in that: " the carrying out the reconstruction of spectrum dimension data according to the binary optical image-forming principle " described in its step (6), refer to carry out according to the binary optical image-forming principle reconstruction of spectrum dimension data, this spectrum dimension data method for reconstructing adopts the linear solution convolution method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310237637.8A CN103345727B (en) | 2013-06-17 | 2013-06-17 | A kind of method for reconstructing of binary optical image spectrum |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310237637.8A CN103345727B (en) | 2013-06-17 | 2013-06-17 | A kind of method for reconstructing of binary optical image spectrum |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103345727A true CN103345727A (en) | 2013-10-09 |
CN103345727B CN103345727B (en) | 2016-12-07 |
Family
ID=49280520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310237637.8A Expired - Fee Related CN103345727B (en) | 2013-06-17 | 2013-06-17 | A kind of method for reconstructing of binary optical image spectrum |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103345727B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105890757A (en) * | 2014-12-15 | 2016-08-24 | 中国科学院国家天文台 | One-dimensional spectrum for extracting multi-fiber spectrum through adoption of deconvolution method |
CN107923238A (en) * | 2015-07-29 | 2018-04-17 | 哈利伯顿能源服务公司 | Use integrated computing element structural remodeling spectrum |
CN110533617A (en) * | 2019-08-30 | 2019-12-03 | Oppo广东移动通信有限公司 | Image processing method and device, storage medium |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006223444A (en) * | 2005-02-16 | 2006-08-31 | Hiroshi Hirata | Electron spin resonance imaging apparatus |
CN102488493A (en) * | 2011-11-15 | 2012-06-13 | 西安电子科技大学 | Small animal living body multi-mode molecule imaging system and imaging method |
-
2013
- 2013-06-17 CN CN201310237637.8A patent/CN103345727B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006223444A (en) * | 2005-02-16 | 2006-08-31 | Hiroshi Hirata | Electron spin resonance imaging apparatus |
CN102488493A (en) * | 2011-11-15 | 2012-06-13 | 西安电子科技大学 | Small animal living body multi-mode molecule imaging system and imaging method |
Non-Patent Citations (2)
Title |
---|
XUHUI WANG等: "High-resolution Holographic Imaging Technology by Microscopic Image Plane Holography", 《PROCEEDINGS OF CONFERENCE 7848 HOLOGRAPHY ,DIFFRACTIVE OPTICS AND APPLICATIONS》, 18 October 2010 (2010-10-18), pages 1 - 8 * |
张萌等: "高光谱数据光谱分辨率对矿物识别的影响分析", 《红外与激光工程》, vol. 35, 15 October 2006 (2006-10-15), pages 493 - 498 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105890757A (en) * | 2014-12-15 | 2016-08-24 | 中国科学院国家天文台 | One-dimensional spectrum for extracting multi-fiber spectrum through adoption of deconvolution method |
CN107923238A (en) * | 2015-07-29 | 2018-04-17 | 哈利伯顿能源服务公司 | Use integrated computing element structural remodeling spectrum |
CN110533617A (en) * | 2019-08-30 | 2019-12-03 | Oppo广东移动通信有限公司 | Image processing method and device, storage medium |
CN110533617B (en) * | 2019-08-30 | 2022-05-27 | Oppo广东移动通信有限公司 | Image processing method and device, and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN103345727B (en) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sahu et al. | Different image fusion techniques–a critical review | |
Guo et al. | High quality multi-focus image fusion using self-similarity and depth information | |
Fugal et al. | Practical methods for automated reconstruction and characterization of particles in digital in-line holograms | |
Fernandes et al. | Pavement pathologies classification using graph-based features | |
Nehme et al. | Learning optimal wavefront shaping for multi-channel imaging | |
CN110992363B (en) | Three-dimensional magnetic flux leakage detection defect reconstruction method based on generation countermeasure network | |
CN103886332A (en) | Method for detecting and identifying defects of metallic meshes | |
US10521918B2 (en) | Method and device for filtering texture, using patch shift | |
CN103345727A (en) | Method for reconstructing binary optical image spectrum | |
Mutahira et al. | A simplified approach using deep neural network for fast and accurate shape from focus | |
Zhu et al. | Spectral dual-channel encoding for image dehazing | |
Hu et al. | Fabric defect detection via un-decimated wavelet decomposition and gumbel distribution model | |
EP2227705A2 (en) | Method and apparatus for providing image data | |
Sun et al. | Double‐flow convolutional neural network for rapid large field of view Fourier ptychographic reconstruction | |
Zhang et al. | Edge detection of three-dimensional objects by manipulating pupil functions in an optical scanning holography system | |
Maalouf | Contribution to fluorescence microscopy, 3D thick samples deconvolution and depth-variant PSF | |
CN108717685B (en) | Method and system for enhancing image resolution | |
Sánchez-Ferreira et al. | Multi-objective differential evolution algorithm for underwater image restoration | |
Rahmat et al. | 3D shape from focus using LULU operators and discrete pulse transform in the presence of noise | |
CN112001956B (en) | CNN-based image denoising method for strong laser far-field focal spot measurement by using schlieren method | |
Drozdova et al. | Radio-astronomical image reconstruction with a conditional denoising diffusion model | |
CN107292859B (en) | Chaotic medium polarization image acquisition method based on optical correlator | |
Proppe et al. | 3D-2D Neural Nets for Phase Retrieval in Noisy Interferometric Imaging | |
Wang et al. | Blind deconvolution subject to sparse representation for fluorescence microscopy | |
Sequeira et al. | PENGUIN: A rapid and efficient image preprocessing tool for multiplexed spatial proteomics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161207 Termination date: 20170617 |