starplot 0.19.1__tar.gz → 0.19.2__tar.gz

{starplot-0.19.1 → starplot-0.19.2}/CITATION.cff

@@ -20,5 +20,5 @@ keywords:
   - graphs
   - plotting
 license: MIT
-version: 0.19.1
-date-released: '2026-01-24'
+version: 0.19.2
+date-released: '2026-01-27'

{starplot-0.19.1 → starplot-0.19.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: starplot
-Version: 0.19.1
+Version: 0.19.2
 Summary: Star charts and maps of the sky
 Keywords: astronomy,stars,charts,maps,constellations,sky,plotting
 Author-email: Steve Berardi <hello@steveberardi.com>

{starplot-0.19.1 → starplot-0.19.2}/docs/changelog.md

@@ -1,10 +1,15 @@
 ## v0.19.x
+
 - Adds a `CollisionHandler` for defining what to do when labels collide with something
 - Introduces catalogs for the Milky Way and constellation borders
 - Makes map plots default to plotting at epoch J2000
 - Replaces the `true_size` kwarg on the DSO plotting function with `where_true_size` which is a list of filters
 - Fixes a bug with auto label offsets from markers
 - [**v0.19.1**] Fixes issue with constellation boundaries
+- [**v0.19.2**]
+    - Adds a `border` field to the constellation model
+    - Modifies the `line` function to accept a Shapely geometry OR a list of coordinates
+    - Fixes bug with plotting true apparent size of solar system objects
 
 ---
 

{starplot-0.19.1 → starplot-0.19.2}/docs/examples/map-canis-major.md

@@ -7,7 +7,7 @@ title: Map of Canis Major
 
 ![map-canis-major](/images/examples/map_canis_major.png)
 
-In this example, we set the plot's `clip_path` to the boundary of Canis Major, so only objects in that constellation are plotted:
+In this example, we set the plot's `clip_path` to the border of Canis Major, so only objects in that constellation are plotted:
 
 ```python
 --8<-- "examples/map_canis_major.py"

{starplot-0.19.1 → starplot-0.19.2}/examples/map_canis_major.py

@@ -1,3 +1,5 @@
+from shapely import Polygon
+
 from starplot import MapPlot, Miller, Constellation, _
 from starplot.styles import PlotStyle, extensions
 
@@ -14,13 +16,16 @@ p = MapPlot(
     dec_max=-10.9,
     style=style,
     resolution=3400,
-    clip_path=canis_major.boundary,
+    clip_path=Polygon(canis_major.border.coords),
     scale=1.2,
 )
 p.constellations(
     where=[_.iau_id == "cma"],
 )
-p.constellation_borders()
+p.line(
+    geometry=canis_major.border,
+    style=p.style.constellation_borders,
+)
 p.open_clusters(where=[_.magnitude < 9], where_true_size=[False])
 p.stars(where=[_.magnitude < 9], where_labels=[_.magnitude < 4], bayer_labels=True)
 p.constellation_labels()

{starplot-0.19.1 → starplot-0.19.2}/examples/optic_iss_transit.py

@@ -26,7 +26,7 @@ p = OpticPlot(
     dec=moon.dec,
     observer=observer,
     optic=Binoculars(
-        magnification=15,
+        magnification=20,
         fov=65,
     ),
     style=style,

{starplot-0.19.1 → starplot-0.19.2}/examples/optic_solar_eclipse.py

@@ -5,7 +5,7 @@ from starplot import Moon, Binoculars, Observer
 from starplot.styles import PlotStyle, extensions
 
 # time of partial eclipse. total eclipse started at 15:13:46
-dt = datetime(2024, 4, 8, 14, 40, 47, 0, tzinfo=ZoneInfo("US/Eastern"))
+dt = datetime(2024, 4, 8, 14, 45, 0, 0, tzinfo=ZoneInfo("US/Eastern"))
 
 
 observer = Observer(
@@ -18,7 +18,7 @@ m = Moon.get(dt=observer.dt, lat=observer.lat, lon=observer.lon)
 
 op = m.create_optic(
     observer=observer,
-    optic=Binoculars(magnification=20, fov=65),
+    optic=Binoculars(magnification=30, fov=65),
     style=PlotStyle().extend(
         extensions.GRAYSCALE_DARK,
         extensions.OPTIC,

{starplot-0.19.1 → starplot-0.19.2}/hash_checks/hashlock.yml

@@ -22,6 +22,10 @@ map_coma_berenices_dso_size:
   dhash: ea92ecacb08cc0c0aa86acacb088c682aa8a9292b08ac282
   filename: /starplot/hash_checks/data/map-coma-berenices-dso-size.png
   phash: 92966d6d4d9616c9
+map_constellation_clip_path:
+  dhash: 3202ee96e47819103002ee96e4781910b002ee96e4781910
+  filename: /starplot/hash_checks/data/map-constellation-clip-path.png
+  phash: 99c45b99c4d98ec6
 map_custom_stars:
   dhash: 5e21a1ceac2c84485e21a1ceac2c84485e21a1ceac2c8448
   filename: /starplot/hash_checks/data/map-custom-stars.png
@@ -43,7 +47,7 @@ map_mollweide:
   filename: /starplot/hash_checks/data/map-mollweide.png
   phash: ebcbb698e3309461
 map_moon_phase_waxing_crescent:
-  dhash: 0814284d4d2814088894a84d4da894888884884d4d888488
+  dhash: 0008344a4a3408008088a04a4aa088808088804a4a808880
   filename: /starplot/hash_checks/data/map-moon-phase-waxing-crescent.png
   phash: b38ccc3333cccc33
 map_orion_base:
@@ -51,9 +55,9 @@ map_orion_base:
   filename: /starplot/hash_checks/data/map-orion-base.png
   phash: bf7952269584b594
 map_orion_extra:
-  dhash: 191b1b3f2b2f6d691b1b195b2b236d695c1b1b1e2b2f6d69
+  dhash: 191b1b3f2f2f6d691b1b1b5b2f236d695c1b1b1b2a2f6d69
   filename: /starplot/hash_checks/data/map-orion-extra.png
-  phash: bf61e20ec19ed09a
+  phash: bf61e606c59ed092
 map_plot_custom_clip_path_virgo:
   dhash: 23032913950c040443032507870c0c0203032507070e0c02
   filename: /starplot/hash_checks/data/map-custom-clip-path-virgo.png
@@ -91,9 +95,9 @@ optic_clipping:
   filename: /starplot/hash_checks/data/optic-clipping.png
   phash: 95696a87add86983
 optic_iss_moon_transit:
-  dhash: 70e8d6b2b28cf07071e896b2b68ef07171cc90b2a48ed071
+  dhash: 70d0ceb0f088f47071d08eb4f68ad47171c488b6b4aac471
   filename: /starplot/hash_checks/data/optic-iss-moon-transit.png
-  phash: cc993366cc993366
+  phash: c49d33624c9d3367
 optic_m45_binoculars:
   dhash: 8e17237133330e4d8e17237133330e4d8e17237133330e4d
   filename: /starplot/hash_checks/data/optic-m45-binoculars.png
@@ -115,15 +119,15 @@ optic_m45_scope_gradient:
   filename: /starplot/hash_checks/data/optic-m45-scope-gradient.png
   phash: aa3b84e48693ced3
 optic_moon_phase_full:
-  dhash: 0e2b553333552b0e0e2b553333552b0e0e2b553333552b0e
+  dhash: 0e334d71714d330e0e334d71714d330e0e334d71714d330e
   filename: /starplot/hash_checks/data/optic-moon-phase-full.png
-  phash: ff7aa02395308735
+  phash: ff7aa02a85318735
 optic_moon_phase_new:
-  dhash: 16088e96968e081616088e96968e081616088e96968e0816
+  dhash: 0814084d4d0814080814084d4d0814080814084d4d081408
   filename: /starplot/hash_checks/data/optic-moon-phase-new.png
-  phash: e5719e96690c6663
+  phash: e63399cc6633998c
 optic_moon_phase_waxing_crescent:
-  dhash: 0e33454d59712b0e0e33454d59712b0e0e33454d59712b0e
+  dhash: 0e3345495145330e0e3345495145330e0e3345495145330e
   filename: /starplot/hash_checks/data/optic-moon-phase-waxing-crescent.png
   phash: bb26e4999166c699
 optic_orion_nebula_refractor:
@@ -135,9 +139,9 @@ optic_polaris_binoculars:
   filename: /starplot/hash_checks/data/optic-binoculars-polaris.png
   phash: aa6b84f182f18ed3
 optic_solar_eclipse_binoculars:
-  dhash: 71ccb28ecc8ed47171ccb28ecc8ed47171c8a28ccc8ac071
+  dhash: 71cc928a8a92cc7171cc928a8a92cc7171c0b28a8a92c071
   filename: /starplot/hash_checks/data/optic-binoculars-eclipse.png
-  phash: 91cc4c333bcc6c73
+  phash: 91cc4e333bcc6c33
 optic_wrapping:
   dhash: 70f0ec92d2d8650a71e0d4acacc471aa71e0d4acacc4792a
   filename: /starplot/hash_checks/data/optic-wrapping.png

{starplot-0.19.1 → starplot-0.19.2}/hash_checks/map_checks.py

@@ -5,6 +5,7 @@ from datetime import datetime
 from zoneinfo import ZoneInfo
 
 import numpy as np
+from shapely import Polygon
 
 from starplot import (
     MapPlot,
@@ -20,6 +21,7 @@ from starplot import (
     Mercator,
     Mollweide,
     StereoNorth,
+    StereoSouth,
     CollisionHandler,
     Binoculars,
     Scope,
@@ -713,3 +715,72 @@ def check_map_allow_marker_and_line_collisions():
     p.export(filename)
     p.close_fig()
     return filename
+
+
+def check_map_constellation_clip_path():
+    constellation = Constellation.get(iau_id="and")
+
+    ra, dec = [p for p in constellation.border.coords.xy]
+    extent = (min(ra) - 2, max(min(dec) - 2, -90), max(ra) + 2, min(max(dec) + 2, 90))
+
+    if constellation.dec > 60:
+        proj = StereoNorth
+    elif constellation.dec < -60:
+        proj = StereoSouth
+    else:
+        proj = Miller
+
+    if extent[0] < 0:
+        extent = (extent[0] + 360, extent[1], extent[2] + 360, extent[3])
+
+    center_ra = (extent[0] + extent[2]) / 2
+    if center_ra < 0:
+        center_ra += 360
+    elif center_ra > 360:
+        center_ra -= 360
+
+    p = MapPlot(
+        projection=proj(center_ra=center_ra),
+        ra_min=extent[0],
+        ra_max=extent[2],
+        dec_min=extent[1],
+        dec_max=extent[3],
+        style=styles.PlotStyle().extend(
+            styles.extensions.BLUE_NIGHT,
+            styles.extensions.MAP,
+        ),
+        clip_path=Polygon(constellation.border.coords),
+        resolution=2000,
+        scale=0.8,
+    )
+
+    p.line(
+        geometry=constellation.border,
+        style=p.style.constellation_borders,
+    )
+
+    for hip1, hip2 in constellation.star_hip_lines:
+        star1 = Star.get(hip=hip1)
+        star2 = Star.get(hip=hip2)
+        p.line(
+            coordinates=[
+                (star1.ra, star1.dec),
+                (star2.ra, star2.dec),
+            ],
+            style=p.style.constellation_lines,
+        )
+
+    p.stars(
+        where=[_.hip.isin(constellation.star_hip_ids)],
+        where_labels=[_.magnitude < 4],
+        bayer_labels=True,
+    )
+
+    p.title(constellation.name, style__line_spacing=80)
+
+    p.ax.set_axis_off()  # hide the axis background that's outside the clip path
+
+    filename = DATA_PATH / "map-constellation-clip-path.png"
+    p.export(filename, padding=0.5)
+
+    return filename

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/__init__.py

@@ -2,7 +2,7 @@
 
 """Star charts and maps of the sky"""
 
-__version__ = "0.19.1"
+__version__ = "0.19.2"
 
 from .plots import (
     MapPlot,

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/data/catalogs.py

@@ -328,8 +328,8 @@ OPEN_NGC = Catalog(
 """
 
 CONSTELLATIONS_IAU = Catalog(
-    path=settings.data_path / "constellations.0.3.1.parquet",
-    url="https://github.com/steveberardi/starplot-constellations/releases/download/v0.3.1/constellations.0.3.1.parquet",
+    path=settings.data_path / "constellations.0.3.3.parquet",
+    url="https://github.com/steveberardi/starplot-constellations/releases/download/v0.3.3/constellations.0.3.3.parquet",
 )
 """
 Constellations recognized by IAU, with lines by Sky & Telescope.

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/data/constellations.py

@@ -43,6 +43,7 @@ def table(
 
     return c.mutate(
         boundary=_.boundary.cast("geometry"),  # cast WKB to geometry type
+        border=_.border.cast("geometry"),  # cast WKB to geometry type
         name=getattr(c, name_column),
     )
 

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/geod.py

@@ -1,7 +1,7 @@
 import math
 
 import pyproj
-
+import numpy as np
 
 GEOD = pyproj.Geod("+a=6378137 +f=0.0", sphere=True)
 
@@ -70,11 +70,10 @@ def ellipse(
 
     height = distance_m(height_degrees / 2)  # b
     width = distance_m(width_degrees / 2)  # a
+    step_size = (end_angle - start_angle) / num_pts
 
     points = []
-    for angle_pt in range(
-        start_angle, end_angle + 1, int((end_angle - start_angle) / num_pts)
-    ):
+    for angle_pt in np.arange(start_angle, end_angle + step_size, step_size):
         radians = math.radians(angle_pt)
         radius_a = (height * width) / math.sqrt(
             height**2 * (math.sin(radians)) ** 2

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/models/constellation.py

@@ -16,7 +16,11 @@ class Constellation(CatalogObject, SkyObject):
     """
 
     boundary: Polygon | MultiPolygon
-    """Shapely Polygon of the constellation's boundary. Right ascension coordinates are in degrees (0...360)."""
+    """
+    Shapely Polygon of the constellation's boundary. Right ascension coordinates are in degrees (0...360).
+    
+    If the constellation's boundary crosses the 0-meridian, then this will be a MultiPolygon split at the meridian.
+    """
 
     star_hip_ids: list[int]
     """List of HIP ids for stars that are part of the _lines_ for this constellation."""
@@ -35,6 +39,13 @@ class Constellation(CatalogObject, SkyObject):
     Serpens Caput has the `iau_id` of `ser1` and Serpens Cauda is `ser2`
     """
 
+    border: LineString = None
+    """
+    Shapely LineString of the constellation's border. Right ascension coordinates are in degrees (0...360).
+    
+    Coordinates in this geometry may extend past RA=360, if the border crosses the 0-meridian. This makes it ready to plot with the `line` function.
+    """
+
     def __repr__(self) -> str:
         return f"Constellation(iau_id={self.iau_id}, name={self.name}, ra={self.ra}, dec={self.dec})"
 

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/models/moon.py

@@ -37,7 +37,7 @@ class Moon(SkyObject):
     """Date/time of moon's position"""
 
     apparent_size: float
-    """Apparent size in the sky (degrees)"""
+    """Apparent diameter in the sky (degrees)"""
 
     phase_angle: float
     """Angle of the moon from the Sun (degrees)"""

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/models/planet.py

@@ -72,7 +72,7 @@ class Planet(SkyObject):
     """Date/time of planet's position"""
 
     apparent_size: float = 0
-    """Apparent size (degrees)"""
+    """Apparent diameter (degrees)"""
 
     geometry: Polygon = None
     """Shapely Polygon of the planet's extent. Right ascension coordinates are in degrees (0...360)."""

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/models/sun.py

@@ -22,7 +22,7 @@ class Sun(SkyObject):
     """Date/time of Sun's position"""
 
     apparent_size: float = 0
-    """Apparent size (degrees)"""
+    """Apparent diameter (degrees)"""
 
     geometry: Polygon = None
     """Shapely Polygon of the Sun's extent. Right ascension coordinates are in degrees (0...360)."""

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/plots/base.py

@@ -8,7 +8,7 @@ from matplotlib import pyplot as plt, patheffects
 from matplotlib.axes import Axes
 from matplotlib.figure import Figure
 from matplotlib.lines import Line2D
-from shapely import Polygon
+from shapely import Polygon, LineString
 
 from starplot.coordinates import CoordinateSystem
 from starplot import geod, models, warnings
@@ -334,9 +334,9 @@ class BasePlot(TextPlotterMixin, ABC):
                 polygon_style = style.marker.to_polygon_style()
                 polygon_style.edge_color = None
                 self.circle(
-                    (p.ra, p.dec),
-                    p.apparent_size,
-                    polygon_style,
+                    center=(p.ra, p.dec),
+                    radius_degrees=p.apparent_size / 2,
+                    style=polygon_style,
                     gid="planet-marker",
                 )
                 self._add_legend_handle_marker(legend_label, style.marker)
@@ -406,10 +406,11 @@ class BasePlot(TextPlotterMixin, ABC):
             polygon_style.edge_color = None
 
             self.circle(
-                (s.ra, s.dec),
-                s.apparent_size,
+                center=(s.ra, s.dec),
+                radius_degrees=s.apparent_size / 2,
                 style=polygon_style,
                 gid="sun-marker",
+                num_pts=200,
             )
 
             style.marker.symbol = MarkerSymbolEnum.CIRCLE
@@ -633,14 +634,27 @@ class BasePlot(TextPlotterMixin, ABC):
             )
 
     @use_style(LineStyle)
-    def line(self, coordinates: list[tuple[float, float]], style: LineStyle, **kwargs):
+    def line(
+        self,
+        style: LineStyle,
+        coordinates: list[tuple[float, float]] = None,
+        geometry: LineString = None,
+        **kwargs,
+    ):
         """Plots a line
 
         Args:
             coordinates: List of coordinates, e.g. `[(ra, dec), (ra, dec)]`
+            geometry: A shapely LineString. If this value is passed, then the `coordinates` kwarg will be ignored.
             style: Style of the line
         """
-        x, y = zip(*[self._prepare_coords(*p) for p in coordinates])
+
+        if coordinates is None and geometry is None:
+            raise ValueError("Must pass coordinates or geometry when plotting lines.")
+
+        coords = geometry.coords if geometry is not None else coordinates
+
+        x, y = zip(*[self._prepare_coords(*p) for p in coords])
 
         self.ax.plot(
             x,
@@ -700,9 +714,9 @@ class BasePlot(TextPlotterMixin, ABC):
 
             if show_phase:
                 self._moon_with_phase(
-                    m.phase_description,
-                    (m.ra, m.dec),
-                    m.apparent_size,
+                    moon_phase=m.phase_description,
+                    center=(m.ra, m.dec),
+                    radius_degrees=m.apparent_size / 2,
                     style=polygon_style,
                     dark_side_color=style.marker.edge_color,
                 )
@@ -797,8 +811,8 @@ class BasePlot(TextPlotterMixin, ABC):
         # Plot left side
         self.ellipse(
             center,
-            radius_degrees * 2,
-            radius_degrees * 2,
+            height_degrees=radius_degrees * 2,
+            width_degrees=radius_degrees * 2,
             style=left,
             num_pts=num_pts,
             angle=0,
@@ -808,8 +822,8 @@ class BasePlot(TextPlotterMixin, ABC):
         # Plot right side
         self.ellipse(
             center,
-            radius_degrees * 2,
-            radius_degrees * 2,
+            height_degrees=radius_degrees * 2,
+            width_degrees=radius_degrees * 2,
             style=right,
             num_pts=num_pts,
             angle=180,
@@ -819,8 +833,8 @@ class BasePlot(TextPlotterMixin, ABC):
         # Plot middle
         self.ellipse(
             center,
-            radius_degrees * 2,
-            radius_degrees,
+            height_degrees=radius_degrees * 2,
+            width_degrees=radius_degrees,
             style=middle,
             gid="moon-marker",
         )

{starplot-0.19.1 → starplot-0.19.2}/src/starplot/plots/map.py

@@ -11,7 +11,7 @@ from skyfield.api import wgs84
 import numpy as np
 
 from starplot.coordinates import CoordinateSystem
-from starplot import geod, geometry
+from starplot import geod
 from starplot.plots.base import BasePlot, DPI
 from starplot.mixins import ExtentMaskMixin
 from starplot.models.observer import Observer
@@ -124,9 +124,6 @@ class MapPlot(
         self.dec_min = dec_min
         self.dec_max = dec_max
 
-        if clip_path and clip_path.geom_type == "MultiPolygon":
-            clip_path = geometry.union_at_zero(clip_path.geoms[0], clip_path.geoms[1])
-
         self.clip_path = clip_path
 
         self._geodetic = ccrs.Geodetic()

{starplot-0.19.1 → starplot-0.19.2}/tutorial/tutorial_06.py

@@ -39,7 +39,7 @@ p.stars(
 )
 
 p.nebula(
-    # select DSOs which have no defined magnitude or less than 7
+    # select nebulae which have no defined magnitude or less than 7
     where=[_.magnitude.isnull() | (_.magnitude < 7)],
     where_true_size=[_.size > 1],  # only plot large nebula as their true size
 )

{starplot-0.19.1 → starplot-0.19.2}/tutorial/tutorial_07.py

@@ -33,8 +33,8 @@ p.planets(
     # since we're plotting the planets as their "true size"
     # the 'auto' offset won't work (it's not supported yet!)
     # so we manually set the offsets here:
-    style__label__offset_x=86,
-    style__label__offset_y=-40,
+    style__label__offset_x=68,
+    style__label__offset_y=-30,
     style__label__font_size=56,
     style__marker__color="#fcdb72",
 )