jtsk 0.0.1

checksums.yaml

@@ -0,0 +1,15 @@
+---
+!binary "U0hBMQ==":
+  metadata.gz: !binary |-
+    ZjdkZWMxZmZhYjhmNDMwNWVhZDNlNTMzYzE4M2VlMDFiMzJmZDIzYw==
+  data.tar.gz: !binary |-
+    M2UxOWU1Y2FjYmE2YzM3YzZmMWI3NzZmMDRlYWIzZDMwYzIxNWM4Mg==
+SHA512:
+  metadata.gz: !binary |-
+    MjRiMjAwNTE0OTI3ZGQyNmU1NTU5M2IzMTAwNjFkNWI5NjlkMDkwZDE0OThj
+    MjA1MmMwNjQ5ODQ2ZWRjNWU5MGM2ODEwMTQ3ZTA4NWJiZjI3MDVhMDhjMTk3
+    NmYxNjA3ZDdhNDVkNjcyMTA1YmZmYTI4OTRjNTQwOGMzZDIxYzM=
+  data.tar.gz: !binary |-
+    MmM5NDY1Y2M3NzBmMGY4ZDUzZDg5ZWY1ZTFhYTRiMWFmNDQyYzM4ZmVmYmJk
+    N2JkOGYzNTM4NDczZDI1MmViNTdkYWFhYTg3NDViNTMwYjE4ZmJkZmJkYmMz
+    OTE2OTFhNzQzY2QxYjFlM2IxN2QyZjM2N2Y5NWY3NTZjODkzY2I=

data/.gitignore

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

data/.rspec

@@ -0,0 +1,3 @@
+--color
+--warnings
+--require spec_helper

data/Gemfile

