title: Cartopy oceans cover land
body: <p><em>Edit</em>: seems to be a known issue (<a href="https://github.com/SciTools/cartopy/issues/803#issuecomment-272286945" rel="nofollow noreferrer">https://github.com/SciTools/cartopy/issues/803#issuecomment-272286945</a>) with certain map projections. </p>
  
  <hr>
  
  <p>For some reason, the NaturalEarth oceans cover the continents, as shown in the example below. If I plot <code>'physical', 'land'</code> it behaves as expected (land is filled), but for <code>'physical', 'ocean'</code> I would expect the opposite, but instead everything is filled. Does anyone know why this happens, and/or how it is possible to fill only the oceans?</p>
  
  <p>I also tried <code>ax.add_feature(cartopy.feature.OCEAN)</code>, but then Cartopy crashes (......) </p>
  
  <pre><code>import cartopy
  import cartopy.crs as ccrs
  import cartopy.feature as cfeature
  import matplotlib.pyplot as pl
  
  projection = ccrs.LambertConformal(central_longitude=4.9, central_latitude=52)
  
  ax = pl.subplot(121, projection=projection)
  ax.add_feature(cfeature.NaturalEarthFeature('physical', 'land', '50m', edgecolor='face', facecolor='g'))
  ax.set_extent([-8, 17, 42, 60], ccrs.PlateCarree())
  
  ax = pl.subplot(122, projection=projection)
  ax.add_feature(cfeature.NaturalEarthFeature('physical', 'ocean', '50m', edgecolor='face', facecolor='g'))
  ax.set_extent([-8, 17, 42, 60], ccrs.PlateCarree())
  </code></pre>
  
  <p><a href="https://i.stack.imgur.com/dbuOH.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/dbuOH.png" alt="enter image description here"></a></p>
  
  <p><strong>EDIT</strong></p>
  
  <p>Seems to be a projection issue; with <code>PlateCarree()</code> it works as expected:</p>
  
  <pre><code>import cartopy
  import cartopy.crs as ccrs
  import cartopy.feature as cfeature
  import matplotlib.pyplot as pl
  
  projection = ccrs.LambertConformal(central_longitude=4.9, central_latitude=52)
  
  ax = pl.subplot(221, projection=projection)
  ax.add_feature(cfeature.NaturalEarthFeature('physical', 'land', '50m', edgecolor='face', facecolor='g'))
  ax.set_extent([-8, 17, 42, 60], ccrs.PlateCarree())
  
  ax = pl.subplot(222, projection=projection)
  ax.add_feature(cfeature.NaturalEarthFeature('physical', 'ocean', '50m', edgecolor='face', facecolor='g'))
  ax.set_extent([-8, 17, 42, 60], ccrs.PlateCarree())
  
  projection = ccrs.PlateCarree()
  
  ax = pl.subplot(223, projection=projection)
  ax.add_feature(cfeature.NaturalEarthFeature('physical', 'land', '50m', edgecolor='face', facecolor='g'))
  ax.set_extent([-8, 17, 42, 60], ccrs.PlateCarree())
  
  ax = pl.subplot(224, projection=projection)
  ax.add_feature(cfeature.NaturalEarthFeature('physical', 'ocean', '50m', edgecolor='face', facecolor='g'))
  ax.set_extent([-8, 17, 42, 60], ccrs.PlateCarree())
  </code></pre>
  
  <p><a href="https://i.stack.imgur.com/zNfjw.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/zNfjw.png" alt="enter image description here"></a></p>