astro_chart 0.1.0

data/ext/astro_chart/swephlib.h

@@ -0,0 +1,189 @@
+
+/************************************************************
+
+  Authors: Dieter Koch and Alois Treindl, Astrodienst Zurich
+
+************************************************************/
+/* Copyright (C) 1997 - 2021 Astrodienst AG, Switzerland.  All rights reserved.
+
+  License conditions
+  ------------------
+
+  This file is part of Swiss Ephemeris.
+
+  Swiss Ephemeris is distributed with NO WARRANTY OF ANY KIND.  No author
+  or distributor accepts any responsibility for the consequences of using it,
+  or for whether it serves any particular purpose or works at all, unless he
+  or she says so in writing.  
+
+  Swiss Ephemeris is made available by its authors under a dual licensing
+  system. The software developer, who uses any part of Swiss Ephemeris
+  in his or her software, must choose between one of the two license models,
+  which are
+  a) GNU Affero General Public License (AGPL)
+  b) Swiss Ephemeris Professional License
+
+  The choice must be made before the software developer distributes software
+  containing parts of Swiss Ephemeris to others, and before any public
+  service using the developed software is activated.
+
+  If the developer choses the AGPL software license, he or she must fulfill
+  the conditions of that license, which includes the obligation to place his
+  or her whole software project under the AGPL or a compatible license.
+  See https://www.gnu.org/licenses/agpl-3.0.html
+
+  If the developer choses the Swiss Ephemeris Professional license,
+  he must follow the instructions as found in http://www.astro.com/swisseph/ 
+  and purchase the Swiss Ephemeris Professional Edition from Astrodienst
+  and sign the corresponding license contract.
+
+  The License grants you the right to use, copy, modify and redistribute
+  Swiss Ephemeris, but only under certain conditions described in the License.
+  Among other things, the License requires that the copyright notices and
+  this notice be preserved on all copies.
+
+  Authors of the Swiss Ephemeris: Dieter Koch and Alois Treindl
+
+  The authors of Swiss Ephemeris have no control or influence over any of
+  the derived works, i.e. over software or services created by other
+  programmers which use Swiss Ephemeris functions.
+
+  The names of the authors or of the copyright holder (Astrodienst) must not
+  be used for promoting any software, product or service which uses or contains
+  the Swiss Ephemeris. This copyright notice is the ONLY place where the
+  names of the authors can legally appear, except in cases where they have
+  given special permission in writing.
+
+  The trademarks 'Swiss Ephemeris' and 'Swiss Ephemeris inside' may be used
+  for promoting such software, products or services.
+*/
+
+#define PREC_IAU_1976_CTIES          2.0 	/* J2000 +/- two centuries */
+#define PREC_IAU_2000_CTIES          2.0 	/* J2000 +/- two centuries */
+/* we use P03 for whole ephemeris */
+#define PREC_IAU_2006_CTIES          75.0 	/* J2000 +/- 75 centuries */
+
+/* For reproducing JPL Horizons to 2 mas (SEFLG_JPLHOR): 
+ * The user has to keep the following files up to date which contain
+ * the earth orientation parameters related to the IAU 1980 nutation
+ * theory. 
+ * Download the file 
+ * datacenter.iers.org/eop/-/somos/5Rgv/document/tx13iers.u24/eopc04_08.62-now
+ * and rename it as eop_1962_today.txt. For current data and estimations for
+ * the near future, also download maia.usno.navy.mil/ser7/finals.all and 
+ * rename it as eop_finals.txt */
+#define DPSI_DEPS_IAU1980_FILE_EOPC04   "eop_1962_today.txt"
+#define DPSI_DEPS_IAU1980_FILE_FINALS   "eop_finals.txt"
+#define DPSI_DEPS_IAU1980_TJD0_HORIZONS  2437684.5 
+#define HORIZONS_TJD0_DPSI_DEPS_IAU1980  2437684.5 
+#define DPSI_IAU1980_TJD0	(64.284 / 1000.0)  // arcsec
+#define DEPS_IAU1980_TJD0	(6.151 / 1000.0)   // arcsec
+
+/* The above files must be available in order to reproduce JPL Horizons 
+ * in agreement with IERS Conventions 1996 (1992), p. 22. 
+ * Call swe_calc_ut() with iflag|SEFLG_JPLHOR.  
+ * This options works only, if the files DPSI_DEPS_IAU1980_FILE_EOPC04 
+ * and DPSI_DEPS_IAU1980_FILE_FINALS are in the ephemeris path.
+ *
+ * If the software does not find the earth orientation files 
+ * in the ephemeris path, then SEFLG_JPLHOR will run as 
+ * SEFLG_JPLHOR_APPROX.
+ */
+
+/* coordinate transformation */
+extern void swi_coortrf(double *xpo, double *xpn, double eps);
+
+/* coordinate transformation */
+extern void swi_coortrf2(double *xpo, double *xpn, double sineps, double coseps);
+
+/* cartesian to polar coordinates */
+extern void swi_cartpol(double *x, double *l);
+ 
+/* cartesian to polar coordinates with velocity */
+extern void swi_cartpol_sp(double *x, double *l);
+extern void swi_polcart_sp(double *l, double *x);
+ 
+/* polar to cartesian coordinates */
+extern void swi_polcart(double *l, double *x);
+
+/* GCRS to J2000 */
+extern void swi_bias(double *x, double tjd, int32 iflag, AS_BOOL backward);
+extern void swi_get_eop_time_range(void);
+/* GCRS to FK5 */
+extern void swi_icrs2fk5(double *x, int32 iflag, AS_BOOL backward);
+
+/* precession */
+extern int swi_precess(double *R, double J, int32 iflag, int direction );
+extern void swi_precess_speed(double *xx, double t, int32 iflag, int direction);
+
+extern int32 swi_guess_ephe_flag(void);
+
+/* from sweph.c, light deflection, aberration, etc. */
+extern void swi_deflect_light(double *xx, double dt, int32 iflag);
+extern void swi_aberr_light(double *xx, double *xe, int32 iflag);
+extern int swi_plan_for_osc_elem(int32 iflag, double tjd, double *xx);
+extern int swi_trop_ra2sid_lon(double *xin, double *xout, double *xoutr, int32 iflag);
+extern int swi_trop_ra2sid_lon_sosy(double *xin, double *xout, int32 iflag);
+extern int swi_get_observer(double tjd, int32 iflag, 
+	AS_BOOL do_save, double *xobs, char *serr);
+extern void swi_force_app_pos_etc(void);
+
+/* obliquity of ecliptic */
+extern void swi_check_ecliptic(double tjd, int32 iflag);
+extern double swi_epsiln(double J, int32 iflag);
+extern void swi_ldp_peps(double J, double *dpre, double *deps);
+
+/* nutation */
+extern void swi_check_nutation(double tjd, int32 iflag);
+extern int swi_nutation(double J, int32 iflag, double *nutlo);
+extern void swi_nutate(double *xx, int32 iflag, AS_BOOL backward);
+
+extern void swi_mean_lunar_elements(double tjd, 
+							 double *node, double *dnode, 
+							 double *peri, double *dperi);
+/* */
+extern double swi_mod2PI(double x);
+
+/* evaluation of chebyshew series and derivative */
+extern double swi_echeb(double x, double *coef, int ncf);
+extern double swi_edcheb(double x, double *coef, int ncf);
+
+/* cross product of vectors */
+extern void swi_cross_prod(double *a, double *b, double *x);
+/* dot product of vecotrs */
+extern double swi_dot_prod_unit(double *x, double *y);
+
+extern double swi_angnorm(double x);
+
+/* generation of SWISSEPH file names */
+extern void swi_gen_filename(double tjd, int ipli, char *fname);
+
+/* cyclic redundancy checksum (CRC), 32 bit */
+extern uint32 swi_crc32(unsigned char *buf, int len);
+
+extern int swi_cutstr(char *s, char *cutlist, char *cpos[], int nmax);
+extern char *swi_right_trim(char *s);
+
+extern double swi_kepler(double E, double M, double ecce);
+
+extern char *swi_get_fict_name(int32 ipl, char *s);
+
+extern void swi_FK4_FK5(double *xp, double tjd);
+
+extern char *swi_strcpy(char *to, char *from);
+extern char *swi_strncpy(char *to, char *from, size_t n);
+
+extern double swi_deltat_ephe(double tjd_ut, int32 epheflag);
+
+#ifdef TRACE
+#  define TRACE_COUNT_MAX         10000
+  extern TLS FILE *swi_fp_trace_c;
+  extern TLS FILE *swi_fp_trace_out;
+  extern TLS int32 swi_trace_count;
+  extern void swi_open_trace(char *serr);
+  static const char *fname_trace_c = "swetrace.c";
+  static const char *fname_trace_out = "swetrace.txt";
+#ifdef FORCE_IFLAG
+  static const char *fname_force_flg = "force.flg";
+#endif
+#endif /* TRACE */