@@ -0,0 +1,6 @@
+source 'https://rubygems.org'
+
+gem 'rspec', '3.0.0.beta2'
+
+# Specify your gem's dependencies in jtsk.gemspec
+gemspec

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2014 Josef Šimánek
+
+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,38 @@
+# JTSK  [![Build Status](https://travis-ci.org/simi/jtsk.svg?branch=master)](https://travis-ci.org/simi/jtsk)
+
+Convert from JTSK to WGS48.
+
+Based on [JTSK Converter](https://github.com/josefzamrzla/JTSK_Converter) by [@josefzamrzla](https://github.com/josefzamrzla).
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+    gem 'jtsk'
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install jtsk
+
+## Usage
+
+```ruby
+converter = JTSK::Converter.new
+x = 1043033.89
+y = 738371.58
+result = converter.to_wgs48(x, y) #=> #<struct JTSK::Wgs48Result latitude=50.092696246901404, longitude=14.482746557404647>
+result.latitude #=> 50.092696246901404
+result.longitude #=> 14.482746557404647
+```
+
+## Contributing
+
+1. Fork it ( https://github.com/[my-github-username]/jtsk/fork )
+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 a new Pull Request

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/rake

@@ -0,0 +1,16 @@
+#!/usr/bin/env ruby
+#
+# This file was generated by Bundler.
+#
+# The application 'rake' is installed as part of a gem, and
+# this file is here to facilitate running it.
+#
+
+require 'pathname'
+ENV['BUNDLE_GEMFILE'] ||= File.expand_path("../../Gemfile",
+  Pathname.new(__FILE__).realpath)
+
+require 'rubygems'
+require 'bundler/setup'
+
+load Gem.bin_path('rake', 'rake')

data/jtsk.gemspec

@@ -0,0 +1,23 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'jtsk/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "jtsk"
+  spec.version       = Jtsk::VERSION
+  spec.authors       = ["Josef Šimánek"]
+  spec.email         = ["josef.simanek@gmail.com"]
+  spec.summary       = %q{Convert from JTSK to WGS48.}
+  spec.homepage      = ""
+  spec.license       = "MIT"
+
+  spec.files         = `git ls-files -z`.split("\x0")
+  spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", ">= 1.5"
+  spec.add_development_dependency "rake"
+  spec.add_development_dependency "rspec"
+end

data/lib/jtsk.rb

@@ -0,0 +1,5 @@
+require 'jtsk/converter'
+require 'jtsk/result'
+
+module JTSK
+end

data/lib/jtsk/converter.rb

@@ -0,0 +1,206 @@
+module JTSK
+  class Converter
+    EPS = 1e-4
+
+    def to_wgs48(x, y)
+      delta = 5.0
+      latitude = 49.0
+      longitude = 14.0
+      steps = 0.0
+
+      loop {
+
+        jtsk = self.wgs48_to_jtsk(latitude - delta, longitude - delta)
+        if(jtsk.x && jtsk.y)
+          v1 = self.dist_points(jtsk.x, jtsk.y, x, y)
+        else
+          v1 = 1e32
+        end
+
+        jtsk = self.wgs48_to_jtsk(latitude - delta, longitude + delta)
+        if(jtsk.x && jtsk.y)
+          v2 = self.dist_points(jtsk.x, jtsk.y, x, y)
+        else
+          v2 = 1e32
+        end
+
+        jtsk = self.wgs48_to_jtsk(latitude + delta, longitude - delta)
+        if(jtsk.x && jtsk.y)
+          v3 = self.dist_points(jtsk.x, jtsk.y, x, y)
+        else
+          v3 = 1e32
+        end
+
+        jtsk = self.wgs48_to_jtsk(latitude + delta, longitude + delta)
+        if(jtsk.x && jtsk.y)
+          v4 = self.dist_points(jtsk.x, jtsk.y, x, y)
+        else
+          v4 = 1e32
+        end
+
+        if ((v1 <= v2) && (v1 <= v3) && (v1 <= v4))
+          latitude = latitude - delta / 2.0
+          longitude = longitude - delta / 2.0
+        end
+
+        if ((v2 <= v1) && (v2 <= v3) && (v2 <= v4))
+          latitude = latitude - delta / 2.0
+          longitude = longitude + delta / 2.0
+        end
+
+        if ((v3 <= v1) && (v3 <= v2) && (v3 <= v4))
+          latitude = latitude + delta / 2.0
+          longitude = longitude - delta / 2.0
+        end
+
+        if ((v4 <= v1) && (v4 <= v2) && (v4 <= v3))
+          latitude = latitude + delta / 2.0
+          longitude = longitude + delta / 2.0
+        end
+
+        delta = delta * 0.55;
+        steps = steps + 4.0;
+ 
+        break if (((delta < 0.00001) || steps > 1000.0));
+      } 
+
+      JTSK::Wgs48Result.new(latitude, longitude)
+    end
+
+    def wgs48_to_jtsk(latitude, longitude)
+      if ((latitude < 40.0) || (latitude > 60.0) || (longitude < 5.0) || (longitude > 25.0)) 
+        JTSK::JtskResult.new(0.0, 0.0)
+      else
+        bessel = self.wgs48_to_bessel(latitude, longitude);
+        self.bessel_to_jtsk(bessel.latitude, bessel.longitude);
+      end
+    end
+
+    def wgs48_to_bessel(latitude, longitude, altitude = 0.0)
+      b = deg2rad(latitude)
+      l = deg2rad(longitude)
+      h = altitude;
+
+      x1, y1, z1 = self.blh_to_geo_coords(b, l, h)
+      x2, y2, z2 = self.transform_coords(x1, y1, z1)
+      b, l, h = self.geo_coords_to_blh(x2, y2, z2)
+
+      latitude = rad2deg(b);
+      longitude = rad2deg(l);
+      # altitude = h;
+
+      JTSK::BesselResult.new(latitude, longitude)
+    end
+
+    def bessel_to_jtsk(latitude, longitude)
+      a     = 6377397.15508
+      e     = 0.081696831215303
+      n     = 0.97992470462083
+      rho_0 = 12310230.12797036
+      sinUQ = 0.863499969506341
+      cosUQ = 0.504348889819882
+      sinVQ = 0.420215144586493
+      cosVQ = 0.907424504992097
+      alfa  = 1.000597498371542
+      k_2   = 1.00685001861538
+
+      b = deg2rad(latitude)
+      l= deg2rad(longitude)
+
+      sinB = Math.sin(b)
+      t = (1 - e * sinB) / (1 + e * sinB)
+      t = pow(1 + sinB, 2) / (1 - pow(sinB, 2)) * Math.exp(e * Math.log(t))
+      t = k_2 * Math.exp(alfa * Math.log(t))
+
+      sinU  = (t - 1) / (t + 1)
+      cosU  = Math.sqrt(1 - sinU * sinU)
+      v     = alfa * l
+      sinV  = Math.sin(v)
+      cosV  = Math.cos(v)
+      cosDV = cosVQ * cosV + sinVQ * sinV
+      sinDV = sinVQ * cosV - cosVQ * sinV
+      sinS  = sinUQ * sinU + cosUQ * cosU * cosDV
+      cosS  = Math.sqrt(1 - sinS * sinS)
+      sinD  = sinDV * cosU / cosS
+      cosD  = Math.sqrt(1 - sinD * sinD)
+
+      eps = n * Math.atan(sinD / cosD)
+      rho = rho_0 * Math.exp(-n * Math.log((1 + sinS) / cosS))
+      
+      JTSK::JtskResult.new(rho * Math.cos(eps), rho * Math.sin(eps))
+    end
+
+    def blh_to_geo_coords(b,l,h)
+      # WGS-84 ellipsoid parameters
+      a   = 6378137.0;
+      f_1 = 298.257223563;
+      e2  = 1 - ((1 - 1 / f_1) ** 2)
+      rho = a / Math.sqrt(1 - e2 * (Math.sin(b) ** 2))
+      x = (rho + h) * Math.cos(b) * Math.cos(l)
+      y = (rho + h) * Math.cos(b) * Math.sin(l)
+      z = ((1 - e2) * rho + h) * Math.sin(b);
+
+      return x, y, z;
+    end
+
+    def transform_coords(xs, ys, zs)
+      # coeficients of transformation from WGS-84 to JTSK
+      dx = -570.69; dy = -85.69; dz = -462.84; # shift
+
+      # rotation
+      wx = 4.99821/3600.0 * Math::PI / 180.0
+      wy = 1.58676/3600.0 * Math::PI / 180.0
+      wz = 5.2611/3600.0 * Math::PI / 180.0
+      m  = -3.543e-6; # scale
+
+      xn = dx + (1 + m) * (+xs + wz * ys - wy * zs);
+      yn = dy + (1 + m) * (-wz * xs + ys + wx * zs);
+      zn = dz + (1 + m) * (+wy * xs - wx * ys + zs);
+
+      return xn, yn, zn
+    end
+
+    def geo_coords_to_blh(x, y, z)
+      # Bessel's ellipsoid parameters
+      a   = 6377397.15508
+      f_1 = 299.152812853
+      a_b = f_1 / (f_1-1)
+      p   = Math.sqrt(pow(x, 2) + pow(y, 2))
+      e2  = 1 - pow(1 - 1 / f_1, 2)
+      th  = Math.atan(z * a_b / p)
+      st  = Math.sin(th)
+      ct  = Math.cos(th)
+      t   = (z + e2 * a_b * a * pow(st, 3)) / (p - e2 * a * pow(ct, 3))
+
+      b = Math.atan(t)
+      h = Math.sqrt(1 + t * t) * (p - a / Math.sqrt(1 + (1 - e2) * t * t))
+      l = 2 * Math.atan(y / (p + x))
+
+      return b, l, h
+    end
+
+    def dist_points(x1, y1, x2, y2)
+      dist = Math.hypot(x1 - x2, y1 - y2)
+
+      if (dist < self.class::EPS)
+        return 0
+      end
+
+      return dist
+    end
+
+    private
+
+    def deg2rad(degrees)
+      degrees * (Math::PI / 180.0 )
+    end
+
+    def rad2deg(radians)
+      radians * (180.0 / Math::PI)
+    end
+
+    def pow(a,b)
+      a ** b
+    end
+  end
+end

data/lib/jtsk/result.rb

@@ -0,0 +1,8 @@
+module JTSK
+  class Wgs48Result < Struct.new(:latitude, :longitude)
+  end
+  class BesselResult < Struct.new(:latitude, :longitude)
+  end
+  class JtskResult < Struct.new(:x, :y)
+  end
+end

data/lib/jtsk/version.rb

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

data/spec/jtsk/converter_spec.rb

@@ -0,0 +1,38 @@
+require 'spec_helper'
+
+describe JTSK::Converter do
+  let(:converter) { described_class.new }
+  let(:x) { 1043033.89 }
+  let(:y) { 738371.58 }
+  let(:result) { converter.to_wgs48(x,y) }
+
+  it 'EPS is properly set' do
+    expect(described_class::EPS).to eq(1e-4)
+  end
+
+  it 'can be initialized' do
+    expect(converter).to be_kind_of(described_class)
+  end
+
+  it 'result is JTSK::Result' do
+    expect(result).to be_kind_of(JTSK::Wgs48Result)
+  end
+
+  context 'Club Velbloud' do
+    let(:x) { 1043033.89 }
+    let(:y) { 738371.58 }
+    it 'calucaltes well' do
+      expect(result.latitude).to eq(50.092696246901404)
+      expect(result.longitude).to eq(14.482746557404647)
+    end
+  end
+
+  context 'Indian Bar' do
+    let(:x) { 1043070.32 }
+    let(:y) { 738988.72 }
+    it 'calculates well' do
+      expect(result.latitude).to eq(50.091628136117045)
+      expect(result.longitude).to eq(14.474266228665831)
+    end
+  end
+end

data/spec/spec_helper.rb

@@ -0,0 +1,5 @@
+require 'bundler/setup'
+require 'jtsk'
+
+RSpec.configure do |config|
+end

metadata

@@ -0,0 +1,103 @@
+--- !ruby/object:Gem::Specification
+name: jtsk
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- Josef Šimánek
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2014-05-13 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '1.5'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '1.5'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ! '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+description: 
+email:
+- josef.simanek@gmail.com
+executables:
+- rake
+extensions: []
+extra_rdoc_files: []
+files:
+- .gitignore
+- .rspec
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/rake
+- jtsk.gemspec
+- lib/jtsk.rb
+- lib/jtsk/converter.rb
+- lib/jtsk/result.rb
+- lib/jtsk/version.rb
+- spec/jtsk/converter_spec.rb
+- spec/spec_helper.rb
+homepage: ''
+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.1.10
+signing_key: 
+specification_version: 4
+summary: Convert from JTSK to WGS48.
+test_files:
+- spec/jtsk/converter_spec.rb
+- spec/spec_helper.rb