CN109883549B - Digital micromirror-based bending spectral line correction method - Google Patents
Digital micromirror-based bending spectral line correction method Download PDFInfo
- Publication number
- CN109883549B CN109883549B CN201910162585.XA CN201910162585A CN109883549B CN 109883549 B CN109883549 B CN 109883549B CN 201910162585 A CN201910162585 A CN 201910162585A CN 109883549 B CN109883549 B CN 109883549B
- Authority
- CN
- China
- Prior art keywords
- line
- spectral line
- bending
- wavelength
- spectral
- 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.)
- Active
Links
Images
Landscapes
- Spectrometry And Color Measurement (AREA)
Abstract
The invention relates to a method for correcting a bending spectral line based on a digital micromirror, which comprises the following steps: obtaining three-dimensional full spectrum data; determining the range of each bending spectral line in the spectrogram; acquiring the spectral center position of each line of each curved spectral line and the corresponding signal intensity by adopting a light intensity value weighting method; determining a fittable range; step five, fitting a relational expression of positions and wavelengths of each row of pixels in a fittable range to obtain wavelength values corresponding to all the pixels; calculating the corresponding wavelength value of the pixel of the part which cannot be fitted; and the bending and overturning detection of the digital micromirror is realized. The invention does not need to consider error terms caused by installation and adjustment errors and environmental errors, and can improve the signal-to-noise ratio of the measurement spectral line.
Description
Technical Field
The invention belongs to the technical field of spectral analysis, and particularly relates to a method for correcting a bending spectral line based on a digital micromirror.
Background
In a spectral analysis instrument, it is a very common method to realize the spatial separation of polychromatic light by using a diffraction grating, wherein an incident slit of a spectrometer is perpendicular to a main cross section of the grating, a principal ray at the center point of the slit is superposed with the main cross section of the grating, and principal rays at other points of the slit form included angles with the main cross section of the grating, thereby causing the bending of spectral lines. Spectral resolution, luminous flux and signal-to-noise ratio are affected by the presence of line bending phenomena. The existing method for correcting the spectral line bending starts from the aspect of optical design, the bending spectral line correction method based on data processing is based on the basic optical formula to carry out complex formula derivation calculation, the calculation process is complicated, and the theoretical derivation is easy to be influenced by the environment to cause deviation. In the method for actually controlling the digital micromirror, the Hadamard template and the array (whole column) overturning mode are mainly used, and the Hadamard template mainly aims at a continuous light source; column inversion, although the simplest control method, affects the accuracy and signal-to-noise ratio of the light intensity signal with respect to the line bending phenomenon. At present, no method for detecting the turning of the bending track of the digital micromirror is available, and the detection method can be used for accurately and rapidly correcting the bent spectral line.
The invention discloses an imaging spectrometer system capable of correcting spectral line bending and a correction method thereof (application publication No. CN 104034419), a spectral element of a grating-prism-grating is designed, the spectral line bending of the center wavelength is compensated and corrected by utilizing the characteristic that the dispersion of the grating and the prism is opposite to the spectral line bending direction, and the residual spectral line bending is compensated and corrected by utilizing the distortion and the image plane inclination generated by a collimating objective and an imaging objective, so that the calibration of the spectral line bending of the full working spectral line is realized.
The invention discloses a spectral line bending long-wave infrared plane grating imaging spectral system (application publication number: CN 103048045A), which utilizes an off-axis lens to eliminate spectral line bending generated by a grating.
The invention discloses a method and a device for correcting spectral line bending of a dispersion type spectral imager (application publication No. CN 104316183A), which carry out spectral line bending correction by utilizing a related wavelength Euclidean distance calculation result and an inverse distance weight method.
The articles "spectral aliasing analysis and correction of hadamard transform spectral imager" (zhou brocade, luqun wave, phase in. "spectral aliasing analysis and correction of hadamard transform spectral imager [ J ] photonics, 2005(10): 1518-.
The paper "analysis and correction of spectral line bending for digital micromirror Hadamard spectrometer" ("analysis and correction of spectral line bending for digital micromirror Hadamard spectrometer," [ J ]. Spectroscopy and Spectroscopy, ", 2016,36(02): 555-.
The above prior art has the following disadvantages: 1. the calculation formula derived from the theoretical formula of the optical system cannot be directly calculated by the formula to obtain a final result due to unavoidable installation and adjustment errors and the like; 2. parameters such as a plurality of angles involved in an optical formula are inconvenient to actually measure, so that the practical application is restricted, and therefore, the existing verification methods are all simulation verification. 3. What is detected in the array (column) flip control mode is local spectral information, resulting in a decrease in spectral signal-to-noise ratio.
Disclosure of Invention
The invention aims to provide a method for correcting a bending spectral line based on a digital micromirror, which does not need to consider error items caused by installation and adjustment errors and environmental errors and has high signal-to-noise ratio of a measured spectral line.
In order to solve the above technical problem, the method for correcting the bending spectral line based on the digital micromirror of the present invention comprises:
step one, obtaining three-dimensional full-spectrum data; controlling the digital micro-mirror array to scan the known spectral line wavelength lambda point by point or block by block1,λ2,···λn···,λNThe excitation light source of a certain element carries out full spectrum scanning to obtain full spectrum data;
step two, determining the range of each bending spectral line in a spectrogram; taking a multiple of a wavelength-free position Noise value Noise in a spectrogram as a signal intensity Limit value Limit, wherein the value of a is selected by taking the aim of extracting all spectrograms as far as possible and not extracting Noise as a target; defining the position of the signal intensity more than or equal to the Limit value in the spectrogram as a spectral line position, and further extracting N curved spectral lines;
thirdly, acquiring the spectral center position of each line of the N curved spectral lines and the corresponding signal intensity by adopting a light intensity value weighting method; for the mth row of the nth curved spectral line, the abscissa of the pixel at the center position of the transverse spectrumWherein xnm-iIs the abscissa of the ith pixel element of the mth line of the nth curved spectral line, Inm-iThe signal intensity corresponding to the ith pixel of the mth line of the nth curved spectral line; i isnmThe number of pixel columns contained in the mth row of the nth curved spectral line is set;
step four, sequencing the longitudinal upper limits of the N bending spectral lines from high to low, and sequencing the longitudinal lower limits of the N bending spectral lines from low to high; determining all action fittable ranges at least comprising two curved spectral line row pixels in one row;
step five, according to the theoretically known spectral line wavelength lambda1,λ2,···λn···,λNWith respect to the position of the pixel columnThe method comprises the steps of solving the corresponding relation between the center position of each line spectrum in each curved spectral line and the wavelength in a fitting range, fitting a relation between the position of each line pixel and the wavelength in the fitting range, and obtaining the wavelength value corresponding to all pixels in the fitting range according to the relation between the position of each line pixel and the wavelength;
sixthly, calculating the wavelength value corresponding to the part of pixels which cannot be fitted; finding curved lines n containing non-fittable partsbCalculating the curved spectral line n by using a light intensity value weighting methodbOrdinate of central pixel of longitudinal spectrumWherein y isnb-kRepresenting curved spectral lines nbThe longitudinal coordinate of the central pixel of the transverse spectrum of the nth b-K line, and nb-K is a bending spectral line nbTotal number of lines of (1)nb-kIs the signal strength of the pixel; obtaining a bending spectral line n by utilizing the relation between all line pixels and the wavelength in the fitting range fitted in the step fivebThe wavelength values corresponding to all the pixels of the row where the center position of the longitudinal spectrum is located; using the central position of transverse spectrum as the alignment reference and using the curved spectral line nbThe wavelength value corresponding to the row pixel where the longitudinal spectrum center position is located completely fills the non-fittable row pixels;
step seven, realizing the bending and overturning detection of the digital micromirror; when the element detection is carried out on the sample, all pixels corresponding to the wavelength value to be detected are searched by using the wavelength value corresponding to each pixel obtained in the fifth step and the sixth step, and the digital micromirrors in the pixels are controlled to simultaneously turn over, so that the measurement of the bending spectral line of the element to be detected of the sample is realized.
The invention has the advantages that: 1. correcting a bending spectral line caused by aberration of an optical element in an optical path of the analysis instrument by controlling bending and overturning of the digital micromirror; 2. compared with the method for correcting the bending spectral line through light path design, the method has the advantages that the bending control is carried out according to real-time detection results without considering error items caused by installation and adjustment errors and environmental errors; 3. by bending the scanning lines, the signal-to-noise ratio of the measured lines can be improved compared to scanning in an array.
Drawings
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
FIG. 1 is a schematic diagram of three-dimensional full-spectrum X-Y view angle curved spectral line detection.
FIG. 2 is a partial enlarged view of 306.80nm wavelength bend line A of the Bi element.
FIG. 3 is a schematic diagram showing wavelength values corresponding to pixels at the longitudinal scanning center position of an 306.80nm wavelength bending spectral line of a Bi element.
FIG. 4 is a diagram showing the correspondence between the unadapted part I at the upper end of the 306.80nm wavelength bending spectral line of the Bi element and the line wavelength value at the longitudinal scan center position.
FIG. 5 is a diagram showing the correspondence between the line wavelength value at the line position of the lower unfixed portion II of the 306.80nm wavelength bending spectral line of the Bi element and the longitudinal scan center position.
Fig. 6 is a flow chart of the present invention.
Detailed Description
As shown in fig. 6, the method for correcting the bending spectral line based on the digital micromirror of the present invention is as follows:
step one, obtaining three-dimensional full spectrum data: as shown in fig. 1, taking a digital micromirror with a specification of 1024 × 768 as an example, controlling a digital micromirror array to perform full-spectrum scanning on an excitation light source of a known element in a point-by-point or block-by-block scanning manner to obtain full-spectrum data; the full spectrum data includes the signal strength corresponding to each pixel (1 or more digital micromirrors may be included in a pixel); wherein the element corresponds to a spectral line wavelength lambda1,λ2,···λn···,λN(ii) a The signal intensity is represented by light color, and the darker the color, the stronger the signal intensity;
step two, determining the range of each bending spectral line in the three-dimensional spectrogram obtained in the step one (namely the abscissa range and the ordinate of a pixel corresponding to each bending spectral line): taking a multiple of a wavelength-free position Noise value Noise in a spectrogram as a signal intensity Limit value Limit, namely, the Limit is a multiplied by Noise, and the value of a is selected by taking the aim of extracting all spectrograms as far as possible and not extracting Noise as a target; the position of the signal intensity which is more than or equal to the Limit value in the spectrogram is specified as the spectral line position, and then the spectral line position is obtainedTaking out N bending spectral lines (corresponding wavelength lambda of N bending spectral lines)1,λ2,···,λNKnown) corresponding pixel coordinate ranges;
step three, acquiring the spectral center position and the corresponding signal intensity of each line of the N curved spectral lines: because the spectral imaging has nonuniformity, a light intensity value weighting method is needed for calculating the transverse spectral center position of each line of the curved spectral line; for the mth row of the nth curved spectral line, the abscissa of the pixel at the center position of the transverse spectrumWherein xnm-iIs the abscissa of the ith pixel element of the mth line of the nth curved spectral line, Inm-iThe signal intensity corresponding to the ith pixel of the mth line of the nth curved spectral line; i isnmThe number of pixel columns contained in the mth row of the nth curved spectral line is set;
determining the range of the relational expression between the pixel positions and the wavelengths of each row in the fittable spectrogram; because the heights of the curved spectral lines are different in the longitudinal direction, some lines on a spectrogram only have one curved spectral line pixel position, and fitting cannot be performed; therefore, the longitudinal upper limits of the N bending spectral lines need to be sequenced from high to low, and similarly, the longitudinal lower limits of the N bending spectral lines need to be sequenced from low to high; determining all action fittable ranges including at least two bending spectral line pixel positions in a row; as shown in fig. 1, the ordinate is YT~YBAll rows in the range are fitting ranges;
step five, according to the theoretically known spectral line wavelength lambda1,λ2,···λn···,λNThe wavelength lambda of the spectral line is measured in a fitting range according to the position of the pixel array1,λ2,···λn···,λNAnd respectively assigning pixels to the transverse spectrum centers of each row of the corresponding curved spectral line, and fitting a relation between the positions and the wavelengths of each row of pixels in a fittable range according to the pixel positions: if the fittable range includes J rows of pixels, the relationship between the pixels in the first row to the J-th row and the wavelength is as follows: lambda [ alpha ]1=f1(X),λ2=f2(X)···λj=fj(X)…λJ=fJ(X), so that the corresponding wavelength values of all the pixels in the fitting range can be obtained; as shown in fig. 2, the pixel at which the wavelength of 306.8nm is located is the central pixel of the transverse spectrum of the 306.80nm wavelength bending spectral line of the Bi element;
sixthly, calculating the wavelength value corresponding to the part of pixels which cannot be fitted; finding curved lines n containing non-fittable partsb(ii) a The light intensity value weighting method is also used for calculating the curved spectral line nbOrdinate of central pixel of longitudinal spectrumWherein y isnb-kRepresenting curved spectral lines nbOrdinate, I, of center pixel of transverse spectrum of the nth b-k linenb-kIs the signal strength of the pixel; obtaining a bending spectral line n by utilizing the relation between all line pixels and the wavelength in the fitting range fitted in the step fivebThe wavelength values corresponding to all the pixels of the row where the center position of the longitudinal spectrum is located; using the central position of transverse spectrum as the alignment reference and using the curved spectral line nbThe wavelength value corresponding to the row pixel where the longitudinal spectrum center position is located completely fills the non-fittable row pixels; taking FIG. 1 as an example, the curved spectral line nbIs the curved line N in fig. 1; as shown in fig. 3, the curved line nbLongitudinal spectrum center pixel line pixel (X)O-3,Ynb-O),(XO-2,Ynb-O),(XO-1,Ynb-O),(XO,Ynb-O),(XO+1,Ynb-O),(XO+2,Ynb-O),(XO+3,Ynb-O) The corresponding wavelength values are 306.2nm, 306.4nm,306.6nm, 306.8nm, 307.0nm,307.2nm and 307.4nm respectively; as shown in FIG. 4, the picture element (X)O,Ynb-O) A wavelength value of 306.8nm, is copied and filled into the curved spectral line nbPixel (X) of the upper unfit part InbT1,YnbT1) Pixel (X)O,Ynb-O),(XO+1,Ynb-O) The wavelength values of 306.8nm and 307.0nm are respectively copied and filled in the bending spectral line nbPixel (X) of the upper unfit part InbT1,YnbT1),(XnbT2,YnbT2),(XnbT2+1,YnbT2) (ii) a As shown in FIG. 5, the picture element (X)O-1,Ynb-O),(XO,Ynb-O),(XO+1,Ynb-O) Respectively 306.6nm, 306.8nm and 307.0, are copied and filled in the bending spectral line nbThe lower extreme does not fit the picture element (X) of part IInbB4-1,YnbB4),(XnbB4,YnbB4),(XnbB4+1,YnbB4) Let pixel (X)O,Ynb-O),(XO+1,Ynb-O) 306.8nm, 307.0nm copied to fill the bend line nbThe lower extreme does not fit the picture element (X) of part IInbB3,YnbB3),(XnbB3+1,YnbB3) (ii) a Pixel (X)O,Ynb-O) The value of 306.8nmm of the wavelength is copied and filled into the curved spectral line nbThe lower extreme does not fit the picture element (X) of part IInbB2,YnbB2) (ii) a Finally, the pixel (X)O,Ynb-O) The value of 306.8nmm of the wavelength is copied and filled into the curved spectral line nbThe lower extreme does not fit the picture element (X) of part IInbB1,YnbB1);
Step seven, realizing the bending and overturning detection of the digital micromirror; when the element detection is carried out on the sample, all pixels corresponding to the wavelength value to be detected are searched by using the wavelength value corresponding to each pixel obtained in the fifth step and the sixth step, and the digital micromirrors in the pixels are controlled to simultaneously turn over, so that the measurement of the bending spectral line of the element to be detected of the sample is realized.
Claims (1)
1. A method for correcting bending spectral lines based on a digital micromirror is characterized by comprising the following steps:
step one, obtaining three-dimensional full-spectrum data; controlling the digital micro-mirror array to scan the known spectral line wavelength lambda point by point or block by block1,λ2,…λn…,λNThe excitation light source of a certain element carries out full spectrum scanning to obtain full spectrum data;
step two, determining the range of each bending spectral line in a spectrogram; taking a multiple of a wavelength-free position Noise value Noise in a spectrogram as a signal intensity Limit value Limit, wherein the value of a is selected by taking the aim of extracting all spectrograms as far as possible and not extracting Noise as a target; defining the position of the signal intensity more than or equal to the Limit value in the spectrogram as a spectral line position, and further extracting N curved spectral lines;
thirdly, acquiring the spectral center position of each line of the N curved spectral lines and the corresponding signal intensity by adopting a light intensity value weighting method; for the mth row of the nth curved spectral line, the abscissa of the pixel at the center position of the transverse spectrumWherein xnm-iIs the abscissa of the ith pixel element of the mth line of the nth curved spectral line, Inm-iThe signal intensity corresponding to the ith pixel of the mth line of the nth curved spectral line; i isnmThe number of pixel columns contained in the mth row of the nth curved spectral line is set;
step four, sequencing the longitudinal upper limits of the N bending spectral lines from high to low, and sequencing the longitudinal lower limits of the N bending spectral lines from low to high; determining all action fittable ranges at least comprising two curved spectral line row pixels in one row;
step five, according to the theoretically known spectral line wavelength lambda1,λ2,…λn…,λNThe corresponding relation between the center position of each line spectrum in each bending spectral line and the wavelength in the fitting range is worked out according to the relation between the center position of each line spectrum and the position of the wavelength in each bending spectral line in the fitting range, the relation between the position of each line pixel in the fitting range and the wavelength is fitted, and the wavelength value corresponding to all the pixels in the fitting range can be obtained according to the relation between the position of each line pixel and the wavelength;
sixthly, calculating the wavelength value corresponding to the part of pixels which cannot be fitted; finding curved lines n containing non-fittable partsbCalculating the curved spectral line n by using a light intensity value weighting methodbOrdinate of central pixel of longitudinal spectrumWherein y isnb-kRepresenting curved spectral lines nbThe longitudinal coordinate of the central pixel of the transverse spectrum of the nth b-K line, and nb-K is a bending spectral line nbTotal number of lines of (1)nb-kIs the signal strength of the pixel; by using stepsObtaining a bending spectral line n by the relation between all line pixels and the wavelength in the fitting range fitted in the fifth stepbThe wavelength values corresponding to all the pixels of the row where the center position of the longitudinal spectrum is located; using the central position of transverse spectrum as the alignment reference and using the curved spectral line nbThe wavelength value corresponding to the row pixel where the longitudinal spectrum center position is located completely fills the non-fittable row pixels;
step seven, realizing the bending and overturning detection of the digital micromirror; when the element detection is carried out on the sample, all pixels corresponding to the wavelength value to be detected are searched by using the wavelength value corresponding to each pixel obtained in the fifth step and the sixth step, and the digital micromirrors in the pixels are controlled to simultaneously turn over, so that the measurement of the bending spectral line of the element to be detected of the sample is realized.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910162585.XA CN109883549B (en) | 2019-03-05 | 2019-03-05 | Digital micromirror-based bending spectral line correction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910162585.XA CN109883549B (en) | 2019-03-05 | 2019-03-05 | Digital micromirror-based bending spectral line correction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109883549A CN109883549A (en) | 2019-06-14 |
CN109883549B true CN109883549B (en) | 2021-03-12 |
Family
ID=66930633
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910162585.XA Active CN109883549B (en) | 2019-03-05 | 2019-03-05 | Digital micromirror-based bending spectral line correction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109883549B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110455760B (en) * | 2019-08-14 | 2021-06-25 | 吉林大学 | Dispersion type AFS light source scattering interference deduction method based on DMD |
CN111580352B (en) * | 2020-06-08 | 2021-06-29 | 吉林大学 | Method for measuring and correcting rotating center of rotating table in digital photoetching system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158599A (en) * | 2007-11-13 | 2008-04-09 | 清华大学 | Wavelength calibration method of optical spectrum instruments |
CN108896179A (en) * | 2018-06-29 | 2018-11-27 | 长春理工大学 | DMD space dimension encodes symmetrical Offner dispersion medium-wave infrared optical spectrum imaging device |
CN108956554A (en) * | 2018-05-18 | 2018-12-07 | 吉林大学 | The method for wavelength calibration of Atomic Fluorescence Spectrometer based on digital micromirror array |
CN208420171U (en) * | 2018-08-10 | 2019-01-22 | 中南民族大学 | A kind of identification device of the communication beam orbital angular momentum mode based on digital micro-mirror |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10054483B1 (en) * | 2017-06-01 | 2018-08-21 | Westco Scientific Instruments, Inc | Micromirror spectrophotometer assembly |
-
2019
- 2019-03-05 CN CN201910162585.XA patent/CN109883549B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158599A (en) * | 2007-11-13 | 2008-04-09 | 清华大学 | Wavelength calibration method of optical spectrum instruments |
CN108956554A (en) * | 2018-05-18 | 2018-12-07 | 吉林大学 | The method for wavelength calibration of Atomic Fluorescence Spectrometer based on digital micromirror array |
CN108896179A (en) * | 2018-06-29 | 2018-11-27 | 长春理工大学 | DMD space dimension encodes symmetrical Offner dispersion medium-wave infrared optical spectrum imaging device |
CN208420171U (en) * | 2018-08-10 | 2019-01-22 | 中南民族大学 | A kind of identification device of the communication beam orbital angular momentum mode based on digital micro-mirror |
Also Published As
Publication number | Publication date |
---|---|
CN109883549A (en) | 2019-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7489396B1 (en) | Spectrophotometric camera | |
US7113281B2 (en) | Device for the pixel-by-pixel photoelectric measurement of a planar measured object | |
JP2005265655A (en) | Spectral reflectance measuring device, film thickness measuring device and spectral reflectance measuring method | |
US10198984B2 (en) | Display panel calibration using detector array measurement | |
CN109883549B (en) | Digital micromirror-based bending spectral line correction method | |
US20080284896A1 (en) | Joint Optics and Image Processing Adjustment of Electro-Optic Imaging Systems | |
JPH08233689A (en) | Method and system for testing performance of at least one electrooptic tester | |
CN106769909B (en) | A kind of modified method and device of imaging spectrometer scaling position | |
JP5424108B2 (en) | Raman imaging equipment | |
JPS63266484A (en) | Checking of display characteristic | |
JP3304692B2 (en) | Spectrometer | |
US6222631B1 (en) | Two-dimensional spectral characteristic measuring apparatus | |
CN111457911B (en) | Bionic polarization compass calibration method based on polarization two-dimensional residual error information | |
CN113175883B (en) | Light source normalization processing method of spectrum confocal measurement system | |
KR100724374B1 (en) | Film thickness acquiring method | |
Shi et al. | Segmented mirror coarse phasing with a dispersed fringe sensor: experiments on NGST's Wavefront Control Testbed | |
CN111076816A (en) | Full-field spectrum calibration error correction method for large-aperture static interference imaging spectrometer | |
CN110174351B (en) | Color measuring device and method | |
CN114441474B (en) | Near-infrared spectrometer and control method and system thereof | |
JP2010078418A (en) | Spectral measurement apparatus, calibration apparatus, spectral measurement method and calibration method | |
KR101642354B1 (en) | calibration apparatus and method for CCD of spectrograph | |
JP2007147507A (en) | Spectrometry and optical spectrometer | |
CN104236710A (en) | Spectrum super-resolution method of handheld light source color illumination photometry meter | |
KR102022836B1 (en) | Apparatus for measuring light, system and method thereof | |
Wöhl et al. | Precise reduction of solar spectra obtained with large CCD arrays |
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 |