Dazzle: Oversampled Image Reconstruction and Difference-Imaging Photometry for the Nancy Grace Roman Space Telescope

Michael D Albrow

REVIEW 3 major objections 6 minor 23 references

Dazzle: Oversampled Image Reconstruction and Difference-Imaging Photometry for the Nancy Grace Roman Space Telescope

T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper presents Dazzle, a pipeline for Roman that builds an oversampled reference from dithered undersampled images and claims high recovery of microlensing events in the resulting difference images.

desk verdict Dazzle is a well-built, openly coded difference-imaging pipeline for Roman's undersampled WFI images; the simulations are honest, but the pure-translation assumption and a sign error in the offset-correction algebra need attention before relying on the recovery rates. read the letter →

arxiv 2412.06905 v1 pith:YQ3W4SMB submitted 2024-12-09 astro-ph.IM astro-ph.EPastro-ph.SR

Michael D Albrow This is my paper · ORCID

classification astro-ph.IMastro-ph.EPastro-ph.SR

keywords differenceimagingoversampledimagereconstructionRomanSpaceTelescopemicrolensingdetectionmatchedfilteringPSFphotometrycrowdedstellarfieldsdithering

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper aims to give the Roman Space Telescope's Galactic Bulge Time Domain Survey a difference-imaging route for transients, especially microlensing, that works when the individual images are undersampled and dithered rather than oversampled by seeing. The central move is to treat every pixel of every dithered frame as a direct sample of one global oversampled scene, solve for that scene pixel-by-pixel in a Legendre basis, and subtract it from each frame to make difference images. The author reports that on 192 simulated Roman images this recovers 90% or more of injected microlensing events with W146 < 23 and better than 80% for W146 < 25, and produces PSF-fitted light curves with a photometric scatter lower envelope of $\log_{10}\sigma[\mathrm{ppt}] = -5.19 + 0.29\, W146$. If these results hold on real data, Dazzle would be a practical complement to direct PSF photometry for faint and blended bulge sources.

What carries the argument

The load-bearing object is the per-pixel analytic oversampled reference: each reference-grid pixel $(i,j)$ is represented as $R_{ij}(x,y)=\sum_{l,m}\theta_{ijlm}B_l(x-x_i)B_m(y-y_j)$, where the $B$'s are Legendre polynomials stretched to overlap neighbouring pixels. The coefficients $\theta_{ijlm}$ are obtained from the weighted linear least-squares solution $\theta=(X^T C^{-1} X)^{-1} X^T C^{-1} z$, with $X$ built from the sub-pixel dither offsets and $z$ from the integer-shifted pixel values. This makes the reconstruction local, parallelizable, and free of image interpolation, and it is what turns the difference-image problem into a per-pixel linear algebra problem rather than a kernel-fitting problem. Detection then rides on a three-dimensional Gaussian matched filter over the aligned difference stack, with trial temporal widths to capture different event durations.

What would settle it

Apply Dazzle to a transient-free set of real Roman GBTDS frames and compare difference-image RMS with the photon-noise floor; if the lower envelope of scatter exceeds the paper's $\log_{10}\sigma[\mathrm{ppt}] = -5.19 + 0.29\, W146$ relation, or if injected events are recovered below the stated rates, the pure-translation and known-PSF assumptions are the failing point.

Watch

Extended reading notes

Core claim

The paper's claim is that a high-precision difference image can be formed for Roman without fitting a convolution kernel at all. Because the PSF is stable in time, the only difference between frames is a dither, so each measured pixel $T_{k,ij}$ equals the unknown oversampled scene $R$ sampled at a known sub-pixel phase: $T_{k,ij}=R(x_i+\Delta x_k+\delta x_k,y_j+\Delta y_k+\delta y_k)+\epsilon$. Representing $R$ inside each reference pixel as a product of shifted Legendre polynomials reduces the reconstruction to an independent linear least-squares problem per pixel, with a design matrix shared by all pixels. Difference images are then direct subtraction $D_{ij}=T_{k,ij}-R_{ij}$, with outlier pixels masked iteratively. The paper further derives a first-order analytic correction for small dither offset errors and demonstrates a three-dimensional matched filter over the difference-image stack that detects injected microlensing events at the stated recovery rates, followed by PSF-fitting photometry on the difference images with Nelder-Mead sub-pixel coordinate refinement.

Load-bearing premise

