A paper analogue computer for seven catalogue stars
The Sky You Can Hold
Turn the star disk now and watch a named star cross your horizon. Then manufacture the same geometry on paper, measure what your printer actually made, and take it outside with the computer switched off.
This is a planispheric astrolabe: a hybrid of a solid star disk and an astrolabe-style latitude plate. It is not a replica of a historical museum instrument. Its useful trick is older and stranger: spherical altitude circles become ordinary circles on paper.
The source audit comes before print
Ptolemy source values, in parts: winter tropic 92;8,15, summer tropic 39;4,19, ecliptic radius 65;36,17, ecliptic offset 26;31,58, Rhodes horizon radius 102;4,45, offset 82;35,3. The full recomputation and residuals appear below.
| HIP / name | RAICRS ° | DEICRS ° | Vmag | pmRA | pmDE |
|---|---|---|---|---|---|
| 11767 Polaris | 37.94614689 | +89.26413805 | 1.97 | +44.22 | -11.74 |
| 24436 Rigel | 78.63446353 | -8.20163919 | 0.18 | +1.87 | -0.56 |
| 27989 Betelgeuse | 88.79287161 | +7.40703634 | 0.45 | +27.33 | +10.86 |
| 32349 Sirius | 101.28854105 | -16.71314306 | -1.44 | -546.01 | -1223.08 |
| 65474 Spica | 201.29835230 | -11.16124491 | 0.98 | -42.50 | -31.73 |
| 69673 Arcturus | 213.91811403 | +19.18726997 | -0.05 | -1093.45 | -1999.40 |
| 91262 Vega | 279.23410832 | +38.78299311 | 0.03 | +201.02 | +287.46 |
All seven rows are catalogue epoch J1991.25 in ICRS. These are the exact numbers plotted. The printed precision is source disclosure, not instrument accuracy.
Move the sky
Choose where you are and a time. The seven points below come from the seven Hipparcos rows printed further down, without hidden replacement coordinates. Drag the disk itself, use the time control, or press an arrow key while the canvas is focused.
Computing the sky.
Gold: selected Hipparcos star. Blue: horizon. Thin curves: 10 degree altitude intervals. Dashed rim: the 85 mm mater crop. Stars below 15 degrees are marked low-confidence for the outdoor comparison.
Manufacture three sheets
Sheets 1 and 2 carry the kit's 150 mm ruler, plus a 100.00 mm horizontal line and a 50.00 × 50.00 mm square. Print at 100% / actual size. Turn off fit, shrink and scale-to-page options. Drivers commonly rescale without making the change obvious.
Measure before building
Not calibrated. This page cannot vouch for a sheet it has not seen measured in both axes.
Printing stays disabled until both scale checks pass. Print sheet 2 on transparent film made for your printer.
- Print and measure first. Do not cut a sheet that the verdict refuses.
- Optionally glue sheet 1 to cereal-box card. Sheet 2 must remain transparent: use printer-safe transparency film. Cut both 170 mm circles and pierce the marked centres. Do not cut away the altitude curves.
- Put the transparent plate over the star disk and join the centres with a pin pushed into a folded-tape pad, or use a split pin.
- Cut the quadrant. Tape a straight drinking straw to the sight line. Tie about 25 cm of thread at the vertex and tape a coin to its free end.
- Set the disk's numbered rim to the triangular index on sheet 2 at the printed mean local sidereal angle. That setting is valid only for the UTC instant and longitude printed on sheet 2.
First anchor: Ptolemy's numbers, left imperfect
The earliest extant description confirmed for this build is Ptolemy's Planisphere, surviving through an Arabic textual history. That does not settle who invented stereographic projection. Sidoli and Berggren's critical translation flags garbling and small discrepancies, so the table keeps the transmitted sexagesimal values beside modern evaluation instead of forcing a match.
R = 60 parts; ε = 23;51,20°; φRhodes = 36°
r(δ) = R tan((90° − δ)/2)
horizon radius = R / |sin φ|; centre offset = R cot φ
| quantity | printed source | decimal | modern equation | residual, parts |
|---|
Residual means modern equation minus the printed source value, in Ptolemy's parts. It is evidence of the textual and rounding seam, not observational error.
Second anchor: the seven rows that become ink
These are ESA Hipparcos Main Catalogue fields at catalogue epoch J1991.25 in ICRS. ICRS orientation is consistent with J2000; these are not therefore “J2000 positions.” Johnson V magnitude sizes only three coarse symbol classes. It does not predict tonight's perceived brightness, and some sources vary or are multiple systems.
| HIP / name | RAICRS ° | DEICRS ° | Vmag | pmRA mas/yr | pmDE mas/yr |
|---|
The disk uses exactly the displayed RA and declination. It does not propagate proper motion, precession, nutation or aberration. Even fast Arcturus moves only about 0.02° from 1991.25 to 2026 by the catalogue proper-motion scale, below a 1° paper tick. Refraction is also omitted from the prediction. Near the horizon it is variable and roughly 0.5° under standard conditions, so the field test excludes altitudes below 15°.
Why the paper can compute
The sophisticated dismissal is fair: a decorative star wheel can still look persuasive. This one earns its altitude curves from the complete circle calculation.
- Put the unit celestial sphere at the origin, north at Z = +1, and project from the south pole onto Z = 0. With printed equator radius R, a point on paper is q = R cos δ/(1 + sin δ), x = q cos α, y = q sin α.
- The inverse is X = 2Rx/(R²+ρ²), Y = 2Ry/(R²+ρ²), Z = (R²−ρ²)/(R²+ρ²), where ρ² = x²+y².
- Every circle on the sphere is the intersection with a plane AX+BY+CZ=D. Substitution and collection gives:
(C + D)(x² + y²) − 2R(Ax + By) − R²(C − D) = 0
When C + D is nonzero, completing the square gives an ordinary Euclidean circle. When C + D is zero, the quadratic term vanishes and the image is a straight line. So “circles map to circles” needs the line exception: generalized circles, including a line when the spherical circle passes through the projection point.
For signed latitude φ and altitude h, the plane of constant altitude becomes:
(sin φ + sin h)(x²+y²) − 2R cos φ x − R²(sin φ − sin h) = 0
centre ch = R cos φ/(sin φ + sin h)
radius ah = R cos h/|sin φ + sin h|
At the equator the horizon is the exact limiting line x = 0. Whenever sin φ + sin h = 0, the implementation draws the corresponding line rather than pretending an enormous circle fits the sheet. Rotation turns right ascension α into hour angle H = θ − α. The fixed altitude circles then return altitude. The theorem is the mechanism.
Take one number outside
Never sight the Sun through the straw. On a clear night choose a printed star predicted above 15°. Use the UTC instant printed on sheet 2; align the disk angle printed there with the triangular index. Read its altitude through the transparent plate.
- First sight a distant level horizon. The quadrant should read 0° within 1°. If not, fix the straw or thread before trusting a star.
- Sight the same star, let the weighted thread settle, and read where it crosses the 1° scale. Repeat three times.
- Compare the median sighting with the planispheric reading. The unvalidated design target is 2° and the declared trial pass band is 3°; neither is a measured performance claim. A miss does not translate to a universal clock error because altitude changes with latitude, declination and hour angle.
If the scale check refused, if the horizon zero misses by more than 1°, or if the star is below 15°, the page refuses the 3° comparison. Record the observation if useful, but do not call it a passed instrument test.
What the ring means
The browser computes Earth Rotation Angle directly from the IERS constants 0.7790572732640 revolutions at J2000.0 and 1.00273781191135448 revolutions per UT1 day, then adds the IERS GMST precession polynomial. UTC stands in for UT1, a below-one-second approximation here. Nutation and the equation of the equinoxes are omitted, so the implementation consistently uses mean sidereal time, not apparent sidereal time.
Sheet 2 prints the selected UTC instant, longitude, and the computed mean local sidereal angle. The paper has no calendar converter: changing the observation time requires generating that sheet again. Longitude matters one degree of hour angle per degree of longitude. Civil-zone and daylight-saving conversion is deliberately outside the artifact.
Sources and boundaries
- Nathan Sidoli and J. L. Berggren, “The Arabic Version of Ptolemy's Planisphere,” SCIAMVS 8 (2007), especially 4.3 and 7. Critical edition and translation.
- ESA, The Hipparcos and Tycho Catalogues, SP-1200 (1997), machine-readable Main Catalogue I/239. The exact seven-row query.
- IERS Conventions (2010), Technical Note 36, Earth rotation and Greenwich sidereal time. Technical Note 36.
No planets or solar sighting are included. Stereographic projection preserves angles locally, but not distance or area. Southern latitude instructions reverse the sky-facing east/west sense explicitly on the print. The layout refuses latitudes outside ±75° because near-polar plates require a different crop.
For a live planetary sky rather than this fixed-star paper mechanism, visit The Sky Above You. For why a catalogue epoch matters, visit The Star That Won't Stay North.