code: def filter_beams(prims, prim_with_staff, bounds):
      n_bounds = []
      n_prims = []
      n_prim_with_staff = []
      for i, prim in enumerate(prims):
          if prim.shape[1] >= 2 * prim.shape[0]:
              continue
          else:
              n_bounds.append(bounds[i])
              n_prims.append(prims[i])
              n_prim_with_staff.append(prim_with_staff[i])
      return n_prims, n_prim_with_staff, n_bounds
  
file_path: Mozart/src/main.py
imports: from commonfunctions import *
  from pre_processing import *
  from connected_componentes import *
  from staff import calculate_thickness_spacing, remove_staff_lines, coordinator
  from segmenter import Segmenter
  from fit import predict
  from glob import glob
  import cv2
  import pickle
  from scipy.ndimage import binary_fill_holes
  from skimage.morphology import thin
  import argparse
label_desc: Use image node_id 3 for calling a global function with example usage: filter_beams(prims, prim_with_staff, bounds) and returns: n_prims,  n_prim_with_staff,  n_bounds