The argument assumes that the dithered frames are pure translations of a fixed scene with a known over-sampled PSF and no differential rotation; if real Roman frames carry residual geometric distortion, rotation, or PSF errors, the reference and every difference image inherit correlated errors before any transient search begins.

Editorial extensions

If this is right

For Roman's Galactic Bulge Time Domain Survey, faint and blended microlensing sources that are poorly served by direct PSF photometry can still be detected and measured through difference images.
Difference-image photometry reaches a bright-end floor near 0.03 mag scatter per point, indicating where direct photometry should take over.
The analytic dither-offset correction offers a way to refine small astrometric misalignments without re-solving the reference from scratch.
The three-dimensional matched filter can be tuned by trial temporal kernels to trade detection rate against false positives for different event timescales.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Beyond the paper's tests, the per-pixel Legendre construction should also work for other undersampled space telescopes with a stable PSF, provided the pure-translation and known-PSF assumptions hold; this is an extrapolation, not a paper claim.
The matched-filter detection, demonstrated for microlensing with a Gaussian temporal kernel, should transfer to any transient with a known time profile by swapping the kernel, making supernova or flare searches a natural next target.
The paper notes that Poisson error bars underestimate the observed scatter, which suggests that correlated noise from the reconstructed reference enters the difference images; modelling that covariance could sharpen both detection and photometric uncertainties.
Because the design matrix depends only on the dither pattern, the reference construction can be precomputed and parallelised per pixel, so the pipeline may scale to the full survey area if geometric distortion corrections stay within the linear regime.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Quick take: this is a serious, useful methods paper, not a breakthrough. It takes known ingredients—per-pixel Legendre-basis oversampled references, iterative outlier masking, 3D matched filtering, PSF fitting on difference images—and assembles them into a coherent, working pipeline for Roman's undersampled WFI data. The code is open-source, the validation uses injected synthetic microlensing events with known truth, and the paper explicitly flags that its Poisson error bars understate scatter. That last point is refreshing; too many pipeline papers hide that kind of thing.

The genuinely new bit is the combination for Roman: the oversampled reference built from randomly dithered undersampled frames, with an analytic correction for small dither errors, and the matched-filter detection stack. The recovery rates (90% for W146<23, >80% for W146<25) are plausible for the simulated geometry.

The soft spots, in order of importance. First, the core assumption is that dithers are pure translations of a fixed scene with a known oversampled PSF. Roman's WFI has real geometric distortion, so real dithered images are not pure translations. The paper does not discuss how distortion would be handled, and resampling to fix it changes the effective PSF. So the recovery rates are demonstrated under an idealized geometry and may not transfer directly. This is not a fatal flaw for a methods paper, but it needs to be addressed in the discussion, and validation on more realistic simulated distortion is needed. Second, the offset-correction algebra in Section 3.4 is internally inconsistent: Eq. 9 defines D = T - R, but Eq. 10/12 effectively use R - T. Likely a typo in the write-up, but as printed the sign story doesn't add up. Third, the tuning parameters (N=5, f=1.2, detection thresholds) are chosen on the same simulated data, so the recovery rates are somewhat optimistic. That is a limitation the authors acknowledge indirectly, but it would be good to see a cross-validation or at least a sensitivity test.

Bottom line: this paper deserves a serious referee. It is a solid, reproducible contribution to a real need—difference imaging for Roman's bulge survey. I'd ask the authors to fix the sign inconsistency, add a distortion discussion, and ideally show how the pipeline behaves with modest distortion injected. The Poisson error-bar caveat is fine for now. I would cite this work and bring it to the attention of anyone building Roman time-domain pipelines.

Referee Report

3 major / 6 minor

Summary. This manuscript presents Dazzle, an open-source Python package for constructing oversampled reference images from dithered, undersampled Roman WFI frames and for performing difference-image detection and photometry. The method represents each reference pixel as a low-order Legendre polynomial in the subpixel dither coordinate, solves for the coefficients using a design matrix assembled from the dithered images, iteratively masks outlier pixels, forms difference images, detects microlensing events with a three-dimensional Gaussian matched filter, and fits PSF photometry to the difference images. The validation is entirely simulation-based: SYNTHPOP and RomanISIM images of the Galactic bulge with injected PSPL microlensing events, recovery rates stated as above 90% for W146 < 23 and above 80% for W146 < 25, illustrative light curves, and a photometric scatter diagram. The paper is clearly written and the public code is a strength.

