docs: def _get_proj_enclosing_polygon(self):
          """
          See :meth:`Mesh._get_proj_enclosing_polygon`.
  
          :class:`RectangularMesh` contains an information about relative
          positions of points, so it allows to define the minimum polygon,
          containing the projection of the mesh, which doesn't necessarily
          have to be convex (in contrast to :class:`Mesh` implementation).
  
          :returns:
              Same structure as :meth:`Mesh._get_proj_convex_hull`.
          """
          if self.lons.size < 4:
              # the mesh doesn't contain even a single cell
              return self._get_proj_convex_hull()
  
          proj = geo_utils.OrthographicProjection(
              *geo_utils.get_spherical_bounding_box(self.lons, self.lats))
          if len(self.lons.shape) == 1:  # 1D mesh
              lons = self.lons.reshape(len(self.lons), 1)
              lats = self.lats.reshape(len(self.lats), 1)
          else:  # 2D mesh
              lons = self.lons.T
              lats = self.lats.T
          mesh2d = numpy.array(proj(lons, lats)).T
          lines = iter(mesh2d)
          # we iterate over horizontal stripes, keeping the "previous"
          # line of points. we keep it reversed, such that together
          # with the current line they define the sequence of points
          # around the stripe.
          prev_line = next(lines)[::-1]
          polygons = []
          for i, line in enumerate(lines):
              coords = numpy.concatenate((prev_line, line, prev_line[0:1]))
              # create the shapely polygon object from the stripe
              # coordinates and simplify it (remove redundant points,
              # if there are any lying on the straight line).
              stripe = shapely.geometry.LineString(coords) \
                                       .simplify(self.DIST_TOLERANCE) \
                                       .buffer(self.DIST_TOLERANCE, 2)
              polygons.append(shapely.geometry.Polygon(stripe.exterior))
              prev_line = line[::-1]
          try:
              # create a final polygon as the union of all the stripe ones
              polygon = shapely.ops.cascaded_union(polygons) \
                                   .simplify(self.DIST_TOLERANCE)
          except ValueError:
              # NOTE(larsbutler): In some rare cases, we've observed ValueErrors
              # ("No Shapely geometry can be created from null value") with very
              # specific sets of polygons such that there are two unique
              # and many duplicates of one.
              # This bug is very difficult to reproduce consistently (except on
              # specific platforms) so the work around here is to remove the
              # duplicate polygons. In fact, we only observed this error on our
              # CI/build machine. None of our dev environments or production
              # machines has encountered this error, at least consistently. >:(
              polygons = [shapely.wkt.loads(x) for x in
                          list(set(p.wkt for p in polygons))]
              polygon = shapely.ops.cascaded_union(polygons) \
                                   .simplify(self.DIST_TOLERANCE)
          return proj, polygon