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 = shaperead('worldcities', 'UseGeo', true);
london = cities(strcmpi('London', {cities.Name}));
extent = [-.25 .25];
request.Latlim = london.Lat + extent;
request.Lonlim = london.Lon + 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']})