Define the name that you will use to specify your custom basemap.

name = 'openstreetmap';

Specify the website that provides the map data. The first character of the URL indicates which server to use to get the data. For load balancing, the provider has three servers that you can use: a, b, or c.

url = 'a.tile.openstreetmap.org';

Create an attribution to display on the map that gives credit to the provider of the map data. Web map providers might define specific requirements for the attribution.

copyright = char(uint8(169));
attribution = copyright + "OpenStreetMap contributors";

Add the custom basemap to the list of basemap layers available.

addCustomBasemap(name,url,'Attribution',attribution)

Plot the locations over the map using a geographic bubble chart. You can specify your custom basemap when you create the geographic bubble chart.

pts = gpxread('boston_placenames');
gb = geobubble(pts.Latitude,pts.Longitude,'Basemap','openstreetmap');
gb.BubbleWidthRange = 25;
gb.MapLayout = 'maximized';
gb.ZoomLevel = 14;

Before adding the basemap, specify the location of the map tiles. To do this, specify the URL of the National Map ArcGIS REST Services Directory. Then, add the path to the map tiles from the USGS Topo basemap service.

url = "https://basemap.nationalmap.gov/ArcGIS/rest/services";
fullurl = url + "/USGSTopo/MapServer/tile/${z}/${y}/${x}";

Specify a name for the basemap and attribution text to display with it.

nm = 'usgstopo';
att = 'Credit: US Geological Survey';

Finally, add the USGS Topo basemap.

addCustomBasemap(nm,fullurl,'Attribution',att)

Plot the path of a glider over the basemap. To do this, import latitude, longitude, and geoid height values.

trk = gpxread('sample_mixed','FeatureType','track');
lat = trk.Latitude;
lon = trk.Longitude;
h = trk.Elevation;

Display the path in 2-D using geographic axes. Specify the basemap using the geobasemap function and the name of the basemap. Use the basemap name given when you created it. Call hold on before plotting the line to prevent the basemap from resetting.

geoplot(lat,lon)
geobasemap('usgstopo')
hold on
geoplot(lat,lon,'r')

Display the path in 3-D using a geographic globe. Specify the basemap using the 'Basemap' name-value pair argument. By default, the view is directly above the path. Tilt the view by holding Ctrl and dragging.

uif = uifigure;
g = geoglobe(uif,'Basemap','usgstopo');
hold(g,'on')
geoplot3(g,lat,lon,h,'r')

Define the name that you will use to specify your custom basemap.

name = 'opentopomap';

Specify the website that provides the map data. The first character of the URL indicates which server to use to get the data. For load balancing, the provider has three servers that you can use: a, b, or c.

url = 'a.tile.opentopomap.org';

Create an attribution to display on the map that gives credit to the provider of the map data. Web map providers might define specific requirements for the attribution.

copyright = char(uint8(169));
attribution = [ ...
      "map data:  " + copyright + "OpenStreetMap contributors,SRTM", ...
      "map style: " + copyright + "OpenTopoMap (CC-BY-SA)"];

Define the name that will appear in the Layer Manager to represent your custom basemap.

displayName = 'Open Topo Map';

Add the custom basemap to the list of available basemap layers.

addCustomBasemap(name,url,'Attribution',attribution, ...
                          'DisplayName',displayName)

Open a web map. Specify the custom basemap using the name you defined when you added it. Then, read glider data into the workspace and plot it over the basemap.

webmap opentopomap
trk = gpxread('sample_mixed.gpx','FeatureType','track');
wmline(trk,'LineWidth',2)

View glider path on basemaps from the USGS National Map. Construct basemap URLs by replacing BASEMAP with the name of the USGS basemap.

Read in the glider path track data.

trk = gpxread('sample_mixed','FeatureType','track');

Specify the custom basemap URL. The USGS National Map supports several tiled web maps. For this example, insert the word "BASEMAP" into the URL string. In this way, you can replace the word BASEMAP with the name of one of the maps supported by the USGS National Map.

baseURL = "https://basemap.nationalmap.gov/ArcGIS/rest/services";
usgsURL = baseURL + "/BASEMAP/MapServer/tile/${z}/${y}/${x}";

Specify a list of the names of USGS basemaps that you want to use. These names will be inserted into the URL in place of "BASEMAP".

basemaps = ["USGSImageryOnly" "USGSImageryTopo" ...
       "USGSTopo" "USGSShadedReliefOnly" "USGSHydroCached"];

Specify a list of display names that you can use with each map. Use display names that are the same as those used by webmap so that webmap does not contain duplicate maps.

displayNames = ["USGS Imagery" "USGS Topographic Imagery" ...
       "USGS Shaded Topographic Map" "USGS Shaded Relief" ...
       "USGS Hydrography"];
maxZoomLevel = 16;

Create a map attribution to give credit to the provider of the map data.

attribution = 'Credit: U.S. Geological Survey';

Create a loop in which you use each basemap with the geoplot function, plotting the glider data on each of the USGS basemaps.

for k =1:length(basemaps)
   basemap = basemaps(k);
   name = lower(basemap);
   url = replace(usgsURL,"BASEMAP",basemap);
   displayName = displayNames(k);
   addCustomBasemap(name,url,'Attribution',attribution, ...
         'DisplayName',displayName,'MaxZoomLevel',maxZoomLevel)
   figure
   geoplot(trk.Latitude,trk.Longitude,'r','LineWidth',5);
   geobasemap(basemap)
   title(displayName)
end