geo-triangle 0.0.1

data/.document

@@ -0,0 +1,5 @@
+lib/**/*.rb
+bin/*
+- 
+features/**/*.feature
+LICENSE.txt

data/.rspec

@@ -0,0 +1 @@
+--color

data/Gemfile

@@ -0,0 +1,13 @@
+source "http://rubygems.org"
+# Add dependencies required to use your gem here.
+# Example:
+#   gem "activesupport", ">= 2.3.5"
+
+# Add dependencies to develop your gem here.
+# Include everything needed to run rake, tests, features, etc.
+group :development do
+  gem "rspec", "~> 2.3.0"
+  gem "bundler", "~> 1.0.0"
+  gem "jeweler", "~> 1.6.0"
+  gem "rcov", ">= 0"
+end

data/Gemfile.lock

@@ -0,0 +1,28 @@
+GEM
+  remote: http://rubygems.org/
+  specs:
+    diff-lcs (1.1.2)
+    git (1.2.5)
+    jeweler (1.6.0)
+      bundler (~> 1.0.0)
+      git (>= 1.2.5)
+      rake
+    rake (0.8.7)
+    rcov (0.9.9)
+    rspec (2.3.0)
+      rspec-core (~> 2.3.0)
+      rspec-expectations (~> 2.3.0)
+      rspec-mocks (~> 2.3.0)
+    rspec-core (2.3.1)
+    rspec-expectations (2.3.0)
+      diff-lcs (~> 1.1.2)
+    rspec-mocks (2.3.0)
+
+PLATFORMS
+  ruby
+
+DEPENDENCIES
+  bundler (~> 1.0.0)
+  jeweler (~> 1.6.0)
+  rcov
+  rspec (~> 2.3.0)

data/LICENSE.txt

@@ -0,0 +1,20 @@
+Copyright (c) 2011 parrot-studio
+
+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.rdoc

@@ -0,0 +1,71 @@
+= geo-triangle
+
+== Introduction
+
+Author:: parrot_studio < parrot at users.sourceforge.jp >
+License:: The MIT License
+
+三点の座標から中心（重心）を求める。
+
+あくまで緯度経度を数値的に処理しているだけで、距離の補正等は考慮していない。
+そのため、あまりに離れすぎた三点を指定すると、値が妥当でない場合があり得る。
+
+== Install
+ gem install geo-triangle
+
+== Usage
+ # require 'rubygems' # if ruby 1.8.x 
+ require 'geo-triangle'
+ 
+ # 緯度(lat), 経度(lon)の組み合わせを3つ渡す
+ GeoTriangle.center(
+   [36.383304,139.073632], [36.391026,139.075107], [36.378536,139.046445]) #=> [36.384289, 139.065061]
+ 
+ # 順序が正しければ引数は適当でもいい
+ GeoTriangle.center(
+   36.383304,139.073632,36.391026,139.075107,36.378536,139.046445) 
+ 
+ # newしてから構築するやり方
+ gt1 = GeoTriangle.new
+ gt1.set_coordinate([36.383304,139.073632])
+ gt1.set_coordinate([36.391026,139.075107])
+ gt1.set_coordinate([36.378536,139.046445])
+ gt1.center #=> [36.384289, 139.065061]
+ 
+ # 4個以上は無視される
+ gt2 = GeoTriangle.new
+ gt2.set_coordinate([36.383304,139.073632])
+ gt2.set_coordinate([36.391026,139.075107])
+ gt2.set_coordinate([36.378536,139.046445])
+ gt2.set_coordinate([37.378536,136.046445]) # 4つ目以降は無視
+ gt2.center #=> [36.384289, 139.065061]
+ 
+ # 同じ座標はまとめられる
+ gt3 = GeoTriangle.new
+ gt3.set_coordinate([36.383304,139.073632])
+ gt3.set_coordinate([36.391026,139.075107])
+ gt3.set_coordinate([36.391026,139.075107]) # 上と同じ座標
+ # 座標が2つと見なされ、座標の数が足りないのでnilを返す
+ gt3.center #=> nil
+
+== Extra Module
+ # require 'rubygems' # if 1.8.x 
+ require 'geo-triangle'
+ 
+ # 重心ではなく外心を求めたい場合
+ cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[5,6])
+ cc.center #=> [2.333333, 3.333333]
+ 
+ # 鈍角三角形の場合、外心は三角形の外になるため、実装はしたが使えない
+
+== Release Note
+
+=== ver 0.0.1
+
+* release for test
+
+== Copyright
+
+Copyright (c) 2011 parrot-studio. See LICENSE.txt for
+further details.
+

