flightgear-phi 2016.4.2+dfsg1-1 / usr / share / games / flightgear / Phi / 3rdparty / leaflet / Leaflet.Geodesic.min.js

if(typeof Number.prototype.toRadians=="undefined"){Number.prototype.toRadians=function(){return this*Math.PI/180}}if(typeof Number.prototype.toDegrees=="undefined"){Number.prototype.toDegrees=function(){return this*180/Math.PI}}var INTERSECT_LNG=179.999;L.Geodesic=L.MultiPolyline.extend({options:{color:"blue",steps:10,dash:1,wrap:true},initialize:function(latlngs,options){this.options=this._merge_options(this.options,options);this.datum={};this.datum.ellipsoid={a:6378137,b:6356752.3142,f:1/298.257223563};L.MultiPolyline.prototype.initialize.call(this,latlngs,this.options)},setLatLngs:function(latlngs){this._latlngs=this.options.dash<1?this._generate_GeodesicDashed(latlngs):this._generate_Geodesic(latlngs);L.MultiPolyline.prototype.setLatLngs.call(this,this._latlngs)},getStats:function(){var obj={distance:0,points:0,polygons:this._latlngs.length},poly,points;for(poly=0;poly<this._latlngs.length;poly++){obj.points+=this._latlngs[poly].length;for(points=0;points<this._latlngs[poly].length-1;points++){obj.distance+=this._vincenty_inverse(this._latlngs[poly][points],this._latlngs[poly][points+1]).distance}}return obj},createCircle:function(center,radius){var _geo=[],_geocnt=0;var prev={lat:0,lng:0,brg:0};var s;_geo[_geocnt]=[];var direct=this._vincenty_direct(L.latLng(center),0,radius,this.options.wrap);prev=L.latLng(direct.lat,direct.lng);_geo[_geocnt].push(prev);for(s=1;s<=this.options.steps;){direct=this._vincenty_direct(L.latLng(center),360/this.options.steps*s,radius,this.options.wrap);var gp=L.latLng(direct.lat,direct.lng);if(Math.abs(gp.lng-prev.lng)>180){var inverse=this._vincenty_inverse(prev,gp);var sec=this._intersection(prev,inverse.initialBearing,{lat:-89,lng:gp.lng-prev.lng>0?-INTERSECT_LNG:INTERSECT_LNG},0);if(sec){_geo[_geocnt].push(L.latLng(sec.lat,sec.lng));_geocnt++;_geo[_geocnt]=[];prev=L.latLng(sec.lat,-sec.lng);_geo[_geocnt].push(prev)}else{_geocnt++;_geo[_geocnt]=[];_geo[_geocnt].push(gp);prev=gp;s++}}else{_geo[_geocnt].push(gp);prev=gp;s++}}this._latlngs=_geo;L.MultiPolyline.prototype.setLatLngs.call(this,this._latlngs)},_generate_Geodesic:function(latlngs){var _geo=[],_geocnt=0,s,poly,points;for(poly=0;poly<latlngs.length;poly++){_geo[_geocnt]=[];for(points=0;points<latlngs[poly].length-1;points++){var inverse=this._vincenty_inverse(L.latLng(latlngs[poly][points]),L.latLng(latlngs[poly][points+1]));var prev=L.latLng(latlngs[poly][points]);_geo[_geocnt].push(prev);for(s=1;s<=this.options.steps;){var direct=this._vincenty_direct(L.latLng(latlngs[poly][points]),inverse.initialBearing,inverse.distance/this.options.steps*s,this.options.wrap);var gp=L.latLng(direct.lat,direct.lng);if(Math.abs(gp.lng-prev.lng)>180){var sec=this._intersection(L.latLng(latlngs[poly][points]),inverse.initialBearing,{lat:-89,lng:gp.lng-prev.lng>0?-INTERSECT_LNG:INTERSECT_LNG},0);if(sec){_geo[_geocnt].push(L.latLng(sec.lat,sec.lng));_geocnt++;_geo[_geocnt]=[];prev=L.latLng(sec.lat,-sec.lng);_geo[_geocnt].push(prev)}else{_geocnt++;_geo[_geocnt]=[];_geo[_geocnt].push(gp);prev=gp;s++}}else{_geo[_geocnt].push(gp);prev=gp;s++}}}_geocnt++}return _geo},_generate_GeodesicDashed:function(latlngs){var _geo=[],_geocnt=0,s,poly,points;for(poly=0;poly<latlngs.length;poly++){_geo[_geocnt]=[];for(points=0;points<latlngs[poly].length-1;points++){var inverse=this._vincenty_inverse(L.latLng(latlngs[poly][points]),L.latLng(latlngs[poly][points+1]));var prev=L.latLng(latlngs[poly][points]);_geo[_geocnt].push(prev);for(s=1;s<=this.options.steps;){var direct=this._vincenty_direct(L.latLng(latlngs[poly][points]),inverse.initialBearing,inverse.distance/this.options.steps*s-inverse.distance/this.options.steps*(1-this.options.dash),this.options.wrap);var gp=L.latLng(direct.lat,direct.lng);if(Math.abs(gp.lng-prev.lng)>180){var sec=this._intersection(L.latLng(latlngs[poly][points]),inverse.initialBearing,{lat:-89,lng:gp.lng-prev.lng>0?-INTERSECT_LNG:INTERSECT_LNG},0);if(sec){_geo[_geocnt].push(L.latLng(sec.lat,sec.lng));_geocnt++;_geo[_geocnt]=[];prev=L.latLng(sec.lat,-sec.lng);_geo[_geocnt].push(prev)}else{_geocnt++;_geo[_geocnt]=[];_geo[_geocnt].push(gp);prev=gp;s++}}else{_geo[_geocnt].push(gp);_geocnt++;var direct2=this._vincenty_direct(L.latLng(latlngs[poly][points]),inverse.initialBearing,inverse.distance/this.options.steps*s,this.options.wrap);_geo[_geocnt]=[];_geo[_geocnt].push(L.latLng(direct2.lat,direct2.lng));s++}}}_geocnt++}return _geo},_vincenty_direct:function(p1,initialBearing,distance,wrap){var φ1=p1.lat.toRadians(),λ1=p1.lng.toRadians();var α1=initialBearing.toRadians();var s=distance;var a=this.datum.ellipsoid.a,b=this.datum.ellipsoid.b,f=this.datum.ellipsoid.f;var sinα1=Math.sin(α1);var cosα1=Math.cos(α1);var tanU1=(1-f)*Math.tan(φ1),cosU1=1/Math.sqrt(1+tanU1*tanU1),sinU1=tanU1*cosU1;var σ1=Math.atan2(tanU1,cosα1);var sinα=cosU1*sinα1;var cosSqα=1-sinα*sinα;var uSq=cosSqα*(a*a-b*b)/(b*b);var A=1+uSq/16384*(4096+uSq*(-768+uSq*(320-175*uSq)));var B=uSq/1024*(256+uSq*(-128+uSq*(74-47*uSq)));var σ=s/(b*A),σʹ,iterations=0;do{var cos2σM=Math.cos(2*σ1+σ);var sinσ=Math.sin(σ);var cosσ=Math.cos(σ);var Δσ=B*sinσ*(cos2σM+B/4*(cosσ*(-1+2*cos2σM*cos2σM)-B/6*cos2σM*(-3+4*sinσ*sinσ)*(-3+4*cos2σM*cos2σM)));σʹ=σ;σ=s/(b*A)+Δσ}while(Math.abs(σ-σʹ)>1e-12&&++iterations);var x=sinU1*sinσ-cosU1*cosσ*cosα1;var φ2=Math.atan2(sinU1*cosσ+cosU1*sinσ*cosα1,(1-f)*Math.sqrt(sinα*sinα+x*x));var λ=Math.atan2(sinσ*sinα1,cosU1*cosσ-sinU1*sinσ*cosα1);var C=f/16*cosSqα*(4+f*(4-3*cosSqα));var L=λ-(1-C)*f*sinα*(σ+C*sinσ*(cos2σM+C*cosσ*(-1+2*cos2σM*cos2σM)));if(wrap)var λ2=(λ1+L+3*Math.PI)%(2*Math.PI)-Math.PI;else var λ2=λ1+L;var revAz=Math.atan2(sinα,-x);return{lat:φ2.toDegrees(),lng:λ2.toDegrees(),finalBearing:revAz.toDegrees()}},_vincenty_inverse:function(p1,p2){var φ1=p1.lat.toRadians(),λ1=p1.lng.toRadians();var φ2=p2.lat.toRadians(),λ2=p2.lng.toRadians();var a=this.datum.ellipsoid.a,b=this.datum.ellipsoid.b,f=this.datum.ellipsoid.f;var L=λ2-λ1;var tanU1=(1-f)*Math.tan(φ1),cosU1=1/Math.sqrt(1+tanU1*tanU1),sinU1=tanU1*cosU1;var tanU2=(1-f)*Math.tan(φ2),cosU2=1/Math.sqrt(1+tanU2*tanU2),sinU2=tanU2*cosU2;var λ=L,λʹ,iterations=0;do{var sinλ=Math.sin(λ),cosλ=Math.cos(λ);var sinSqσ=cosU2*sinλ*(cosU2*sinλ)+(cosU1*sinU2-sinU1*cosU2*cosλ)*(cosU1*sinU2-sinU1*cosU2*cosλ);var sinσ=Math.sqrt(sinSqσ);if(sinσ==0)return 0;var cosσ=sinU1*sinU2+cosU1*cosU2*cosλ;var σ=Math.atan2(sinσ,cosσ);var sinα=cosU1*cosU2*sinλ/sinσ;var cosSqα=1-sinα*sinα;var cos2σM=cosσ-2*sinU1*sinU2/cosSqα;if(isNaN(cos2σM))cos2σM=0;var C=f/16*cosSqα*(4+f*(4-3*cosSqα));λʹ=λ;λ=L+(1-C)*f*sinα*(σ+C*sinσ*(cos2σM+C*cosσ*(-1+2*cos2σM*cos2σM)))}while(Math.abs(λ-λʹ)>1e-12&&++iterations<100);if(iterations>=100){console.log("Formula failed to converge. Altering target position.");return this._vincenty_inverse(p1,{lat:p2.lat,lng:p2.lng-.01})}var uSq=cosSqα*(a*a-b*b)/(b*b);var A=1+uSq/16384*(4096+uSq*(-768+uSq*(320-175*uSq)));var B=uSq/1024*(256+uSq*(-128+uSq*(74-47*uSq)));var Δσ=B*sinσ*(cos2σM+B/4*(cosσ*(-1+2*cos2σM*cos2σM)-B/6*cos2σM*(-3+4*sinσ*sinσ)*(-3+4*cos2σM*cos2σM)));var s=b*A*(σ-Δσ);var fwdAz=Math.atan2(cosU2*sinλ,cosU1*sinU2-sinU1*cosU2*cosλ);var revAz=Math.atan2(cosU1*sinλ,-sinU1*cosU2+cosU1*sinU2*cosλ);s=Number(s.toFixed(3));return{distance:s,initialBearing:fwdAz.toDegrees(),finalBearing:revAz.toDegrees()}},_intersection:function(p1,brng1,p2,brng2){var φ1=p1.lat.toRadians(),λ1=p1.lng.toRadians();var φ2=p2.lat.toRadians(),λ2=p2.lng.toRadians();var θ13=Number(brng1).toRadians(),θ23=Number(brng2).toRadians();var Δφ=φ2-φ1,Δλ=λ2-λ1;var δ12=2*Math.asin(Math.sqrt(Math.sin(Δφ/2)*Math.sin(Δφ/2)+Math.cos(φ1)*Math.cos(φ2)*Math.sin(Δλ/2)*Math.sin(Δλ/2)));if(δ12==0)return null;var θ1=Math.acos((Math.sin(φ2)-Math.sin(φ1)*Math.cos(δ12))/(Math.sin(δ12)*Math.cos(φ1)));if(isNaN(θ1))θ1=0;var θ2=Math.acos((Math.sin(φ1)-Math.sin(φ2)*Math.cos(δ12))/(Math.sin(δ12)*Math.cos(φ2)));if(Math.sin(λ2-λ1)>0){var θ12=θ1;var θ21=2*Math.PI-θ2}else{var θ12=2*Math.PI-θ1;var θ21=θ2}var α1=(θ13-θ12+Math.PI)%(2*Math.PI)-Math.PI;var α2=(θ21-θ23+Math.PI)%(2*Math.PI)-Math.PI;if(Math.sin(α1)==0&&Math.sin(α2)==0)return null;if(Math.sin(α1)*Math.sin(α2)<0)return null;var α3=Math.acos(-Math.cos(α1)*Math.cos(α2)+Math.sin(α1)*Math.sin(α2)*Math.cos(δ12));var δ13=Math.atan2(Math.sin(δ12)*Math.sin(α1)*Math.sin(α2),Math.cos(α2)+Math.cos(α1)*Math.cos(α3));var φ3=Math.asin(Math.sin(φ1)*Math.cos(δ13)+Math.cos(φ1)*Math.sin(δ13)*Math.cos(θ13));var Δλ13=Math.atan2(Math.sin(θ13)*Math.sin(δ13)*Math.cos(φ1),Math.cos(δ13)-Math.sin(φ1)*Math.sin(φ3));var λ3=λ1+Δλ13;λ3=(λ3+3*Math.PI)%(2*Math.PI)-Math.PI;return{lat:φ3.toDegrees(),lng:λ3.toDegrees()}},_merge_options:function(obj1,obj2){var obj3={};for(var attrname in obj1){obj3[attrname]=obj1[attrname]}for(var attrname in obj2){obj3[attrname]=obj2[attrname]}return obj3}});L.geodesic=function(latlngs,options){return new L.Geodesic(latlngs,options)};(function(){var orig_L_GeoJSON_geometryToLayer=L.GeoJSON.geometryToLayer;L.GeoJSON.geometryToLayer=function(geojson,pointToLayer,coordsToLatLng,vectorOptions){if(geojson.properties&&geojson.properties.geodesic){var geometry=geojson.type==="Feature"?geojson.geometry:geojson,coords=geometry.coordinates,props=geojson.properties,latlngs;coordsToLatLng=coordsToLatLng||this.coordsToLatLng;if(props.geodesic_steps)vectorOptions=L.extend({steps:props.geodesic_steps},vectorOptions);if(props.geodesic_wrap)vectorOptions=L.extend({wrap:props.geodesic_wrap},vectorOptions);switch(geometry.type){case"LineString":latlngs=this.coordsToLatLngs(coords,0,coordsToLatLng);return new L.Geodesic([latlngs],vectorOptions);case"MultiLineString":latlngs=this.coordsToLatLngs(coords,1,coordsToLatLng);return new L.Geodesic(latlngs,vectorOptions);default:console.log("Not yet supported drawing GeoJSON "+geometry.type+" as a geodesic: Drawing as non-geodesic.")}}return orig_L_GeoJSON_geometryToLayer.apply(this,arguments)}})();