Ruby - 日・月の出・入・南中計算 gem の作成!

Updated:


※この記事は8年以上前に投稿されたもので、情報が古い可能性があります。

当方、 「日の出・日の入りの計算―天体の出没時刻の求め方」を参考に日・月の出・入・南中を計算する Ruby スクリプトを作成しておりましたが、あらゆる面で流用したくなったために、今回 RubyGems ライブラリにし公開することとしました。

以下では、今回作成した gem の簡単な利用方法をご紹介します。

0. 前提条件Permalink

  • Ruby 2.3.1-p112 での作業を想定。
  • 自作した gem ライブラリの名称は “mk_sunmoon” で、計算対象年月日は 0000-01-01 〜 9999-01-01。
  • 当ライブラリの計算可能項目
    • 日の出(時刻、方位角)
    • 日南中(時刻、高度)
    • 日の入(時刻、方位角)
    • 月の出(時刻、方位角)
    • 月南中(時刻、高度)
    • 月の入(時刻、方位角)

1. インストールPermalink

$ sudo gem install mk_sunmoon

2. Ruby スクリプトの作成例Permalink

File: sun_moon.rb

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
#! /usr/local/bin/ruby
# coding: utf-8
# 日の出・日の入の時刻・方位、日の南中の時刻・高度、
# 月の出・月の入の時刻・方位、月の南中の時刻・高度  の算出
#
# date          name            version
# 2016.06.13    mk-mode         1.00 新規作成
#
# Copyright(C) 2016 mk-mode.com All Rights Reserved.
#---------------------------------------------------------------------------------
# [引数] 第1: 日付 (Format: 99999999, Default: システム日付)
#              計算対象の日付(グレゴリオ暦)(JST)を半角8桁数字で指定。
#        第2: 緯度 (Format: [+-]999.99999999, Default: 35.47222222)
#              (度で指定。北緯は +, 南緯は -, 小数点以下は省略可)
#        第3: 経度 (Format: [+-]999.99999999, Default: 133.05055556)
#              (度で指定。東経は +, 西経は - ,小数点以下は省略可)
#        第4: 標高 (Format: [+]9999.99999999, Default: 0.0)
#              (メートルで指定。小数点以下は省略可)
#---------------------------------------------------------------------------------
#+
require 'mk_sunmoon'

class SunMoon
  USAGE =<<-EOS
  [USAGE] ./sun_moon.rb [date[, latitude[, longitude[, altitude]]]]
  引数
  第1: 日付 (Format: 99999999, Default: システム日付)
        計算対象の日付(グレゴリオ暦)(JST)を半角8桁数字で指定。
  第2: 緯度 (Format: [+-]999.99999999, Default: 35.47222222)
        (度で指定。北緯は +, 南緯は -, 小数点以下は省略可)
  第3: 経度 (Format: [+-]999.99999999, Default: 133.05055556)
        (度で指定。東経は +, 西経は - ,小数点以下は省略可)
  第4: 標高 (Format: [+]9999.99999999, Default: 0.0)
        (メートルで指定。小数点以下は省略可)
EOS
  def initialize
    check_args(ARGV)
    @obj = MkSunmoon.new(*ARGV)
  end

  def exec
    begin
      sign_lat = @obj.lat < 0 ? "S" : "N"
      sign_lon = @obj.lon < 0 ? "W" : "E"
      sr = @obj.sunrise
      ss = @obj.sunset
      sm = @obj.sun_mp
      mr = @obj.moonrise
      ms = @obj.moonset
      mm = @obj.moon_mp
      str  = sprintf("  %04d-%02d-%02d", @obj.year, @obj.month, @obj.day)
      str << " [#{@obj.lat}#{sign_lat}, #{@obj.lon}#{sign_lon}, #{@obj.alt}m]\n"
      str << "  日の出: #{sr[0]} (方位角: "
      str << (sr[1] == "---" ? "--------°)\n" : sprintf("%8.4f°)\n", sr[1]))
      str << "  日南中: #{sm[0]} (  高度: "
      str << (sm[1] == "---" ? "--------°)\n" : sprintf("%8.4f°)\n", sm[1]))
      str << "  日の入: #{ss[0]} (方位角: "
      str << (ss[1] == "---" ? "--------°)\n" : sprintf("%8.4f°)\n", ss[1]))
      str << "  月の出: #{mr[0]} (方位角: "
      str << (mr[1] == "---" ? "--------°)\n" : sprintf("%8.4f°)\n", mr[1]))
      str << "  月南中: #{mm[0]} (  高度: "
      str << (mm[1] == "---" ? "--------°)\n" : sprintf("%8.4f°)\n", mm[1]))
      str << "  月の入: #{ms[0]} (方位角: "
      str << (ms[1] == "---" ? "--------°)\n" : sprintf("%8.4f°)\n", ms[1]))
      puts str
    rescue => e
      $stderr.puts "[#{e.class}] #{e.message}"
      e.backtrace.each { |tr| $stderr.puts "\t#{tr}" }
      exit 1
    end
  end

  private

  def check_args(args)
    date, lat, lon, alt = args

    begin
      if (date && date !~ /^\d{8}$/) ||
         (lat && lat !~ /^[+-]?\d{,2}(\.\d*)?$/) ||
         (lon && lon !~ /^[+-]?\d{,3}(\.\d*)?$/) ||
         (alt && alt !~ /^+?\d{,4}(\.\d*)?$/)
        puts USAGE
        exit 0
      end
    rescue => e
      raise
    end
  end
end

SunMoon.new.exec if __FILE__ == $0

3. サンプルスクリプトの実行Permalink

$ ./sun_moon.rb 20160613 35.47222222 133.05055556 0
  2016-06-13 [35.47222222N, 133.05055556E, 0.0m]
  日の出: 04:51:52 (方位角:  60.3626°)
  日南中: 12:07:52 (  高度:  77.7587°)
  日の入: 19:23:59 (方位角: 299.6796°)
  月の出: 12:48:14 (方位角:  89.5617°)
  月南中: 18:58:52 (  高度:  54.1050°)
  月の入: 00:32:24 (方位角: 272.8283°)

4. gem ライブラリPermalink


以前作成した Ruby スクリプトを何かと応用したかったので gem ライブラリ化した次第です。

以上。





 

Sponsored Link

 

Comments