data/lib/astro_chart/aspects.rb

@@ -0,0 +1,40 @@
+module AstroChart
+  module Aspects
+    # Major aspects: angle => [name, max_orb]
+    MAJOR = {
+      60  => ["六分相", 6],   # Sextile
+      90  => ["四分相", 8],   # Square
+      120 => ["三分相", 8],   # Trine
+      180 => ["對分相", 10],  # Opposition
+    }.freeze
+
+    CONJUNCTION_ORB = 15
+
+    # Calculate the aspect between two ecliptic positions.
+    # Returns [aspect_name, orb] or [nil, nil].
+    def self.calculate(pos1, pos2)
+      return [nil, nil] if pos1.nil? || pos2.nil?
+
+      p1 = pos1 % 360.0
+      p2 = pos2 % 360.0
+
+      diff = (p1 - p2).abs % 360.0
+      diff = 360.0 - diff if diff > 180.0
+
+      # Conjunction (widest orb)
+      if diff <= CONJUNCTION_ORB
+        return ["合相", diff.round(2)]
+      end
+
+      # Check major aspects
+      MAJOR.each do |angle, (name, max_orb)|
+        orb = (diff - angle).abs
+        if orb <= max_orb
+          return [name, orb.round(2)]
+        end
+      end
+
+      [nil, nil]
+    end
+  end
+end