data/Rakefile

@@ -0,0 +1,50 @@
+# encoding: utf-8
+
+require 'rubygems'
+require 'bundler'
+begin
+  Bundler.setup(:default, :development)
+rescue Bundler::BundlerError => e
+  $stderr.puts e.message
+  $stderr.puts "Run `bundle install` to install missing gems"
+  exit e.status_code
+end
+require 'rake'
+
+require 'jeweler'
+Jeweler::Tasks.new do |gem|
+  # gem is a Gem::Specification... see http://docs.rubygems.org/read/chapter/20 for more options
+  gem.name = "geo-triangle"
+  gem.homepage = "http://github.com/parrot-studio/geo-triangle"
+  gem.license = "MIT"
+  gem.summary = %Q{calculate geo-triangle's gravity center}
+  gem.description = %Q{calculate center of 3 geo-coordinates(lat, lon)}
+  gem.email = "parrot@users.sourceforge.jp"
+  gem.authors = ["parrot-studio"]
+  # dependencies defined in Gemfile
+end
+Jeweler::RubygemsDotOrgTasks.new
+
+require 'rspec/core'
+require 'rspec/core/rake_task'
+RSpec::Core::RakeTask.new(:spec) do |spec|
+  spec.rspec_opts = ["-c", "-r ./spec/spec_helper.rb"]
+  spec.pattern = FileList['spec/**/*_spec.rb']
+end
+
+RSpec::Core::RakeTask.new(:rcov) do |spec|
+  spec.pattern = 'spec/**/*_spec.rb'
+  spec.rcov = true
+end
+
+task :default => :spec
+
+require 'rake/rdoctask'
+Rake::RDocTask.new do |rdoc|
+  version = File.exist?('VERSION') ? File.read('VERSION') : ""
+
+  rdoc.rdoc_dir = 'rdoc'
+  rdoc.title = "geo-triangle #{version}"
+  rdoc.rdoc_files.include('README*')
+  rdoc.rdoc_files.include('lib/**/*.rb')
+end

data/VERSION

@@ -0,0 +1 @@
+0.0.1

data/lib/geo-triangle.rb

@@ -0,0 +1,60 @@
+# coding: utf-8
+require 'bigdecimal'
+require 'gravity_center'
+require 'geo_triangle_ext/linear_function'
+require 'geo_triangle_ext/circum_center'
+
+class GeoTriangle
+
+  class << self
+
+    def center(*val)
+      gt = self.new
+      v = [val].flatten
+      gt.set_coordinate(v[0], v[1])
+      gt.set_coordinate(v[2], v[3])
+      gt.set_coordinate(v[4], v[5])
+      gt.center
+    end
+
+  end
+
+  def exist?(lat, lon)
+    return false unless lat
+    return false unless lon
+    return false if lat.abs > 90.0
+    return false if lon.abs > 180.0
+    true
+  end
+
+  def set_coordinate(*val)
+    v = [val].flatten
+    lat = v[0]
+    lon = v[1]
+    coordinates << [lat, lon] if exist?(lat, lon)
+    self
+  end
+
+  def coordinates
+    @coordinates ||= []
+    @coordinates = @coordinates.uniq
+    @coordinates = @coordinates.take(3) if @coordinates.size > 3
+    @coordinates
+  end
+
+  def valid_coordinates?
+    return false unless coordinates.size == 3
+    return false unless coordinates.all?{|c| c.size == 2}
+    true
+  end
+
+  def center
+    return unless valid_coordinates?
+    gc = GravityCenter.create(*coordinates)
+    lat = gc.center_x.round(6).to_f
+    lon = gc.center_y.round(6).to_f
+    exist?(lat, lon) ? [lat, lon] : nil
+  end
+
+end
+

data/lib/geo_triangle_ext/circum_center.rb

@@ -0,0 +1,84 @@
+# coding: utf-8
+module GeoTriangleExt
+
+  class CircumCenter
+
+    class << self
+
+      def create(*co)
+        cc = self.new
+        co.take(3).each{|c| cc.coordinates << c}
+        cc
+      end
+
+    end
+
+    def set_coordinate(x, y)
+      coordinates << [x, y]
+    end
+
+    def coordinates
+      @coordinates ||= []
+      @coordinates = @coordinates.uniq
+      @coordinates = @coordinates.take(3) if @coordinates.size > 3
+      @coordinates
+    end
+
+    def valid_coordinates?
+      return false unless coordinates.size == 3
+      return false unless coordinates.all?{|c| c.size == 2}
+      true
+    end
+
+    def functions
+      @functions ||= create_functions
+      @functions
+    end
+
+    def valid_functions?
+      return false if functions.size < 2
+      return false if functions.combination(2).any?{|f1, f2|
+        f1.orthogonal_slope == f2.orthogonal_slope}
+      true
+    end
+
+    def center
+      return unless valid_functions?
+      @center_point ||= calc_center
+      @center_point
+    end
+
+    private
+
+    def create_functions
+      return [] unless valid_coordinates?
+      ret = []
+      coordinates.combination(2) do |a, b|
+        f = LinearFunction.create(*([a, b].flatten))
+        next unless f
+        ret << f if f.valid?
+      end
+      ret
+    end
+
+    def calc_center
+      return unless valid_functions?
+
+      l1 = functions.first
+      l2 = functions.last
+
+      a1 = l1.orthogonal_slope
+      b1 = l1.orthogonal_intercept
+      a2 = l2.orthogonal_slope
+      b2 = l2.orthogonal_intercept
+
+      x = (b2 - b1) / (a1 - a2)
+      y = (a1 * b2 - b1 * a2) / (a1 - a2)
+
+      [x, y]
+    end
+
+  end
+
+end
+

data/lib/geo_triangle_ext/linear_function.rb

@@ -0,0 +1,77 @@
+# coding: utf-8
+module GeoTriangleExt
+
+  class LinearFunction
+
+    attr_writer :ax, :ay, :bx, :by
+
+    class << self
+
+      def create(ax, ay, bx, by)
+        lf = self.new
+        lf.ax = ax
+        lf.ay = ay
+        lf.bx = bx
+        lf.by = by
+        lf
+      end
+
+    end
+
+    def ax
+      @bax ||= BigDecimal(@ax.to_s)
+      @bax
+
+    end
+
+    def ay
+      @bay ||= BigDecimal(@ay.to_s)
+      @bay
+    end
+
+    def bx
+      @bbx ||= BigDecimal(@bx.to_s)
+      @bbx
+    end
+
+    def by
+      @bby ||= BigDecimal(@by.to_s)
+      @bby
+    end
+
+    def valid?
+      return false if ax == bx
+      return false if ay == by
+      true
+    end
+
+    def middle_point
+      x = (ax + bx) / 2.0
+      y = (ay + by) / 2.0
+      [x, y]
+    end
+
+    def slope
+      return unless valid?
+      (ay - by) / (ax - bx)
+    end
+
+    def intercept
+      return unless valid?
+      (ax * by - ay * bx) / (ax - bx)
+    end
+
+    def orthogonal_slope
+      return unless valid?
+      (bx - ax) / (ay - by)
+    end
+
+    def orthogonal_intercept
+      return unless valid?
+      (ax*ax + ay*ay - bx*bx - by*by) / (2 * (ay - by))
+    end
+
+  end
+
+end
+

data/lib/gravity_center.rb

@@ -0,0 +1,75 @@
+# coding: utf-8
+require 'bigdecimal'
+
+class GravityCenter
+
+  attr_reader :x1, :y1, :x2, :y2, :x3, :y3
+
+  def self.create(*val)
+    gc = self.new
+    gc.instance_eval do
+      v = [val].flatten
+      @x1 = BigDecimal(v[0].to_s)
+      @y1 = BigDecimal(v[1].to_s)
+      @x2 = BigDecimal(v[2].to_s)
+      @y2 = BigDecimal(v[3].to_s)
+      @x3 = BigDecimal(v[4].to_s)
+      @y3 = BigDecimal(v[5].to_s)
+    end
+    gc
+  end
+
+  def sx1
+    @sx1 ||= x1 - x3
+    @sx1
+  end
+
+  def sy1
+    @sy1 ||= y1 - y3
+    @sy1
+  end
+
+  def sx2
+    @sx2 ||= x2 - x3
+    @sx2
+  end
+
+  def sy2
+    @sy2 ||= y2 - y3
+    @sy2
+  end
+
+  def center_sx
+    @center_sx ||= calc_center_sx
+    @center_sx
+  end
+
+  def center_x
+    center_sx + x3
+  end
+
+  def center_sy
+    @center_sy ||= calc_center_sy
+    @center_sy
+  end
+
+  def center_y
+    center_sy + y3
+  end
+
+  def center
+    [center_x, center_y]
+  end
+
+  private
+
+  def calc_center_sx
+    (sx1 + sx2) / 3
+  end
+
+  def calc_center_sy
+    (sy1 + sy2) / 3
+  end
+
+end
+

data/spec/geo-triangle_spec.rb

@@ -0,0 +1,130 @@
+# coding: utf-8
+describe GeoTriangle do
+
+  describe "#exist?" do
+
+    before do
+      @gt = GeoTriangle.new
+    end
+
+    context "lat or lon is nil" do
+      it {@gt.exist?(40.2, nil).should be_false}
+      it {@gt.exist?(nil, 135.0).should be_false}
+    end
+
+    context "over value" do
+      it {@gt.exist?(90.1, 135.0).should be_false}
+      it {@gt.exist?(42.0, 180.1).should be_false}
+    end
+
+    context "valid lat and lon" do
+      it {@gt.exist?(36.383304, 139.073632).should be_true}
+    end
+
+  end
+
+  describe "#set_coordinate" do
+
+    before do
+      @gt = GeoTriangle.new
+    end
+
+    context "exist coordinate" do
+      before do
+        @gt.set_coordinate([36.383304, 139.073632])
+        @gt.set_coordinate(36.391026, 139.075107)
+      end
+      it {@gt.coordinates.should == [[36.383304, 139.073632], [36.391026, 139.075107]]}
+    end
+
+    context "not exist coordinate" do
+      before do
+        @gt.set_coordinate([136.383304, 139.073632])
+        @gt.set_coordinate(36.391026, 239.075107)
+      end
+      it {@gt.coordinates.should be_empty}
+    end
+
+    context "not coordinate pair" do
+      before do
+        @gt.set_coordinate(36.383304)
+      end
+      it {@gt.coordinates.should be_empty}
+    end
+
+  end
+
+  describe "#valid_coordinates?" do
+
+    before do
+      @gt = GeoTriangle.new
+    end
+
+    context "coordinates size not enough" do
+      before do
+        @gt.set_coordinate([36.383304, 139.073632])
+        @gt.set_coordinate(36.391026, 139.075107)
+      end
+      it {@gt.coordinates.size.should == 2}
+      it {@gt.valid_coordinates?.should be_false}
+    end
+
+    context "invalid coordinate format" do
+      before do
+        @gt.set_coordinate([36.383304, 139.073632])
+        @gt.set_coordinate(36.391026, 139.075107)
+        @gt.coordinates << [135.0]
+      end
+      it {@gt.coordinates.size.should == 3}
+      it {@gt.valid_coordinates?.should be_false}
+    end
+
+    context "valid coordinates" do
+      before do
+        @gt.set_coordinate([36.383304,139.073632])
+        @gt.set_coordinate([36.391026,139.075107])
+        @gt.set_coordinate([36.378536,139.046445])
+        @gt.set_coordinate([37.378536,138.046445])
+      end
+      it {@gt.coordinates.size.should == 3}
+      it {@gt.valid_coordinates?.should be_true}
+    end
+
+  end
+
+  describe "#center" do
+
+    before do
+      @gt = GeoTriangle.new
+    end
+
+    context "valid coordinates center" do
+      before do
+        @gt.set_coordinate([36.383304,139.073632])
+        @gt.set_coordinate([36.391026,139.075107])
+        @gt.set_coordinate([36.378536,139.046445])
+      end
+      it {@gt.center.should == [36.384289, 139.065061]}
+    end
+
+    context "invalid coordinates" do
+      before do
+        @gt.set_coordinate([36.383304,139.073632])
+        @gt.set_coordinate([36.391026,139.075107])
+        @gt.set_coordinate([36.391026,139.075107])
+      end
+      it {@gt.coordinates.size.should == 2}
+      it {@gt.valid_coordinates?.should be_false}
+      it {@gt.center.should be_nil}
+    end
+
+  end
+
+  describe ".center" do
+    subject {GeoTriangle.center(
+      [36.383304,139.073632], [36.391026,139.075107], [36.378536,139.046445])}
+    it {should == [36.384289, 139.065061]}
+  end
+
+end
+

data/spec/geo_triangle_ext/circum_center_spec.rb

@@ -0,0 +1,178 @@
+# coding: utf-8
+describe GeoTriangleExt::CircumCenter do
+
+  describe "#set_coordinate" do
+    before do
+      @cc = GeoTriangleExt::CircumCenter.new
+      @cc.set_coordinate(1, 2)
+    end
+    it {@cc.coordinates.should == [[1,2]]}
+  end
+
+  describe "#coordinates" do
+
+    context "default" do
+      it {GeoTriangleExt::CircumCenter.new.coordinates.should be_empty}
+    end
+
+    context "has only 3 coordinates" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.new
+        @cc.set_coordinate(1, 2)
+        @cc.set_coordinate(3, 4)
+        @cc.set_coordinate(5, 6)
+        @cc.set_coordinate(7, 8)
+      end
+      it {@cc.coordinates.size.should == 3}
+      it {@cc.coordinates.should == [[1,2],[3,4],[5,6]]}
+    end
+
+  end
+
+  describe "#valid_coordinates?" do
+
+    context "coordinates size not enough" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.new
+        @cc.set_coordinate(1, 2)
+        @cc.set_coordinate(3, 4)
+      end
+      it {@cc.coordinates.size.should == 2}
+      it {@cc.valid_coordinates?.should be_false}
+    end
+
+    context "invalid coordinate format" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.new
+        @cc.set_coordinate(1, 2)
+        @cc.set_coordinate(3, 4)
+        @cc.coordinates << 5
+      end
+      it {@cc.coordinates.size.should == 3}
+      it {@cc.valid_coordinates?.should be_false}
+    end
+
+    context "valid coordinates" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.new
+        @cc.set_coordinate(1, 2)
+        @cc.set_coordinate(3, 4)
+        @cc.set_coordinate(5, 6)
+        @cc.set_coordinate(7, 8)
+      end
+      it {@cc.coordinates.size.should == 3}
+      it {@cc.valid_coordinates?.should be_true}
+    end
+
+  end
+
+  describe ".create" do
+
+    context "create object" do
+      subject {GeoTriangleExt::CircumCenter.create([1,2], [3,4], [5,6])}
+      it {should be_instance_of(GeoTriangleExt::CircumCenter)}
+      its(:coordinates) {should == [[1,2],[3,4],[5,6]]}
+    end
+
+    context "take 3 coordinates if given over 3" do
+      subject {GeoTriangleExt::CircumCenter.create([1,2], [3,4], [5,6],[7,8])}
+      it {should be_instance_of(GeoTriangleExt::CircumCenter)}
+      its(:coordinates) {should == [[1,2],[3,4],[5,6]]}
+    end
+
+  end
+
+  describe "#functions" do
+
+    context "3 functions from triangle coordinates" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[5,6])
+      end
+      it {@cc.functions.size.should == 3}
+    end
+
+    context "not created if not valid coordinates" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.new
+        @cc.set_coordinate(1, 2)
+        @cc.set_coordinate(3, 4)
+      end
+      it {@cc.functions.should be_empty}
+    end
+
+    context "2 functions if 2 coordinates x or y is same" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[5,4])
+      end
+      it {@cc.functions.size.should == 2}
+    end
+
+    context "1 functions if right triangle" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[0,0])
+      end
+      it {@cc.functions.size.should == 1}
+    end
+
+  end
+
+  describe "#valid_functions?" do
+
+    context "valid triangle coordinates" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[5,6])
+      end
+      it {@cc.functions.size.should == 3}
+      it {@cc.valid_functions?.should be_true}
+    end
+
+    context "2 functions is valid" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[5,4])
+      end
+      it {@cc.functions.size.should == 2}
+      it {@cc.valid_functions?.should be_true}
+    end
+
+    context "1 functions is not valid" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[0,0])
+      end
+      it {@cc.functions.size.should == 1}
+      it {@cc.valid_functions?.should be_false}
+    end
+
+    context "3 coordinates on same line" do
+      before do
+        @cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[4,-4])
+      end
+      it {@cc.functions.size.should == 3}
+      it {@cc.valid_functions?.should be_false}
+    end
+
+  end
+
+  describe "#center" do
+
+    context "triangle coordinates" do
+      before do
+        cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[5,6])
+        @x, @y = cc.center
+      end
+      it {@x.round(5).to_f.should == 3.25}
+      it {@y.round(5).to_f.should == 3.125}
+    end
+
+    context "2 coordinates x or y is same" do
+      before do
+        cc = GeoTriangleExt::CircumCenter.create([0,4],[2,0],[5,4])
+        @x, @y = cc.center
+      end
+      it {@x.round(5).to_f.should == 2.5}
+      it {@y.round(5).to_f.should == 2.75}
+    end
+
+  end
+
+end
+

