.. note::
:class: sphx-glr-download-link-note
Click :ref:`here <sphx_glr_download_auto_tutorials_01_introductory_viz_earth_coordinates.py>` to download the full example code
.. rst-class:: sphx-glr-example-title
.. _sphx_glr_auto_tutorials_01_introductory_viz_earth_coordinates.py:
============================
Earth Coordinate Conversion
============================
In this tutorial, we will show how to place actors on specific locations
on the surface of the Earth using a new function.
.. code-block:: default
from fury import window, actor, utils, io
from fury.data import read_viz_textures, fetch_viz_textures
import math
import numpy as np
import itertools
Create a new scene, and load in the image of the Earth using
``fetch_viz_textures`` and ``read_viz_textures``. We will use a 16k
resolution texture for maximum detail.
.. code-block:: default
scene = window.Scene()
fetch_viz_textures()
earth_file = read_viz_textures("1_earth_16k.jpg")
earth_image = io.load_image(earth_file)
earth_actor = actor.texture_on_sphere(earth_image)
scene.add(earth_actor)
.. rst-class:: sphx-glr-script-out
Out:
.. code-block:: none
Dataset is already in place. If you want to fetch it again please first remove the folder /Users/koudoro/.fury/textures
/Users/koudoro/miniconda3/envs/fury-env-37/lib/python3.7/site-packages/PIL/Image.py:2850: DecompressionBombWarning: Image size (131220000 pixels) exceeds limit of 89478485 pixels, could be decompression bomb DOS attack.
DecompressionBombWarning,
Rotate the Earth to make sure the texture is correctly oriented. Change it's
scale using ``actor.SetScale()``.
.. code-block:: default
utils.rotate(earth_actor, (-90, 1, 0, 0))
utils.rotate(earth_actor, (180, 0, 1, 0))
earth_actor.SetScale(2, 2, 2)
Define the function to convert geographical coordinates of a location in
latitude and longitude degrees to coordinates on the ``earth_actor`` surface.
In this function, convert to radians, then to spherical coordinates, and
lastly, to cartesian coordinates.
.. code-block:: default
def latlong_coordinates(lat, lon):
# Convert latitude and longitude to spherical coordinates
degrees_to_radians = math.pi/180.0
# phi = 90 - latitude
phi = (90-lat)*degrees_to_radians
# theta = longitude
theta = lon*degrees_to_radians*-1
# now convert to cartesian
x = np.sin(phi)*np.cos(theta)
y = np.sin(phi)*np.sin(theta)
z = np.cos(phi)
# flipping z to y for FURY coordinates
return (x, z, y)
Use this new function to place some sphere actors on several big cities
around the Earth.
.. code-block:: default
locationone = latlong_coordinates(40.730610, -73.935242) # new york city, us
locationtwo = latlong_coordinates(39.916668, 116.383331) # beijing, china
locationthree = latlong_coordinates(48.864716, 2.349014) # paris, france
Set the centers, radii, and colors of these spheres, and create a new
``sphere_actor`` for each location to add to the scene.
.. code-block:: default
centers = np.array([[*locationone], [*locationtwo], [*locationthree]])
colors = np.random.rand(3, 3)
radii = np.array([0.005, 0.005, 0.005])
sphere_actor = actor.sphere(centers, colors, radii)
scene.add(sphere_actor)
Create some text actors to add to the scene indicating each location and its
geographical coordinates.
.. code-block:: default
nyc_actor = actor.text_3d("New York City, New York\n40.7128° N, 74.0060° W",
(locationone[0]-0.04, locationone[1],
locationone[2]+0.07),
window.colors.white, 0.01)
paris_actor = actor.text_3d("Paris, France\n48.8566° N, 2.3522° E",
(locationthree[0]-0.04, locationthree[1],
locationthree[2]-0.07),
window.colors.white, 0.01)
beijing_actor = actor.text_3d("Beijing, China\n39.9042° N, 116.4074° E",
(locationtwo[0]-0.06, locationtwo[1],
locationtwo[2]-0.07),
window.colors.white, 0.01)
utils.rotate(paris_actor, (85, 0, 1, 0))
utils.rotate(beijing_actor, (180, 0, 1, 0))
utils.rotate(nyc_actor, (5, 1, 0, 0))
Create a ShowManager object, which acts as the interface between the scene,
the window and the interactor.
.. code-block:: default
showm = window.ShowManager(scene,
size=(900, 768), reset_camera=False,
order_transparent=True)
Let's create a ``timer_callback`` function to add some animation to the
Earth. Change the camera position and angle to fly over and zoom in on each
new location.
.. code-block:: default
counter = itertools.count()
def timer_callback(_obj, _event):
cnt = next(counter)
showm.render()
if cnt == 0:
scene.set_camera(position=(1.5, 3.5, 7.0))
if cnt < 200 and cnt > 25:
scene.zoom(1.015)
scene.pitch(0.01)
if cnt == 200:
scene.add(nyc_actor)
if cnt > 250 and cnt < 350:
scene.zoom(0.985)
if cnt > 350 and cnt < 425:
scene.azimuth(1)
if cnt > 425 and cnt < 525:
scene.zoom(1.015)
scene.pitch(0.011)
if cnt == 525:
scene.add(paris_actor)
if cnt > 575 and cnt < 700:
scene.zoom(0.985)
if cnt > 700 and cnt < 820:
scene.azimuth(1)
if cnt > 820 and cnt < 930:
scene.zoom(1.015)
if cnt == 930:
scene.add(beijing_actor)
if cnt == 1000:
showm.exit()
Initialize the ShowManager object, add the timer_callback, and watch the
new animation take place!
.. code-block:: default
showm.initialize()
showm.add_timer_callback(True, 25, timer_callback)
showm.start()
window.record(showm.scene, size=(900, 768),
out_path="viz_earth_coordinates.png")
.. image:: /auto_tutorials/01_introductory/images/sphx_glr_viz_earth_coordinates_001.png
:class: sphx-glr-single-img
.. rst-class:: sphx-glr-timing
**Total running time of the script:** ( 0 minutes 32.530 seconds)
.. _sphx_glr_download_auto_tutorials_01_introductory_viz_earth_coordinates.py:
.. only :: html
.. container:: sphx-glr-footer
:class: sphx-glr-footer-example
.. container:: sphx-glr-download
:download:`Download Python source code: viz_earth_coordinates.py <viz_earth_coordinates.py>`
.. container:: sphx-glr-download
:download:`Download Jupyter notebook: viz_earth_coordinates.ipynb <viz_earth_coordinates.ipynb>`
.. only:: html
.. rst-class:: sphx-glr-signature
`Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_