Significance. If the central claims hold, Dazzle would be a practical component of a Roman GBTDS transient pipeline, especially for faint and blended sources where direct PSF fitting is expected to struggle. The paper's strengths are its analytic, interpolation-free reference model, its injection-recovery validation against an external benchmark (the injected PSPL light curves), and the fact that the implementation and reproduction scripts are publicly available. However, two equation-level errors and the idealized translation-only geometry assumed in Eq. (2) currently limit the strength of the Roman-readiness claim; these issues are fixable but require new or corrected analysis.

major comments (3)

[Sec. 3.4, Eqs. (9)-(16)] There is an internal sign inconsistency in the offset-correction derivation. Equation (9) defines the difference image as D = T - R, but Eq. (10) writes the corrected difference image as D' = R(x+beta, y+gamma) - T, which is -D before the shift terms. Expanding Eq. (10) gives D' = -D + beta dR/dx + gamma dR/dy, yet Eq. (12) states D' = D + beta dR/dx + gamma dR/dy. If D is kept as T - R, the corrected residual is D - beta dR/dx - gamma dR/dy, so the minimisation in Eq. (13) has the wrong sign for the shift terms. As written, the normal equations (14)-(16) will return offset corrections of the wrong sign. This section should be rewritten with a single sign convention and tested on simulated frames with deliberately misregistered dither offsets.
[Sec. 3.6, Eqs. (17)-(18), and Fig. 7] Equation (18) is not the variance of the flux estimate in Eq. (17). For constant pixel noise sigma, Eq. (17) reduces to F = sum(P D)/sum(P^2), whose variance is sigma^2 / sum(P^2); Eq. (18) instead evaluates to sigma^2, independent of the PSF normalization. The correct variance of the weighted least-squares estimate is (sum P^2/sigma^2)^(-1). Because the paper uses these error bars to judge photometric scatter and then states that the Poisson error bars underestimate the scatter, this formula needs to be corrected and the comparison in Fig. 7 reassessed. If the PSF grid is normalized in a special way, that normalization must be stated explicitly.
[Secs. 3.1-3.5, Eq. (2)] The reference construction assumes a single global translation (integer plus subpixel) between each frame and the reference grid, with no differential rotation (Sec. 3.1). Roman's WFI has strong geometric distortion over its 4k x 4k field, so real dithered frames are not pure translations of a fixed scene. A field-dependent mapping violates Eq. (2), and the per-pixel Legendre fit can absorb slowly varying distortion into the reference, leaving correlated residuals that scale with local intensity gradients. The simulations do not appear to include or correct for this effect, and no test with differential rotation is reported. To support the claim that Dazzle is ready for Roman GBTDS data, the recovery rates in Sec. 3.5 should be recomputed on simulated images with a realistic distortion model, or with images resampled to a distortion-free frame and the ePSF appropriately transformed. At a minimum, the quoted recovery rates should be explicitly stated as applying to the idealized translation-only case.

minor comments (6)

[Secs. 2 and 3.6] The text contains duplicated words: 'the of the input stars' near the Fig. 2 discussion and 'by by' in Sec. 3.6; these should be corrected.
[Sec. 3.5] The recovery statistic is defined only as 'within 3 pixels' and the false-positive rate is described as 'a few' without a number. Reporting completeness and false-positive curves as a function of magnitude, u0, spatial kernel, temporal kernel, and detection threshold would make the detection claim quantitative and reproducible.
[Sec. 3.2] The choices N = 5 and f = 1.2 are justified only as working well 'from experimentation'; a short sensitivity test showing reference residuals or recovery rates versus N and f would make this tuning reproducible.
[Sec. 3.4] In Eqs. (15)-(16), the vector b should be written as b_k to make clear that the normal equations are solved separately for each image k, while the matrix A is common to all frames.
[Sec. 3.6] The statement that the Poisson error bars underestimate the scatter should be quantified, for example by plotting the ratio of observed scatter to formal error versus magnitude, rather than left as a qualitative remark.
[Sec. 1] The introduction suggests difference imaging will be superior to direct PSF fitting for faint, blended sources, but no side-by-side comparison on the same simulated frames is shown; such a comparison would substantiate the motivation.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the pipeline's detection and photometry claims are validated against externally injected synthetic events and models, not derived from fitted parameters.

full rationale