data/lib/astro_chart/chart.rb

@@ -0,0 +1,78 @@
+module AstroChart
+  class Chart
+    attr_reader :birth_date, :birth_time, :latitude, :longitude, :timezone
+
+    # birth_date: "YYYY-MM-DD"
+    # birth_time: "HH:MM"
+    # latitude / longitude: Float
+    # timezone: IANA timezone string (e.g. "Asia/Taipei")
+    def initialize(birth_date:, birth_time:, latitude:, longitude:, timezone:)
+      @birth_date = birth_date
+      @birth_time = birth_time
+      @latitude   = latitude.to_f
+      @longitude  = longitude.to_f
+      @timezone   = timezone
+    end
+
+    # Generate complete chart data.
+    # Returns a Hash with string keys, compatible with the existing JSON API.
+    def generate
+      jd = TimeConversion.to_julian_day(birth_date, birth_time, timezone)
+
+      cusps, ascendant = Houses.calculate(jd, latitude, longitude)
+      asc_zodiac = Zodiac.sign_name(ascendant)
+
+      positions = Planets.calculate_positions(jd)
+      planet_details = Planets.build_details(positions, cusps)
+
+      kp = Planets.key_points_data(positions, cusps, ascendant)
+
+      # Merge aspects into planet details
+      aspect_map = {
+        "太陽"  => "sun_aspects",
+        "月亮"  => "moon_aspects",
+        "土星"  => "saturn_aspects",
+        "金星"  => "venus_aspects",
+        "北交點" => "north_node_aspects",
+        "南交點" => "south_node_aspects",
+      }
+
+      planet_details.each do |planet|
+        key = aspect_map[planet["planet"]]
+        planet["aspects"] = kp[key] if key
+      end
+
+      # Append ruler points
+      planet_details.concat(kp["additional_points"])
+
+      houses_data = cusps.each_with_index.map do |deg, i|
+        {
+          "house_number" => i + 1,
+          "degree"       => deg.round(4),
+          "zodiac"       => Zodiac.sign_name(deg),
+        }
+      end
+
+      {
+        "input" => {
+          "birth_date"  => birth_date,
+          "birth_time"  => birth_time,
+          "coordinates" => {
+            "latitude"  => latitude,
+            "longitude" => longitude,
+          },
+          "timezone"    => timezone,
+        },
+        "chart" => {
+          "ascendant" => {
+            "zodiac"       => asc_zodiac,
+            "degree"       => (ascendant % 30).round(4),
+            "total_degree" => ascendant.round(4),
+          },
+          "planets" => planet_details,
+          "houses"  => houses_data,
+        },
+      }
+    end
+  end
+end

