mk_coord 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 057fe547daa1b16ae2cc6ee1d096bb05624f53c1
+  data.tar.gz: 25945cca78c92dca05268b5be3802f56bc53c972
+SHA512:
+  metadata.gz: 16b4b5547a73286dc8d56bc17fa5df76c59f01853a747472266385bba6aebc1b118136e026f79c3330700abbab03ce270777d40f14e85e9501038980fa05e228
+  data.tar.gz: f21b48b716863f1358164e0b4ec02fafdd2005a8319e223fff02c6ec73091aea3c0c68e7cca0efa8edb5c2e826e00a32c856a6c4b20e23aad34c1260f612de37

data/.gitignore

@@ -0,0 +1,9 @@
+/.bundle/
+/.yardoc
+/Gemfile.lock
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/

data/.rspec

@@ -0,0 +1,2 @@
+--format documentation
+--color

data/.travis.yml

@@ -0,0 +1,5 @@
+sudo: false
+language: ruby
+rvm:
+  - 2.3.1
+before_install: gem install bundler -v 1.12.5

data/Gemfile

@@ -0,0 +1,9 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in mk_coord.gemspec
+gemspec
+
+group :development do
+  gem "guard"
+  gem "guard-rspec", "~> 4.7.0"
+end

data/Guardfile

@@ -0,0 +1,42 @@
+# A sample Guardfile
+# More info at https://github.com/guard/guard#readme
+
+## Uncomment and set this to only include directories you want to watch
+# directories %w(app lib config test spec features) \
+#  .select{|d| Dir.exists?(d) ? d : UI.warning("Directory #{d} does not exist")}
+
+## Note: if you are using the `directories` clause above and you are not
+## watching the project directory ('.'), then you will want to move
+## the Guardfile to a watched dir and symlink it back, e.g.
+#
+#  $ mkdir config
+#  $ mv Guardfile config/
+#  $ ln -s config/Guardfile .
+#
+# and, you'll have to watch "config/Guardfile" instead of "Guardfile"
+
+# Note: The cmd option is now required due to the increasing number of ways
+#       rspec may be run, below are examples of the most common uses.
+#  * bundler: 'bundle exec rspec'
+#  * bundler binstubs: 'bin/rspec'
+#  * spring: 'bin/rspec' (This will use spring if running and you have
+#                          installed the spring binstubs per the docs)
+#  * zeus: 'zeus rspec' (requires the server to be started separately)
+#  * 'just' rspec: 'rspec'
+
+guard :rspec, cmd: "bundle exec rspec" do
+  require "guard/rspec/dsl"
+  dsl = Guard::RSpec::Dsl.new(self)
+
+  # Feel free to open issues for suggestions and improvements
+
+  # RSpec files
+  rspec = dsl.rspec
+  watch(rspec.spec_helper) { rspec.spec_dir }
+  watch(rspec.spec_support) { rspec.spec_dir }
+  watch(rspec.spec_files)
+
+  # Ruby files
+  ruby = dsl.ruby
+  dsl.watch_spec_files_for(ruby.lib_files)
+end

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2016 Masaru Koizumi
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.

data/README.md