data/spec/geo_triangle_ext/linear_function_spec.rb

@@ -0,0 +1,97 @@
+# coding: utf-8
+describe GeoTriangleExt::LinearFunction do
+
+  describe "#valid?" do
+
+    context "ax == bx" do
+      before do
+        @lf = GeoTriangleExt::LinearFunction.new
+        @lf.ax = 1
+        @lf.ay = 2
+        @lf.bx = 1
+        @lf.by = 5
+      end
+      it {@lf.valid?.should be_false}
+    end
+
+    context "bx == by" do
+      before do
+        @lf = GeoTriangleExt::LinearFunction.new
+        @lf.ax = 1
+        @lf.ay = 2
+        @lf.bx = 7
+        @lf.by = 2
+      end
+      it {@lf.valid?.should be_false}
+    end
+
+    context "unrelated" do
+      before do
+        @lf = GeoTriangleExt::LinearFunction.new
+        @lf.ax = 1
+        @lf.ay = 2
+        @lf.bx = 7
+        @lf.by = 0
+      end
+      it {@lf.valid?.should be_true}
+    end
+
+  end
+
+  describe ".create" do
+    subject{GeoTriangleExt::LinearFunction.create(1, 2, 7, 0)}
+    it {should be_instance_of(GeoTriangleExt::LinearFunction)}
+    its(:ax) {should == 1.0}
+    its(:ay) {should == 2.0}
+    its(:bx) {should == 7.0}
+    its(:by) {should == 0.0}
+  end
+
+  describe "#middle_point" do
+
+    context "positive number only" do
+      before do
+        @lf = GeoTriangleExt::LinearFunction.create(1, 2, 7, 1)
+      end
+      it {@lf.middle_point.should == [4.0, 1.5]}
+    end
+
+    context "include negative number" do
+      before do
+        @lf = GeoTriangleExt::LinearFunction.create(1, 2, -7, -1)
+      end
+      it {@lf.middle_point.should == [-3.0, 0.5]}
+    end
+
+  end
+
+  describe "#slope" do
+    before do
+      @lf = GeoTriangleExt::LinearFunction.create(0, 4, 2, 0)
+    end
+    it {@lf.slope.should == -2.0}
+  end
+
+  describe "#intercept" do
+    before do
+      @lf = GeoTriangleExt::LinearFunction.create(0, 4, 2, 0)
+    end
+    it {@lf.intercept.should == 4.0}
+  end
+
+  describe "#orthogonal_slope" do
+    before do
+      @lf = GeoTriangleExt::LinearFunction.create(0, 4, 2, 0)
+    end
+    it {@lf.orthogonal_slope.should == 0.5}
+  end
+
+  describe "#orthogonal_intercept" do
+    before do
+      @lf = GeoTriangleExt::LinearFunction.create(0, 4, 2, 0)
+    end
+    it {@lf.orthogonal_intercept.should == 1.5}
+  end
+
+end
+

