CN105023281B - Asterism based on point spread function wavefront modification is as centroid computing method - Google Patents
Asterism based on point spread function wavefront modification is as centroid computing method Download PDFInfo
- Publication number
- CN105023281B CN105023281B CN201510272682.6A CN201510272682A CN105023281B CN 105023281 B CN105023281 B CN 105023281B CN 201510272682 A CN201510272682 A CN 201510272682A CN 105023281 B CN105023281 B CN 105023281B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- asterism
- mtd
- mtr
- 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
Abstract
The present invention relates to a kind of asterism based on point spread function wavefront modification as centroid computing method, including：Gather multiframe star chart；In the star chart, every asterism pixel as where peak value is found, and a rectangular area is selected centered on asterism pixel as where peak value, this rectangular area includes asterism picture overwhelming majority energy, this region is determined that rectangular area ranks pixel count is odd number by measurement accuracy demand；In every asterism in the rectangular area as where, the k rank squares of every asterism picture are calculated with the method for photon weighted average；According to k ranks square and one and the poor directly proportional deviation of the k ranks square and peak pixel coordinate, the real centroid coordinate of asterism picture is obtained.
Description
Technical field
The present invention relates to astronomy and field of space technology, more particularly to a kind of star based on point spread function wavefront modification
Point is as centroid computing method.
Background technology
Uranometry measures the position of celestial body, voluntarily and parallax by measuring coordinate of the celestial body on celestial sphere, and it is to day
Literature and astrophysics have basic meaning, and the basis of geodesic survey, terrestrial positioning and navigation.Star sensor is that astronomy is led
Critical equipment equipment in boat system, high-precision data can be provided for Spacecraft Attitude Control.It is by measuring fixed star
Spacecraft three-axis attitude is exported in the transformation relation of star sensor coordinate system and geocentric coordinate system.The center coordination precision of asterism picture
Directly determine position of heavenly body measurement accuracy and spacecraft attitude measurement accuracy.
Currently, in uranometry field and star sensor field, center coordination is more using directly calculating center of gravity method person solution
The methods of analysing Function Fitting.Directly calculating gravity model appoach has the advantages that amount of calculation is small, calculating speed is fast, easily operates, but equally
Also there is the defects of measurement accuracy is low.Analytic Functions Fitting method has artificially been carried out it is assumed that leading to the shape of point spread function
It is difficult further to improve to cause final measurement accuracy.Therefore, calculating speed can how be accelerated, and can considers point spread function
True form, turn into professional improve method for positioning mass center target.
The content of the invention
The defects of it is an object of the invention to overcome existing method for positioning mass center measurement accuracy low, so as to provide a kind of energy
Effectively improve the method for positioning mass center of precision.
To achieve these goals, the invention provides a kind of asterism based on point spread function wavefront modification as barycenter meter
Calculation method, including：
Step 101, collection multiframe star chart；
Step 102, in the star chart, find every asterism as pixel where peak value, and with asterism as where peak value picture
A rectangular area is selected centered on element, this rectangular area includes asterism picture overwhelming majority energy, and this region is by measurement essence
Degree demand determines that rectangular area ranks pixel count is odd number；
Step 103, in every asterism obtained by step 102 in the rectangular area as where, with photon weighted average
Method calculates the k rank squares of every asterism picture；
Step 104, according to the k ranks square obtained by step 103 and one and the difference of the k ranks square and peak pixel coordinate
Directly proportional deviation, obtain the real centroid coordinate of asterism picture.
In above-mentioned technical proposal, in a step 101, in addition to the multiframe star chart to being gathered pre-processes, described pre-
Processing includes deducting dark noise, background noise, and carries out flat field amendment.
In above-mentioned technical proposal, in a step 101, the star chart gathered is true star chart, or utilizes optics in laboratory
The star chart of Imaging Simulation.
In above-mentioned technical proposal, in step 103, the k rank squares of every asterism picture are calculated with the method for photon weighted average
Calculation formula it is as follows：
Wherein, k is natural number, (rg, cg) it is the k rank squares that photon weighted average calculates, r represents row, and c represents row；Detection
Device array is m × n, and (i, j) is the ranks index of pel array, IijThe gray value arranged for the i-th row, jth.
In above-mentioned technical proposal, at step 104, the real centroid coordinate of asterism picture is calculated using equation below：
Wherein, l, N are natural number, (ra, ca) it is the real centroid coordinate that calculates after photon weighted average amendment, (krl,
kcl) be ranks direction the wavefront modification factor, (pr, pc) be ranks direction peak pixel coordinate.
In above-mentioned technical proposal, at step 104, the real centroid coordinate of asterism picture is obtained by least square fitting,
Specifically include：
Step 104-1, the star chart generation such as minor function vector by being gathered：
Wherein,
mRGil=(RGi-PRi)l,
mCGil=(CGi-kPCi)l,
Wherein, l, N are natural number, and the star chart gathered has M frames, have N number of asterism picture on every frame star chart；(RAi,CAi) for the
The real centroid coordinate of i asterism picture, (RGi,CGi) for the k rank squares of i-th asterism picture, (kRil,kCil) it is i-th asterism picture
The wavefront modification factor, (mRGil,mCGil) for photon weighted average k ranks square and the peak pixel coordinate of i-th asterism picture
Difference；(dRAi,dCAi) where i-th asterism picture two field picture relative to the translation of the first two field picture, spin matrix RiFor i-th
Two field picture where asterism picture relative to the first two field picture spin matrix；Wherein spin matrix RiFor
θiThe anglec of rotation for a two field picture relative to the first two field picture；
Step 104-2, generating parameter vector to be fitted by the function vector is：
Step 104-3, each variable in parameter vector to be fitted is substituted into following iterative solution, obtains every star
The real centroid coordinate and wavefront modification factor of point picture, translation vector and rotating vector：
xnew=xold+δxi；
Wherein, δ xi=-(Jij TJij)-1Jij TFj；J=1,2 ..., 2M+3N-1.
The advantage of the invention is that：
The method of the present invention, which has, calculates the advantages of simple, calculating speed is fast and result precision is high.
Brief description of the drawings
Fig. 1 (a) is the schematic diagram of wave front aberration；
Fig. 1 (b) is the schematic diagram of point spread function；
Fig. 2 is the asterism of the present invention as the flow chart of centroid computing method.
Fig. 3 is experimental verification measurement result figure of the present invention.
Embodiment
In conjunction with accompanying drawing, the invention will be further described.
In astronomical telescope imaging or star sensor imaging process, due to wave front aberration (referring to Fig. 1 (a)) presence,
So that the point spread function of fixed star is not preferable Airy disk, but asymmetrical Airy disk in irregular shape is (referring to Fig. 1
(b), wherein RMS=λ/20).The gravity model appoach for center coordination calculates point spread function using photon weighted average in the prior art
Several barycenter, imply a precondition：Assuming that the distribution of point spread function is preferable Airy disk or Gaussian Profile, do not have
In view of the influence of wavefront deviation.Therefore center coordination precision is difficult further to improve.
The present invention takes into full account influence of the optical system wavefront distortion to point spread function shape, proposes a kind of based on point expansion
The asterism of function wavefront modification is dissipated as centroid computing method, and on the basis of photon weighted average, wavefront is obtained by measured data
Modifying factor.
With reference to figure 2, the asterism of the invention based on point spread function wavefront modification as centroid computing method specifically include with
Lower step：
Step 1), utilize imaging sensor collection multiframe star chart；The star chart gathered can be true star chart, or test
Room utilizes the star chart of optical imagery simulation；
Step 2), the star chart obtained to step 1) pre-process, and deduct dark noise, background noise, and carry out flat field and repair
Just, the star chart after being pre-processed；
Step 3), using the star chart after the pretreatment obtained by step 2), find every asterism pixel as where peak value,
And a rectangular area is selected centered on pixel where peak value, this rectangular area includes asterism picture overwhelming majority energy, square
Shape region ranks pixel count is odd number；Wherein, the size of rectangular area determines truncated error, and the selection in this region needs basis
Specific measurement accuracy determines.
Step 4), using rectangular area of the every asterism as where obtained by step 3), with the side of photon weighted average
Method calculates the k rank squares of every asterism picture；
Step 5), the real centroid coordinate of hypothesis asterism picture are that the k ranks square that step 4) obtains adds a deviation, and this is inclined
Difference is poor directly proportional with k ranks square and peak pixel coordinate, then obtains the real centroid seat of asterism picture using least square fitting
Mark.
Each step in the inventive method is described further below.
The star chart frame number of collection is set to M in step 1), and the number of asterism picture is N, in order that the picture that asterism picture is covered
Prime number amount is as more as possible, typically using defocus or optimization Optical System Design, makes asterism as over-sampling.
In step 4), if detector array is classified as m × n, (i, j) be pel array ranks index, IijFor the i-th row, jth
The gray value of row.Equation below can be used to calculate photon weighted average k rank squares：
Wherein, k is natural number, (rg, cg) it is the k rank squares that photon weighted average calculates, r represents row, and c represents row.
In step 5), it is assumed that the real centroid coordinate of asterism picture is that k ranks square adds a deviation, and this deviation and k
Rank square and peak pixel coordinate it is poor directly proportional, be formulated, i.e.,：
Wherein, l, N are natural number, (ra, ca) it is the real centroid coordinate that calculates after photon weighted average amendment, (krl,
kcl) be ranks direction the wavefront modification factor, (pr, pc) be ranks direction peak pixel coordinate.
Next the specific steps of the real centroid coordinate of fitting asterism picture are provided.
Typically, for L variable xiThe group of functions of (i=1,2 ..., L)：
Fi(x1,x2,…,xL)=0, i=1,2 ..., L (5)
Wherein, xiThe vector formed for all variables, FiThe vector formed for all functions.In xiInfinitesimal area in,
Group of functions has following Taylor expansion：
Known Jacobian matrix is：
Ignore the quadratic component and higher order term in formula (6), and make Fi(xi+δxi)=0, convolution (6) can obtain with formula (7)
To a system of linear equations：
Jijδxi=-Fj (8)
Then the solution of infinitely small variable quantity is obtained, i.e.,
δxi=-(Jij TJij)-1Jij TFj (9)
Therefore, variable xiThere is iterative solution：
xnew=xold+δxi (10)
Specific to the M frame star charts of the collection in step 1), N number of asterism picture, just like minor function vector：
Wherein,
mRGil=(RGi-PRi)l (12)
mCGil=(CGi-PCi)l (13)
Wherein, (RAi,CAi) it is real centroid coordinate, (RGi,CGi) it is k rank squares, (kRil,kCil) it is the wavefront modification factor,
(mRGil,mCGil) for the difference of photon weighted average k ranks square and peak pixel coordinate, (PRi, PCi) be ranks direction peak value picture
Plain coordinate.
Here it is considered that actual conditions, introduced in formula (11) each two field picture relative to the translation of the first two field picture and
Rotation, i.e. (dRAi,dCAi) and Ri, wherein spin matrix RiFor
It would therefore be desirable to the parameter vector of fitting is
Therefore, can be with using formula (10) (after substituting into formula (10), j span is j=1,2 ..., 2M+3N-1)
Obtain the unknown parameters such as the real centroid coordinate and wavefront modification factor, the translation vector and rotating vector of every asterism picture.
Experimental verification is carried out using the method for the present invention, the measurement result obtained by the method for the present invention is existing with using
K ranks moments method, the measurement result of Gauss curve fitting method in technology are compared.As shown in figure 3, it can be seen that from measurement result
The standard deviation of the method method more of the prior art of the present invention is smaller, and center coordination precision is significantly improved.
It should be noted last that the above embodiments are merely illustrative of the technical solutions of the present invention and it is unrestricted.Although ginseng
The present invention is described in detail according to embodiment, it will be understood by those within the art that, to the technical side of the present invention
Case is modified or equivalent substitution, and without departure from the spirit and scope of technical solution of the present invention, it all should cover in the present invention
Right among.
Claims (6)
1. a kind of asterism based on point spread function wavefront modification is as centroid computing method, including：
Step 101, collection multiframe star chart；
Step 102, in the star chart, find every asterism as pixel where peak value, and using asterism as where peak value pixel as
Center selects a rectangular area, and this rectangular area includes asterism picture overwhelming majority energy, and this rectangular area is by measurement essence
Degree demand determines that rectangular area ranks pixel count is odd number；
Step 103, in every asterism obtained by step 102 in the rectangular area as where, with the method for photon weighted average
Calculate the k rank squares of every asterism picture；
Step 104, according to the k ranks square obtained by step 103 and one with the difference of the k ranks square and peak pixel coordinate into just
The deviation of ratio, obtain the real centroid coordinate of asterism picture.
2. the asterism according to claim 1 based on point spread function wavefront modification exists as centroid computing method, its feature
In in a step 101, in addition to the multiframe star chart to being gathered pre-processes, and the pretreatment includes removing dark noise, sheet
Back noise, and carry out flat field amendment.
3. the asterism according to claim 1 or 2 based on point spread function wavefront modification is as centroid computing method, its feature
It is, in a step 101, the star chart gathered is true star chart, or utilizes the star chart of optical imagery simulation in laboratory.
4. the asterism according to claim 1 or 2 based on point spread function wavefront modification is as centroid computing method, its feature
It is, in step 103, the calculation formula that the k rank squares of every asterism picture are calculated with the method for photon weighted average is as follows：
<mrow>
<msub>
<mi>r</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>k</mi>
</msup>
<mi>i</mi>
<mo>/</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>k</mi>
</msup>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mi>g</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>k</mi>
</msup>
<mi>j</mi>
<mo>/</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>k</mi>
</msup>
</mrow>
Wherein, k is natural number, (rg, cg) it is the k rank squares that photon weighted average calculates, r represents row, and c represents row；Detector array
Be classified as m × n, (i, j) be pel array ranks index, IijThe gray value arranged for the i-th row, jth.
5. the asterism according to claim 1 or 2 based on point spread function wavefront modification is as centroid computing method, its feature
It is, at step 104, the real centroid coordinate of asterism picture is calculated using equation below：
<mrow>
<msub>
<mi>r</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<msub>
<mi>r</mi>
<mi>g</mi>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>k</mi>
<mrow>
<mi>r</mi>
<mi>l</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>l</mi>
</msup>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<msub>
<mi>c</mi>
<mi>g</mi>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>k</mi>
<mrow>
<mi>r</mi>
<mi>l</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>l</mi>
</msup>
</mrow>
Wherein, l, N are natural number, (ra, ca) it is the real centroid coordinate that calculates after photon weighted average amendment, (krl, kcl)
The respectively wavefront modification factor in ranks direction, (pr, pc) be ranks direction peak pixel coordinate, (rg, cg) it is photon weight
The k rank squares that average computation goes out.
6. the asterism according to claim 5 based on point spread function wavefront modification exists as centroid computing method, its feature
In at step 104, obtaining the real centroid coordinate of asterism picture by least square fitting, specifically include：
Step 104-1, the star chart generation such as minor function vector by being gathered：
<mrow>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mrow>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mrow>
<mi>C</mi>
<mi>i</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>RA</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>dRA</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>RG</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>kR</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>mRG</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>CA</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>dCA</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>CG</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>kC</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
<msub>
<mi>mCG</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mn>2</mn>
<mi>M</mi>
<mi>N</mi>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>2</mn>
<mi>M</mi>
<mi>N</mi>
</mrow>
Wherein,
mRGil=(RGi-PRi)l,
mCGil=(CGi-kPCi)l,
Wherein, l, N are natural number, and the star chart gathered has M frames, have N number of asterism picture on every frame star chart；(RAi, CAi) it is i-th
The real centroid coordinate of asterism picture, (RGi, CGi) for the k rank squares of i-th asterism picture, (PRi, PCi) for the row of i-th asterism picture
The peak pixel coordinate of column direction；(kRil, kCil) for the wavefront modification factor of i-th asterism picture, (mRGil, mCGil) it is i-th
The photon weighted average k ranks square of asterism picture and the difference of peak pixel coordinate；(dRAi, dCAi) it is two field picture where i-th asterism picture
Relative to the translation of the first two field picture, spin matrix RiIt is rotation of i-th asterism picture place two field picture relative to the first two field picture
Torque battle array；Wherein spin matrix RiFor
<mrow>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
θiIt is the anglec of rotation of i-th asterism picture place two field picture relative to the first two field picture；
Step 104-2, generating parameter vector to be fitted by the function vector is：
<mrow>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>R</mi>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>CA</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>kR</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>kC</mi>
<mrow>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>dRA</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>dCA</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>N</mi>
<mo>+</mo>
<mn>3</mn>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
<mo>&times;</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
Step 104-3, each variable in parameter vector to be fitted is substituted into following iterative solution, obtains every asterism picture
The real centroid coordinate and wavefront modification factor, translation vector and rotating vector：
xnew=xold+δxi；
xoldFor the x before iterationi, xnewFor the x after iterationi；
Wherein, δ xi=-(Jij TJij)-1Jij TFi；
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510272682.6A CN105023281B (en) | 2015-05-25 | 2015-05-25 | Asterism based on point spread function wavefront modification is as centroid computing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510272682.6A CN105023281B (en) | 2015-05-25 | 2015-05-25 | Asterism based on point spread function wavefront modification is as centroid computing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105023281A CN105023281A (en) | 2015-11-04 |
| CN105023281B true CN105023281B (en) | 2018-02-16 |
Family
ID=54413222
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510272682.6A Active CN105023281B (en) | 2015-05-25 | 2015-05-25 | Asterism based on point spread function wavefront modification is as centroid computing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105023281B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106570902B (en) * | 2016-11-04 | 2019-09-24 | 中国科学院国家空间科学中心 | The opposite centroid computing method of PSF obtained based on detector pixel response spectrum |
| CN106530351B (en) * | 2016-11-24 | 2019-05-17 | 中国科学院国家空间科学中心 | A kind of method for positioning mass center obtained based on image sensor pixel internal quantum efficiency |
| CN108507564B (en) * | 2018-06-29 | 2020-06-16 | 中国科学院国家空间科学中心 | Star sensor centroid positioning method based on point spread function fitting |
| CN109443381B (en) * | 2018-10-17 | 2022-05-20 | 北京遥感设备研究所 | Star sensor centroid precision self-adaptive compensation method |
| CN114964677B (en) * | 2022-05-27 | 2023-05-02 | 中国科学院西安光学精密机械研究所 | Device and method for correcting MTF test error caused by airflow disturbance and environmental vibration |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103913166A (en) * | 2014-03-12 | 2014-07-09 | 哈尔滨工程大学 | Star extraction method based on energy distribution |
| CN104316049A (en) * | 2014-10-28 | 2015-01-28 | 中国科学院长春光学精密机械与物理研究所 | High-precision and low-signal-to-noise-ratio elliptic star spot subdivision location method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7260456B2 (en) * | 2004-01-05 | 2007-08-21 | The Boeing Company | Pixel-frequency slews and filters for star data measurements |
-
2015
- 2015-05-25 CN CN201510272682.6A patent/CN105023281B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103913166A (en) * | 2014-03-12 | 2014-07-09 | 哈尔滨工程大学 | Star extraction method based on energy distribution |
| CN104316049A (en) * | 2014-10-28 | 2015-01-28 | 中国科学院长春光学精密机械与物理研究所 | High-precision and low-signal-to-noise-ratio elliptic star spot subdivision location method |
Non-Patent Citations (2)
| Title |
|---|
| Micro-pixel accuracy centroid displacement estimation and detector calibration;Chengxing Zhai 等;《THE ROYAL SOCIETY》;20110803;摘要, 正文第7页第1段，第4.4节第1-3段，图2 * |
| 星点中心高精度质心定位算法;原玉磊 等;《测绘科学技术学报》;20121231;第29卷(第2期);摘要，正文第1节，第2节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105023281A (en) | 2015-11-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105023281B (en) | Asterism based on point spread function wavefront modification is as centroid computing method | |
| US8433515B2 (en) | Method for measuring precision of star sensor and system using the same | |
| Zhang | Star identification | |
| CN104462776B (en) | A kind of low orbit earth observation satellite is to moon absolute radiation calibration method | |
| Wei et al. | Star sensor calibration based on integrated modelling with intrinsic and extrinsic parameters | |
| CN104573251A (en) | Method for determining full-field-of-view apparent spectral radiance of satellite-borne optical remote sensor | |
| CN105528500A (en) | Imaging simulation method and system for decimeter-scale satellite-borne TDI CCD stereoscopic mapping camera | |
| CN106949907B (en) | A kind of quick system detection method of side of a ship window star | |
| CN110849331B (en) | Monocular vision measurement and ground test method based on three-dimensional point cloud database model | |
| CN102135435B (en) | Error correction method and device for digital sun sensor | |
| Del Genio et al. | Constraints on Saturn’s tropospheric general circulation from Cassini ISS images | |
| CN111581884A (en) | Satellite-borne multi-resolution target optical radiation image simulation method | |
| CN103234556A (en) | In-orbit calibration method of star sensor lens distortion and satellite attitude determination method based on starlight vector correction | |
| CN105444778A (en) | Star sensor in-orbit attitude determination error obtaining method based on imaging geometric inversion | |
| CN105758400A (en) | Extracting method for fixed star sensitive east and west parameters of stationary satellite imaging navigation and registration | |
| CN107451957A (en) | A kind of spaceborne TDI CMOS camera imagings emulation mode and equipment | |
| CN103591966A (en) | Star simulator test platform and test calibration method | |
| CN105486315B (en) | Remote sensing satellite is to moon absolute calibration attitude adjusting method | |
| CN105547286B (en) | A kind of compound three visual fields star sensor star map simulation method | |
| CN103743488B (en) | Infrared imaging simulation method for globe limb background characteristics of remote sensing satellite | |
| Hayashi et al. | Comparison of potential field solutions for Carrington Rotation 2144 | |
| CN108072387A (en) | A kind of base-line check method and system of the low in-orbit deviation of precision sensor | |
| CN107329190B (en) | Imaging test method for fixed star sensitivity of static meteorological satellite | |
| Feng et al. | Spherical cap harmonic analysis of regional magnetic anomalies based on CHAMP satellite data | |
| CN111366162B (en) | Small celestial body detector pose estimation method based on solar panel projection and template matching |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1 Patentee after: NATIONAL SPACE SCIENCE CENTER, CAS Address before: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1 Patentee before: Space Science & Applied Research Centre, Chinese Academy of Sciences |
|
| CP01 | Change in the name or title of a patent holder |