data/lib/astro_chart/ephemeris.rb

@@ -0,0 +1,36 @@
+module AstroChart
+  # Thin Ruby wrapper around the C extension (AstroChart::Ext).
+  # All public callers should go through this module so we can
+  # swap the backend without touching the rest of the code.
+  module Ephemeris
+    PLANETS = {
+      "太陽"   => Ext::SUN,
+      "月亮"   => Ext::MOON,
+      "水星"   => Ext::MERCURY,
+      "金星"   => Ext::VENUS,
+      "火星"   => Ext::MARS,
+      "木星"   => Ext::JUPITER,
+      "土星"   => Ext::SATURN,
+      "天王星"  => Ext::URANUS,
+      "海王星"  => Ext::NEPTUNE,
+      "冥王星"  => Ext::PLUTO,
+      "北交點"  => Ext::TRUE_NODE,
+    }.freeze
+
+    # Convert date/time to Julian Day number.
+    def self.julday(year, month, day, hour)
+      Ext.julday(year, month, day, hour)
+    end
+
+    # Calculate planet ecliptic longitude (degrees 0-360).
+    def self.calc_ut(jd, planet_id)
+      Ext.calc_ut(jd, planet_id)
+    end
+
+    # Calculate house cusps + ascendant.
+    # Returns { "cusps" => [12 floats], "ascendant" => float, "mc" => float }
+    def self.houses(jd, latitude, longitude, system = "P")
+      Ext.houses(jd, latitude, longitude, system.ord)
+    end
+  end
+end

data/lib/astro_chart/houses.rb

@@ -0,0 +1,30 @@
+module AstroChart
+  module Houses
+    # Calculate house cusps and ascendant from Julian Day + coordinates.
+    # Returns [cusps_array(12), ascendant_degree].
+    def self.calculate(jd, latitude, longitude, system = "P")
+      data = Ephemeris.houses(jd, latitude, longitude, system)
+      [data["cusps"], data["ascendant"]]
+    end
+
+    # Determine which house (1-12) a planet falls in based on house cusps.
+    def self.find_house(planet_pos, cusps)
+      return nil if planet_pos.nil? || cusps.nil? || cusps.empty?
+
+      n = cusps.length
+      base = cusps[0]
+
+      # Normalize house cusps relative to first house cusp
+      norm_cusps = cusps.map { |h| h >= base ? h : h + 360.0 }
+      norm_pos = planet_pos >= base ? planet_pos : planet_pos + 360.0
+
+      n.times do |i|
+        start_deg = norm_cusps[i]
+        end_deg = i < n - 1 ? norm_cusps[i + 1] : norm_cusps[0] + 360.0
+        return i + 1 if norm_pos >= start_deg && norm_pos < end_deg
+      end
+
+      nil
+    end
+  end
+end

