New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
DM-11456 revise jointcal doc and write example for new photometry model #46
Conversation
e8bcbb7
to
838868c
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
A collection mostly of wording suggestions, but a couple of places clarification would be helpful.
\begin{document} | ||
|
||
\maketitle | ||
|
||
\abstract{The jointcal package aims at optimizing simultaneously |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Suggest: "The jointcal package simultaneously optimizes... images. In principle and often in practice, this approach produces..."
simultaneous astrometry functionality is fulfilled by ``SCAMP''. The | ||
code we describe here has similar aims, but follows a slightly | ||
different route. It is meant to be used within the LSST software | ||
stack framework. } |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is "stack framework" a standard phrasing? Or is it just "LSST software stack?"
doc/jointcal.tex
Outdated
plotting tools. | ||
optimal in a statistical sense. | ||
|
||
This package performs a least-square fit to a set of images. Since |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
least-squares
doc/jointcal.tex
Outdated
This package performs a least-square fit to a set of images. Since | ||
it aims at statistical optimality, we maximize the likelihood of the | ||
measurements with respect to all unknown parameters required to describe | ||
the data. These parameters consists mostly in two sets: the position |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Suggest: These parameters are mostly in two sets: the position on the sky of the objects in common between images, and...
doc/jointcal.tex
Outdated
ours: it does not optimize over the position of common objects but | ||
rather minimizes the distance between pairs of transformed measurements | ||
of the same object. This approach is not a maximum likelihood optimization, | ||
and hence is likely sub-optimal statistics wise. We do not know how |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
is likely statistically sub-optimal.
doc/jointcal.tex
Outdated
mapping scaling, which is very well determined from the beginning. | ||
However, small changes of scaling can lead to the $\chi^{2}$ increasing | ||
between iterations. This is why we provide the \ClName{AstrometryModel::freezeErrorScales} | ||
which allows one to uses for there on the \textit{current} state of |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
"to uses for there": needs to be rehashed.
doc/jointcal.tex
Outdated
is in practice obtained using the routine \ClName{Gtransfo::transformPosAndErrors} | ||
which is attached to the projector. The derivative of the reference | ||
residual $R_{j}^{r}$ with respect to the \ClName{FittedStar} position | ||
$F_{j}$ (see eq. \ref{eq:ref_res}), is just the $2\times2$ matrix |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
equation reference was broken in my build; please verify
doc/jointcal.tex
Outdated
and eventually to offset parameters in the \RoutineName{AstrometryFit::OffsetParams}. | ||
|
||
Since the indexing of parameters is done dynamically, it is straightforward | ||
to only fit a subset of parameters. this is why the routine \RoutineName{AstrometryFit::AssignIndices} |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
capitalize This
doc/jointcal.tex
Outdated
|
||
\section{Fitting the transformations between a set of images} | ||
|
||
Some applications require to determine transformations between images |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
suggest: require determination of transformations
doc/jointcal.tex
Outdated
one just avoids to call \RoutineName{Associations::CollectRefStars}. | ||
\item One chooses polynomial mappings for all Calexp but one which will | ||
serve as a reference and have a fixed identity mapping. The distortion | ||
model \ClName{SimplePolyModel} allows to do that. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
allows that.
1054173
to
1bb0d83
Compare
Move images into separate directory fix image link in original TeX
Exported new version of TeX from LyX
19b0fdb
to
33e1581
Compare
I imported the text into LyX to edit it, moved some sections around, and generalized most of the math involving the minimization. Primary changes are to the abstract and sections 3.1, 3.2, 3.6, 3.7, 3.8., plus moving the Associations section after the math. Many other sections need work (especially now that photometry is implemented), but the above changes will help me implement the new photometry model, plus they I feel they make the math overall much clearer.