japan_plane_rectangular 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 857c201c090e63bc80b3ace5062fd573568fed57
+  data.tar.gz: 804f6dbfe9b1f445699fae03c35842e13513d51b
+SHA512:
+  metadata.gz: c719c3d4a18995056c415cedce862f048348fb9707a5562fb853c8cfcd883a172632fd9cae8083d840b8102aad5465e05627cad03bd0e81e34ff1bad62cb5301
+  data.tar.gz: 2bd05010662da314f90e062f2a9a7b1e574536eb17ca87a7f363a3c965cbe1d0b2a44321567b035c70cb1b899e900d52a4e01d0c34aeb6a3bbd2b115911fb784

data/.gitignore

@@ -0,0 +1,50 @@
+*.gem
+*.rbc
+/.config
+/coverage/
+/InstalledFiles
+/pkg/
+/spec/reports/
+/spec/examples.txt
+/test/tmp/
+/test/version_tmp/
+/tmp/
+
+# Used by dotenv library to load environment variables.
+# .env
+
+## Specific to RubyMotion:
+.dat*
+.repl_history
+build/
+*.bridgesupport
+build-iPhoneOS/
+build-iPhoneSimulator/
+
+## Specific to RubyMotion (use of CocoaPods):
+#
+# We recommend against adding the Pods directory to your .gitignore. However
+# you should judge for yourself, the pros and cons are mentioned at:
+# https://guides.cocoapods.org/using/using-cocoapods.html#should-i-check-the-pods-directory-into-source-control
+#
+# vendor/Pods/
+
+## Documentation cache and generated files:
+/.yardoc/
+/_yardoc/
+/doc/
+/rdoc/
+
+## Environment normalization:
+/.bundle/
+/vendor/bundle
+/lib/bundler/man/
+
+# for a library or gem, you might want to ignore these files since the code is
+# intended to run in multiple environments; otherwise, check them in:
+# Gemfile.lock
+# .ruby-version
+# .ruby-gemset
+
+# unless supporting rvm < 1.11.0 or doing something fancy, ignore this:
+.rvmrc

data/Gemfile

@@ -0,0 +1,6 @@
+source "https://rubygems.org"
+
+git_source(:github) {|repo_name| "https://github.com/#{repo_name}" }
+
+# Specify your gem's dependencies in japan_plane_rectangular.gemspec
+gemspec

data/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2018 dash14
+
+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,52 @@
+# Japan Plane Rectangular
+
+Conversion between world geodetic system and Japan plane rectangular coordinate system.
+
+世界測地系と日本平面直角座標系との変換
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'japan_plane_rectangular'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install japan_plane_rectangular
+
+## Usage
+
+```ruby
+require 'japan_plane_rectangular'
+
+# TokyoTower
+latlon = [35.658596, 139.745403]
+
+# Lat, Lon -> X, Y
+xy = JapanPlaneRectangular.to_xy(latlon, 9)
+p xy # [-37873.418395058005, -7961.358270500216]
+
+# X, Y -> Lat, Lon
+latlon = JapanPlaneRectangular.to_latlon(xy, 9)
+p latlon # [35.658596, 139.745403]
+
+# Get nearest zone number by Lat, Lon
+zone = JapanPlaneRectangular.nearest_zone(latlon)
+p zone # 9
+```
+
+## Development
+
+After checking out the repo, run `bin/setup` to install dependencies. 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/dash14/ruby-japan-plane-rectangular.

data/Rakefile

@@ -0,0 +1,10 @@
+require "bundler/gem_tasks"
+require "rake/testtask"
+
+Rake::TestTask.new(:test) do |t|
+  t.libs << "test"
+  t.libs << "lib"
+  t.test_files = FileList["test/**/*_test.rb"]
+end
+
+task :default => :test

data/bin/console

@@ -0,0 +1,14 @@
+#!/usr/bin/env ruby
+
+require "bundler/setup"
+require "japan_plane_rectangular"
+
+# 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(__FILE__)

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/japan_plane_rectangular.gemspec

