code: def get_modis_tile_list(ds):
      
      from demcoreg import modis_grid
      modis_dict = {}
      for key in modis_grid.modis_dict:
          modis_dict[key] = ogr.CreateGeometryFromWkt(modis_grid.modis_dict[key])
      geom = geolib.ds_geom(ds)
      geom_dup = geolib.geom_dup(geom)
      ct = osr.CoordinateTransformation(geom_dup.GetSpatialReference(), geolib.wgs_srs)
      geom_dup.Transform(ct)
      tile_list = []
      for key, val in list(modis_dict.items()):
          if geom_dup.Intersects(val):
              tile_list.append(key)
      return tile_list
docs: Helper function to identify MODIS tiles that intersect input geometry
  
      modis_gird.py contains dictionary of tile boundaries (tile name and WKT polygon ring from bbox)
  
      See: https://modis-land.gsfc.nasa.gov/MODLAND_grid.html