code: def _indiameter_from_fibers(self):
          r
          Np = self.num_pores()
          indiam = np.zeros(Np, dtype=float)
          incen = np.zeros([Np, 3], dtype=float)
          hull_pores = np.unique(self._hull_image)
          (Lx, Ly, Lz) = np.shape(self._hull_image)
          (indx, indy, indz) = np.indices([Lx, Ly, Lz])
          indx = indx.flatten()
          indy = indy.flatten()
          indz = indz.flatten()
          for i, pore in enumerate(hull_pores):
              logger.info("Processing pore: "+str(i)+" of "+str(len(hull_pores)))
              dt_pore = self._dt_image*(self._hull_image == pore)
              indiam[pore] = dt_pore.max()*2
              max_ind = np.argmax(dt_pore)
              incen[pore, 0] = indx[max_ind]
              incen[pore, 1] = indy[max_ind]
              incen[pore, 2] = indz[max_ind]
          indiam *= self.network.resolution
          incen *= self.network.resolution
          return (indiam, incen)
docs: r"""
          Calculate an indiameter by distance transforming sections of the
          fiber image. By definition the maximum value will be the largest radius
          of an inscribed sphere inside the fibrous hull