mars_geo 0.0.1

data/.gitignore

@@ -0,0 +1,17 @@
+*.gem
+*.rbc
+.bundle
+.config
+.yardoc
+Gemfile.lock
+InstalledFiles
+_yardoc
+coverage
+doc/
+lib/bundler/man
+pkg
+rdoc
+spec/reports
+test/tmp
+test/version_tmp
+tmp

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in mars_geo.gemspec
+gemspec

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2013 jonny
+
+MIT License
+
+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,51 @@
+# MarsGeo
+
+这是一个将真实经纬度转换成火星地址的Gem. 我们在做手机App开发时经常会用到地图，通常手机获得的经纬度放到
+Google地图的时候会有几百米的偏差，这是因为某些原因国内地图都做了偏移加密，这个俗称火星坐标系。MarsGeo是
+根据网上的一段Java代码翻译过来的，目的是把真实经纬度转换成火星坐标系好和Google地图匹配，至于准确性如何
+我没有全面测试过，不过我自己使用的情况还是可以滴。
+
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+    gem 'mars_geo'
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install mars_geo
+
+## Usage
+
+在要使用的文件中 include MarsGeo Moudle
+
+```
+class HomeController < ApplicationController
+    include MarsGeo
+
+    def index
+        point = Point.new(31.108949, 121.333383)
+        offset_point = point.offset
+        p offset_point.lat 
+        p offset_point.lng 
+    end
+
+end
+
+```
+
+
+
+
+## Contributing
+
+1. Fork it
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Add some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create new Pull Request

data/Rakefile

@@ -0,0 +1 @@
+require "bundler/gem_tasks"

data/lib/mars_geo/converter.rb

