arxiv: v2 [math.ho] 5 May 2018

Transcription

1 Longest Straight Line Paths on Water or Land on the Earth Rohan Chabukswar 1 and Kushal Mukherjee 1 United Tehnologies Researh Center Ireland IBM Researh India May 8, 018 arxiv: v [math.ho] 5 May 018 Abstrat There has been some interest reently in determining the longest distane one an sail for on the earth without hitting land, as well as in the onverse problem of determining the longest distane one ould drive for on the earth without enountering a major body of water. In its basi form, this is an optimisation problem, rendered haoti by the presene of islands and lakes, and indeed the fratal nature of the oasts. In this paper we present a methodology for alulating the two paths using the branh-and-bound algorithm. 1 Introdution On Deember 9, 01, Reddit user kepleronlyknows posted a map Figure 1) to r/mapporn [1] showing the longest straight line that an be sailed for on the earth without hitting land, from Pakistan to eastern Russia. This was generated a lot of interest and led to subsequent attempts to prove and disprove the user, along with disussions on solving the onverse problem of determining the longest distane that an be driven on land without hitting a major water body. ChabukR@utr.ut.om KushMukh@in.ibm.om Figure 1: kepleronlyknows s original map posted to r/mapporn 1

2 In their basi forms, the two are optimisation problems in two dimensions however, the objetive funtions are disontinuous due to the presene of islands in the ase of the longest sailable path), and lakes in the ase of the longest drivable path). Globally optimal solutions would require an exhaustive searh over the two dimensions, a not inonsiderable obstale in terms of time and memory. For a 1-arminute resolution km at the equator), whih is the resolution provided by the ETOPO1 Global Relief model of Earth s surfae by National Oeani and Atmospheri Administration, there would be 33,80,000 great irles to onsider to find the global optimum, and eah great irle would have 1600 individual points to proess a staggering 5,038,88,000,000 points to verify. 1.1 Previous Work Apart from the disovery of the longest sailable path by kepleronlyknows, Reddit user Groke also posted [] a path from Norway to Antartia in a straight line without hitting land, whih is shorter than the path from Pakistan to Eastern Russia. YouTube user David Cooke laimed [3] to have disovered an even longer path The Cooke Passage ) from Quebe to British Columbia, whih was later proved [] to not be a straight line in a blog post at GeoGarage, a nautial harts platform. Guy Bruneau of IT/GIS Consulting servies alulated [5] a path from Eastern China to Western Liberia as being the longest distane you an travel between two points in straight line without rossing any oean or any major water bodies. However, the path rosses through the Dead Sea whih an be onsidered to be a major water body), and hene does not satisfy the onstraints originally set out. Problem Statement Thus, in an effort to prove or disprove kepleronlyknows, and to find the longest straight line path on land without rossing a major water body, the authors set out to systematially arrive at the globally optimal solutions to the problem. The land/water data for the whole earth required to set up the problem was not available readily. In the end, the data set used was the ETOPO1 Global Relief Model [6], and onsider any negative altitude to be water and positive altitude to be land. Although this is not stritly true in general, in the absene of atual data it is the losest approximation to real relief that is freely available..1 ETOPO1 Global Relief Model ETOPO1 is a 1 ar-minute global relief model of Earth s surfae that integrates land topography and oean bathymetry. It was built from numerous global and regional data sets, and is available in Ie Surfae top of Antarti and Greenland ie sheets) and Bedrok base of the ie sheets) versions. The version suited for this problem was the Ie Surfae one, sine a sailable path should not hit an ie layer, and tehnially ie ould be driven upon.. Optimisation The real rux of the problems lay in solving the optimisation problem globally. After an initial aborted attempt to do a brute-fore searh aross the paths, the authors settled on utilising the branh-and-bound algorithm. 3 Problem Formulation The first step in formulating the problem was to implement a method for enumerating the different straight line paths on the Earth. To do this, it is important to note that all straight line paths lie on a great irle. Eah great irle an have as many paths as there are rossings of the great irle with the land-water boundary, but sine the interest lies in only the longest paths either land or water) in eah great irle, enumerating all the great irles would be enough. This an be done in a systemati way by observing that any great irle exept the equator itself) intersets the equator at two antipodal points, and at a speifi angle. Exepting the trivial ase of the equator, the enumeration an thus be hosen by hoosing the longitude of the intersetion with the equator the origin ), and the angle between the

