In [13]:

    

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
%matplotlib inline


    

In [14]:

    

def make_map(projection=ccrs.PlateCarree(), extent=[-42, 0, -32, 0.5]):
    subplot_kw = dict(projection=projection)
    fig, ax = plt.subplots(subplot_kw=subplot_kw)
    ax.set_extent(extent)
    gl = ax.gridlines(draw_labels=True)
    gl.xlabels_top = gl.ylabels_right = False
    gl.xformatter = LONGITUDE_FORMATTER
    gl.yformatter = LATITUDE_FORMATTER
    return fig, ax


    

In [15]:

    

from owslib.wms import WebMapService
url = "http://nowcoast.noaa.gov/wms/com.esri.wms.Esrimap/analyses"
wms = WebMapService(url)

layer = 'NCEP_RAS_ANAL_RTG_SST'
fig, ax = make_map(projection=ccrs.PlateCarree())
ax.add_wms(wms, layer)
_ = ax.set_title(layer)


    

    






<matplotlib.figure.Figure at 0x7f64e821c9d0>

In [16]:

    

url = "http://nowcoast.noaa.gov/wms/com.esri.wms.Esrimap/analyses"
layer = 'NCEP_RAS_ANAL_RTG_SST'

ax = plt.axes(projection=ccrs.PlateCarree())
ax.add_wms(url, layer)
ax.set_extent((-15, 25, 35, 60))


    

    






<matplotlib.figure.Figure at 0x7f64e013b050>

In [17]:

    

url = 'http://map1c.vis.earthdata.nasa.gov/wmts-geo/wmts.cgi'
layer = 'VIIRS_CityLights_2012'

ax = plt.axes(projection=ccrs.PlateCarree())
ax.add_wmts(url, layer)
ax.set_extent((-15, 25, 35, 60))

plt.title('Suomi NPP Earth at night April/October 2012')
plt.show()


    

In [18]:

    

url='http://coast.noaa.gov/arcgis/services/CCAP/All_Change_0106/MapServer/WMSServer?request=GetCapabilities&service=WMS'
wms = WebMapService(url)


    

In [19]:

    

for k,v in wms.contents.iteritems():
    print v.name,v.title


    

    




1 pacific_coast_0106_all_change_wm
0 southeast_0106_all_change_wm
3 gulf_coast_0106_all_change_wm
2 northeast_0106_all_change_wm
4 great_lakes_0106_all_change_wm

In [20]:

    

layer = '2'


    

In [ ]:

    

ax = plt.axes(projection=ccrs.PlateCarree())
ax.add_wms(url, layer, crs='CRS:84')
ax.set_extent((-71.5, 41.5, -71.0, 42))


    

In [22]:

    

url='http://services.arcgisonline.com/arcgis/rest/services/Ocean_Basemap/MapServer/WMTS/1.0.0/WMTSCapabilities.xml'
layer='Ocean_Basemap'
ax = plt.axes(projection=ccrs.PlateCarree())
ax.add_wmts(url, layer)
ax.set_extent((-15, 25, 35, 60))


    

    




---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-22-7856174d5762> in <module>()
      2 layer='Ocean_Basemap'
      3 ax = plt.axes(projection=ccrs.PlateCarree())
----> 4 ax.add_wmts(url, layer)
      5 ax.set_extent((-15, 25, 35, 60))

/home/usgs/miniconda/envs/ioos/lib/python2.7/site-packages/cartopy/mpl/geoaxes.pyc in add_wmts(self, wmts, layer_name, **kwargs)
   1516         """
   1517         from cartopy.io.ogc_clients import WMTSRasterSource
-> 1518         wmts = WMTSRasterSource(wmts, layer_name)
   1519         return self.add_raster(wmts, **kwargs)
   1520 

/home/usgs/miniconda/envs/ioos/lib/python2.7/site-packages/cartopy/io/ogc_clients.pyc in __init__(self, wmts, layer_name)
    191                 hasattr(wmts, 'contents') and
    192                 hasattr(wmts, 'gettile')):
--> 193             wmts = owslib.wmts.WebMapTileService(wmts)
    194 
    195         try:

/home/usgs/miniconda/envs/ioos/lib/python2.7/site-packages/OWSLib-0.8.13-py2.7.egg/owslib/wmts.pyc in __init__(self, url, version, xml, username, password, parse_remote_metadata, vendor_kwargs)
    171 
    172         # build metadata objects
--> 173         self._buildMetadata(parse_remote_metadata)
    174 
    175     def _getcapproperty(self):

/home/usgs/miniconda/envs/ioos/lib/python2.7/site-packages/OWSLib-0.8.13-py2.7.egg/owslib/wmts.pyc in _buildMetadata(self, parse_remote_metadata)
    191         # serviceProvider metadata
    192         serviceprov = self._capabilities.find(_SERVICE_PROVIDER_TAG)
--> 193         self.provider = ServiceProvider(serviceprov)
    194 
    195         # serviceOperations metadata

/home/usgs/miniconda/envs/ioos/lib/python2.7/site-packages/OWSLib-0.8.13-py2.7.egg/owslib/ows.pyc in __init__(self, infoset, namespace)
     78     def __init__(self, infoset,namespace=DEFAULT_OWS_NAMESPACE):
     79         self._root = infoset
---> 80         val = self._root.find(util.nspath('ProviderName', namespace))
     81         self.name = util.testXMLValue(val)
     82         self.contact = ServiceContact(infoset, namespace)

AttributeError: 'NoneType' object has no attribute 'find'





    

In [ ]:

    




    

In [ ]:

    

http://coast.noaa.gov/arcgis/services/CCAP/All_Change_0106/MapServer/WMSServer?request=GetCapabilities&service=WMS