vmap0 = wmsfind('vmap0.tiles','SearchField','serverurl');
vmap0 = wmsupdate(vmap0);
layer = refine(vmap0,'basic');
request = WMSMapRequest(layer);
request.Transparent = true;
imageLength = 720;
request = boundImageSize(request,imageLength);
globalImage = getMap(request.Server, request.RequestURL);

figure
worldmap('world')
geoshow(globalImage,request.RasterReference);
title(['VMAP0 ' layer.LayerTitle ' Layer'])

vmap0 = wmsfind('vmap0.tiles','SearchField','serverurl');
vmap0 = wmsupdate(vmap0);
layers = [ refine(vmap0, 'rail'); refine(vmap0, 'river'); ...
       refine(vmap0, 'priroad'); refine(vmap0, 'secroad'); ...
       refine(vmap0, 'ctylabel'); refine(vmap0, 'basic')];
request = WMSMapRequest(layers);

cities = readgeotable('worldcities.shp');
row = cities.Name == 'London';
london = cities(row,:);

extent = [-.25 .25];
request.Latlim = london.Shape.Latitude + extent;
request.Lonlim = london.Shape.Longitude + extent;
 
request.Transparent = true;
imageLength = 1024;
request = boundImageSize(request,imageLength);
londonImage = getMap(request.Server,request.RequestURL);

figure
worldmap(londonImage,request.RasterReference)
geoshow(londonImage,request.RasterReference)
title({'Region Surrounding London, England', ...
   ['with Primary and Secondary Roads, ', ...
    'Rivers, Rails, City Label, and Basic Layers']})