@@ -0,0 +1,167 @@
+#encoding: utf-8
+
+module MarsGeo
+  class Converter
+    attr_accessor :casm_rr,:casm_t1,:casm_t2,:casm_x1,:casm_y1, :casm_x2,:casm_y2,:casm_f
+
+    def initialize
+      self.casm_rr = 0
+    end
+
+    def yj_sin2(x)
+      ff = 0
+
+      if x<0
+        x = -x
+        ff = 1
+      end
+
+      cc = (x / 6.28318530717959).to_i
+      tt = x - cc * 6.28318530717959
+
+      if tt > 3.1415926535897932 
+        tt = tt - 3.1415926535897932
+        if ff == 1
+          ff = 0
+        elsif ff == 0
+          ff = 1
+        end
+      end
+      x = tt
+      ss = x
+      s2 = x
+      tt = tt * tt
+      s2 = s2 * tt
+      ss = ss - s2 * 0.166666666666667
+      s2 = s2 * tt
+      ss = ss + s2 * 8.33333333333333E-03
+      s2 = s2 * tt
+      ss = ss - s2 * 1.98412698412698E-04
+      s2 = s2 * tt
+      ss = ss + s2 * 2.75573192239859E-06
+      s2 = s2 * tt
+      ss = ss - s2 * 2.50521083854417E-08
+      ss = -ss if (ff == 1)
+
+      return ss;
+    end
+
+
+    def transform_yj5(x, y)
+      tt = 300 + 1 * x + 2 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.sqrt(x * x))
+      tt = tt + (20 * yj_sin2(18.849555921538764 * x) + 20 * yj_sin2(6.283185307179588 * x)) * 0.6667
+      tt = tt + (20 * yj_sin2(3.141592653589794 * x) + 40 * yj_sin2(1.047197551196598 * x)) * 0.6667
+      tt = tt + (150 * yj_sin2(0.2617993877991495 * x) + 300 * yj_sin2(0.1047197551196598 * x)) * 0.6667
+      tt
+    end
+
+    def transform_yjy5(x,y)
+      tt = -100 + 2 * x + 3 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.sqrt(x * x))
+      tt = tt + (20 * yj_sin2(18.849555921538764 * x) + 20 * yj_sin2(6.283185307179588 * x)) * 0.6667
+      tt = tt + (20 * yj_sin2(3.141592653589794 * y) + 40 * yj_sin2(1.047197551196598 * y)) * 0.6667
+      tt = tt + (160 * yj_sin2(0.2617993877991495 * y) + 320 * yj_sin2(0.1047197551196598 * y)) * 0.6667
+    end
+
+    def transform_jy5(x, xx)
+      a = 6378245
+      e = 0.00669342
+      n = Math.sqrt(1 - e * yj_sin2(x * 0.0174532925199433) * yj_sin2(x * 0.0174532925199433))
+      n = (xx * 180) / (a / n * Math.cos(x * 0.0174532925199433) * 3.1415926)
+    end
+
+    def transform_jyj5(x,yy)
+      a = 6378245
+      e = 0.00669342
+      mm = 1 - e * yj_sin2(x * 0.0174532925199433) * yj_sin2(x * 0.0174532925199433)
+      m = (a * (1 - e)) / (mm * Math.sqrt(mm))
+      return (yy * 180) / (m * 3.1415926)
+    end
+
+    def random_yj
+      casm_a = 314159269
+      casm_c = 453806245
+      self.casm_rr = casm_a * self.casm_rr + casm_c
+      t =  (self.casm_rr / 2).to_i
+      self.casm_rr = self.casm_rr - t * 2
+      self.casm_rr = self.casm_rr.to_f / 2
+    end
+
+    def ini_casm(w_time, w_lng,  w_lat)
+      self.casm_t1 = w_time
+      self.casm_t2 = w_time
+      tt = (w_time / 0.357).to_i
+      self.casm_rr = w_time - tt * 0.357
+      self.casm_rr = 0.3 if (w_time == 0)
+      self.casm_x1 = w_lng
+      self.casm_y1 = w_lat
+      self.casm_x2 = w_lng
+      self.casm_y2 = w_lat
+      self.casm_f = 3
+    end
+
+
+
+    def wgtochina_lb(wg_flag,wg_lat,wg_lng,wg_heit,wg_week,wg_time)
+
+      point = {lat:0.0,lng:0.0}
+
+      return nil if (wg_heit > 5000) 
+      x_l = wg_lng;
+      x_l = x_l / 3686400.0;
+      y_l = wg_lat;
+      y_l = y_l / 3686400.0;
+
+      return nil if (x_l < 72.004)
+      return nil if (x_l > 137.8347)
+      return nil if (y_l < 0.8293)
+      return nil if (y_l > 55.8271)
+
+      if (wg_flag == 0) 
+        ini_casm(wg_time, wg_lng, wg_lat)
+        point[:lat] = wg_lng
+        point[:lng] = wg_lat
+        return point;
+      end
+
+      casm_t2 = wg_time
+      casm_t1,casm_f = 0,0
+      t1_t2 = (casm_t2 - casm_t1) / 1000.0
+      if (t1_t2 <= 0) 
+        casm_t1 = casm_t2
+        casm_f = casm_f + 1
+        casm_x1 = casm_x2
+        casm_f = casm_f + 1
+        casm_y1 = casm_y2
+        casm_f = casm_f + 1
+      else 
+        if t1_t2 > 120
+            if casm_f == 3
+                casm_f = 0
+                casm_x2 = wg_lng
+                casm_y2 = wg_lat
+                x1_x2 = casm_x2 - casm_x1
+                y1_y2 = casm_y2 - casm_y1
+                casm_v = Math.sqrt(x1_x2 * x1_x2 + y1_y2 * y1_y2) / t1_t2
+                return point if (casm_v > 3185) 
+            end
+            casm_t1 = casm_t2
+            casm_f = casm_f + 1
+            casm_x1 = casm_x2
+            casm_f = casm_f + 1
+            casm_y1 = casm_y2
+            casm_f = casm_f + 1
+        end
+      end
+      x_add = transform_yj5(x_l - 105, y_l - 35)
+      y_add = transform_yjy5(x_l - 105, y_l - 35)
+      h_add = wg_heit
+
+
+      x_add = x_add + h_add * 0.001 + yj_sin2(wg_time * 0.0174532925199433) + random_yj()
+      y_add = y_add + h_add * 0.001 + yj_sin2(wg_time * 0.0174532925199433) + random_yj()
+      point[:lng] =  ((x_l + transform_jy5(y_l, x_add)) * 3686400).to_i
+      point[:lat] = ((y_l + transform_jyj5(y_l, y_add)) * 3686400).to_i
+      return point
+    end
+  end
+end

data/lib/mars_geo/point.rb