data/lib/astro_chart/planets.rb

@@ -0,0 +1,122 @@
+module AstroChart
+  module Planets
+    # Calculate raw ecliptic longitudes for all planets + nodes.
+    # Returns Hash: { "太陽" => 123.45, "月亮" => 67.89, ... }
+    def self.calculate_positions(jd)
+      positions = {}
+
+      Ephemeris::PLANETS.each do |name, planet_id|
+        positions[name] = Ephemeris.calc_ut(jd, planet_id)
+      end
+
+      # South node = North node + 180
+      positions["南交點"] = (positions["北交點"] + 180.0) % 360.0
+
+      positions
+    end
+
+    # Build detailed planet info list with zodiac, house, degree.
+    def self.build_details(positions, cusps)
+      positions.map do |planet, pos|
+        house = Houses.find_house(pos, cusps)
+        next if house.nil?
+
+        {
+          "planet"       => planet,
+          "zodiac"       => Zodiac.sign_name(pos),
+          "house"        => house,
+          "degree"       => (pos % 30).round(4),
+          "total_degree" => pos.round(4),
+        }
+      end.compact
+    end
+
+    # Calculate aspects for a given point against all planet positions.
+    def self.aspects_for(point_pos, point_name, positions)
+      return [] if point_pos.nil?
+
+      results = []
+      positions.each do |planet, pos|
+        next if planet == point_name
+
+        aspect_type, orb = Aspects.calculate(point_pos, pos)
+        if aspect_type
+          results << {
+            "planet"      => planet,
+            "aspect_type" => aspect_type,
+            "orb"         => orb,
+          }
+        end
+      end
+      results
+    end
+
+    # Calculate key points data: aspects for major planets + ruler info.
+    def self.key_points_data(positions, cusps, ascendant)
+      sun_aspects         = aspects_for(positions["太陽"],  "太陽",  positions)
+      moon_aspects        = aspects_for(positions["月亮"],  "月亮",  positions)
+      saturn_aspects      = aspects_for(positions["土星"],  "土星",  positions)
+      venus_aspects       = aspects_for(positions["金星"],  "金星",  positions)
+      north_node_aspects  = aspects_for(positions["北交點"], "北交點", positions)
+      south_node_aspects  = aspects_for(positions["南交點"], "南交點", positions)
+
+      additional_points = []
+
+      # North Node ruler
+      nn_pos = positions["北交點"]
+      if nn_pos
+        nn_zodiac = Zodiac.sign_name(nn_pos)
+        nn_ruler  = Zodiac.ruler(nn_zodiac)
+        nn_ruler_pos = positions[nn_ruler]
+        if nn_ruler_pos
+          additional_points << build_ruler_point("北交點定位星", nn_ruler, nn_ruler_pos, cusps, positions)
+        end
+      end
+
+      # South Node ruler
+      sn_pos = positions["南交點"]
+      if sn_pos
+        sn_zodiac = Zodiac.sign_name(sn_pos)
+        sn_ruler  = Zodiac.ruler(sn_zodiac)
+        sn_ruler_pos = positions[sn_ruler]
+        if sn_ruler_pos
+          additional_points << build_ruler_point("南交點定位星", sn_ruler, sn_ruler_pos, cusps, positions)
+        end
+      end
+
+      # Ascendant ruler
+      if ascendant
+        asc_zodiac = Zodiac.sign_name(ascendant)
+        asc_ruler  = Zodiac.ruler(asc_zodiac)
+        asc_ruler_pos = positions[asc_ruler]
+        if asc_ruler_pos
+          additional_points << build_ruler_point("上升星座定位星", asc_ruler, asc_ruler_pos, cusps, positions)
+        end
+      end
+
+      {
+        "sun_aspects"        => sun_aspects,
+        "moon_aspects"       => moon_aspects,
+        "saturn_aspects"     => saturn_aspects,
+        "venus_aspects"      => venus_aspects,
+        "north_node_aspects" => north_node_aspects,
+        "south_node_aspects" => south_node_aspects,
+        "additional_points"  => additional_points,
+      }
+    end
+
+    private_class_method
+
+    def self.build_ruler_point(label, ruling_planet, pos, cusps, positions)
+      {
+        "planet"        => label,
+        "ruling_planet" => ruling_planet,
+        "zodiac"        => Zodiac.sign_name(pos),
+        "house"         => Houses.find_house(pos, cusps),
+        "degree"        => (pos % 30).round(4),
+        "total_degree"  => pos.round(4),
+        "aspects"       => aspects_for(pos, ruling_planet, positions),
+      }
+    end
+  end
+end

