starplot 0.21.0__tar.gz → 0.21.1__tar.gz

{starplot-0.21.0 → starplot-0.21.1}/CITATION.cff

@@ -20,5 +20,5 @@ keywords:
   - graphs
   - plotting
 license: MIT
-version: 0.21.0
+version: 0.21.1
 date-released: '2026-09-27'

{starplot-0.21.0 → starplot-0.21.1}/Makefile

@@ -88,7 +88,7 @@ build-star-designations:
 	uv run $(DOTENV) data/scripts/star_designations.py
 
 version:
-	uv run $(DOTENV) python -c 'import starplot as sp; print(sp.__version__)'
+	@uv run $(DOTENV) python -c 'import starplot as sp; print(sp.__version__)'
 
 setup:
 	uv run $(DOTENV) starplot setup

{starplot-0.21.0 → starplot-0.21.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: starplot
-Version: 0.21.0
+Version: 0.21.1
 Summary: Star charts and maps of the sky
 Keywords: astronomy,stars,charts,maps,constellations,sky,plotting
 Author-email: Steve Berardi <steve@starplot.dev>
@@ -111,11 +111,6 @@ p.export("starchart.svg")
 [https://starplot.dev](https://starplot.dev)
 
 
-## Demo
-For a demo of Starplot's zenith charts, check out: 
-
-[Sky Atlas - Star Chart Creator](https://skyatlas.app/star-charts/)
-
 ## Getting Help / Updates
 
 - Chat with other starplotters on our [Discord server](https://discord.gg/WewJJjshFu)

{starplot-0.21.0 → starplot-0.21.1}/README.md

@@ -71,11 +71,6 @@ p.export("starchart.svg")
 [https://starplot.dev](https://starplot.dev)
 
 
-## Demo
-For a demo of Starplot's zenith charts, check out: 
-
-[Sky Atlas - Star Chart Creator](https://skyatlas.app/star-charts/)
-
 ## Getting Help / Updates
 
 - Chat with other starplotters on our [Discord server](https://discord.gg/WewJJjshFu)

{starplot-0.21.0 → starplot-0.21.1}/docs/changelog.md

@@ -8,6 +8,7 @@
 - Refactors styling framework with more consistent naming
 - Removed fonts from PyPI package, to reduce our filesize
 - Adds support for Python 3.14
+- [**v0.21.1**] Fixes a couple bugs with a few azimuthal projections
 
 ---
 

{starplot-0.21.0 → starplot-0.21.1}/examples/map_milky_way_stars.py

@@ -17,7 +17,7 @@ p = MapPlot(
     scale=0.36,
 )
 p.stars(
-    where=[_.magnitude < 11],
+    where=[_.magnitude < 12],
     where_labels=False,
     opacity_fn=lambda s: 0.9 if s.magnitude < 9 else 0.6,
     catalog=BIG_SKY,

{starplot-0.21.0 → starplot-0.21.1}/src/starplot/__init__.py

@@ -2,7 +2,7 @@
 
 """Star charts and maps of the sky"""
 
-__version__ = "0.21.0"
+__version__ = "0.21.1"
 
 import contextlib
 

{starplot-0.21.0 → starplot-0.21.1}/src/starplot/svg/canvas.py

@@ -16,7 +16,10 @@ from starplot import geometry as _geometry
 from starplot.config import settings
 from starplot.projections import (
     CoordinateReferenceSystem,
+    Equidistant,
+    LambertAzEqArea,
     ProjectionBase,
+    Stereographic,
     latlon_bounds_to_projection,
 )
 from starplot.styles import (
@@ -150,6 +153,53 @@ class Canvas:
             and abs(self.bounds[1] - self.bounds[3]) >= 179
         )
 
+    def _extent_contains_projection_singularity(self) -> bool:
+        """
+        Returns True if the extent contains the antipode of an azimuthal
+        projection's center, where the projection is either undefined
+        (stereographic) or stretched into a ring around the whole disc
+        (equidistant, Lambert equal-area).
+        """
+        if not isinstance(
+            self.projection, (Stereographic, Equidistant, LambertAzEqArea)
+        ):
+            return False
+
+        ra_min, dec_min, ra_max, dec_max = self.bounds
+        antipode_ra = (self.projection.center_ra + 180) % 360
+        antipode_dec = -self.projection.center_dec
+
+        if not dec_min <= antipode_dec <= dec_max:
+            return False
+
+        ra_span = ra_max - ra_min
+        return ra_span >= 360 or (antipode_ra - ra_min) % 360 <= ra_span
+
+    def _radec_bounds_from_projected(self) -> tuple[float, float, float, float]:
+        """
+        Converts the projected bounds back to RA/DEC bounds, making sure they include
+        any celestial pole that's inside the projected bounds.
+        """
+        ra_min, dec_min, ra_max, dec_max = self.tx.transform_bounds(
+            *self.projected_bounds, direction="INVERSE"
+        )
+
+        # transform_bounds only samples the edges of the projected box, so it can miss a
+        # pole that's inside the box (e.g. the south pole in an equatorial Stereographic
+        # plot), which leaves things like meridian gridlines stopping short of that pole
+        for pole_dec in (90, -90):
+            px, py = self.tx.transform(0, pole_dec)
+            if not (np.isfinite(px) and np.isfinite(py)):
+                continue
+            if self.minx <= px <= self.maxx and self.miny <= py <= self.maxy:
+                ra_min, ra_max = -180.0, 180.0
+                if pole_dec > 0:
+                    dec_max = 90.0
+                else:
+                    dec_min = -90.0
+
+        return ra_min, dec_min, ra_max, dec_max
+
     def _init_bounds(self):
         """
         Calculates true bounds from user-bounds, which can change slightly as a result of the map projection used.
@@ -163,9 +213,7 @@ class Canvas:
                 *self.clip_path.bounds, densify_pts=1_000
             )
             self.projected_bounds = self.minx, self.miny, self.maxx, self.maxy
-            self.bounds = self.tx.transform_bounds(
-                *self.projected_bounds, direction="INVERSE"
-            )
+            self.bounds = self._radec_bounds_from_projected()
         elif self._is_global():
             self.minx, self.miny, self.maxx, self.maxy = self.projection.global_bounds
             self.projected_bounds = self.minx, self.miny, self.maxx, self.maxy
@@ -184,10 +232,22 @@ class Canvas:
                 curved=self.projection.curved,
                 transformer=self.tx,
             )
+
+            if self._extent_contains_projection_singularity():
+                cx, cy = self.tx.transform(
+                    self.projection.center_ra, self.projection.center_dec
+                )
+                radius = max(
+                    abs(self.minx - cx),
+                    abs(self.maxx - cx),
+                    abs(self.miny - cy),
+                    abs(self.maxy - cy),
+                )
+                self.minx, self.maxx = cx - radius, cx + radius
+                self.miny, self.maxy = cy - radius, cy + radius
+
             self.projected_bounds = self.minx, self.miny, self.maxx, self.maxy
-            self.bounds = self.tx.transform_bounds(
-                *self.projected_bounds, direction="INVERSE"
-            )
+            self.bounds = self._radec_bounds_from_projected()
 
         self._refresh_figure_dimensions()
 
@@ -239,9 +299,7 @@ class Canvas:
             self.maxy = lerp(maxy, miny, ay0)
 
             self.projected_bounds = self.minx, self.miny, self.maxx, self.maxy
-            self.bounds = self.tx.transform_bounds(
-                *self.projected_bounds, direction="INVERSE"
-            )
+            self.bounds = self._radec_bounds_from_projected()
 
             self._refresh_figure_dimensions()
 

{starplot-0.21.0 → starplot-0.21.1}/tests/svg/test_canvas.py

@@ -14,6 +14,8 @@ from starplot.projections import (
     Equidistant,
     Orthographic,
     PlateCarree,
+    Stereographic,
+    StereoNorth,
 )
 from starplot.styles import (
     GradientStyle,
@@ -666,6 +668,45 @@ class TestInitBounds:
         assert canvas.minx != canvas.maxx
         assert math.isfinite(canvas.minx) and math.isfinite(canvas.maxx)
 
+    def test_bounds_include_pole_inside_clip_path(self):
+        # GIVEN an equatorial Stereographic canvas whose circular clip path
+        # (radius 110 degrees) contains both celestial poles
+        clip_path = _geometry.circle(center=(180, 0), diameter_degrees=220, num_pts=200)
+
+        # WHEN it initializes its bounds
+        canvas = _canvas(Stereographic(), bounds=[0, -90, 360, 90], clip_path=clip_path)
+
+        # THEN the RA/DEC bounds reach both poles
+        assert canvas.bounds[1] == -90
+        assert canvas.bounds[3] == 90
+
+    def test_extent_containing_stereographic_antipode_is_centered(self):
+        # GIVEN a Stereographic projection centered at dec 40, and an extent
+        # that contains its antipode (RA 280, dec -40)
+        projection = Stereographic(center_ra=100, center_dec=40)
+
+        # WHEN constructing the canvas
+        canvas = _canvas(projection, bounds=[0, -60, 360, 90])
+
+        # THEN the projection's center lands in the middle of the plot, and a
+        # point 100 degrees from center is still inside the plot
+        assert canvas._extent_contains_projection_singularity()
+        assert canvas._to_axes(100, 40) == pytest.approx((0.5, 0.5), abs=1e-6)
+        x, y = canvas._to_axes(100 + 180, 50)
+        assert 0 <= x <= 1 and 0 <= y <= 1
+
+    def test_extent_away_from_center_is_not_recentered(self):
+        # GIVEN a StereoNorth projection (centered on the north pole) with an
+        # extent that doesn't include the pole or its antipode
+        canvas = _canvas(StereoNorth(), bounds=[161.25, 47, 213, 65])
+
+        # WHEN it initializes its bounds
+        # THEN the bounds are not re-centered on the pole, so the pole stays
+        # outside the plot instead of being padded into view
+        assert not canvas._extent_contains_projection_singularity()
+        _, y = canvas._to_axes(0, 90)
+        assert not (0 <= y <= 1)
+
 
 class TestRefreshFigureDimensions:
     def test_wide_bounds_fill_width_and_derive_height_from_aspect_ratio(self):

starplot-0.21.0/examples/galaxy_custom_marker.py

@@ -1,177 +0,0 @@
-import numpy as np
-import matplotlib.path as mpath
-import matplotlib.patches as mpatches
-import cartopy.crs as ccrs
-from astropy.coordinates import SkyCoord
-
-from starplot import MapPlot, LambertAzEqArea, _
-from starplot.styles import PlotStyle, extensions
-
-group_blue = {
-    "M81": {"RA": 148.888, "DEC": 69.065, "distance": 3.70},
-    "M82": {"RA": 148.971, "DEC": 69.679, "distance": 3.61},
-    "NGC3077": {"RA": 150.829, "DEC": 68.734, "distance": 3.85},
-    "NGC2403": {"RA": 114.213, "DEC": 65.602, "distance": 3.19},
-    "NGC4236": {"RA": 184.175, "DEC": 69.462, "distance": 4.41},
-    "NGC4244": {"RA": 184.372, "DEC": 37.807, "distance": 4.31},
-    "NGC4449": {"RA": 187.046, "DEC": 44.093, "distance": 4.27},
-}
-gal_ra = np.array([g["RA"] for g in group_blue.values()])
-gal_dec = np.array([g["DEC"] for g in group_blue.values()])
-gal_names = list(group_blue.keys())
-
-group_silver = {
-    "NGC4214": {"RA": 183.912, "DEC": 36.327, "distance": 2.94},
-    "HolmII": {"RA": 124.767, "DEC": 70.714, "distance": 3.39},
-    "NGC4395": {"RA": 186.458, "DEC": 33.5461, "distance": 4.61},
-    "M94": {"RA": 192.723, "DEC": 41.1194, "distance": 4.66},
-    "NGC2366": {"RA": 112.228, "DEC": 69.2053, "distance": 3.19},
-    "DDO82": {"RA": 157.646, "DEC": 70.6194, "distance": 4.00},
-    "DDO168": {"RA": 198.619, "DEC": 45.9194, "distance": 4.33},
-    "DO165": {"RA": 196.612, "DEC": 67.7042, "distance": 4.57},
-    "IC3687": {"RA": 190.563, "DEC": 38.5019, "distance": 4.57},
-    "HolmI": {"RA": 145.135, "DEC": 71.1864, "distance": 3.84},
-}
-other_gal_ra = np.array([g["RA"] for g in group_silver.values()])
-other_gal_dec = np.array([g["DEC"] for g in group_silver.values()])
-other_gal_names = list(group_silver.keys())
-
-
-def get_hurricane():
-    u = np.array(
-        [
-            [2.444, 7.553],
-            [0.513, 7.046],
-            [-1.243, 5.433],
-            [-2.353, 2.975],
-            [-2.578, 0.092],
-            [-2.075, -1.795],
-            [-0.336, -2.870],
-            [2.609, -2.016],
-        ]
-    )
-    u[:, 0] -= 0.098
-    codes = [1] + [2] * (len(u) - 2) + [2]
-    u = np.append(u, -u[::-1], axis=0)
-    codes += codes
-    # Scale to angular size
-    scale = 1.15 / np.max(np.abs(u))
-    return mpath.Path(scale * u, codes, closed=False)
-
-
-hurricane = get_hurricane()
-
-center_coord = SkyCoord(ra=160, dec=58, unit="deg")
-center_ra = center_coord.ra.deg
-center_dec = center_coord.dec.deg
-
-style = PlotStyle().extend(
-    extensions.BLUE_LIGHT,
-    extensions.MAP,
-)
-p = MapPlot(
-    projection=LambertAzEqArea(center_ra=center_ra, center_dec=center_dec),
-    ra_min=120,
-    ra_max=200,
-    dec_min=33,
-    dec_max=78,
-    style=style,
-    resolution=4000,
-    autoscale=True,
-)
-
-p.stars(
-    where=[_.magnitude < 8],
-    bayer_labels=True,
-    where_labels=[_.magnitude < 1.75],
-    opacity_fn=lambda star: max(0.01, 1.0 - 0.1 * star.magnitude),
-)
-
-for ra, dec, name in zip(gal_ra, gal_dec, gal_names):
-    # Build hurricane polygon in sky coordinates (degrees)
-    ra_offsets = hurricane.vertices[:, 0] / np.cos(
-        np.radians(dec)
-    )  # Adjust angular size for declination
-    dec_offsets = hurricane.vertices[:, 1]
-    hurricane_lon = ra + ra_offsets
-    hurricane_lat = dec + dec_offsets
-    # Convert RA to the longitude convention starplot/cartopy use (RA increases eastward)
-    hurricane_lon_cartopy = np.mod(360.0 - hurricane_lon, 360.0)
-    # Transform each vertex into the map projection
-    projected_vertices = np.array(
-        [
-            p.ax.projection.transform_point(lon, lat, ccrs.PlateCarree())
-            for lon, lat in zip(hurricane_lon_cartopy, hurricane_lat)
-        ]
-    )
-
-    patch = mpatches.PathPatch(
-        mpath.Path(projected_vertices, hurricane.codes),
-        facecolor="lightblue",
-        edgecolor="navy",
-        linewidth=5.0,
-        zorder=30,
-    )
-    p.ax.add_patch(patch)
-    p.marker(
-        ra=ra,
-        dec=dec,
-        style={
-            "marker": {
-                "symbol": "circle",
-                "size": 10,
-                "fill": "darkblue",
-                "stroke": "navy",
-                "stroke_width": 3.0,
-            }
-        },
-    )
-    p.text(
-        text=name,
-        ra=ra,
-        dec=dec,
-        style={
-            "fill": "navy",
-            "font_size": 28,
-            "font_weight": 700,
-            "offset_x": 35,
-            "offset_y": 35,
-        },
-    )
-
-for ra, dec, name in zip(other_gal_ra, other_gal_dec, other_gal_names):
-    p.marker(
-        ra=ra,
-        dec=dec,
-        style={
-            "marker": {
-                "symbol": "ellipse",
-                "size": 30,
-                "fill": "orange",
-                "stroke": "black",
-                "stroke_width": 2.0,
-            }
-        },
-    )
-    p.text(
-        text=name,
-        ra=ra,
-        dec=dec,
-        style={
-            "fill": "black",
-            "font_size": 28,
-            "offset_x": 15.0,
-            "offset_y": 40.0,
-        },
-    )
-
-p.gridlines()
-p.constellations(style__stroke="#78d78e", style__opacity=0.25)
-p.constellation_borders(style__stroke="#035019", style__opacity=0.25)
-p.constellation_labels(style__font_size=28, style__fill="#035019")
-p.milky_way()
-p.ecliptic()
-p.celestial_equator()
-p.ax.tick_params(axis="both", labelsize=28)
-
-p.export("galaxy_custom_marker.png", padding=0.3, transparent=True)

starplot-0.21.0/examples/horizon_sgr.py

@@ -1,39 +0,0 @@
-from datetime import datetime
-from zoneinfo import ZoneInfo
-
-from starplot import HorizonPlot, Observer, _
-from starplot.styles import PlotStyle, extensions
-
-style = PlotStyle().extend(
-    extensions.BLUE_MEDIUM,
-    extensions.HORIZON,
-)
-
-dt = datetime(2024, 8, 30, 21, 0, 0, 0, tzinfo=ZoneInfo("US/Pacific"))
-
-observer = Observer(
-    dt=dt,
-    lat=36.606111,  # Lone Pine, California
-    lon=-118.079444,
-)
-
-p = HorizonPlot(
-    altitude=(0, 60),
-    azimuth=(135, 225),
-    observer=observer,
-    style=style,
-    resolution=4000,
-    scale=0.9,
-)
-
-p.constellations()
-p.milky_way()
-
-p.stars(where=[_.magnitude < 5], where_labels=[_.magnitude < 2])
-p.messier(where=[_.magnitude < 11], where_true_size=[False])
-
-p.constellation_labels()
-p.horizon(labels={180: "SOUTH"})
-p.gridlines()
-
-p.export("horizon_sgr.png")