@@ -0,0 +1,22 @@
+module MarsGeo
+  class Point
+    attr_accessor :lat,:lng
+
+
+    def initialize(lat,lng)
+      self.lat, self.lng = lat, lng
+    end
+
+
+
+    def offset
+        temp_lat = self.lat * 3686400
+        temp_lng = self.lng * 3686400
+        converter = Converter.new
+        hash_point = converter.wgtochina_lb(1,temp_lat.to_i,temp_lng.to_i,0,0,0)
+        temp_lat = hash_point[:lat] / 3686400.0
+        temp_lng = hash_point[:lng] / 3686400.0
+        Point.new(temp_lat,temp_lng)
+    end
+  end
+end

data/lib/mars_geo/version.rb

@@ -0,0 +1,3 @@
+module MarsGeo
+  VERSION = "0.0.1"
+end

data/lib/mars_geo.rb

@@ -0,0 +1,6 @@
+require "mars_geo/version"
+
+module MarsGeo
+  autoload :Converter, 'mars_geo/converter'
+  autoload :Point, 'mars_geo/point'
+end

data/mars_geo.gemspec

@@ -0,0 +1,19 @@
+# -*- encoding: utf-8 -*-
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'mars_geo/version'
+
+Gem::Specification.new do |gem|
+  gem.name          = "mars_geo"
+  gem.version       = MarsGeo::VERSION
+  gem.authors       = ["Jonny"]
+  gem.email         = ["mars131@gmail.com"]
+  gem.description   = %q{Mars geograpic lib}
+  gem.summary       = %q{Implementation of china geograpic encrypt which offset the position for display correctly in chinese map service }
+  gem.homepage      = "https://github.com/jonnyzheng/mars_geo"
+
+  gem.files         = `git ls-files`.split($/)
+  gem.executables   = gem.files.grep(%r{^bin/}).map{ |f| File.basename(f) }
+  gem.test_files    = gem.files.grep(%r{^(test|spec|features)/})
+  gem.require_paths = ["lib"]
+end

data/spec/converter_spec.rb

@@ -0,0 +1,46 @@
+require 'mars_geo'
+
+include MarsGeo
+
+describe MarsGeo::Converter do
+
+  it 'method yj_sin2 should correct' do 
+    converter = Converter.new
+    value = converter.yj_sin2(31.5566)
+    value.should eq(0.14020995787539087)
+  end
+
+  it 'method transform_yj5 should correct' do 
+    converter = Converter.new
+    value = converter.transform_yj5(31.5566,121.77887615)
+    value.should eq(1116.0853927686755)
+  end
+
+  it 'method transform_yjy5 should correct' do 
+    converter = Converter.new
+    value = converter.transform_yjy5(31.5566,121.77887615)
+    value.should eq(3767.9839409629326)
+  end
+
+
+  it 'method transform_jy5 should correct' do
+    converter = Converter.new
+    value = converter.transform_jy5(31.5566,121.77887615)
+    value.should eq(0.0012826033424097633)
+  end
+
+  it 'method transform_jyj5 should correct' do
+    converter = Converter.new
+    value = converter.transform_jyj5(31.5566,121.77887615)
+    value.should eq(0.0010982842446974356)
+  end
+
+  it 'method wgtochina_lb should correct' do
+      point = Point.new(31.208949, 121.533383)
+      offset_point = point.offset
+      offset_point.lat.should eq(31.2068310546875)
+      offset_point.lng.should eq(121.53770209418403)
+  end
+
+
+end

metadata

@@ -0,0 +1,59 @@
+--- !ruby/object:Gem::Specification
+name: mars_geo
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+  prerelease: 
+platform: ruby
+authors:
+- Jonny
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2013-06-08 00:00:00.000000000 Z
+dependencies: []
+description: Mars geograpic lib
+email:
+- mars131@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- .gitignore
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- lib/mars_geo.rb
+- lib/mars_geo/converter.rb
+- lib/mars_geo/point.rb
+- lib/mars_geo/version.rb
+- mars_geo.gemspec
+- spec/converter_spec.rb
+homepage: https://github.com/jonnyzheng/mars_geo
+licenses: []
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 1.8.24
+signing_key: 
+specification_version: 3
+summary: Implementation of china geograpic encrypt which offset the position for display
+  correctly in chinese map service
+test_files:
+- spec/converter_spec.rb
+has_rdoc: