def plot_pore(self, pores, fig=None, axis_bounds=None,
                  include_points=False):
        r"""
        Plot all throats around a given pore or list of pores in 3D

        Parameters
        ----------
        pores : list or array containing pore indices tp include in figure

        fig : matplotlib figure object to place plot in

        axis_bounds : list of [xmin, xmax, ymin, ymax, zmin, zmax] values
            to limit axes to

        include_points : bool
            Determines whether to scatter pore and throat centroids
        """
        from mpl_toolkits.mplot3d.art3d import Poly3DCollection
        if len(pores) > 0:
            net = self.network
            net_pores = net.map_pores(pores=pores, origin=self)
            centroids = self['pore.centroid'][pores]
            coords = net['pore.coords'][net_pores]
            net_throats = net.find_neighbor_throats(pores=net_pores)
            throats = self.map_throats(throats=net_throats, origin=net)
            tcentroids = self["throat.centroid"][throats]
            # Can't create volume from one throat
            if 1 <= len(throats):
                verts = self['throat.vertices'][throats]
                normals = self['throat.normal'][throats]
                # Get verts in hull order
                ordered_verts = []
                for i in range(len(verts)):
                    vert_2D = self._rotate_and_chop(verts[i], normals[i],
                                                    [0, 0, 1])
                    hull = ConvexHull(vert_2D, qhull_options='QJ Pp')
                    ordered_verts.append(verts[i][hull.vertices])
                offsets = self['throat.offset_vertices'][throats]
                ordered_offs = []
                for i in range(len(offsets)):
                    offs_2D = self._rotate_and_chop(offsets[i], normals[i],
                                                    [0, 0, 1])
                    offs_hull = ConvexHull(offs_2D, qhull_options='QJ Pp')
                    ordered_offs.append(offsets[i][offs_hull.vertices])
                # Get domain extents for setting axis
                if axis_bounds is None:
                    [xmin, xmax, ymin, ymax, zmin, zmax] = \
                        self.vertex_dimension(net_pores, parm='minmax')
                else:
                    [xmin, xmax, ymin, ymax, zmin, zmax] = axis_bounds
                if fig is None:
                    fig = plt.figure()
                ax = fig.gca(projection='3d')
                outer_items = Poly3DCollection(ordered_verts, linewidths=1,
                                               alpha=0.2, zsort='min')
                outer_face_colours = [(1, 0, 0, 0.01)]
                outer_items.set_facecolor(outer_face_colours)
                ax.add_collection(outer_items)
                inner_items = Poly3DCollection(ordered_offs, linewidths=1,
                                               alpha=0.2, zsort='min')
                inner_face_colours = [(0, 0, 1, 0.01)]
                inner_items.set_facecolor(inner_face_colours)
                ax.add_collection(inner_items)
                ax.set_xlim(xmin, xmax)
                ax.set_ylim(ymin, ymax)
                ax.set_zlim(zmin, zmax)
                if include_points:
                    ax.scatter(centroids[:, 0],
                               centroids[:, 1],
                               centroids[:, 2],
                               c='y')
                    ax.scatter(tcentroids[:, 0],
                               tcentroids[:, 1],
                               tcentroids[:, 2],
                               c='r')
                    ax.scatter(coords[:, 0], coords[:, 1], coords[:, 2], c='b')
                ax.ticklabel_format(style='sci', scilimits=(0, 0))
            else:
                self.plot_throat(throats, fig)
        else:
            logger.error('Please provide pore indices')
        return fig
r"""
        Plot all throats around a given pore or list of pores in 3D

        Parameters
        ----------
        pores : list or array containing pore indices tp include in figure

        fig : matplotlib figure object to place plot in

        axis_bounds : list of [xmin, xmax, ymin, ymax, zmin, zmax] values
            to limit axes to

        include_points : bool
            Determines whether to scatter pore and throat centroids