code: def dt_filter_rel_annual_idx(dt_list, min_rel_dt=(1,1), max_rel_dt=(12,31)):
      
      dt_list = np.array(dt_list)
      years = get_unique_years(dt_list)
      from collections import OrderedDict
      out = OrderedDict() 
      for year in years:
          #If within the same year
          if min_rel_dt[0] < max_rel_dt[1]:
              dt1 = datetime(year, min_rel_dt[0], min_rel_dt[1])
              dt2 = datetime(year, max_rel_dt[0], max_rel_dt[1])
          #Or if our relative values include Jan 1
          else:
              dt1 = datetime(year, min_rel_dt[0], min_rel_dt[1])
              dt2 = datetime(year+1, max_rel_dt[0], max_rel_dt[1])
          idx = np.logical_and((dt_list >= dt1), (dt_list <= dt2))
          if np.any(idx):
              out[year] = idx
      return out
docs: Return dictionary containing indices of timestamps that fall within relative month/day bounds of each year