data/spec/gravity_center_spec.rb

@@ -0,0 +1,60 @@
+# coding: utf-8
+describe GravityCenter do
+
+  describe ".create" do
+
+    context "3 coordinates set" do
+      subject {GravityCenter.create([1,2],[3,4],[5,6])}
+      it {should be_instance_of(GravityCenter)}
+      its(:x1) {should == 1.0}
+      its(:y1) {should == 2.0}
+      its(:x2) {should == 3.0}
+      its(:y2) {should == 4.0}
+      its(:x3) {should == 5.0}
+      its(:y3) {should == 6.0}
+    end
+
+    context "6 numbers set" do
+      subject {GravityCenter.create(10,20,30,40,50,60)}
+      it {should be_instance_of(GravityCenter)}
+      its(:x1) {should == 10.0}
+      its(:y1) {should == 20.0}
+      its(:x2) {should == 30.0}
+      its(:y2) {should == 40.0}
+      its(:x3) {should == 50.0}
+      its(:y3) {should == 60.0}
+    end
+
+  end
+
+  describe "#sx1, #sx2, #sy1, #sy2 is shifted coordinates for (x3,y3) to (0,0)" do
+    subject {GravityCenter.create([0,4],[2,0],[5,6])}
+    its(:sx1) {should == -5.0}
+    its(:sy1) {should == -2.0}
+    its(:sx2) {should == -3.0}
+    its(:sy2) {should == -6.0}
+  end
+
+  describe "#center_sx, #center_sy is gravity center of shifted triangle" do
+    subject {GravityCenter.create([0,4],[2,0],[5,6])}
+    it {subject.center_sx.round(6).to_f.should == -2.666667}
+    it {subject.center_sy.round(6).to_f.should == -2.666667}
+  end
+
+  describe "#center_x, #center_y is gravity center of original triangle" do
+    subject {GravityCenter.create([0,4],[2,0],[5,6])}
+    it {subject.center_x.round(6).to_f.should == 2.333333}
+    it {subject.center_y.round(6).to_f.should == 3.333333}
+  end
+
+  describe "#center is gravity center coordinate of triangle" do
+    before do
+      gc = GravityCenter.create([0,4],[2,0],[5,6])
+      @x, @y = gc.center
+    end
+    it {@x.round(6).to_f.should == 2.333333}
+    it {@y.round(6).to_f.should == 3.333333}
+  end
+
+end
+