3 great irle and the equator at that point the heading ), either of whih ould be hosen between 0 and 360. This ounts eah great irle four times, however, so the atual hoie an be restrited to between 0 and 180 for eah. Eah point of the great irle an then be pinpointed by the angle of the vetor in the plane of the great irle the referene vetor an be arbitrarily set to the origin). The domain of this angle φ) is between 0 and 360. Any great irle defined by the origin δ) and the heading α) an thus be denoted by GC δ, α). To math the ETOPO1 data set, whih is given as altitude in meters at eah latitude and longitude, eah point φ) on eah great irle GC δ, α) needs to be translated into latitude β) and longitude λ), whih an be done using basi spherial trigonometry: β = sin 1 sin α sin φ), 1) λ = δ + Tan 1 os α sin φ, os φ), ) 1 where the Tan y, x) is the four-quadrant artangent funtion. Traversing eah great irle defined by β, λ), the relief enountered rδ,α) φ), φ 180, 180 ] is given by rδ,α) φ) = etopo1 sin 1 sin α sin φ), δ + Tan 1 os α sin φ, os φ) 3) For example, as shown in Figure, GC 163 0, ) passes through both the Challenger Deep N, 1 10 E) at an angular distane of 10 from the origin, and Mount Everest N, E) at an angular distane of 78 0 from the origin. Figure : Great irle passing through the lowest and highest points on the Earth. The relief enountered r163 0,151 0 ) φ), φ 180, 180 ] is shown in Figure 3a. Given the relief along eah great irle, the assumption of land and water being given by positive and negative heights above sea level an be applied to lassify eah point of the great irle being either land or water: if rδ,α) φ) > 0 1 landδ,α) φ) = 0 ) if rδ,α) φ) < 0. indeterminate if rδ,α) φ) = 0 The indeterminay is resolved onservatively by assuming rδ,α) φ) = 0 to be land when searhing for the longest water path and water when searhing for the longest land path. For the example great irle, land163 0,151 0 ) φ), φ 180, 180 ] is shown in Figure 3b. One the lassifiation is obtained, it is a simple matter to determine the longest land and water paths in the great irle, whih form the objetive funtions that need to be maximised. 3

4 a) Atual relief b) Land/Water lassifiation Figure 3: Relief along the great irle passing through the lowest and highest points on the Earth. Branh-And-Bound Branh-and-bound is a method to solve ombinatorial optimisation problems, espeially when it is possible to get bounds on the objetive funtion for a subset in the spae of solutions. For the problem addressed in this paper, we ontinue to denote eah great irle by the pair of angles that uniquely haraterise it as GC δ, α). A set of solutions refers to a set of great irle trajetories that interset the equator between latitudes δ min and δ max and with heading in the range α min to α max. We denote the set as [δ min, δ max ] [α min, α max ]. Eah great irle GC δ, α) has an assoiated length l δ, α) that is defined as the maximum distane that an be traversed along the great irle without hitting land..1 Relaxation Tehniques The next task is to obtain an upper bound on the length of longest ontinuous setion of sea l δ, α) along all great irles suh that δ [δ min, δ max ], α [α min, α max ]. Theorem 1. Considering GC δ min+δ max ), αmin+αmax as representative of the set of great irles uh that δ [δ min, δ max ], α [α min, α max ], the maximum angular separation between the representative great irle and any other great irle in the set is always bounded by δ max δ min + α max α max. Proof. Let θ = αmax+αmin, θ = α max α min, and φ = δ max δ min. Using ) spherial trigonometry, )) the maximum angular separation between them is os 1 os θ) os φ + os θ) sin φ. This is non-trivial to use, however, the maximum deviation ours when θ = 90, and has a theoretial upper limit of θ + φ, as θ, φ 0. Even for θ, φ = 5, the limit is aurate to.6%. Sine we know that θ + φ θ + φ, the angular separation is always bounded by this value, whih is, δ max δ min + α max α max. Consider the great irle GC δ min+δ max ) as representative of the set of great irles. Using, αmin+αmax Theorem 1, the maximum separation between the representative great irle and any other great irle in the set is always bounded by R δ max δ min + α max α max ), where R is the radius of the earth.

