Пример #1
0
    def gshhs_line(self, outline_color='k', domain=None, resolution='low', **kwargs):
        # domain is a shapely geometry (Polygon or MultiPolygon)
        import cartopy.gshhs as gshhs
#        import cartopy.spherical as spherical
        from matplotlib.collections import PatchCollection, LineCollection

        paths = []

        projection = self.projection

        if domain is None:
            domain = self.map_domain(ccrs.PlateCarree())

        for points in gshhs.read_gshhc(gshhs.fnames[resolution], poly=False, domain=domain):
            paths.extend(patch.geos_to_path(shapely.geometry.LineString(points)))

#            slinestring = shapely.geometry.LineString(points)
#            projected = projection.project_geometry(slinestring)            
#            paths.extend(patch.geos_to_path(projected))

        collection = PatchCollection([mpatches.PathPatch(pth) for pth in paths],
                             edgecolor=outline_color, facecolor='none',
                             transform=ccrs.PlateCarree(),
                             **kwargs
                             )

        self.add_collection(collection, autolim=False)
Пример #2
0
    def gshhs(self, outline_color='k', land_fill='green', ocean_fill='None', resolution='coarse', domain=None):
        import cartopy.gshhs as gshhs
        from matplotlib.collections import PatchCollection

        if ocean_fill is not None and land_fill is None:
            land_fill = self.outline_patch.get_facecolor()

        if domain is None:
            domain = self.ll_boundary_poly()

        paths = []
        for points in gshhs.read_gshhc(gshhs.fnames[resolution], poly=True, domain=domain):
            # XXX Sometimes we only want to do lines...
            poly = shapely.geometry.Polygon(points[::-1, :])
            projected = self.projection.project_polygon(poly)

            paths.extend(patch.geos_to_path(projected))

        if ocean_fill is not None:
            self.outline_patch.set_facecolor(ocean_fill)

        collection = PatchCollection([mpatches.PathPatch(pth) for pth in paths],
                                     edgecolor=outline_color,
                                     facecolor=land_fill,
                                     zorder=2,
                                     )
        # XXX Should it update the limits??? (AND HOW???)
        self.add_collection(collection)