CN114329872A - Method and system for calculating and simulating ellipticity stability of optical remote sensing camera - Google Patents
Method and system for calculating and simulating ellipticity stability of optical remote sensing camera Download PDFInfo
- Publication number
- CN114329872A CN114329872A CN202011048679.3A CN202011048679A CN114329872A CN 114329872 A CN114329872 A CN 114329872A CN 202011048679 A CN202011048679 A CN 202011048679A CN 114329872 A CN114329872 A CN 114329872A
- Authority
- CN
- China
- Prior art keywords
- ellipticity
- optical
- stability
- ellipsometry
- remote sensing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
The method and the system for calculating and simulating the ellipticity stability of the optical remote sensing camera can comprehensively consider the influence of various error items on the ellipticity stability, realize the evaluation of related optical indexes and provide data support for astronomical research.
Description
Technical Field
The invention relates to the field of integrated simulation of optical remote sensing cameras, in particular to a method and a system for calculating and simulating the ellipticity stability of an optical remote sensing camera.
Background
With the continuous development of the optical remote sensing imaging technology, the space optical remote sensing system becomes more and more complex, and develops towards the direction of large caliber, large view field and high resolution, the difficulty of corresponding processing, manufacturing, supporting and adjusting is increasing, and the observation performance is seriously influenced by the interference of emission, gravity environment and temperature factors in the actual imaging process. In the process, whether the ellipsometry and the ellipsometry stability which are important indexes for astronomy research influence meet the requirements or not needs to be simulated and calculated.
Disclosure of Invention
In view of the above, it is necessary to provide a method for simulating ellipsometry stability of an optical remote sensing camera, which can perform simulation calculation on ellipsometry stability as an index required for astronomical observation.
A method for calculating and simulating the ellipticity stability of an optical remote sensing camera comprises the following steps:
establishing an optical model of the optical remote sensing camera;
determining a suitable sampling interval for the point spread function;
calculating a point spread function of a certain field of view in the optical model containing diffraction effect according to the suitable sampling interval;
calculating the ellipticity and the ellipticity component value of the point spread function;
and obtaining the stability of the ellipse according to the ellipse and the ellipse component value.
In some embodiments, in the step of establishing the optical model of the optical remote sensing camera, the optical model includes an ideal optical system and an optical system considering error factors including, but not limited to, an optical design residual, a processing surface shape error, a system setup residual, an error caused by a gravity environment change, and a thermal environment change.
In some embodiments, in the step of determining the suitable sampling interval of the point spread function, specifically according to the sampling theorem, a sampling interval range that is not under-sampled is determined, and then the suitable sampling interval is determined by calculating the ellipticity through PSFs taken at different sampling intervals.
In some embodiments, the step of calculating a point spread function of a field of view containing diffraction effects in the optical model according to the suitable sampling interval specifically includes:
and calculating a point spread function PSF (point spread function) of a certain field of view of the optical model, including a diffraction effect, by using optical software CodeV through Fourier transform, tracking the coordinate XY of a chief ray on an image surface, and calculating the distance centroid of the maximum value of the PSF deviating from the chief ray.
In some embodiments, the step of calculating the ellipticity and the ellipticity component value of the point spread function includes:
obtaining the coordinate of the maximum PSF value according to the coordinate XY of the principal ray on the image surface and the centroid of the distance of the maximum PSF value deviating from the principal ray;
the coordinates of the other points are deduced by the suitable sampling interval, and further, within a circle of radius 0.5 ", the ellipsometry e of the point spread function PSF and the two components of the ellipsometry e1, e2 are calculated for each field of view according to the following formula,
wherein x and y are each value pair in PSFThe coordinates of the response are determined,
for PSF centroid coordinates, e is the ellipsometry value e1, e2 is the two components of the ellipsometry.
In some embodiments, the step of obtaining the stability of the ellipse according to the ellipse and the ellipse component values specifically includes the following steps:
fitting the ellipsoids and ellipsometric component values;
obtaining the ellipticity and the ellipticity component of the middle field of view through interpolation;
and obtaining the real value of the middle view field through optical software CodeV, and subtracting the real value from the data obtained by interpolation to obtain the ellipticity stability.
In some embodiments, the step of fitting the ellipsoids and the ellipsometry component values comprises:
the ellipsoids and ellipsometric stability were fitted with a binary cubic polynomial by MATLAB software.
In some embodiments, the step of obtaining the ellipticity and the ellipticity component of the intermediate field of view by interpolation specifically includes:
interpolation is performed in the software MATLAB by interpolation at the middle position of every four adjacent grid points to obtain the ellipticity e ' and the ellipticity component values e1 ', e2 ' of the middle field of view.
In addition, the invention also provides a system for computing and simulating the ellipticity stability of the optical remote sensing camera, which comprises the following components:
the optical model building unit is used for building an optical model of the optical remote sensing camera;
a sampling unit determining a suitable sampling interval of the point spread function;
a point spread function unit, which calculates the point spread function of a certain field of view in the optical model containing diffraction effect according to the suitable sampling interval;
the first calculating unit is used for calculating the ellipticity and the ellipticity component value of the point spread function;
and the second calculating unit is used for obtaining the stability of the ellipse according to the ellipse and the ellipse component value.
The method and the system for calculating and simulating the ellipticity stability of the optical remote sensing camera can comprehensively consider the influence of various error items on the ellipticity stability, realize the evaluation of related optical indexes and provide data support for astronomical research.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the application and together with the description serve to explain the application and not to limit the application. In the drawings:
fig. 1 is a flowchart of steps of a method for computing and simulating ellipticity stability of an optical remote sensing camera according to embodiment 1 of the present invention.
Fig. 2 is a schematic diagram illustrating an influence of a PSF sampling interval of the optical remote sensing camera provided in embodiment 1 of the present invention on ellipticity.
Fig. 3 is a schematic structural diagram of a simulation system for calculating the ellipsometry stability of an optical remote sensing camera according to embodiment 2 of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the technical solutions of the present application will be described in detail and completely with reference to the following specific embodiments of the present application and the accompanying drawings. It should be apparent that the described embodiments are only some of the embodiments of the present application, and not all of the embodiments.
Referring to fig. 1, the method for computing and simulating the ellipsometry stability of an optical remote sensing camera provided by the present application includes the following steps:
step S110: and establishing an optical model of the optical remote sensing camera.
Specifically, the optical model of the optical remote sensing camera includes an ideal optical system and an optical system in consideration of various error factors including, but not limited to: optical design residual error, processing surface shape error, system installation and adjustment residual error, error caused by gravity environment change, error caused by thermal environment change and the like.
Step S120: an appropriate sampling interval for the point spread function is determined.
Please give the concrete steps.
It can be appreciated that the above steps to determine the appropriate sampling interval of the Point Spread Function (PSF) ensure that under-sampling does not occur and the accuracy of the simulation computation is not affected.
Step S130: and calculating a point spread function of a certain field of view in the optical model containing diffraction effect according to the suitable sampling interval.
Specifically, an optical software CodeV is used for calculating a point spread function PSF of a certain field of view of the optical model, wherein the point spread function PSF comprises a diffraction effect, the coordinate XY of a tracing principal ray on an image plane, and the distance centroid of the maximum PSF deviating from the principal ray.
Step S140: and calculating the ellipticity and the ellipticity component value of the point spread function.
In the step of calculating the ellipticity and the ellipticity component value of the point spread function, the method specifically comprises the following steps:
obtaining the coordinate of the maximum PSF value according to the coordinate XY of the principal ray on the image surface and the centroid of the distance of the maximum PSF value deviating from the principal ray;
the coordinates of the other points are deduced by the suitable sampling interval, and further, within a circle of radius 0.5 ", the ellipsometry e of the point spread function PSF and the two components of the ellipsometry e1, e2 are calculated for each field of view according to the following formula,
wherein x and y are coordinates corresponding to each numerical value in the PSF,
for PSF centroid coordinates, e is the ellipsometry value e1, e2 is the two components of the ellipsometry.
It is understood that the coordinates of other points are obtained by linear superposition of the coordinates of the maximum value and the sampling interval, and the distance between the data in the PSF is the value of the sampling interval.
Step S150: and obtaining the stability of the ellipse according to the ellipse and the ellipse component value.
The step of obtaining the stability of the ellipse according to the ellipse and the ellipse component values specifically comprises the following steps:
step S151: fitting the ellipsoids and ellipsometric component values.
Further, by MATLAB software, fitting the ellipse ratio and the ellipse ratio stability by using a binary cubic polynomial avoids the high-order steep change or mutation characteristic in the fitting form with high-order.
Step S152: and obtaining the ellipticity and the ellipticity component of the middle field of view by interpolation.
Specifically, interpolation is performed at the middle position of every four adjacent grid points in the software MATLAB by interpolation, resulting in the ellipticity e ' and the ellipticity component values e1 ', e2 ' of the middle field of view.
Step S153: and obtaining the real value of the middle view field through optical software CodeV, and subtracting the real value from the data obtained by interpolation to obtain the ellipticity stability.
Fig. 2 is a schematic diagram illustrating an influence of a PSF sampling interval of the optical remote sensing camera 1 on ellipticity according to an embodiment of the present invention. In fig. 2, ellipse values at different PSF sampling intervals are calculated, and a suitable PSF sampling interval can be determined according to a certain calculation.
The method for calculating and simulating the ellipticity stability of the optical remote sensing camera comprises the steps of establishing an optical model of the optical remote sensing camera, determining a suitable sampling interval of a point spread function, calculating the point spread function of which a certain field of view contains a diffraction effect in the optical model according to the suitable sampling interval, calculating the ellipticity and the ellipticity component value of the point spread function, and obtaining the ellipticity stability according to the ellipticity and the ellipticity component value.
Example 2
Referring to fig. 3, a schematic structural diagram of a simulation system for calculating the ellipsometry stability of an optical remote sensing camera according to embodiment 2 of the present invention includes: an optical model construction unit 110, a sampling unit 120, a point spread function unit 130, a first calculation unit 140, and a second calculation unit 150. Wherein:
the optical model construction unit 110 is used for establishing an optical model of the optical remote sensing camera. Specifically, the optical model of the optical remote sensing camera includes an ideal optical system and an optical system in consideration of various error factors including, but not limited to: optical design residual error, processing surface shape error, system installation and adjustment residual error, error caused by gravity environment change, error caused by thermal environment change and the like.
The sampling unit 120 is used to determine the appropriate sampling interval for the point spread function.
The point spread function unit 130 is configured to calculate a point spread function of a field of view in the optical model that includes diffraction effects according to the suitable sampling interval.
Specifically, an optical software CodeV is used for calculating a point spread function PSF of a certain field of view of the optical model, wherein the point spread function PSF comprises a diffraction effect, the coordinate XY of a tracing principal ray on an image plane, and the distance centroid of the maximum PSF deviating from the principal ray.
The first calculating unit 140 is configured to calculate the ellipticity and the ellipticity component value of the point spread function, and specifically includes:
obtaining the coordinate of the maximum PSF value according to the coordinate XY of the principal ray on the image surface and the centroid of the distance of the maximum PSF value deviating from the principal ray;
the coordinates of the other points are deduced by the suitable sampling interval, and further, within a circle of radius 0.5 ", the ellipsometry e of the point spread function PSF and the two components of the ellipsometry e1, e2 are calculated for each field of view according to the following formula,
wherein x and y are coordinates corresponding to each numerical value in the PSF,
for PSF centroid coordinates, e is the ellipsometry value e1, e2 is the two components of the ellipsometry.
It is understood that the coordinates of other points are obtained by linear superposition of the coordinates of the maximum value and the sampling interval, and the distance between the data in the PSF is the value of the sampling interval.
The second calculating unit 150 is used for obtaining the stability of the ellipse according to the ellipse and the ellipse component values.
The method specifically comprises the following steps:
fitting the ellipsoids and ellipsometric component values. Further, by MATLAB software, fitting the ellipse ratio and the ellipse ratio stability by using a binary cubic polynomial avoids the high-order steep change or mutation characteristic in the fitting form with high-order.
And obtaining the ellipticity and the ellipticity component of the middle field of view by interpolation. Specifically, interpolation is performed at the middle position of every four adjacent grid points in the software MATLAB by interpolation, resulting in the ellipticity e ' and the ellipticity component values e1 ', e2 ' of the middle field of view.
And obtaining the real value of the middle view field through optical software CodeV, and subtracting the real value from the data obtained by interpolation to obtain the ellipticity stability.
The system for calculating and simulating the ellipticity stability of the optical remote sensing camera comprises an optical model for establishing the optical remote sensing camera, a suitable sampling interval of a point spread function is determined, the point spread function with a diffraction effect in a certain view field in the optical model is calculated according to the suitable sampling interval, the ellipticity and the ellipticity component value of the point spread function are calculated, and the ellipticity stability is obtained according to the ellipticity and the ellipticity component value.
The above are merely examples of the present application and are not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.
Claims (9)
1. A method for calculating and simulating the ellipticity stability of an optical remote sensing camera is characterized by comprising the following steps:
establishing an optical model of the optical remote sensing camera;
determining a suitable sampling interval for the point spread function;
calculating a point spread function of a certain field of view in the optical model containing diffraction effect according to the suitable sampling interval;
calculating the ellipticity and the ellipticity component value of the point spread function;
and obtaining the stability of the ellipse according to the ellipse and the ellipse component value.
2. The method for computational simulation of ellipsometry stability of an optical remote sensing camera according to claim 1, wherein in the step of establishing an optical model of the optical remote sensing camera, the optical model includes an ideal optical system and an optical system after considering error factors including but not limited to optical design residual, processing surface shape error, system setup residual, error caused by gravity environment change and thermal environment change.
3. The method for computation simulation of ellipsometry stability of an optical remote sensing camera according to claim 2, wherein in the step of determining a suitable sampling interval of the point spread function, the method specifically comprises:
according to the sampling theorem, the sampling interval range which can not be undersampled is determined, and then the ellipse ratio is calculated through the PSFs which are acquired at different sampling intervals, so that the proper sampling interval is determined.
4. The method for computation simulation of ellipsometry stability of an optical remote sensing camera according to claim 3, wherein in the step of computing a point spread function containing diffraction effect in a certain field of view in said optical model according to said suitable sampling interval, specifically comprises:
and calculating a point spread function PSF (point spread function) of a certain field of view of the optical model, including a diffraction effect, by using optical software CodeV through Fourier transform, tracking the coordinate XY of a chief ray on an image surface, and calculating the distance centroid of the maximum value of the PSF deviating from the chief ray.
5. The method for computing simulation of ellipsometry stability of optical remote sensing camera according to claim 4, wherein in the step of computing ellipsometry and ellipsometry component values of said point spread function, the method specifically comprises:
obtaining the coordinate of the maximum PSF value according to the coordinate XY of the principal ray on the image surface and the centroid of the distance of the maximum PSF value deviating from the principal ray;
the coordinates of the other points are calculated by said suitable sampling interval, and further within a circle of radius 0.5 ", the ellipsometry e of the point spread function PSF and the two components of the ellipsometry e1, e2 are calculated for each field of view according to the following formula,
wherein x and y are coordinates corresponding to each numerical value in the PSF,
for PSF centroid coordinates, e is the ellipsometry value e1, e2 is the two components of the ellipsometry.
6. The method for ellipsometric stability calculation and simulation of optical remote sensing camera of claim 5, wherein the step of obtaining the ellipsometric stability according to the ellipsometric and ellipsometric component values specifically comprises the following steps:
fitting the ellipsoids and ellipsometric component values;
obtaining the ellipticity and the ellipticity component of the middle field of view through interpolation;
and obtaining the real value of the middle view field through optical software CodeV, and subtracting the real value from the data obtained by interpolation to obtain the ellipticity stability.
7. The method for computation simulation of ellipsometry stability of optical remote sensing camera according to claim 6, wherein in said step of fitting said ellipsometry and ellipsometry component values, specifically:
the ellipsoids and ellipsometric stability were fitted with a binary cubic polynomial by MATLAB software.
8. The method for computing and simulating ellipsometry stability of an optical remote sensing camera according to claim 7, wherein in the step of obtaining ellipsometry and ellipsometry components of the intermediate field of view through interpolation, the method specifically comprises:
interpolation is performed in the software MATLAB by interpolation at the middle position of every four adjacent grid points to obtain the ellipticity e ' and the ellipticity component values e1 ', e2 ' of the middle field of view.
9. A simulation system for calculating the ellipticity stability of an optical remote sensing camera is characterized by comprising the following components:
the optical model building unit is used for building an optical model of the optical remote sensing camera;
a sampling unit determining a suitable sampling interval of the point spread function;
a point spread function unit, which calculates the point spread function of a certain field of view in the optical model containing diffraction effect according to the suitable sampling interval;
the first calculating unit is used for calculating the ellipticity and the ellipticity component value of the point spread function;
and the second calculating unit is used for obtaining the stability of the ellipse according to the ellipse and the ellipse component value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011048679.3A CN114329872A (en) | 2020-09-29 | 2020-09-29 | Method and system for calculating and simulating ellipticity stability of optical remote sensing camera |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011048679.3A CN114329872A (en) | 2020-09-29 | 2020-09-29 | Method and system for calculating and simulating ellipticity stability of optical remote sensing camera |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114329872A true CN114329872A (en) | 2022-04-12 |
Family
ID=81011437
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011048679.3A Pending CN114329872A (en) | 2020-09-29 | 2020-09-29 | Method and system for calculating and simulating ellipticity stability of optical remote sensing camera |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114329872A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117036962A (en) * | 2023-10-08 | 2023-11-10 | 中国科学院空天信息创新研究院 | Remote sensing image change detection method, device, equipment and storage medium |
-
2020
- 2020-09-29 CN CN202011048679.3A patent/CN114329872A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117036962A (en) * | 2023-10-08 | 2023-11-10 | 中国科学院空天信息创新研究院 | Remote sensing image change detection method, device, equipment and storage medium |
| CN117036962B (en) * | 2023-10-08 | 2024-02-06 | 中国科学院空天信息创新研究院 | Remote sensing image change detection method, device, equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102183214B (en) | Light detection method for large-aperture aspherical mirror structure | |
| JP7406853B2 (en) | Computer-implemented normal vector-based point cloud smoothing filtering method | |
| Wieneke | Improvements for volume self-calibration | |
| CN112967342B (en) | High-precision three-dimensional reconstruction method and system, computer equipment and storage medium | |
| CN110163905A (en) | A kind of high-precision ellipse fitting method | |
| CN113393577B (en) | Oblique photography terrain reconstruction method | |
| EP4127457A1 (en) | System and method for wind flow turbulence measurement by lidar in a complex terrain | |
| CN114943774B (en) | Temperature compensation calculation method for infrared binocular camera | |
| CN115775303A (en) | High-reflectivity object three-dimensional reconstruction method based on deep learning and illumination model | |
| CN114329872A (en) | Method and system for calculating and simulating ellipticity stability of optical remote sensing camera | |
| Zhang et al. | Improving the extracting precision of stripe center for structured light measurement | |
| Zeng et al. | An error similarity-based error-compensation method for measurement in the nonuniform temperature field | |
| CN115272110A (en) | Projector distortion correction method and device in structured light three-dimensional reconstruction | |
| JP2018179577A (en) | Position measuring device | |
| CN107255458B (en) | Resolving method of vertical projection grating measurement simulation system | |
| CN112985258A (en) | Calibration method and measurement method of three-dimensional measurement system | |
| CN110108271B (en) | Star deflection compensation method caused by pneumatic optical effect | |
| Li et al. | Posed relationship calibration with parallel mirror reflection for stereo deflectometry | |
| Fan et al. | Near-field photometric stereo using a ring-light imaging device | |
| CN114119747B (en) | Three-dimensional flow field flow display method based on PMD wave front detection | |
| CN109342008A (en) | Model in wind tunnel angle of attack one camera video measuring method based on homography matrix | |
| CN115127473A (en) | Large-caliber optical element surface shape detection device and detection image distortion correction method | |
| Cai et al. | High-precision projection moiré measurement method based on virtual reality bridging | |
| Hou et al. | High-precision visual imaging model and calibration method for multi-depth-of-field targets | |
| Jing et al. | Path integration guided with a quality map for shape reconstruction in the fringe reflection technique |
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 |