@@ -0,0 +1,25 @@
+
+lib = File.expand_path("../lib", __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require "japan_plane_rectangular/version"
+
+Gem::Specification.new do |spec|
+  spec.name          = "japan_plane_rectangular"
+  spec.version       = JapanPlaneRectangular::VERSION
+  spec.authors       = ["dash14"]
+  spec.email         = ["dash14.ack@gmail.com"]
+
+  spec.summary       = %q{Conversion between WGS and Japan plane rectangular CS}
+  spec.homepage      = "https://github.com/dash14/ruby-japan-plane-rectangular"
+
+  spec.files         = `git ls-files -z`.split("\x0").reject do |f|
+    f.match(%r{^(test|spec|features)/})
+  end
+  spec.bindir        = "exe"
+  spec.executables   = spec.files.grep(%r{^exe/}) { |f| File.basename(f) }
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", "~> 1.16"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency "minitest", "~> 5.0"
+end

data/lib/japan_plane_rectangular.rb

@@ -0,0 +1,9 @@
+# frozen_string_literal: true
+
+require 'japan_plane_rectangular/version'
+require 'japan_plane_rectangular/functions'
+
+# Japan Plane Rectangular
+module JapanPlaneRectangular
+  extend JapanPlaneRectangular::Functions
+end

data/lib/japan_plane_rectangular/constants.rb

@@ -0,0 +1,37 @@
+
+module JapanPlaneRectangular
+  module Constants
+    # Positions of zone origins
+    ORIGINS = [
+      [33.00000, 129.5000],
+      [33.00000, 131.00000],
+      [36.00000, 132.166666666666669],
+      [33.00000, 133.50000],
+      [36.00000, 134.333333333333334],
+      [36.00000, 136.00000],
+      [36.00000, 137.166666666666666],
+      [36.00000, 138.5],
+      [36.00000, 139.833333333333334],
+      [40.00000, 140.833333333333334],
+      [44.00000, 140.25000],
+      [44.00000, 142.25000],
+      [44.00000, 144.25000],
+      [26.00000, 142.00000],
+      [26.00000, 127.50000],
+      [26.00000, 124.00000],
+      [26.00000, 131.00000],
+      [20.00000, 136.00000],
+      [26.00000, 154.000000]
+    ]
+
+    M0 = 0.9999
+
+    # Equatorial radius of earth
+    GRS80_ER = 6378137.0
+
+    # Inverse flattening
+    GRS80_IF = 298.257222101
+
+    ECCENTRICITY = Math.sqrt(2.0 * GRS80_IF - 1.0) / GRS80_IF
+  end
+end

data/lib/japan_plane_rectangular/functions.rb

@@ -0,0 +1,189 @@
+# frozen_string_literal: true
+
+require 'japan_plane_rectangular/constants'
+
+module JapanPlaneRectangular
+  module Functions
+    include JapanPlaneRectangular::Constants
+
+    # Find nearest zone number by [Latitude, Longitude]
+    # @param [Array<Float>] latlon [Latitude, Longitude]
+    # @return Integer zone number (1~19)
+    def nearest_zone(latlon)
+      min_index = 0
+      min_val = Float::MAX
+      ORIGINS.each_with_index do |origin, i|
+        xy1 = to_xy(latlon, 7)
+        xy2 = to_xy(origin, 7)
+        d = distance(xy1, xy2)
+        if d < min_val
+          min_index = i
+          min_val = d
+        end
+      end
+      min_index + 1 # zone number
+    end
+
+    # Convert to XY from latlon in zone
+    # @param [Array<Float>] latlon [latitude, longitude]
+    # @param [Integer] zone zone number (1~19)
+    # @return [Array<Float>] [x, y]
+    def to_xy(point, zone)
+      raise ArgumentError, 'Invalid zone number' unless 1 <= zone && zone <= 19
+      origin_point = ORIGINS[zone - 1]
+
+      phi0 = to_radian(origin_point[0])
+      lamda0 = to_radian(origin_point[1])
+
+      phi1 = to_radian(point[0])
+      lamda1 = to_radian(point[1])
+
+      s0 = meridian_arc_length(phi0)
+      s1 = meridian_arc_length(phi1)
+
+      ut = GRS80_ER / Math.sqrt(1.0 - ECCENTRICITY**2.0 * Math.sin(phi1)**2.0)
+      conp = Math.cos(phi1)
+      t1 = Math.tan(phi1)
+      dlamda = lamda1 - lamda0
+      eta2 = (ECCENTRICITY**2.0 / (1.0 - ECCENTRICITY**2.0)) * conp**2.0
+
+      v1 = 5.0 - t1**2.0 + 9.0 * eta2 + 4.0 * eta2**2.0
+      v2 = -61.0 + 58.0 * t1**2.0 - t1**4.0 - 270.0 * eta2 + 330.0 * t1**2.0 * eta2
+      v3 = -1385.0 + 3111.0 * t1**2.0 - 543.0 * t1**4.0 + t1**6.0
+
+      x = ((s1 - s0) + ut * conp**2.0 * t1 * dlamda**2.0 / 2.0 +
+          ut * conp**4.0 * t1 * v1 * dlamda**4.0 / 24.0 -
+          ut * conp**6.0 * t1 * v2 * dlamda**6.0 / 720.0 -
+          ut * conp**8.0 * t1 * v3 * dlamda**8.0 / 40320.0) * M0
+
+      v1 = -1.0 + t1**2.0 - eta2
+      v2 = -5.0 + 18.0 * t1**2.0 - t1**4.0 - 14.0 * eta2 + 58.0 * t1**2.0 * eta2
+      v3 = -61.0 + 479.0 * t1**2.0 - 179.0 * t1**4.0 + t1**6.0
+
+      y = (ut * conp * dlamda -
+          ut * conp**3.0 * v1 * dlamda**3.0 / 6.0 -
+          ut * conp**5.0 * v2 * dlamda**5.0 / 120.0 -
+          ut * conp**7.0 * v3 * dlamda**7.0 / 5040.0) * M0
+
+      [x, y]
+    end
+
+    # Convert to latlon from XY in zone
+    # @param [Array<Float>] xy [x, y]
+    # @param [Integer] zone zone number (1~19)
+    # @return [Array<Float>] [latitude, longitude]
+    def to_latlon(xy, zone)
+      raise ArgumentError, 'Invalid zone number' unless 1 <= zone && zone <= 19
+      origin_point = ORIGINS[zone - 1]
+      (x, y) = xy
+
+      phi0 = to_radian(origin_point[0])
+      lamda0 = to_radian(origin_point[1])
+
+      phi1 = perpendicular(x, phi0)
+
+      ut = GRS80_ER / Math.sqrt(1.0 - ECCENTRICITY**2.0 * Math.sin(phi1)**2.0)
+      conp = Math.cos(phi1)
+      t1 = Math.tan(phi1)
+      eta2 = (ECCENTRICITY**2.0 / (1.0 - ECCENTRICITY**2.0)) * conp**2.0
+
+      yy = y / M0
+      v1 = 1.0 + eta2
+      v2 = 5.0 + 3.0 * t1**2.0 + 6.0 * eta2 - 6.0 * t1**2.0 * eta2 - 3.0 * eta2**2.0 - 9.0 * t1**2.0 * eta2**2.0
+      v3 = 61.0 + 90.0 * t1**2.0 + 45.0 * t1**4.0 + 107.0 * eta2 - 162.0 * t1**2.0 * eta2 - 45.0 * t1**4.0 * eta2
+      v4 = 1385.0 + 3633.0 * t1**2.0 + 4095.0 * t1**4.0 + 1575.0 * t1**6.0
+
+      phir = -(v1 / (2.0 * ut**2.0)) * yy**2.0
+      phir += (v2 / (24.0 * ut**4.0)) * yy**4.0
+      phir -= (v3 / (720.0 * ut**6.0)) * yy**6.0
+      phir += (v4 / (40320.0 * ut**8.0)) * yy**8.0
+      phir *= t1
+      phir += phi1
+      phir = to_degree(phir)
+
+      v1 = ut * conp
+      v2 = 1.0 + 2.0 * t1**2.0 + eta2
+      v3 = 5.0 + 28.0 * t1**2.0 + 24.0 * t1**4.0 + 6.0 * eta2 + 8.0 * t1**2.0 * eta2
+      v4 = 61.0 + 662.0 * t1**2.0 + 1320.0 * t1**4.0 + 720.0 * t1**6.0
+
+      lamdar = (1.0 / v1) * yy
+      lamdar -= (v2 / (6.0 * ut**2.0 * v1)) * yy**3.0
+      lamdar += (v3 / (120.0 * ut**4.0 * v1)) * yy**5.0
+      lamdar -= (v4 / (5040.0 * ut**6.0 * v1)) * yy**7.0
+      lamdar += lamda0
+
+      lamdar = to_degree(lamdar)
+
+      [phir, lamdar]
+    end
+
+    private
+
+    def distance(p1, p2)
+      (x1, y1) = p1
+      (x2, y2) = p2
+      Math.sqrt((x2 - x1)**2.0 + (y2 - y1)**2.0)
+    end
+
+    def to_radian(degree)
+      degree * Math::PI / 180.0
+    end
+
+    def to_degree(radian)
+      radian * 180.0 / Math::PI
+    end
+
+    def perpendicular(x, phi0)
+      s0 = meridian_arc_length(phi0)
+      m = s0 + x / M0
+      cnt = 0
+      phin = phi0
+      e2 = ECCENTRICITY**2.0
+      phi0 = phin
+
+      loop do
+        cnt += 1
+        phi0 = phin
+        sn = meridian_arc_length(phin)
+        v1 = 2.0 * (sn - m) * ((1.0 - e2 * Math.sin(phin)**2.0)**1.5)
+        v2 = 3.0 * e2 * (sn - m) * Math.sin(phin) * Math.cos(phin) * Math.sqrt(1.0 - e2 * Math.sin(phin)**2.0) - 2.0 * GRS80_ER * (1.0 - e2)
+        phin += v1 / v2
+        break if ((phin - phi0).abs < 0.00000000000001) || cnt > 100
+      end
+      phin
+    end
+
+    def meridian_arc_length(lat_rad)
+      e0 = ECCENTRICITY
+      e2 = e0**2.0
+      e4 = e0**4.0
+      e6 = e0**6.0
+      e8 = e0**8.0
+      e10 = e0**10.0
+      e12 = e0**12.0
+      e14 = e0**14.0
+      e16 = e0**16.0
+
+      a = 1.0 + 3.0 / 4.0 * e2 + 45.0 / 64.0 * e4 + 175.0 / 256.0 * e6 + 11025.0 / 16384.0 * e8 + 43659.0 / 65536.0 * e10 + 693693.0 / 1048576.0 * e12 + 19324305.0 / 29360128.0 * e14 + 4927697775.0 / 7516192768.0 * e16
+      b = 3.0 / 4.0 * e2 + 15.0 / 16.0 * e4 + 525.0 / 512.0 * e6 + 2205.0 / 2048.0 * e8 + 72765.0 / 65536.0 * e10 + 297297.0 / 262144.0 * e12 + 135270135.0 / 117440512.0 * e14 + 547521975.0 / 469762048.0 * e16
+      c = 15.0 / 64.0 * e4 + 105.0 / 256.0 * e6 + 2205.0 / 4096.0 * e8 + 10395.0 / 16384.0 * e10 + 1486485.0 / 2097152.0 * e12 + 45090045.0 / 58720256.0 * e14 + 766530765.0 / 939524096.0 * e16
+      d = 35.0 / 512.0 * e6 + 315.0 / 2048.0 * e8 + 31185.0 / 131072.0 * e10 + 165165.0 / 524288.0 * e12 + 45090045.0 / 117440512.0 * e14 + 209053845.0 / 469762048.0 * e16
+      e = 315.0 / 16384.0 * e8 + 3465.0 / 65536.0 * e10 + 99099.0 / 1048576.0 * e12 + 4099095.0 / 29360128.0 * e14 + 348423075.0 / 1879048192.0 * e16
+      f = 693.0 / 131072 * e10 + 9009.0 / 524288.0 * e12 + 4099095.0 / 117440512.0 * e14 + 26801775.0 / 469762048.0 * e16
+      g = 3003 / 2097152.0 * e12 + 315315.0 / 58720256.0 * e14 + 11486475.0 / 939524096.0 * e16
+      h = 45045.0 / 117440512.0 * e14 + 765765.0 / 469762048.0 * e16
+      i = 765765.0 / 7516192768.0 * e16
+
+      meridian = GRS80_ER * (1.0 - e2) * (a * lat_rad -
+                  b * Math.sin(lat_rad * 2.0) / 2.0 +
+                  c * Math.sin(lat_rad * 4.0) / 4.0 -
+                  d * Math.sin(lat_rad * 6.0) / 6.0 +
+                  e * Math.sin(lat_rad * 8.0) / 8.0 -
+                  f * Math.sin(lat_rad * 10.0) / 10.0 +
+                  g * Math.sin(lat_rad * 12.0) / 12.0 -
+                  h * Math.sin(lat_rad * 14.0) / 14.0 +
+                  i * Math.sin(lat_rad * 16.0) / 16.0)
+      meridian
+    end
+  end
+end

data/lib/japan_plane_rectangular/version.rb

@@ -0,0 +1,5 @@
+# frozen_string_literal: true
+
+module JapanPlaneRectangular
+  VERSION = '0.1.0'
+end

metadata

@@ -0,0 +1,97 @@
+--- !ruby/object:Gem::Specification
+name: japan_plane_rectangular
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- dash14
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2018-02-17 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.16'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.16'
+- !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: minitest
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '5.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '5.0'
+description: 
+email:
+- dash14.ack@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- Gemfile
+- LICENSE
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- japan_plane_rectangular.gemspec
+- lib/japan_plane_rectangular.rb
+- lib/japan_plane_rectangular/constants.rb
+- lib/japan_plane_rectangular/functions.rb
+- lib/japan_plane_rectangular/version.rb
+homepage: https://github.com/dash14/ruby-japan-plane-rectangular
+licenses: []
+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.13
+signing_key: 
+specification_version: 4
+summary: Conversion between WGS and Japan plane rectangular CS
+test_files: []