skymap 0.1.0 → 0.2.0

data/lib/skymap/catalog/constellation_boundaries.rb

@@ -0,0 +1,43 @@
+# frozen_string_literal: true
+
+require "csv"
+
+module Skymap
+  module Catalog
+    class ConstellationBoundaries
+      include Enumerable
+
+      PATH = File.expand_path(
+        "../../../data/constellations/iau/boundaries.csv",
+        __dir__
+      )
+
+      def each
+        CSV.foreach(PATH, headers: true).chunk_while do |previous, row|
+          previous["border"] == row["border"]
+        end.each do |rows|
+          yield boundary_from(rows)
+        end
+      end
+
+      private
+
+      def boundary_from(rows)
+        ConstellationBoundary.new(
+          constellations: rows.first["constellations"].split,
+          points: rows.map { |row| coordinates_from(row) }
+        )
+      end
+
+      def coordinates_from(row)
+        Astronoby::Coordinates::Equatorial.new(
+          right_ascension: Astronoby::Angle.from_hours(
+            row["right_ascension"].to_f
+          ),
+          declination: Astronoby::Angle.from_degrees(row["declination"].to_f),
+          epoch: Astronoby::JulianDate::J2000
+        )
+      end
+    end
+  end
+end

data/lib/skymap/catalog/constellation_lines.rb

@@ -8,7 +8,7 @@ module Skymap
       include Enumerable
 
       PATH = File.expand_path(
-        "../../../data/constellations/iau.csv",
+        "../../../data/constellations/iau/lines.csv",
         __dir__
       )
 

data/lib/skymap/constellation_boundary.rb

@@ -0,0 +1,5 @@
+# frozen_string_literal: true
+
+module Skymap
+  ConstellationBoundary = Data.define(:constellations, :points)
+end

data/lib/skymap/context.rb

@@ -2,7 +2,7 @@
 
 module Skymap
   class Context
-    attr_reader :canvas
+    attr_reader :observer, :instant, :canvas
 
     def initialize(observer:, instant:, canvas:)
       @observer = observer
@@ -13,7 +13,10 @@ module Skymap
     end
 
     def horizontal(star)
-      coordinates = star.equatorial_coordinates
+      horizontal_coordinates(star.equatorial_coordinates)
+    end
+
+    def horizontal_coordinates(coordinates)
       key = [
         coordinates.right_ascension,
         coordinates.declination,

data/lib/skymap/horizon.rb

@@ -0,0 +1,32 @@
+# frozen_string_literal: true
+
+module Skymap
+  module Horizon
+    module_function
+
+    def direction(horizontals)
+      vectors = horizontals.map do |horizontal|
+        [
+          horizontal.altitude.cos * horizontal.azimuth.cos,
+          horizontal.altitude.cos * horizontal.azimuth.sin,
+          horizontal.altitude.sin
+        ]
+      end
+      vectors.transpose.map(&:sum)
+    end
+
+    def crossing(above, below)
+      above.zip(below).map do |from, to|
+        above.last * to - below.last * from
+      end
+    end
+
+    def project(context, direction)
+      x, y, = direction
+      context.project(
+        altitude: Astronoby::Angle.zero,
+        azimuth: Astronoby::Angle.from_radians(Math.atan2(y, x))
+      )
+    end
+  end
+end

data/lib/skymap/layers/constellation_boundaries.rb

@@ -0,0 +1,67 @@
+# frozen_string_literal: true
+
+module Skymap
+  module Layers
+    class ConstellationBoundaries
+      LINE_WIDTH_RATIO = 0.0017
+
+      def initialize(boundaries:)
+        @boundaries = boundaries
+      end
+
+      def elements(context)
+        width = LINE_WIDTH_RATIO * context.canvas.radius
+
+        @boundaries.flat_map do |boundary|
+          visible_runs(context, boundary.points).map do |points|
+            Renderer::Boundary.new(points: points, width: width)
+          end
+        end
+      end
+
+      private
+
+      def visible_runs(context, coordinates)
+        horizontals = coordinates.map do |point|
+          context.horizontal_coordinates(point)
+        end
+        path = [position(context, horizontals.first)] +
+          horizontals.each_cons(2).flat_map do |from, to|
+            crossing(context, from, to) + [position(context, to)]
+          end
+
+        path.chunk { |point| !point.nil? }.filter_map do |visible, points|
+          points if visible && points.size > 1
+        end
+      end
+
+      def position(context, horizontal)
+        return if below?(horizontal)
+
+        context.project(
+          altitude: horizontal.altitude,
+          azimuth: horizontal.azimuth
+        )
+      end
+
+      def crossing(context, from, to)
+        return [] if below?(from) == below?(to)
+
+        above, below = below?(from) ? [to, from] : [from, to]
+        [
+          Horizon.project(
+            context,
+            Horizon.crossing(
+              Horizon.direction([above]),
+              Horizon.direction([below])
+            )
+          )
+        ]
+      end
+
+      def below?(horizontal)
+        horizontal.altitude.negative?
+      end
+    end
+  end
+end

data/lib/skymap/layers/constellation_lines.rb

@@ -18,15 +18,17 @@ module Skymap
           next if ends.flatten.include?(nil)
 
           directions = ends.map do |end_stars|
-            direction(end_stars.map { |star| context.horizontal(star) })
+            Horizon.direction(
+              end_stars.map { |star| context.horizontal(star) }
+            )
           end
           next if directions.all? { |direction| direction.last.negative? }
 
           from, to = ends.each_with_index.map do |end_stars, index|
             if directions[index].last.negative?
-              project_direction(
+              Horizon.project(
                 context,
-                horizon_crossing(directions[1 - index], directions[index])
+                Horizon.crossing(directions[1 - index], directions[index])
               )
             else
               Point.midpoint(end_stars.map { |star| context.position(star) })
@@ -39,33 +41,6 @@ module Skymap
           )
         end
       end
-
-      private
-
-      def direction(horizontals)
-        vectors = horizontals.map do |horizontal|
-          [
-            horizontal.altitude.cos * horizontal.azimuth.cos,
-            horizontal.altitude.cos * horizontal.azimuth.sin,
-            horizontal.altitude.sin
-          ]
-        end
-        vectors.transpose.map(&:sum)
-      end
-
-      def horizon_crossing(above, below)
-        above.zip(below).map do |from, to|
-          above.last * to - below.last * from
-        end
-      end
-
-      def project_direction(context, direction)
-        x, y, = direction
-        context.project(
-          altitude: Astronoby::Angle.zero,
-          azimuth: Astronoby::Angle.from_radians(Math.atan2(y, x))
-        )
-      end
     end
   end
 end

data/lib/skymap/layers/moon.rb

@@ -0,0 +1,40 @@
+# frozen_string_literal: true
+
+module Skymap
+  module Layers
+    class Moon
+      MOON_RADIUS_RATIO = 0.03
+
+      def initialize(ephem:)
+        @ephem = ephem
+      end
+
+      def elements(context)
+        moon = Astronoby::Moon.at(context.instant, ephem: @ephem)
+        horizontal = moon.observed_by(context.observer).horizontal
+        return [] if horizontal.altitude.negative?
+
+        [
+          Renderer::Moon.new(
+            center: context.project(
+              altitude: horizontal.altitude,
+              azimuth: horizontal.azimuth
+            ),
+            radius: MOON_RADIUS_RATIO * context.canvas.radius,
+            illuminated_fraction: moon.illuminated_fraction,
+            rotation: rotation(moon, context.observer)
+          )
+        ]
+      end
+
+      private
+
+      def rotation(moon, observer)
+        from_zenith = moon.bright_limb_position_angle.degrees -
+          moon.parallactic_angle(observer: observer).degrees
+
+        -90 - from_zenith
+      end
+    end
+  end
+end

data/lib/skymap/renderer/boundary.rb

@@ -0,0 +1,7 @@
+# frozen_string_literal: true
+
+module Skymap
+  module Renderer
+    Boundary = Data.define(:points, :width)
+  end
+end

data/lib/skymap/renderer/moon.rb

@@ -0,0 +1,7 @@
+# frozen_string_literal: true
+
+module Skymap
+  module Renderer
+    Moon = Data.define(:center, :radius, :illuminated_fraction, :rotation)
+  end
+end

data/lib/skymap/renderer/svg.rb

@@ -10,6 +10,11 @@ module Skymap
       HORIZON_WIDTH = 3
       STAR_COLOR = "#ffffff"
       LINE_COLOR = "#46689c"
+      BOUNDARY_COLOR = "#2b3f5e"
+      BOUNDARY_DASH_RATIO = 4
+      MOON_COLOR = "#f1ecd6"
+      MOON_SHADOW_COLOR = "#232a35"
+      MOON_SHADOW_OPACITY = 0.6
       LABEL_COLOR = HORIZON_COLOR
       LABEL_FONT = "sans-serif"
 
@@ -35,9 +40,11 @@ module Skymap
 
         elements.each do |element|
           case element
+          when Boundary then draw_boundary(svg, element)
           when Line then draw_line(svg, element)
           when Dot then draw_dot(svg, element)
           when Label then draw_label(svg, element)
+          when Moon then draw_moon(svg, element)
           else raise ArgumentError, "unknown element: #{element.class}"
           end
         end
@@ -59,6 +66,20 @@ module Skymap
         )
       end
 
+      def draw_boundary(svg, boundary)
+        dash = (boundary.width * BOUNDARY_DASH_RATIO).round(2)
+
+        svg.polyline(
+          points: boundary.points.map { |point|
+            "#{point.x.round(2)},#{point.y.round(2)}"
+          }.join(" "),
+          fill: "none",
+          stroke: BOUNDARY_COLOR,
+          stroke_width: boundary.width.round(2),
+          stroke_dasharray: "#{dash} #{dash}"
+        )
+      end
+
       def draw_dot(svg, dot)
         svg.circle(
           cx: dot.center.x.round(2),
@@ -68,6 +89,32 @@ module Skymap
         )
       end
 
+      def draw_moon(svg, moon)
+        x = moon.center.x.round(2)
+        y = moon.center.y.round(2)
+        radius = moon.radius.round(2)
+        fraction = moon.illuminated_fraction
+        terminator = (radius * (1 - 2 * fraction).abs).round(2)
+        bulge = (fraction > 0.5) ? 1 : 0
+        top = "#{x},#{(y - radius).round(2)}"
+        bottom = "#{x},#{(y + radius).round(2)}"
+        lit = "M #{top} A #{radius},#{radius} 0 0 1 #{bottom} " \
+          "A #{terminator},#{radius} 0 0 #{bulge} #{top} Z"
+
+        svg.build do
+          g(transform: "rotate(#{moon.rotation.round(2)} #{x} #{y})") do
+            circle(
+              cx: x,
+              cy: y,
+              r: radius,
+              fill: MOON_SHADOW_COLOR,
+              fill_opacity: MOON_SHADOW_OPACITY
+            )
+            path(d: lit, fill: MOON_COLOR)
+          end
+        end
+      end
+
       def draw_label(svg, label)
         svg.text(
           label.text,

data/lib/skymap/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module Skymap
-  VERSION = "0.1.0"
+  VERSION = "0.2.0"
 end

data/lib/skymap.rb

@@ -3,19 +3,26 @@
 require "astronoby"
 
 require_relative "skymap/canvas"
+require_relative "skymap/catalog/constellation_boundaries"
 require_relative "skymap/catalog/constellation_lines"
 require_relative "skymap/catalog/stars"
 require_relative "skymap/chart"
+require_relative "skymap/constellation_boundary"
 require_relative "skymap/constellation_line"
 require_relative "skymap/context"
+require_relative "skymap/horizon"
 require_relative "skymap/layers/cardinal_directions"
+require_relative "skymap/layers/constellation_boundaries"
 require_relative "skymap/layers/constellation_lines"
+require_relative "skymap/layers/moon"
 require_relative "skymap/layers/stars"
 require_relative "skymap/point"
 require_relative "skymap/projection/stereographic"
+require_relative "skymap/renderer/boundary"
 require_relative "skymap/renderer/dot"
 require_relative "skymap/renderer/label"
 require_relative "skymap/renderer/line"
+require_relative "skymap/renderer/moon"
 require_relative "skymap/renderer/svg"
 require_relative "skymap/star"
 require_relative "skymap/star_size"

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: skymap
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.2.0
 platform: ruby
 authors:
 - Rémy Hannequin
@@ -145,23 +145,31 @@ files:
 - CHANGELOG.md
 - LICENSE.txt
 - README.md
-- data/constellations/iau.csv
+- data/constellations/iau/boundaries.csv
+- data/constellations/iau/lines.csv
 - data/stars.csv
 - lib/skymap.rb
 - lib/skymap/canvas.rb
+- lib/skymap/catalog/constellation_boundaries.rb
 - lib/skymap/catalog/constellation_lines.rb
 - lib/skymap/catalog/stars.rb
 - lib/skymap/chart.rb
+- lib/skymap/constellation_boundary.rb
 - lib/skymap/constellation_line.rb
 - lib/skymap/context.rb
+- lib/skymap/horizon.rb
 - lib/skymap/layers/cardinal_directions.rb
+- lib/skymap/layers/constellation_boundaries.rb
 - lib/skymap/layers/constellation_lines.rb
+- lib/skymap/layers/moon.rb
 - lib/skymap/layers/stars.rb
 - lib/skymap/point.rb
 - lib/skymap/projection/stereographic.rb
+- lib/skymap/renderer/boundary.rb
 - lib/skymap/renderer/dot.rb
 - lib/skymap/renderer/label.rb
 - lib/skymap/renderer/line.rb
+- lib/skymap/renderer/moon.rb
 - lib/skymap/renderer/svg.rb
 - lib/skymap/star.rb
 - lib/skymap/star_size.rb