data/spec/spec_helper.rb

@@ -0,0 +1,12 @@
+$LOAD_PATH.unshift(File.join(File.dirname(__FILE__), '..', 'lib'))
+$LOAD_PATH.unshift(File.dirname(__FILE__))
+require 'rspec'
+require 'geo-triangle'
+
+# Requires supporting files with custom matchers and macros, etc,
+# in ./support/ and its subdirectories.
+Dir["#{File.dirname(__FILE__)}/support/**/*.rb"].each {|f| require f}
+
+RSpec.configure do |config|
+  
+end

metadata

@@ -0,0 +1,119 @@
+--- !ruby/object:Gem::Specification 
+name: geo-triangle
+version: !ruby/object:Gem::Version 
+  prerelease: 
+  version: 0.0.1
+platform: ruby
+authors: 
+- parrot-studio
+autorequire: 
+bindir: bin
+cert_chain: []
+
+date: 2011-05-19 00:00:00 +09:00
+default_executable: 
+dependencies: 
+- !ruby/object:Gem::Dependency 
+  name: rspec
+  requirement: &id001 !ruby/object:Gem::Requirement 
+    none: false
+    requirements: 
+    - - ~>
+      - !ruby/object:Gem::Version 
+        version: 2.3.0
+  type: :development
+  prerelease: false
+  version_requirements: *id001
+- !ruby/object:Gem::Dependency 
+  name: bundler
+  requirement: &id002 !ruby/object:Gem::Requirement 
+    none: false
+    requirements: 
+    - - ~>
+      - !ruby/object:Gem::Version 
+        version: 1.0.0
+  type: :development
+  prerelease: false
+  version_requirements: *id002
+- !ruby/object:Gem::Dependency 
+  name: jeweler
+  requirement: &id003 !ruby/object:Gem::Requirement 
+    none: false
+    requirements: 
+    - - ~>
+      - !ruby/object:Gem::Version 
+        version: 1.6.0
+  type: :development
+  prerelease: false
+  version_requirements: *id003
+- !ruby/object:Gem::Dependency 
+  name: rcov
+  requirement: &id004 !ruby/object:Gem::Requirement 
+    none: false
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        version: "0"
+  type: :development
+  prerelease: false
+  version_requirements: *id004
+description: calculate center of 3 geo-coordinates(lat, lon)
+email: parrot@users.sourceforge.jp
+executables: []
+
+extensions: []
+
+extra_rdoc_files: 
+- LICENSE.txt
+- README.rdoc
+files: 
+- .document
+- .rspec
+- Gemfile
+- Gemfile.lock
+- LICENSE.txt
+- README.rdoc
+- Rakefile
+- VERSION
+- lib/geo-triangle.rb
+- lib/geo_triangle_ext/circum_center.rb
+- lib/geo_triangle_ext/linear_function.rb
+- lib/gravity_center.rb
+- spec/geo-triangle_spec.rb
+- spec/geo_triangle_ext/circum_center_spec.rb
+- spec/geo_triangle_ext/linear_function_spec.rb
+- spec/gravity_center_spec.rb
+- spec/spec_helper.rb
+has_rdoc: true
+homepage: http://github.com/parrot-studio/geo-triangle
+licenses: 
+- MIT
+post_install_message: 
+rdoc_options: []
+
+require_paths: 
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement 
+  none: false
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      hash: 449437378648933997
+      segments: 
+      - 0
+      version: "0"
+required_rubygems_version: !ruby/object:Gem::Requirement 
+  none: false
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      version: "0"
+requirements: []
+
+rubyforge_project: 
+rubygems_version: 1.6.2
+signing_key: 
+specification_version: 3
+summary: calculate geo-triangle's gravity center
+test_files: []
+