Module: Mercator

Defined in:
lib/mercator.rb,
lib/mercator/datum.rb,
lib/mercator/version.rb,
lib/mercator/projection.rb

Defined Under Namespace

Classes: Datum, GRS80, Projection, Sweref99, Sweref991200, Sweref991330, Sweref991415, Sweref991500, Sweref991545, Sweref991630, Sweref991715, Sweref991800, Sweref991845, Sweref992015, Sweref992145, Sweref992315, Sweref99TM, WGS84

Constant Summary collapse

VERSION = 
"1.0.1"

Class Method Summary collapse

Class Method Details

.grid_to_spherical(north, east, datum, projection) ⇒ Object

Convert from grid coordinates to spherical coordinates using the specified datum and projection.



58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
 
# File 'lib/mercator.rb', line 58

def self.grid_to_spherical(north, east, datum, projection)
  e²  = datum.flattening * (2 - datum.flattening)
  n   = datum.flattening / (2 - datum.flattening)
  â   = datum.semi_major_axis/(1+n) * (1 + n**2/4 + n**4/64)
  λ₀  = projection.central_meridian * PI / 180.0

  ξ   = (north - projection.false_northing) / projection.scale / â
  η   = (east - projection.false_easting) / projection.scale / â

  δ₁  = n/2 - 2*n**2/3 + 37*n**3/96 - n**4/360
  δ₂  = n**2/48 + n**3/15 - 437*n**4/1440
  δ₃  = 17*n**3/480 - 37*n**4/840
  δ₄  = 4397*n**4/161280

  ξ′  = ξ - δ₁*sin(2*ξ)*cosh(2*η) - δ₂*sin(4*ξ)*cosh(4*η) -
        δ₃*sin(6*ξ)*cosh(6*η) - δ₄*sin(8*ξ)*cosh(8*η)
  η′  = η - δ₁*cos(2*ξ)*sinh(2*η) - δ₂*cos(4*ξ)*sinh(4*η) -
        δ₃*cos(6*ξ)*sinh(6*η) - δ₄*cos(8*ξ)*sinh(8*η)

  φ⋆  = asin(sin(ξ′) / cosh(η′))
  δ   = atan(sinh(η′) / cos(ξ′))

  a⋆  = e² + e²**2 + e²**3 + e²**4
  b⋆  = -1.0/6*(7*e²**2 + 17*e²**3 + 30*e²**4)
  c⋆  = (224*e²**3 + 889*e²**4)/120
  d⋆  = -1.0/1260*4279*e²**4

  λ   = λ₀ + δ
  φ   = φ⋆ + sin(φ⋆)*cos(φ⋆) * (a⋆ + b⋆*sin(φ⋆)**2 + c⋆*sin(φ⋆)**4 + d⋆*sin(φ⋆)**6)

  return [λ * 180.0 / PI, φ * 180.0 / PI]
end

.spherical_to_grid(latitude, longitude, datum, projection) ⇒ Object

Convert from spherical coordinates to grid coordinates using the specified datum and projection.



22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
 
# File 'lib/mercator.rb', line 22

def self.spherical_to_grid(latitude, longitude, datum, projection)
  φ   = latitude * PI / 180.0
  λ   = longitude * PI / 180.0
  λ₀  = projection.central_meridian.to_rad

  e²  = datum.flattening * (2 - datum.flattening)
  n   = datum.flattening / (2 - datum.flattening)
  â   = datum.semi_major_axis/(1+n) * (1 + n**2/4 + n**4/64)

  a   = e²
  b   = (5*e²**2 - e²**3)/6
  c   = (104*e²**3 - 45*e²**4)/120
  d   = 1237*e²**4/1260

  φ′  = φ - sin(φ)*cos(φ)*(a + b*sin(φ)**2 + c*sin(φ)**4 + d*sin(φ)**6)
  δ   = λ - λ₀

  ξ   = atan(tan(φ′) / cos(δ))
  η   = atanh(cos(φ′) * sin(δ))

  β₁  = n/2 - 2*n**2/3 + 5*n**3/16 + 41*n**4/180
  β₂  = 13*n**2/48 - 3*n**3/5 + 557*n**4/1440
  β₃  = 61*n**3/240 - 103*n**4/180
  β₄  = 49561*n**4/161280

  x   = projection.scale * â * (ξ + β₁*sin(2*ξ)*cosh(2*η) + β₂*sin(4*ξ)*cosh(4*η) +
        β₃*sin(6*ξ)*cosh(6*η) + β₄*sin(8*ξ)*cosh(8*η) ) + projection.false_northing
  y   = projection.scale * â * (η + β₁*cos(2*ξ)*sinh(2*η) + β₂*cos(4*ξ)*sinh(4*η) +
        β₃*cos(6*ξ)*sinh(6*η) + β₄*cos(8*ξ)*sinh(8*η) ) + projection.false_easting

  return [x, y]
end