code: def Plane(center=(0, 0, 0), direction=(0, 0, 1), i_size=1, j_size=1,
            i_resolution=10, j_resolution=10):
      
      planeSource = vtk.vtkPlaneSource()
      planeSource.SetXResolution(i_resolution)
      planeSource.SetYResolution(j_resolution)
      planeSource.Update()
  
      surf = PolyData(planeSource.GetOutput())
  
      surf.points[:, 0] *= i_size
      surf.points[:, 1] *= j_size
      surf.rotate_y(-90)
      translate(surf, center, direction)
      return surf
docs: Create a plane
  
      Parameters
      ----------
      center : list or np.ndarray
          Location of the centroid in [x, y, z]
  
      direction : list or np.ndarray
          Direction cylinder points to  in [x, y, z]
  
      i_size : float
          Size of the plane in the i direction.
  
      j_size : float
          Size of the plane in the i direction.
  
      i_resolution : int
          Number of points on the plane in the i direction.
  
      j_resolution : int
          Number of points on the plane in the j direction.
  
      Returns
      -------
      plane : vtki.PolyData
          Plane mesh