Planetary surveys

With ipyaladin, you can also display planetary surfaces.

from ipyaladin import Aladin
from astropy.table import Table

First, let’s create an Aladin Lite widget with the Mars Viking MDIM21 survey.

mars = Aladin(
    target="159.2135528 -58.6241989",
    survey="https://alasky.u-strasbg.fr/Planets/Mars_Viking_MDIM21",
    fov=10,
)
mars

The target does not return a SkyCoord object anymore, since we don’t represent the sky here. It is a couple of Longitude and Latitude.

mars.target

The WCS are also defined on planetary surfaces (in Marmo et.al. 2018). This is why you see a CTYPE starting with MA for Mars.

mars.wcs

Let’s change the center of the view to Olympus Mons.

mars.target = "Olympus Mons"

Any name recognized by the Gazetter of Planetary Nomenclature will work.

We can also add tables:

longitudes = [
    226.2,
    70.5,
    250.4,
    -59.2,
    147.21,
    316.0,
    32.53,
    -112.58,
    298.0,
    30.0,
    70.5,
    280.0,
    87.0,
    117.5,
    350.0,
]
latitudes = [
    18.65,
    -42.4,
    40.5,
    -13.9,
    25.02,
    -49.7,
    70.0,
    1.57,
    25.0,
    19.79,
    -42.4,
    45.0,
    12.9,
    46.7,
    -45.0,
]
names = [
    "Olympus Mons",
    "Hellas Planitia",
    "Alba Mons",
    "Valles Marineris",
    "Elysium Mons",
    "Argyre Basin",
    "Vastitas Borealis",
    "Tharsis Montes",
    "Outflow channels",
    "Arabia Terra",
    "Hellas Basin",
    "Tempe Terra",
    "Isidis Basin",
    "Utopia Basin",
    "Noachis Terra",
]

table = Table([longitudes, latitudes, names], names=("Longitude", "Latitude", "Name"))
mars.add_table(
    table, color="#67d38d", source_size=15, shape="cross", name="Mars_features"
)

Unzoom and look for the other points.

Surveys of Earth

earth = Aladin(
    survey="CDS/P/Earth/DEM/elevation",
    fov=100,
)
earth

earth.target

earth.wcs

earth.target = "Strasbourg"  # ipyaladin's home

Any name recognized by OpenStreetMaps will work.