data/lib/astro_chart/time_conversion.rb

@@ -0,0 +1,30 @@
+require "tzinfo"
+
+module AstroChart
+  module TimeConversion
+    # Convert a local date/time + timezone to Julian Day.
+    #
+    # birth_date: "YYYY-MM-DD"
+    # birth_time: "HH:MM"
+    # timezone:   IANA timezone string, e.g. "Asia/Taipei"
+    def self.to_julian_day(birth_date, birth_time, timezone)
+      tz = TZInfo::Timezone.get(timezone)
+
+      parts = birth_date.split("-").map(&:to_i)
+      year, month, day = parts
+
+      time_parts = birth_time.split(":").map(&:to_i)
+      hour, minute = time_parts
+
+      # Build local time and convert to UTC
+      local_time = Time.new(year, month, day, hour, minute, 0, tz.observed_utc_offset)
+
+      # TZInfo gives us the proper UTC offset; convert to UTC
+      utc_time = tz.local_to_utc(Time.new(year, month, day, hour, minute, 0))
+
+      ut_hour = utc_time.hour + utc_time.min / 60.0
+
+      Ephemeris.julday(utc_time.year, utc_time.month, utc_time.day, ut_hour)
+    end
+  end
+end

data/lib/astro_chart/version.rb

@@ -0,0 +1,3 @@
+module AstroChart
+  VERSION = "0.1.0"
+end

data/lib/astro_chart/zodiac.rb

@@ -0,0 +1,37 @@
+module AstroChart
+  module Zodiac
+    SIGNS = [
+      "牡羊座", "金牛座", "雙子座", "巨蟹座",
+      "獅子座", "處女座", "天秤座", "天蠍座",
+      "射手座", "摩羯座", "水瓶座", "雙魚座"
+    ].freeze
+
+    RULERS = {
+      "牡羊座" => "火星",
+      "金牛座" => "金星",
+      "雙子座" => "水星",
+      "巨蟹座" => "月亮",
+      "獅子座" => "太陽",
+      "處女座" => "水星",
+      "天秤座" => "金星",
+      "天蠍座" => "冥王星",
+      "射手座" => "木星",
+      "摩羯座" => "土星",
+      "水瓶座" => "天王星",
+      "雙魚座" => "海王星",
+    }.freeze
+
+    # Return zodiac sign name for a given ecliptic longitude.
+    def self.sign_name(degree)
+      index = (degree % 360).floor / 30
+      SIGNS[index]
+    end
+
+    # Return the ruling planet of a zodiac sign.
+    def self.ruler(zodiac)
+      r = RULERS[zodiac]
+      raise ArgumentError, "無效的星座名稱: #{zodiac}" if r.nil?
+      r
+    end
+  end
+end

data/lib/astro_chart.rb

@@ -0,0 +1,12 @@
+require_relative "astro_chart/version"
+require_relative "astro_chart/astro_chart_ext"
+require_relative "astro_chart/ephemeris"
+require_relative "astro_chart/zodiac"
+require_relative "astro_chart/aspects"
+require_relative "astro_chart/houses"
+require_relative "astro_chart/time_conversion"
+require_relative "astro_chart/planets"
+require_relative "astro_chart/chart"
+
+module AstroChart
+end