Basemaps

Basemaps#

Here we will show the three options we have for basemaps. These can either be created as backgrounds for plotting your gridded / point data on top of, or used by themselves.

Either set the parameters simple_basemap, imagery_basemap,or modis_basemap to True in either function maps.plot_grd() or maps.basemap().

Alternatives, use the equivalent functions on existing figures: maps.add_simple_basemap(), maps.add_imagery(), or maps.add_modis().

Import the modules

[1]:

%load_ext autoreload
%autoreload 2

import os

from polartoolkit import fetch, maps, regions

[2]:

# set default to southern hemisphere for this notebook
os.environ["POLARTOOLKIT_HEMISPHERE"] = "south"

Basemap of grounded vs floating ice#

[3]:

# by default uses Depoorter et al. 2013 shapefiles
fig = maps.basemap(
    simple_basemap=True,
    region=regions.getz_ice_shelf,
)

fig.show(dpi=200)
../_images/how_to_basemaps_5_0.png 
[4]:

# change the version of the simple basemap to use MEaSUREs v2 shapefiles
# specify different colors
fig = maps.basemap(
    simple_basemap=True,
    simple_basemap_version="depoorter-2013",
    simple_basemap_grounded_color="purple",
    simple_basemap_floating_color="orange",
    simple_basemap_pen="4p,red",
    region=regions.getz_ice_shelf,
)

fig.show(dpi=200)
../_images/how_to_basemaps_6_0.png 
[5]:

# equivalent for the northern hemisphere only shows grounded ice
fig = maps.basemap(
    simple_basemap=True,
    hemisphere="north",
    region=regions.north_greenland,
)

fig.show(dpi=200)
../_images/how_to_basemaps_7_0.png

Add to existing figure#

This will cover any portions of the existing figure for grounded or floating ice. You can supply options kwarg transparency to make the simple basemap transparent.

[6]:

# create figure
fig = maps.plot_grd(
    fetch.etopo(hemisphere="north"),
    region=regions.north_greenland,
    cbar_label="Elevation (m)",
    hemisphere="north",
)

# add the basemap
fig.add_simple_basemap(
    grounded_color="red",
    transparency=60,
)
fig.show(dpi=200)
../_images/how_to_basemaps_9_0.png

Satellite imagery basemap#

[7]:

fig = maps.basemap(
    imagery_basemap=True,
    region=regions.getz_ice_shelf,
)

fig.show(dpi=200)
../_images/how_to_basemaps_11_0.png

Add to existing figure#

[8]:

# create figure
fig = maps.plot_grd(
    fetch.etopo(),
    region=regions.getz_ice_shelf,
    cbar_label="Elevation (m)",
)

# add the LIMA imagery, but make it semi-transparent to be able to see grid values
# behind
fig.add_imagery(
    transparency=50,
)
fig.show(dpi=200)
../_images/how_to_basemaps_13_0.png

MODIS imagery basemap#

[9]:

# by default for Antarctica is 750 m MODIS resolution
fig = maps.basemap(
    modis_basemap=True,
    region=regions.getz_ice_shelf,
    # modis_version="125m", # uncomment this to use 125 m resolution MODIS imagery
)
fig.show(dpi=200)
../_images/how_to_basemaps_15_0.png 
[10]:

# equivalent for the northern hemisphere, by default uses 500 m MODIS resolution
fig = maps.basemap(
    modis_basemap=True,
    hemisphere="north",
    region=regions.north_greenland,
)

fig.show(dpi=200)
../_images/how_to_basemaps_16_0.png

Add to existing figure#

[11]:

# create figure
fig = maps.plot_grd(
    fetch.etopo(),
    region=regions.getz_ice_shelf,
    cbar_label="Elevation (m)",
)

# add the MODIS imagery, but make it semi-transparent to be able to see grid values
# behind
fig.add_modis(
    transparency=50,
)
fig.show(dpi=200)
../_images/how_to_basemaps_18_0.png 
[ ]:
Contents