5 The novelty in our approah lies in a method for evaluation the upper bound on the longest ontinuous setion of sea. Let there be a great irle GC δ o, α o ) where δ o [δ min, δ max ] and α o [α min, α max ] whih has the longest ontinuous setion of sea. That implies, l δ o, α o ) = max l δ, α). 5) δ [δ min,δ max] α [α min,α max] l δ, α) max l δ, α), δ [δ min, δ max ], α [α min, α max ]. 6) δ [δ min,δ max] α [α min,α max] Further, onsider the globe where the sea has eroded into the land by a distane ɛ. In this new globe ɛ-eroded globe), the land masses have shrunk by a distane of epsilon from all sides. All piees of land, that are within a distane of ɛ from the sea, are overed by sea on the ɛ-globe. We define l ɛ δ, α) as the maximum distane that an be traversed along the great irle without hitting land on the ɛ-eroded globe. Clearly, l ɛ δ, α) max δ [δ min,δ max] α [α min,α max] l δ, α), δ [δ min, δ max ], α [α min, α max ]. 7) Let s onsider the ase where ɛ R δ max δ min + α max α max ). That is, the extent of erosion is larger that the maximum distane between any great irle in the set and the representative great irle. Sine, the distane between the great irle δ o, α o ) and the representative great irle is less that ɛ, the longest setion over sea) of the great irle δ o, α o ) would erode piees of land that lie on the representative great irle. Thus, the longest setion over sea) of the representative great irle on the ɛ-eroded globe would be at least as large as l δ o, α o ). δmin + δ max l ɛ, α ) min + α max l δ, α). 8) Hene the upper bound of l δ, α), δ [δ min, δ max ], α [α min, α max ] an be evaluated as the length of the longest streth of the representative great irle on the ɛ-eroded globe suh that ɛ R δmax δmin + αmax αmax ). The erosion should tehnially be arried out in three dimensions assuming that the angular landwater boundaries ontinue till the entre). However, this would be omputationally intensive, to store and alulate at the resolution required. Sine the ETOPO1 database is equivalent to the equiretangular projetion, in the two-dimensional ontext, the eroding element would have to be ontinuously deformed as per the Tissot s indiatrix for equiretangular projetion, as shown in Figure. Continuously deforming the eroding element takes the problem to using non-standard algorithms for erosion. Instead, for speed and omplexity, we hoose to use a single over-estimated erosion radius for the whole map. Over-estimating the erosion radius redues the relaxation, and hene might potentially affet the omputation required, but does not affet the results. The fator for over-estimation, however, needs to be hosen with are. It is ertain, due to the existene of Antartia for water), and the absene of land near the North Pole, that the great irle with maximum distane won t extend beyond 78 8 North or South. Thus the erosion only needs to be over-estimated till that latitude, giving us a reasonable fator of 5.. Branh-And-Bound Algorithm Algorithm 1 depits the standard branh-and-bound algorithm that has been developed to address the problem of finding the longest straight line path over sea. The suessive branhing of the set of trajetories into four subsets by dividing the origin and heading into two subsets eah) based on their upper bounds enables the algorithm to quikly onverge at the solution. 5

6 Figure : The equiretangular projetion with Tissot s indiatrix of deformation, from [7]. Algorithm 1 Algorithm to obtain the great irle with the longest setion over sea 1: Initialise Nodes with one element: δ min = 0, δ max = π, α min = 0, α max = π, ɛ = πr, B = πr). : Remove element δmin, δ max, αmin, α max, ɛ, B ) with largest B from Nodes. 3: while ɛ ɛ threshold do : Add δmin, δ max +δ min, αmin, α max +α min, ɛ = ɛ, B = l ɛ δ max +δ min 5: Add, δmax, αmin, α max +α min, ɛ = ɛ, B = l ɛ 6: Add δmin, δ max +δ min, α max +α min, αmax, ɛ = ɛ, B = l ɛ δ max 7: Add +δ min, δmax, α max +α min, αmax, ɛ = ɛ, B = l ɛ 3δ min +δ max, 3α min +α max δ min +3δ max, 3α min +α max 3δ min +δ max, α min +3α max δ min +3δ max, α min +3α max 8: Remove element δmin, δ max, αmin, α max, ɛ, B ) with largest B from Nodes. 9: Return δ min, δ max, α min, α max, ɛ, B ). )) to Nodes. )) to Nodes. )) to Nodes. )) to Nodes. 6

7 5 5.1 Results Erosion a) Water Mass Erosion b) Land Mass Erosion Figure 5: Erosion steps 3 7 for water 5a) and land 5b) masses respetively. The erosion steps 3 7 for both water and land an be seen in Figures 5a and 5b respetively. The nuanes of the erosion are better seen by zooming in on the British Isles, as has been done in Figures 6a and 6b. 5. Longest Paths The algorithm returned the longest path in about 10 minutes of omputation for water path, and 5 minutes of omputation for land path on a standard laptop. 7

8 a) Water Mass Erosion b) Land Mass Erosion Figure 6: Erosion steps 3 7 for water 6a) and land 6b) masses respetively, around the British Isles. 8

