code: def ps_disk(radius):
      r
      rad = int(sp.ceil(radius))
      other = sp.ones((2 * rad + 1, 2 * rad + 1), dtype=bool)
      other[rad, rad] = False
      disk = spim.distance_transform_edt(other) < radius
      return disk
docs: r"""
      Creates circular disk structuring element for morphological operations
  
      Parameters
      ----------
      radius : float or int
          The desired radius of the structuring element
  
      Returns
      -------
      strel : 2D-array
          A 2D numpy bool array of the structring element