docs: def within(self, region):
          """
          :param region: a shapely polygon
          :returns: a filtered SiteCollection of sites within the region
          """
          mask = numpy.array([
              geometry.Point(rec['lon'], rec['lat']).within(region)
              for rec in self.array])
          return self.filter(mask)