starplot 0.11.4__py2.py3-none-any.whl → 0.12.0__py2.py3-none-any.whl

starplot/__init__.py

@@ -1,6 +1,6 @@
 """Star charts and maps"""
 
-__version__ = "0.11.4"
+__version__ = "0.12.0"
 
 from .base import BasePlot  # noqa: F401
 from .map import MapPlot, Projection  # noqa: F401

starplot/base.py

@@ -22,6 +22,7 @@ from starplot.styles import (
     LabelStyle,
     LegendLocationEnum,
     LegendStyle,
+    LineStyle,
     MarkerSymbolEnum,
     PathStyle,
     PolygonStyle,
@@ -84,7 +85,7 @@ class BasePlot(ABC):
 
         self.text_border = patheffects.withStroke(
             linewidth=self.style.text_border_width,
-            foreground=self.style.background_color.as_hex(),
+            foreground=self.style.text_border_color.as_hex(),
         )
         self._size_multiplier = self.resolution / 3000
         self.timescale = load.timescale().from_datetime(self.dt)
@@ -594,6 +595,25 @@ class BasePlot(ABC):
             num_pts=num_pts,
         )
 
+    @use_style(LineStyle)
+    def line(self, coordinates: list[tuple[float, float]], style: LineStyle):
+        """Plots a line
+
+        Args:
+            coordinates: List of coordinates, e.g. `[(ra, dec), (ra, dec)]`
+            style: Style of the line
+        """
+        x, y = zip(*[self._prepare_coords(*p) for p in coordinates])
+
+        self.ax.plot(
+            x,
+            y,
+            clip_on=True,
+            clip_path=self._background_clip_path,
+            **style.matplot_kwargs(self._size_multiplier),
+            **self._plot_kwargs(),
+        )
+
     @use_style(ObjectStyle, "moon")
     def moon(
         self,

starplot/data/constellations.py

@@ -112,12 +112,27 @@ CONSTELLATIONS_ABBREVIATIONS = {k.lower(): k.upper() for k in properties.keys()}
 """Constellation 3-letter abbreviations"""
 
 
-def load():
+def _load_deprecated():
     with _load.open("constellations_hip.fab") as f:
         consdata = stellarium.parse_constellations(f)
     return consdata
 
 
+def load(**kwargs):
+    import geopandas as gpd
+    import numpy as np
+    from starplot.data import DataFiles
+
+    cons = gpd.read_file(
+        DataFiles.CONSTELLATIONS.value,
+        engine="pyogrio",
+        use_arrow=True,
+        **kwargs,
+    )
+    cons = cons.replace({np.nan: None})
+    return cons
+
+
 def get(constellation_id: str):
     return properties.get(constellation_id)
 

starplot/data/dsos.py

@@ -81,7 +81,7 @@ ONGC_TYPE = {
 
 ONGC_TYPE_MAP = {v: k.value for k, v in ONGC_TYPE.items()}
 
-DEFAULT_DSO_TYPES = [
+BASIC_DSO_TYPES = [
     # Star Clusters ----------
     DsoType.OPEN_CLUSTER,
     DsoType.GLOBULAR_CLUSTER,

starplot/data/prep/constellations.py

@@ -0,0 +1,108 @@
+import geopandas as gpd
+import pandas as pd
+
+from shapely.geometry import Polygon, MultiPolygon
+
+from starplot.data import constellations
+from starplot.data.prep.utils import RAW_DATA_PATH, DATA_LIBRARY
+
+
+DATA_PATH = RAW_DATA_PATH / "iau"
+CRS = "+ellps=sphere +f=0 +proj=latlong +axis=wnu +a=6378137 +no_defs"
+
+# - Three letter (index)
+# - Full name
+# - Centroid
+# - Min bounding box of lines
+# - Polygon of bounds
+# - List of HIP ids for lines
+
+# constellation_dict = {
+#     "id": [],
+#     "name": [],
+#     "centroid": [],
+#     "lines_bounding_box": [],
+#     "lines_hip_ids": [],
+#     "borders": [],
+# }
+
+
+def parse_ra(ra_str):
+    """Parses RA from border file HH MM SS to 0...360 degree float"""
+    h, m, s = ra_str.strip().split(" ")
+    return round(15 * (float(h) + float(m) / 60 + float(s) / 3600), 4)
+
+
+def parse_dec(dec_str):
+    """Parses DEC from ONGC CSV from HH:MM:SS to -90...90 degree float"""
+    return round(float(dec_str), 4)
+
+
+def parse_borders(lines):
+    coords = []
+    for line in lines:
+        if "|" not in line:
+            continue
+        ra_str, dec_str, _ = line.split("|")
+        ra = parse_ra(ra_str)
+        dec = parse_dec(dec_str)
+        coords.append((ra, dec))
+    return coords
+
+
+def build_constellations():
+    constellation_records = []
+    con_lines = constellations.lines()
+
+    for cid, props in constellations.properties.items():
+        constellation_dict = {
+            "id": cid.lower(),
+            "iau_id": cid.lower(),
+            "name": props[0],
+            "center_ra": props[1] * 15,
+            "center_dec": props[2],
+            "lines_hip_ids": ",".join(
+                "-".join([str(h) for h in hips]) for hips in con_lines[cid.lower()]
+            ),
+        }
+        # print(cid)
+
+        if cid == "Ser":
+            ser1_coords = []
+            ser2_coords = []
+            with open(DATA_PATH / "ser1.txt", "r") as ser1:
+                ser1_coords = parse_borders(ser1.readlines())
+
+            with open(DATA_PATH / "ser2.txt", "r") as ser2:
+                ser2_coords = parse_borders(ser2.readlines())
+
+            constellation_dict["geometry"] = MultiPolygon(
+                [Polygon(ser1_coords), Polygon(ser2_coords)]
+            )
+
+        else:
+            with open(DATA_PATH / f"{cid.lower()}.txt", "r") as borderfile:
+                coords = parse_borders(borderfile.readlines())
+                constellation_dict["geometry"] = Polygon(coords)
+
+        constellation_records.append(constellation_dict)
+
+    return constellation_records
+
+
+constellation_records = build_constellations()
+df = pd.DataFrame.from_records(constellation_records)
+
+gdf = gpd.GeoDataFrame(
+    df,
+    geometry=df["geometry"],
+    crs=CRS,
+)
+gdf = gdf.set_index("id")
+gdf.to_file(
+    DATA_LIBRARY / "constellations.gpkg", driver="GPKG", engine="pyogrio", index=True
+)
+
+print(gdf.loc["uma"])
+
+print("Total Constellations: " + str(len(constellation_records)))

starplot/data/prep/dsos.py

@@ -0,0 +1,299 @@
+import os
+from pathlib import Path
+
+import geopandas as gpd
+import pandas as pd
+import numpy as np
+
+from shapely.geometry import Polygon, MultiPolygon, Point
+
+from starplot import geod
+from starplot.data.prep.utils import RAW_DATA_PATH, DATA_LIBRARY, zip_file
+
+DATA_PATH = RAW_DATA_PATH / "ongc" / "outlines"
+CRS = "+ellps=sphere +f=0 +proj=latlong +axis=wnu +a=6378137 +no_defs"
+IGNORE_OUTLINES = [
+    "IC0424",  # seems too big?
+]
+MIN_SIZE = 0.1
+
+
+def _size(d):
+    """Returns size (in sq degrees) of minimum bounding rectangle of a DSO"""
+    size = None
+    geometry_types = d["geometry"].geom_type
+
+    if "Polygon" in geometry_types and "MultiPolygon" not in geometry_types:
+        size = d.geometry.envelope.area
+
+    elif "MultiPolygon" in geometry_types:
+        size = sum([p.envelope.area for p in d.geometry.geoms])
+
+    elif d.maj_ax and not np.isnan(d.min_ax):
+        size = (d.maj_ax / 60) * (d.min_ax / 60)
+
+    elif d.maj_ax:
+        size = (d.maj_ax / 60) ** 2
+
+    if size:
+        size = round(size, 8)
+
+    return size
+
+
+def read_csv():
+    df = pd.read_csv(
+        "raw/ongc/NGC.csv",
+        sep=";",
+    )
+    df_addendum = pd.read_csv(
+        "raw/ongc/addendum.csv",
+        sep=";",
+    )
+    df = pd.concat([df, df_addendum])
+
+    df["ra_degrees"] = df.apply(parse_ra, axis=1)
+    df["dec_degrees"] = df.apply(parse_dec, axis=1)
+    df["M"] = df.apply(parse_m, axis=1)
+    df["IC"] = df.apply(parse_ic, axis=1)
+    df["NGC"] = df.apply(parse_ngc, axis=1)
+
+    df = df.rename(
+        columns={
+            "Name": "name",
+            "Type": "type",
+            "M": "m",
+            "NGC": "ngc",
+            "IC": "ic",
+            "MajAx": "maj_ax",
+            "MinAx": "min_ax",
+            "PosAng": "angle",
+            "B-Mag": "mag_b",
+            "V-Mag": "mag_v",
+            "NED notes": "ned_notes",
+            "Common names": "common_names",
+            "Const": "constellation",
+        }
+    )
+
+    df = df.drop(
+        columns=[
+            "Cstar U-Mag",
+            "Cstar B-Mag",
+            "Cstar V-Mag",
+            "Cstar Names",
+            "ned_notes",
+            "constellation",
+            "OpenNGC notes",
+            "Sources",
+            "Hubble",
+            "Identifiers",
+            "J-Mag",
+            "H-Mag",
+            "K-Mag",
+            "RA",
+            "Dec",
+        ],
+        axis=1,
+    )
+
+    gdf = gpd.GeoDataFrame(
+        df,
+        geometry=gpd.points_from_xy(df.ra_degrees, df.dec_degrees),
+        crs=CRS,
+    )
+
+    gdf["geometry"] = gdf.apply(create_ellipse, axis=1)
+    # gdf["geometry"] = gdf.apply(create_point, axis=1)
+
+    return gdf
+
+
+def create_point(d):
+    return Point([d.ra_degrees, d.dec_degrees])
+
+
+def parse_m(row):
+    """Parses messier number"""
+    if not np.isnan(row.M):
+        return str(int(row.M)).lstrip("0")
+
+    return None
+
+
+def parse_ic(row):
+    """Parses IC number if name starts with IC"""
+    if row.Name.startswith("IC"):
+        return row.Name[2:].lstrip("0")
+
+    if str(row.IC) != "nan":
+        return str(row.IC).lstrip("0")
+
+    return None
+
+
+def parse_ngc(row):
+    """Parses NGC number if name starts with NGC"""
+    if row.Name.startswith("NGC"):
+        return row.Name[3:].lstrip("0")
+
+    if str(row.NGC) != "nan":
+        return str(row.NGC).lstrip("0")
+
+    return None
+
+
+def parse_ra(row):
+    """Parses RA from ONGC CSV from HH:MM:SS to 0...360 degree float"""
+    if row.Type == "NonEx":
+        print(f"Non Existent object, ignoring... {row.Name}")
+        return
+
+    ra = row.RA
+    h, m, s = ra.split(":")
+    return round(15 * (float(h) + float(m) / 60 + float(s) / 3600), 4)
+
+
+def parse_dec(row):
+    """Parses DEC from ONGC CSV from HH:MM:SS to -90...90 degree float"""
+    if row.Type == "NonEx":
+        print(f"Non Existent object, ignoring... {row.Name}")
+        return
+
+    dec = row.Dec
+    if dec[0] == "-":
+        mx = -1
+    else:
+        mx = 1
+    h, m, s = dec[1:].split(":")
+    return round(mx * (float(h) + float(m) / 60 + float(s) / 3600), 4)
+
+
+def create_ellipse(d):
+    maj_ax, min_ax, angle = d.maj_ax, d.min_ax, d.angle
+
+    if np.isnan(maj_ax):
+        return d.geometry
+
+    if np.isnan(angle):
+        angle = 0
+
+    maj_ax_degrees = (maj_ax / 60) / 2
+
+    if np.isnan(min_ax):
+        min_ax_degrees = maj_ax_degrees
+    else:
+        min_ax_degrees = (min_ax / 60) / 2
+
+    points = geod.ellipse(
+        (d.ra_degrees / 15, d.dec_degrees),
+        min_ax_degrees * 2,
+        maj_ax_degrees * 2,
+        angle,
+        num_pts=100,
+    )
+
+    points360 = []
+    for lon, dec in points:
+        if lon < 0:
+            lon += 360
+        points360.append([round(lon, 4), round(dec, 4)])
+
+    return Polygon(points360)
+
+
+def parse_designation_from_filename(filename):
+    designation, level = filename.split("_")
+
+    if designation.startswith("IC"):
+        ic = designation[2:]
+    else:
+        ic = None
+
+    if designation.startswith("NGC"):
+        ngc = designation[3:]
+    else:
+        ngc = None
+
+    return designation, ic, ngc, level
+
+
+def walk_files(path=DATA_PATH):
+    for dirpath, dirnames, filenames in os.walk(path):
+        for filename in sorted(filenames):
+            yield Path(os.path.join(dirpath, filename))
+
+
+gdf = read_csv()
+gdf = gdf.set_index("name")
+
+outlines = {}
+
+for f in walk_files():
+    designation, ic, ngc, level = parse_designation_from_filename(f.name)
+    name, _ = f.name.split("_")
+
+    if level == "lv3" or name in IGNORE_OUTLINES:
+        continue
+    # if designation in d['designation']:
+    #     continue
+
+    dso_df = pd.read_csv(f, sep="\t")
+    polygons = []
+    current_poly = []
+
+    for i, row in dso_df.iterrows():
+        cont_flag = row["Cont_Flag"]
+        ra = row["RAJ2000"]
+        dec = row["DEJ2000"]
+        current_poly.append([ra, dec])
+
+        if cont_flag == "*":
+            # a * indicates this is the last point in the current polygon
+            polygons.append(current_poly)
+            current_poly = []
+
+    if len(polygons) > 1:
+        dso_geom = MultiPolygon([Polygon(p) for p in polygons])
+    else:
+        dso_geom = Polygon(polygons[0])
+
+    if dso_geom.area > MIN_SIZE:
+        outlines[designation] = dso_geom
+
+    if not ic and not ngc:
+        print(designation)
+        centroid = dso_geom.centroid
+        gdf.loc[name, "geometry"] = dso_geom
+        gdf.loc[name, "ra_degrees"] = round(centroid.x, 4)
+        gdf.loc[name, "dec_degrees"] = round(centroid.y, 4)
+        gdf.loc[name, "type"] = "Neb"
+    else:
+        if gdf.loc[name].empty:
+            print(f"NGC/IC object not found: {name}")
+        elif dso_geom.area > MIN_SIZE:
+            gdf.loc[name, "geometry"] = dso_geom
+
+
+# add size column
+gdf["size_deg2"] = gdf.apply(_size, axis=1)
+
+# print(gdf.loc["NGC1976"])
+
+print(gdf.loc["NGC6720"])  # ring nebula
+print(gdf.loc["Mel022"])  # M45
+# print("INDEX")
+# print(gdf.sindex.size)
+
+# gdf.to_file(BUILD_PATH / "ongc.gpkg", driver="GPKG", crs=CRS, index=True)
+
+gdf.set_crs(CRS, inplace=True)
+gdf.to_file(DATA_LIBRARY / "ongc.gpkg", driver="GPKG", engine="pyogrio", index=True)
+# crs=CRS, engine="pyogrio",
+
+print("Total nebula outlines: " + str(len(outlines)))
+# result = gpd.read_file(HERE.parent / "build" / "ngc.gpkg")
+# print(result)
+
+# Zip it up!
+zip_file(DATA_LIBRARY / "ongc.gpkg")

starplot/data/prep/utils.py

@@ -0,0 +1,16 @@
+import zipfile
+from pathlib import Path
+
+
+HERE = Path(__file__).resolve().parent
+ROOT = HERE.parent.parent.parent.parent
+RAW_DATA_PATH = ROOT / "raw"
+BUILD_PATH = ROOT / "build"
+DATA_LIBRARY = HERE.parent / "library"
+
+CRS = "+ellps=sphere +f=0 +proj=latlong +axis=wnu +a=6378137 +no_defs"
+
+
+def zip_file(filename):
+    zipped = zipfile.ZipFile(f"{filename}.zip", "w", zipfile.ZIP_DEFLATED)
+    zipped.write(filename)