9 5..1 Longest Sailable Straight Line Path on Earth The algorithm returned the path shown in Figure 7. Although it does not look like a straight line on the map, the algorithm using great irles ensures that it is. Figure 7: Longest Sailable Straight Line Path on Earth. The fat that it is indeed a straight line an be seen from Figures 8a, 8b, and 8, whih are shown on a globe with a perspetive from above the path. Figure 8a shows the path originating in Sonmiani, Las Bela, Balohistan, Pakistan N, E), threading the needle between Afria and Madagasar, between Antartia and Tiera del Fuego in South Ameria Figure 8b), and ending in Karaginsky Distrit, Kamhatka Krai, Russia N, E) Figure 8). The path overs an astounding total angular distane of , for a distane of kilometres. This path is visually the same one as found by kepleronlyknows, thus proving his assertion. a) Start of Path b) Threading the Needle ) End of Path Figure 8: Longest Sailable Straight Line Path on Earth. 5.. Longest Drivable Straight Line Path on Earth The algorithm returned the path shown in Figure 9. Again, although it does not look like a straight line on the map, the algorithm using great irles ensures that it is. 9

10 Figure 9: Longest Drivable Straight Line Path on Earth. The fat that it is indeed a straight line an be seen from Figure 10, whih is shown on a globe with a perspetive from above the path. Figure 9 shows the path originating near Jinjiang, Quanzhou, Fujian, China 330 N, E), weaving through China and Mongolia for a bit, passing though Kazakhstan and Russia to further weave through Belarus and Ukraine, and passing through Poland, Czeh Republi, Germany, Austria, Liehtenstein, Switzerland, Frane, and Spain, to end near Sagres, Portugal 37 0 N W), traversing a total of 15 ountries. Although not as long as the longest sailable path, the longest drivable path overs a still-respetable total angular distane of , for a distane of 11 1 kilometres. Figure 10: Longest Drivable Straight Line Path on Earth. 6 Limitations and Future Work As mentioned in Setion, our approah assumes the partition between land and water to be determined by the height with respet to the mean sea level. This unfortunately disounts, among other features, highland rivers, and low-lying plains. We also annot onsider bridges, although the probability of a bridge lying at the preise loation and diretion of any given path is negligible. Seondly, the ETOPO1 database has a resolution of 1-arminute, about 1.8 kilometres at the equator. The optimisation problem, therefore, has by neessity used this as the lower-limit of auray for analysis. Both of these fators are the result of the input dataset. 10

11 Thirdly, an impliit assumption is that the Earth is perfetly spherial. This is untrue beause of the ellipsoidal shape of the earth, as well as due to geoid undulations. Future work would involve running the analysis on a higher resolution dataset, or on one that provides the shape of all the land and water bodies in higher detail. It would also require us to stop using spherial trigonometry, or at least, onsider deviations due to the non-spherial nature of the geoid, using loal urvatures. Consideration all of these issues, and more, does seem to be an overkill for what is essentially a rereational problem. 7 Conlusion We proposed an innovative approah for relaxation of an optimisation problem for utilising the branhand-bound algorithm. On the way, we managed to prove that kepleronlyknows was right about the longest sailable straight line path on the Earth, and we found the analogous longest drivable straight line path. The problem was approahed as a purely mathematial exerise. The authors do not reommend sailing or driving along the found paths. Referenes [1] kepleronlyknows Reddit). The longest straight line: You an sail almost 0,000 miles in a straight line from Pakistan to the Kamhatka Peninsula, Russia. Last retrieved May 8, 018, from an_sail_almost/. [] Groke Reddit). You an just barely sail north from Norway to Antartia in a straight line without hitting land. Album with lose-ups in omments. Retrieved May 8, 018, from om/r/mapporn/omments/8ktzt/you_an_just_barely_sail_north_from_norway_to/. [3] Cooke, David Reddit). A new reord has been disovered: the longest straight line over the oean without touhing land. You an sail from eastern Canada all the way around the world to western Canada in a straight line. Retrieved May 8, 018, from omments/36v0p/a_new_reord_has_been_disovered_the_longest/. [] GeoGarage. The Cooke passage : a new world s longest straight-line sail? Retrieved May 8, 018, from [5] Guy Bruneau IT/GIS Consulting Servies). Longest distane on land. Retrieved May 8, 018, from [6] Amante, C. and B.W. Eakins, 009. ETOPO1 1 Ar-Minute Global Relief Model: Proedures, Data Soures and Analysis. NOAA Tehnial Memorandum NESDIS NGDC-. National Geophysial Data Center, NOAA. doi:10.789/v5c876m [May 8, 018]. [7] Justin Kunimune Wikimedia). File:Plate Carrée with Tissot s Indiatries of Distortion.svg Last retrieved May 8, 018, from Tissot%7s_Indiatries_of_Distortion.svg. 11