.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/geodesic.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_auto_examples_geodesic.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_auto_examples_geodesic.py:


Geodesic Calculation
====================

Example of using the geodesic module and comparision
between Haversine and Vincenty implementation

.. GENERATED FROM PYTHON SOURCE LINES 9-154




.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Haversine distance: 2670.135807683966
    Vincenty distance:  2674.0325533379682

    Haversine initial bearing: 43.44379502125787
    Vincenty initial bearing:  43.5279748443246

    Haversine destination point: (48.912637935409954, 1.7753897932233385)
    Vincenty destination point:  (48.912613664862455, 1.7771301984818817)

    Haversine destination point: (48.912637935409954, 1.7753897932233385)
    Vincenty destination point:  (48.912613664862455, 1.7771301984818817)


.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <div style="width:100%;"><div style="position:relative;width:100%;height:0;padding-bottom:60%;"><span style="color:#565656">Make this Notebook Trusted to load map: File -> Trust Notebook</span><iframe srcdoc="&lt;!DOCTYPE html&gt;
    &lt;html&gt;
    &lt;head&gt;
    
        &lt;meta http-equiv=&quot;content-type&quot; content=&quot;text/html; charset=UTF-8&quot; /&gt;
    
            &lt;script&gt;
                L_NO_TOUCH = false;
                L_DISABLE_3D = false;
            &lt;/script&gt;
    
        &lt;style&gt;html, body {width: 100%;height: 100%;margin: 0;padding: 0;}&lt;/style&gt;
        &lt;style&gt;#map {position:absolute;top:0;bottom:0;right:0;left:0;}&lt;/style&gt;
        &lt;script src=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.js&quot;&gt;&lt;/script&gt;
        &lt;script src=&quot;https://code.jquery.com/jquery-3.7.1.min.js&quot;&gt;&lt;/script&gt;
        &lt;script src=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/js/bootstrap.bundle.min.js&quot;&gt;&lt;/script&gt;
        &lt;script src=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.js&quot;&gt;&lt;/script&gt;
        &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/leaflet@1.9.3/dist/leaflet.css&quot;/&gt;
        &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/bootstrap@5.2.2/dist/css/bootstrap.min.css&quot;/&gt;
        &lt;link rel=&quot;stylesheet&quot; href=&quot;https://netdna.bootstrapcdn.com/bootstrap/3.0.0/css/bootstrap-glyphicons.css&quot;/&gt;
        &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6.2.0/css/all.min.css&quot;/&gt;
        &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdnjs.cloudflare.com/ajax/libs/Leaflet.awesome-markers/2.0.2/leaflet.awesome-markers.css&quot;/&gt;
        &lt;link rel=&quot;stylesheet&quot; href=&quot;https://cdn.jsdelivr.net/gh/python-visualization/folium/folium/templates/leaflet.awesome.rotate.min.css&quot;/&gt;
    
                &lt;meta name=&quot;viewport&quot; content=&quot;width=device-width,
                    initial-scale=1.0, maximum-scale=1.0, user-scalable=no&quot; /&gt;
                &lt;style&gt;
                    #map_0d318ca7a1ec94c542a91ca2dff96030 {
                        position: relative;
                        width: 100.0%;
                        height: 100.0%;
                        left: 0.0%;
                        top: 0.0%;
                    }
                    .leaflet-container { font-size: 1rem; }
                &lt;/style&gt;
        
    &lt;/head&gt;
    &lt;body&gt;
    
    
                &lt;div class=&quot;folium-map&quot; id=&quot;map_0d318ca7a1ec94c542a91ca2dff96030&quot; &gt;&lt;/div&gt;
        
    &lt;/body&gt;
    &lt;script&gt;
    
    
                var map_0d318ca7a1ec94c542a91ca2dff96030 = L.map(
                    &quot;map_0d318ca7a1ec94c542a91ca2dff96030&quot;,
                    {
                        center: [48.596289, 2.35164],
                        crs: L.CRS.EPSG3857,
                        ...{
      &quot;zoom&quot;: 10,
      &quot;zoomControl&quot;: true,
      &quot;preferCanvas&quot;: false,
    }

                    }
                );

            

        
    
                var tile_layer_ff6fa41179bb85ae5507d2b2f3320df2 = L.tileLayer(
                    &quot;https://tile.openstreetmap.org/{z}/{x}/{y}.png&quot;,
                    {
      &quot;minZoom&quot;: 0,
      &quot;maxZoom&quot;: 19,
      &quot;maxNativeZoom&quot;: 19,
      &quot;noWrap&quot;: false,
      &quot;attribution&quot;: &quot;\u0026copy; \u003ca href=\&quot;https://www.openstreetmap.org/copyright\&quot;\u003eOpenStreetMap\u003c/a\u003e contributors&quot;,
      &quot;subdomains&quot;: &quot;abc&quot;,
      &quot;detectRetina&quot;: false,
      &quot;tms&quot;: false,
      &quot;opacity&quot;: 1,
    }

                );
        
    
                tile_layer_ff6fa41179bb85ae5507d2b2f3320df2.addTo(map_0d318ca7a1ec94c542a91ca2dff96030);
        
    
                var poly_line_7725ff6a7578b4ccbddbe1e84acbfce1 = L.polyline(
                    [[48.596289, 2.35164], [48.912637935409954, 1.7753897932233385]],
                    {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;red&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: false, &quot;fillColor&quot;: &quot;red&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;noClip&quot;: false, &quot;opacity&quot;: 1.0, &quot;smoothFactor&quot;: 1.0, &quot;stroke&quot;: true, &quot;weight&quot;: 2}
                ).addTo(map_0d318ca7a1ec94c542a91ca2dff96030);
        
    
                var poly_line_f7d20cf790c2778a42f2e90d784a6de2 = L.polyline(
                    [[48.596289, 2.35164], [48.912613664862455, 1.7771301984818817]],
                    {&quot;bubblingMouseEvents&quot;: true, &quot;color&quot;: &quot;blue&quot;, &quot;dashArray&quot;: null, &quot;dashOffset&quot;: null, &quot;fill&quot;: false, &quot;fillColor&quot;: &quot;blue&quot;, &quot;fillOpacity&quot;: 0.2, &quot;fillRule&quot;: &quot;evenodd&quot;, &quot;lineCap&quot;: &quot;round&quot;, &quot;lineJoin&quot;: &quot;round&quot;, &quot;noClip&quot;: false, &quot;opacity&quot;: 1.0, &quot;smoothFactor&quot;: 1.0, &quot;stroke&quot;: true, &quot;weight&quot;: 2}
                ).addTo(map_0d318ca7a1ec94c542a91ca2dff96030);
        
    
                var marker_489b50604cbb3a04ee9609fd1ba9136c = L.marker(
                    [48.596289, 2.35164],
                    {
    }
                ).addTo(map_0d318ca7a1ec94c542a91ca2dff96030);
        
    
                var marker_2f5b807af93f0c95bcdfbac513263a56 = L.marker(
                    [48.912637935409954, 1.7753897932233385],
                    {
    }
                ).addTo(map_0d318ca7a1ec94c542a91ca2dff96030);
        
    
                var icon_c436d18eb22351f97bb74c289c1f50e9 = L.AwesomeMarkers.icon(
                    {
      &quot;markerColor&quot;: &quot;red&quot;,
      &quot;iconColor&quot;: &quot;white&quot;,
      &quot;icon&quot;: &quot;info-sign&quot;,
      &quot;prefix&quot;: &quot;glyphicon&quot;,
      &quot;extraClasses&quot;: &quot;fa-rotate-0&quot;,
    }
                );
        
    
            var popup_ccbcce920f8940edd0189140fe056d24 = L.popup({
      &quot;maxWidth&quot;: &quot;100%&quot;,
    });

        
            
                    var html_f26b1284a38f67aad5242f3585026a99 = $(`&lt;div id=&quot;html_f26b1284a38f67aad5242f3585026a99&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;&lt;i&gt;Haversine&lt;/i&gt;&lt;/div&gt;`)[0];
                    popup_ccbcce920f8940edd0189140fe056d24.setContent(html_f26b1284a38f67aad5242f3585026a99);
            
        

            marker_2f5b807af93f0c95bcdfbac513263a56.bindPopup(popup_ccbcce920f8940edd0189140fe056d24)
            ;

        
    
    
                marker_2f5b807af93f0c95bcdfbac513263a56.bindTooltip(
                    `&lt;div&gt;
                         Haversine
                     &lt;/div&gt;`,
                    {
      &quot;sticky&quot;: true,
    }
                );
        
    
                    marker_2f5b807af93f0c95bcdfbac513263a56.setIcon(icon_c436d18eb22351f97bb74c289c1f50e9);
            
    
                var marker_185cda25470a5a092872a1dac0ff9107 = L.marker(
                    [48.912613664862455, 1.7771301984818817],
                    {
    }
                ).addTo(map_0d318ca7a1ec94c542a91ca2dff96030);
        
    
                var icon_00e4f7efb64b37da0c4e8de27aba31dc = L.AwesomeMarkers.icon(
                    {
      &quot;markerColor&quot;: &quot;blue&quot;,
      &quot;iconColor&quot;: &quot;white&quot;,
      &quot;icon&quot;: &quot;info-sign&quot;,
      &quot;prefix&quot;: &quot;glyphicon&quot;,
      &quot;extraClasses&quot;: &quot;fa-rotate-0&quot;,
    }
                );
        
    
            var popup_8b363eb52d92484bece33185fd13df04 = L.popup({
      &quot;maxWidth&quot;: &quot;100%&quot;,
    });

        
            
                    var html_5dd593e8f9bf07fc353b1b4f99457614 = $(`&lt;div id=&quot;html_5dd593e8f9bf07fc353b1b4f99457614&quot; style=&quot;width: 100.0%; height: 100.0%;&quot;&gt;&lt;i&gt;Vincenty&lt;/i&gt;&lt;/div&gt;`)[0];
                    popup_8b363eb52d92484bece33185fd13df04.setContent(html_5dd593e8f9bf07fc353b1b4f99457614);
            
        

            marker_185cda25470a5a092872a1dac0ff9107.bindPopup(popup_8b363eb52d92484bece33185fd13df04)
            ;

        
    
    
                marker_185cda25470a5a092872a1dac0ff9107.bindTooltip(
                    `&lt;div&gt;
                         Vincenty
                     &lt;/div&gt;`,
                    {
      &quot;sticky&quot;: true,
    }
                );
        
    
                    marker_185cda25470a5a092872a1dac0ff9107.setIcon(icon_00e4f7efb64b37da0c4e8de27aba31dc);
            
    &lt;/script&gt;
    &lt;/html&gt;" style="position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>
    </div>
    <br />
    <br />