@@ -0,0 +1,81 @@
+# MkCoord
+
+## Introduction
+
+This is the gem library which converts coordinates.
+
+## Computable items
+
+* Rectangular position:
+  - Equator coordinates -> Ecliptic coordinates
+  - Ecliptic coordinates -> Equator coordinates
+* Polar position:
+  - Equator coordinates -> Ecliptic coordinates
+  - Ecliptic coordinates -> Equator coordinates
+* Rectangular coordinates -> Polar coordinates
+* Polar coordinates -> Rectangular coordinates
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'mk_coord'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install mk_coord
+
+## Usage
+
+    require 'mk_coord'
+    
+    eps = 0.409
+    pos_rect = [
+       0.99443659220700997281,
+      -0.03816291768957833647,
+      -0.01655177670960058384
+    ]
+    pos_pol  = [
+       6.24482770879939,
+      -0.016630599800372015
+    ]
+    r = 0.9951686008947793
+    p pos_rect
+    p pos_pol
+    p r
+    puts "---"
+    p MkCoord.rect_eq2ec(pos_rect, eps)
+    p MkCoord.rect_ec2eq(pos_rect, eps)
+    p MkCoord.pol_eq2ec(pos_pol, eps)
+    p MkCoord.pol_ec2eq(pos_pol, eps)
+    p MkCoord.rect2pol(pos_rect)
+    p MkCoord.pol2rect(pos_pol, r)
+
+About the units:
+
+* The unit of epsilon(eps) should be radian.
+* The unit of polar postion(pos_pol) should be radian.
+* Everything is good for the unit of rectangular postion(pos_pol).
+* Everything is good for the unit of distance(r).
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake spec` to run the tests. You can also run `bin/console` for an interactive prompt that will allow you to experiment.
+
+To install this gem onto your local machine, run `bundle exec rake install`. To release a new version, update the version number in `version.rb`, and then run `bundle exec rake release`, which will create a git tag for the version, push git commits and tags, and push the `.gem` file to [rubygems.org](https://rubygems.org).
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/komasaru/mk_coord.
+
+
+## License
+
+The gem is available as open source under the terms of the [MIT License](http://opensource.org/licenses/MIT).
+

data/Rakefile

@@ -0,0 +1,6 @@
+require "bundler/gem_tasks"
+require "rspec/core/rake_task"
+
+RSpec::Core::RakeTask.new(:spec)
+
+task :default => :spec

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "mk_coord"
+
+# You can add fixtures and/or initialization code here to make experimenting
+# with your gem easier. You can also use a different console, if you like.
+
+# (If you use this, don't forget to add pry to your Gemfile!)
+# require "pry"
+# Pry.start
+
+require "irb"
+IRB.start

data/bin/setup

@@ -0,0 +1,8 @@
+#!/usr/bin/env bash
+set -euo pipefail
+IFS=$'\n\t'
+set -vx
+
+bundle install
+
+# Do any other automated setup that you need to do here

data/exe/mk_coord

@@ -0,0 +1,26 @@
+#!/usr/bin/env ruby
+
+require "mk_coord"
+
+eps = 23.43929 * MkCoord::Const::PI / 180.0
+pos_rect = [
+   0.99443659220700997281,
+  -0.03816291768957833647,
+  -0.01655177670960058384
+]
+pos_pol  = [
+   6.24482770879939,
+  -0.016630599800372015
+]
+r = 0.9951686008947793
+p pos_rect
+p pos_pol
+p r
+puts "---"
+p MkCoord.rect_eq2ec(pos_rect, eps)
+p MkCoord.rect_ec2eq(pos_rect, eps)
+p MkCoord.pol_eq2ec(pos_pol, eps)
+p MkCoord.pol_ec2eq(pos_pol, eps)
+p MkCoord.rect2pol(pos_rect)
+p MkCoord.pol2rect(pos_pol, r)
+

data/lib/mk_coord/const.rb

@@ -0,0 +1,11 @@
+module MkCoord
+  module Const
+    PI  = 3.141592653589793238462
+    PI2 = PI * 2
+    R_UNIT = [
+      [1.0, 0.0, 0.0],
+      [0.0, 1.0, 0.0],
+      [0.0, 0.0, 1.0]
+    ].freeze
+  end
+end

data/lib/mk_coord/matrix.rb

@@ -0,0 +1,102 @@
+module MkCoord
+  module Matrix
+    module_function
+
+    #=========================================================================
+    # 回転行列（x 軸中心）
+    #
+    # @param:  phi  (回転量(rad))
+    # @param:  r    (回転行列)
+    # @return: r_a  (回転後)
+    #=========================================================================
+    def r_x(phi, r = R_UNIT)
+      r_a = Array.new
+
+      begin
+        s, c = Math.sin(phi), Math.cos(phi)
+        a10 =   c * r[1][0] + s * r[2][0]
+        a11 =   c * r[1][1] + s * r[2][1]
+        a12 =   c * r[1][2] + s * r[2][2]
+        a20 = - s * r[1][0] + c * r[2][0]
+        a21 = - s * r[1][1] + c * r[2][1]
+        a22 = - s * r[1][2] + c * r[2][2]
+        r_a << r[0]
+        r_a << [a10, a11, a12]
+        r_a << [a20, a21, a22]
+        return r_a
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # 回転行列（y 軸中心）
+    #
+    # @param:  theta  (回転量(rad))
+    # @param:  r      (回転行列)
+    # @return: r_a    (回転後)
+    #=========================================================================
+    def r_y(theta, r = R_UNIT)
+      r_a = Array.new
+
+      begin
+        s, c = Math.sin(theta), Math.cos(theta)
+        a00 = c * r[0][0] - s * r[2][0]
+        a01 = c * r[0][1] - s * r[2][1]
+        a02 = c * r[0][2] - s * r[2][2]
+        a20 = s * r[0][0] + c * r[2][0]
+        a21 = s * r[0][1] + c * r[2][1]
+        a22 = s * r[0][2] + c * r[2][2]
+        r_a << [a00, a01, a02]
+        r_a << r[1]
+        r_a << [a20, a21, a22]
+        return r_a
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # 回転行列（z 軸中心）
+    #
+    # @param:  psi  (回転量(rad))
+    # @param:  r    (回転行列)
+    # @return: r_a  (回転後)
+    #=========================================================================
+    def r_z(psi, r = R_UNIT)
+      r_a = Array.new
+
+      begin
+        s, c = Math.sin(psi), Math.cos(psi)
+        a00 =   c * r[0][0] + s * r[1][0];
+        a01 =   c * r[0][1] + s * r[1][1];
+        a02 =   c * r[0][2] + s * r[1][2];
+        a10 = - s * r[0][0] + c * r[1][0];
+        a11 = - s * r[0][1] + c * r[1][1];
+        a12 = - s * r[0][2] + c * r[1][2];
+        r_a << [a00, a01, a02]
+        r_a << [a10, a11, a12]
+        r_a << r[2]
+        return r_a
+      rescue => e
+        raise
+      end
+    end
+
+    #=========================================================================
+    # 座標回転
+    #
+    # @param:  r      (回転行列)
+    # @param:  pos    (回転前直角座標)
+    # @return: pos_r  (回転後直角座標)
+    #=========================================================================
+    def rotate(r, pos)
+      return (0..2).map do |i|
+        (0..2).inject(0) { |sum, j| sum + r[i][j] * pos[j] }
+      end
+    rescue => e
+      raise
+    end
+  end
+end
+

data/lib/mk_coord/trigonometric.rb

@@ -0,0 +1,84 @@
+module MkCoord
+  module Trigonometric
+    module_function
+
+    #=========================================================================
+    # λ の計算
+    #
+    # * λ = arctan((sinδ sinε + cosδ sinα cosε ) / cosδ cosα)
+    #
+    # @param:  alpha   (unit: rad)
+    # @param:  delta   (unit: rad)
+    # @param:  eps     (unit: rad)
+    # @return: lambda  (unit: rad)
+    #=========================================================================
+    def compute_lambda(alpha, delta, eps)
+      a = Math.sin(delta) * Math.sin(eps) \
+        + Math.cos(delta) * Math.sin(alpha) * Math.cos(eps)
+      b = Math.cos(delta) * Math.cos(alpha)
+      lambda = Math.atan2(a, b)
+      lambda %= PI2 if lambda < 0
+      return lambda
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # β の計算
+    #
+    # * β = arcsisn((sinδ cosε - cosδ sinα sinε )
+    #
+    # @param:  alpha   (unit: rad)
+    # @param:  delta   (unit: rad)
+    # @param:  eps     (unit: rad)
+    # @return: lambda  (unit: rad)
+    #=========================================================================
+    def compute_beta(alpha, delta, eps)
+      a = Math.sin(delta) * Math.cos(eps) \
+        - Math.cos(delta) * Math.sin(alpha) * Math.sin(eps)
+      return Math.asin(a)
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # α の計算
+    #
+    # * α = arctan((-sinβ sinε + cosβ sinλ cosε ) / cosβ cosλ)
+    #
+    # @param:  lambda  (unit: rad)
+    # @param:  beta    (unit: rad)
+    # @param:  eps     (unit: rad)
+    # @return: alpha   (unit: rad)
+    #=========================================================================
+    def compute_alpha(lambda, beta, eps)
+      a = -Math.sin(beta) * Math.sin(eps) \
+        +  Math.cos(beta) * Math.sin(lambda) * Math.cos(eps)
+      b =  Math.cos(beta) * Math.cos(lambda)
+      alpha = Math.atan2(a, b)
+      alpha %= PI2 if a < 0
+      return alpha
+    rescue => e
+      raise
+    end
+
+    #=========================================================================
+    # δ の計算
+    #
+    # * δ = arcsisn((sinβ cosε + cosβ sinλ sinε )
+    #
+    # @param:  lambda  (unit: rad)
+    # @param:  beta    (unit: rad)
+    # @param:  eps     (unit: rad)
+    # @return: delta   (unit: rad)
+    #=========================================================================
+    def compute_delta(lambda, beta, eps)
+      a = Math.sin(beta) * Math.cos(eps) \
+        + Math.cos(beta) * Math.sin(lambda) * Math.sin(eps)
+      return Math.asin(a)
+    rescue => e
+      raise
+    end
+  end
+end
+

data/lib/mk_coord/version.rb

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

data/lib/mk_coord.rb

@@ -0,0 +1,59 @@
+require "mk_coord/version"
+require "mk_coord/const"
+require "mk_coord/matrix"
+require "mk_coord/trigonometric"
+include MkCoord::Const
+include MkCoord::Matrix
+include MkCoord::Trigonometric
+
+module MkCoord
+  def self.rect_eq2ec(rect_eq, eps)
+    mtx = r_x(eps)
+    return rotate(mtx, rect_eq)
+  rescue => e
+    raise
+  end
+
+  def self.rect_ec2eq(rect_ec, eps)
+    mtx = r_x(-eps)
+    return rotate(mtx, rect_ec)
+  rescue => e
+    raise
+  end
+
+  def self.pol_eq2ec(pol_eq, eps)
+    lambda = compute_lambda(*pol_eq, eps)
+    beta   = compute_beta(*pol_eq, eps)
+    return [lambda, beta]
+  rescue => e
+    raise
+  end
+
+  def self.pol_ec2eq(pol_ec, eps)
+    alpha = compute_alpha(*pol_ec, eps)
+    delta = compute_delta(*pol_ec, eps)
+    return [alpha, delta]
+  rescue => e
+    raise
+  end
+
+  def self.rect2pol(rect)
+    x, y, z = rect
+    r = Math.sqrt(x * x + y * y)
+    lambda = Math.atan2(y, x)
+    phi    = Math.atan2(z, r)
+    lambda %= PI2 if lambda < 0
+    return [lambda, phi]
+  rescue => e
+    raise
+  end
+
+  def self.pol2rect(pol, r)
+    lambda, phi = pol
+    mtx = r_y(lambda)
+    mtx = r_z(phi, mtx)
+    return rotate(mtx, [r, 0.0, 0.0])
+  rescue => e
+    raise
+  end
+end

data/mk_coord.gemspec

@@ -0,0 +1,33 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'mk_coord/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "mk_coord"
+  spec.version       = MkCoord::VERSION
+  spec.authors       = ["komasaru"]
+  spec.email         = ["masaru@mk-mode.com"]
+
+  spec.summary       = %q{Coordinate conversion library.}
+  spec.description   = %q{This is a coordinate conversion library.}
+  spec.homepage      = "https://github.com/komasaru/mk_coord"
+  spec.license       = "MIT"
+
+  # Prevent pushing this gem to RubyGems.org. To allow pushes either set the 'allowed_push_host'
+  # to allow pushing to a single host or delete this section to allow pushing to any host.
+  #if spec.respond_to?(:metadata)
+  #  spec.metadata['allowed_push_host'] = "TODO: Set to 'http://mygemserver.com'"
+  #else
+  #  raise "RubyGems 2.0 or newer is required to protect against public gem pushes."
+  #end
+
+  spec.files         = `git ls-files -z`.split("\x0").reject { |f| f.match(%r{^(test|spec|features)/}) }
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", "~> 1.12"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "rspec", "~> 3.0"
+end

metadata

@@ -0,0 +1,104 @@
+--- !ruby/object:Gem::Specification
+name: mk_coord
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- komasaru
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2016-07-31 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.12'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.12'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.0'
+description: This is a coordinate conversion library.
+email:
+- masaru@mk-mode.com
+executables:
+- mk_coord
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- ".rspec"
+- ".travis.yml"
+- Gemfile
+- Guardfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- exe/mk_coord
+- lib/mk_coord.rb
+- lib/mk_coord/const.rb
+- lib/mk_coord/matrix.rb
+- lib/mk_coord/trigonometric.rb
+- lib/mk_coord/version.rb
+- mk_coord.gemspec
+homepage: https://github.com/komasaru/mk_coord
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.6.6
+signing_key: 
+specification_version: 4
+summary: Coordinate conversion library.
+test_files: []