Get attributes from features in one layer that intersect features from another layer

getAttributesFromntersectingFeatures.py

from osgeo import ogr

def ogrIntersection(tiles_layer, site_fname):

    # Get layer for the site
    driver  = ogr.GetDriverByName("GeoJSON")
    src_GeoJSON = driver.Open(site_fname)
    site_layer = src_GeoJSON.GetLayer()

    # List with tiles for every feature in site GeoJSON
    tiles = []

    # Find overlapping features
    for feat1 in tiles_layer:
        geom1 = feat1.GetGeometryRef()
        for feat2 in site_layer:
            geom2 = feat2.GetGeometryRef()
            if geom2.Intersects(geom1):
                # Get field 0 containing the tile indicex
                tiles.append(feat1.GetField(0))

    return(tiles);

# MODIS tiles KML
fname = '/home/glopez/Projects/TATSSI/data/kmz/modis_sin.kml'
driver  = ogr.GetDriverByName('KML')
src_kml = driver.Open(fname)
tiles_layer = src_kml.GetLayer()

# One of the WorldPeatland sites
site_fname = '/data/world_peatlands/src/WorldPeatland/sites/kampar.geojson'

tiles_for_site = ogrIntersection(tiles_layer, site_fname)
print(tiles_for_site)