I examined the derivation chain for circularity of the types enumerated. The central construction (Sec. 3.2) solves for each reference-grid pixel's Legendre coefficients from the dithered pixel measurements via the design matrix solution in Eq. 7, minimizing Eq. 8. The coefficients are determined by the input images and assumed dither offsets; they are not defined in terms of any target detection or photometry result. The claimed recovery rates ('90% or better for sources with W146 < 23 and better than 80% for W146 < 25', Sec. 3.5) are explicitly measured against injected synthetic microlensing events with known parameters, i.e., an external benchmark. The photometry validation (Sec. 3.6, Figs. 5-6) compares the measured difference-image light curves to the known injected PSPL model, and the RMS-scatter comparison in Fig. 7 uses independent noise expectations from Wilson et al. (2023). No fitted parameter is renamed as a prediction. The self-citation to Albrow et al. (2009) is used only to justify the general need for accurate sub-pixel coordinates in PSF photometry, a standard external result, not to load-bearing support Dazzle's specific claims. Likewise, Anderson & King (2000) is cited for the definition of an effective PSF, and Gilliland et al. (2000) for the general per-pixel representation approach, neither of which forces the paper's conclusions. The acknowledged assumption of pure translational dithering with no differential rotation (Sec. 3.1) is a real limitation for Roman's geometric distortion, but that is a correctness/robustness concern, not a circularity: the derivation does not assume what it claims to predict. Therefore no circular step is present; the appropriate score is 0.

Assumptions & free parameters 6 free parameters · 6 assumptions · 0 invented entities

No invented physical entities. The central claim depends on standard image-processing assumptions plus hand-chosen tuning values (N, f, kernel widths, thresholds), all selected using the same simulations that produce the results. These enumerated choices are the main hidden degrees of freedom.

free parameters (6)

Legendre basis order N = 5
Chosen by experimentation for smooth over-sampled representation and clean difference images (Section 3.2).
extension factor f = 1.2
Chosen by experimentation to smooth across pixel boundaries (Section 3.2).
spatial Gaussian kernel sigma = 1.0 image units
Set to PSF width for the 3D matched filter (Section 3.5, Fig. 4 caption).
temporal Gaussian kernel sigmas = 2, 4, 8, 16 image units
Trial widths to capture varying microlensing timescales; injected events have FWHM in flux of about 3.5 units (Section 3.5).
detection threshold = not specified (configurable)
The threshold for local-maximum peaks depends on temporal kernel width and is left configurable; recovery rates depend on it (Section 3.5).
photometric scatter lower-envelope fit = log10 sigma[ppt] = -5.19 + 0.29 * W146
Linear fit to the lower envelope of RMS scatter for W146 > 23; descriptive calibration rather than a derived law (Section 3.5, Fig. 7).

assumptions (6)

domain assumption The dithered frames are related by pure translations with no differential rotation.
Section 3.1 states rotation-free linear translations; if violated by astrometric residuals, the shared pixel-space reference is biased.
domain assumption A known over-sampled effective PSF can be evaluated at arbitrary image locations.
Section 3.1, used in PSF fitting and in constructing the reference; real ePSF uncertainties are not modeled.
domain assumption The dither pattern well-samples the sub-pixel space so that 25 coefficients per pixel are constrained.
Section 3.2 requires at least 25 images and in practice 75 with random dithers; poor sub-pixel coverage would leave the reference ill-conditioned.
domain assumption Pixel noise is Gaussian with known variances and a diagonal covariance matrix.
Sections 3.2-3.3; the paper acknowledges outlier pixels and masks them iteratively, but cosmic rays and saturation are treated as exceptions.
ad hoc to paper A per-pixel Legendre basis of order 5 with extension factor 1.2 spans the local scene.
Chosen from experimentation on the same simulated images; no uniqueness or convergence argument is given.
standard math The first-order Taylor expansion in the offset correction is valid for small dither errors.
Section 3.4 linearizes the reference shift; the paper notes it works only for small offsets and a small number of discrepant images.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Dazzle: Oversampled Image Reconstruction and Difference-Imaging Photometry for the Nancy Grace Roman Space Telescope." pith.science (2026). https://pith.science/paper/YQ3W4SMB

@misc{pith2026241206905,
  author       = {Pith},
  title        = {Pith review of: Dazzle: Oversampled Image Reconstruction and Difference-Imaging Photometry for the Nancy Grace Roman Space Telescope},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YQ3W4SMB}},
  note         = {Machine review of arXiv:2412.06905}
}

read the original abstract