|

.. code-block:: Python


    import folium

    from pyBADA import geodesic as geo


    # input parameters
    LAT_init = 48.596289
    LON_init = 2.351640
    LAT_final = 48.613721
    LON_final = 2.376616

    distance = 54972.271  # [m]
    bearing = 310

    # ===========================================================
    # calculation of distance between 2 latitude/longitude points
    # ===========================================================

    # Haversine formula
    distance_haversine = geo.Haversine.distance(
        LAT_init=LAT_init,
        LON_init=LON_init,
        LAT_final=LAT_final,
        LON_final=LON_final,
    )
    print("Haversine distance:", distance_haversine)

    # Vincenty formula
    distance_vincenty = geo.Vincenty.distance(
        LAT_init=LAT_init,
        LON_init=LON_init,
        LAT_final=LAT_final,
        LON_final=LON_final,
    )
    print("Vincenty distance: ", distance_vincenty)
    print("")

    # ===========================================================
    # calculation of bearing between 2 latitude/longitude points
    # ===========================================================

    # Haversine formula
    bearing_haversine = geo.Haversine.bearing(
        LAT_init=LAT_init,
        LON_init=LON_init,
        LAT_final=LAT_final,
        LON_final=LON_final,
    )
    print("Haversine initial bearing:", bearing_haversine)

    # Vincenty formula for initial bearing
    bearing_vincenty = geo.Vincenty.bearing_initial(
        LAT_init=LAT_init,
        LON_init=LON_init,
        LAT_final=LAT_final,
        LON_final=LON_final,
    )
    print("Vincenty initial bearing: ", bearing_vincenty)

    # Vincenty formula for final bearing
    bearing_vincenty = geo.Vincenty.bearing_final(
        LAT_init=LAT_init,
        LON_init=LON_init,
        LAT_final=LAT_final,
        LON_final=LON_final,
    )
    # print("Vincenty final bearing:   ", bearing_vincenty)
    print("")

    # =====================================================================================
    # calculation of destination point based in initial point, distance and initial bearing
    # =====================================================================================

    # Haversine formula
    destPoint_haversine = geo.Haversine.destinationPoint(
        LAT_init=LAT_init, LON_init=LON_init, distance=distance, bearing=bearing
    )
    print("Haversine destination point:", destPoint_haversine)

    # Vincenty formula
    destPoint_vincenty = geo.Vincenty.destinationPoint(
        LAT_init=LAT_init, LON_init=LON_init, distance=distance, bearing=bearing
    )
    print("Vincenty destination point: ", destPoint_vincenty)
    print("")

    # =====================================================================================
    # calculation of destination point based in initial point, distance and initial bearing
    # =====================================================================================

    # Haversine formula
    destPoint_haversine = geo.Haversine.destinationPoint(
        LAT_init=LAT_init, LON_init=LON_init, distance=distance, bearing=bearing
    )
    print("Haversine destination point:", destPoint_haversine)

    # Vincenty formula
    destPoint_vincenty = geo.Vincenty.destinationPoint(
        LAT_init=LAT_init, LON_init=LON_init, distance=distance, bearing=bearing
    )
    print("Vincenty destination point: ", destPoint_vincenty)

    # =====================================================================================
    # drawing destination point on the 2D map
    # =====================================================================================

    initPoint = (LAT_init, LON_init)

    # create a base map
    myMap = folium.Map(location=initPoint, zoom_start=10)

    # construct the lines
    haversineLine = []
    haversineLine.append(initPoint)
    haversineLine.append(destPoint_haversine)

    vincentyLine = []
    vincentyLine.append(initPoint)
    vincentyLine.append(destPoint_vincenty)

    # draw the lines on the map
    folium.PolyLine(haversineLine, color="red", weight=2).add_to(myMap)
    folium.PolyLine(vincentyLine, color="blue", weight=2).add_to(myMap)

    # add markers on the map
    folium.Marker(initPoint).add_to(myMap)
    folium.Marker(
        destPoint_haversine,
        popup="<i>Haversine</i>",
        tooltip="Haversine",
        icon=folium.Icon(color="red"),
    ).add_to(myMap)
    folium.Marker(
        destPoint_vincenty,
        popup="<i>Vincenty</i>",
        tooltip="Vincenty",
        icon=folium.Icon(color="blue"),
    ).add_to(myMap)

    # saving the map in the HTML file. This map then can be viewed in the browser and you can interact with it
    # myMap.save("map_geodesic.html")

    # display the map
    myMap


.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 0.168 seconds)


.. _sphx_glr_download_auto_examples_geodesic.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: geodesic.ipynb <geodesic.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: geodesic.py <geodesic.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: geodesic.zip <geodesic.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_