We present algorithms and software for constructing high-precision difference images to detect and measure transients, such as microlensing events, in crowded stellar fields using the Nancy Grace Roman Space Telescope. Our method generates difference images by subtracting an over-sampled reference, with iterative masking to address outlier pixels. We also provide an analytic correction for small dither offset errors. Microlensing event detection is achieved through a three-dimensional matched-filtering technique, optimized with Gaussian kernels to capture varying event durations, and verified through synthetic tests with high recovery rates. Transient photometry is performed via PSF fitting on difference images, using Nelder-Mead optimization for sub-pixel accuracy. The software, Dazzle, is available as an open-source Python package built on widely used libraries, offering accessible tools for the detection and characterization of transient phenomena in crowded fields.

Figures

Figures reproduced from arXiv: 2412.06905 by the authors.

**Figure 2.** Luminosity function of stars in the SYNTHPOP catalogue used as input to RomanISIM. There are a total of 10.5 million stars in each of the simulated images. A Hess diagram and W146 (AB magnitudes) luminosity function of the of the input stars are shown in [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗

**Figure 3.** Sample 80 x 80 region of a simulated Roman image of the Galactic Bulge. Left panel: raw simulated image. Right panel: image at ten times the original sampling from a representation constructed from a stack of 192 randomly-dithered images. Panels have identical logarithmic intensity scaling and white pixels are saturated. 3.1. Assumptions We assume that our raw data consists of a set of images, Tk(i, j), of a single… view at source ↗

Figures from the paper (4 more)

**Figure 4.** Figure 4: Left panel: Sample 1k × 1k region showing the locations of injected microlensing events (black dots) and detected variables (coloured circles) using a gaussian kernel with σ = 4.0 in the temporal direction and σ = 1.0 in each of the spatial directions. The colour scale…

**Figure 5.** Figure 5: Four examples showing time series of sample 7 x 7 pixel stamp images centered on a star. In each case a PSPL microlensing light curve with parameters u0 = 0.1, t0 = 2160, tE = 400 minutes has been injected directly into the stellar flux before image simulation. The top…

**Figure 6.** Figure 6: Complete light curves from PSF-fitting to the sets of simulated difference images shown in [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]

**Figure 7.** Figure 7: RMS photometric scatter, σ, in parts per thousand as a function of W146 magnitude. Each point corresponds to single light curve, with the scatter normalized to the baseline (unmagnified) flux in each case. The blue crosses are the total expected noise from the simula…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references · 2 canonical work pages

[1]

2000, A&AS, 144, 363, doi: 10.1051/aas:2000214

Alard, C. 2000, A&AS, 144, 363, doi: 10.1051/aas:2000214

doi:10.1051/aas:2000214 2000
[2]

Alard, C., & Lupton, R. H. 1998, ApJ, 503, 325, doi: 10.1086/305984

doi:10.1086/305984 1998
[3]

D., Gilliland, R

Albrow, M. D., Gilliland, R. L., Brown, T. M., et al. 2001, ApJ, 559, 1060, doi: 10.1086/322353

doi:10.1086/322353 2001
[4]

D., Horne, K., Bramich, D

Albrow, M. D., Horne, K., Bramich, D. M., et al. 2009, MNRAS, 397, 2099, doi: 10.1111/j.1365-2966.2009.15098.x

arXiv 2009
[5]

Anderson, J., & King, I. R. 2000, PASP, 112, 1360, doi: 10.1086/316632 Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167, doi: 10.3847/1538-4357/ac7c74

doi:10.1086/316632 2000
[6]

A., Abe, F., Dodd, R

Bond, I. A., Abe, F., Dodd, R. J., et al. 2001, MNRAS, 327, 868, doi: 10.1046/j.1365-8711.2001.04776.x

arXiv 2001
[7]

M., Horne, K., Albrow, M

Bramich, D. M., Horne, K., Albrow, M. D., et al. 2013, MNRAS, 428, 2275, doi: 10.1093/mnras/sts184

doi:10.1093/mnras/sts184 2013
[8]

L., et al

Bruntt, H., Frandsen, S., Gilliland, R. L., et al. 2001, A&A, 371, 614, doi: 10.1051/0004-6361:20010369

work page doi:10.1051/0004-6361:20010369 2001

Show all 23 references

[9]

A., Clayton, G

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245, doi: 10.1086/167900

1989 doi
[10]

2012, Computational Optimization and Applications, 51, 259

Gao, F., & Han, L. 2012, Computational Optimization and Applications, 51, 259

2012
[11]

L., Brown, T

Gilliland, R. L., Brown, T. M., Guhathakurta, P., et al. 2000, ApJL, 545, L47, doi: 10.1086/317334

2000 doi
[12]

R., Millman, K

Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357

2020
[13]

W., & Casey, A

Hogg, D. W., & Casey, A. R. 2024, arXiv e-prints, arXiv:2403.11011, doi: 10.48550/arXiv.2403.11011 Kl¨ uter, J., Huston, M. J., Aronica, A., et al. 2024, arXiv e-prints, arXiv:2411.18821, doi: 10.48550/arXiv.2411.18821

2024 doi arXiv
[14]

2006, A&A, 453, 635, doi: 10.1051/0004-6361:20053842

Picaud, S. 2006, A&A, 453, 635, doi: 10.1051/0004-6361:20053842

2006 doi
[15]

1965, The Computer Journal, 7, 308 O’Donnell, J

Nelder, J., & Mead, R. 1965, The Computer Journal, 7, 308 O’Donnell, J. E. 1994, ApJ, 422, 158, doi: 10.1086/173713

1965 doi
[16]

T., Gaudi, B

Penny, M. T., Gaudi, B. S., Kerins, E., et al. 2019, ApJS, 241, 3, doi: 10.3847/1538-4365/aafb69

2019 doi
[17]

D., Sivaramakrishnan, A., Lajoie, C.-P., et al

Perrin, M. D., Sivaramakrishnan, A., Lajoie, C.-P., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 9143, Space Telescopes and Instrumentation 2014: Optical, Infrared, and Millimeter Wave, ed. J. Oschmann, Jacobus M., M. Clampi...

2014 doi
[18]

C., Reyl´ e, C., Derri` ere, S., & Picaud, S

Robin, A. C., Reyl´ e, C., Derri` ere, S., & Picaud, S. 2003, A&A, 409, 523, doi: 10.1051/0004-6361:20031117

2003 doi
[19]

C., Casertano, S., Bond, H

Sahu, K. C., Casertano, S., Bond, H. E., et al. 2006, Nature, 443, 534, doi: 10.1038/nature05158

2006 doi
[20]

B., & Crotts, A

Tomaney, A. B., & Crotts, A. P. S. 1996, AJ, 112, 2872, doi: 10.1086/118228

1996 doi
[21]

E., et al

Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261, doi: 10.1038/s41592-019-0686-2

2020 doi
[22]

F., Barclay, T., Powell, B

Wilson, R. F., Barclay, T., Powell, B. P., et al. 2023, ApJS, 269, 5, doi: 10.3847/1538-4365/acf3df

2023 doi
[23]

Wozniak, P. R. 2000, AcA, 50, 421, doi: 10.48550/arXiv.astro-ph/0012143

2000 doi arXiv

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.

[1] [1]

2000, A&AS, 144, 363, doi: 10.1051/aas:2000214

Alard, C. 2000, A&AS, 144, 363, doi: 10.1051/aas:2000214

doi:10.1051/aas:2000214 2000

[2] [2]

Alard, C., & Lupton, R. H. 1998, ApJ, 503, 325, doi: 10.1086/305984

doi:10.1086/305984 1998

[3] [3]

D., Gilliland, R

Albrow, M. D., Gilliland, R. L., Brown, T. M., et al. 2001, ApJ, 559, 1060, doi: 10.1086/322353

doi:10.1086/322353 2001

[4] [4]

D., Horne, K., Bramich, D

Albrow, M. D., Horne, K., Bramich, D. M., et al. 2009, MNRAS, 397, 2099, doi: 10.1111/j.1365-2966.2009.15098.x

arXiv 2009

[5] [5]

Anderson, J., & King, I. R. 2000, PASP, 112, 1360, doi: 10.1086/316632 Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167, doi: 10.3847/1538-4357/ac7c74

doi:10.1086/316632 2000

[6] [6]

A., Abe, F., Dodd, R

Bond, I. A., Abe, F., Dodd, R. J., et al. 2001, MNRAS, 327, 868, doi: 10.1046/j.1365-8711.2001.04776.x

arXiv 2001

[7] [7]

M., Horne, K., Albrow, M

Bramich, D. M., Horne, K., Albrow, M. D., et al. 2013, MNRAS, 428, 2275, doi: 10.1093/mnras/sts184

doi:10.1093/mnras/sts184 2013

[8] [8]

L., et al

Bruntt, H., Frandsen, S., Gilliland, R. L., et al. 2001, A&A, 371, 614, doi: 10.1051/0004-6361:20010369

work page doi:10.1051/0004-6361:20010369 2001