vector.py

from pathlib import Path

from osgeo import ogr


ogr.UseExceptions()


def get_features(vector_path: str | Path) -> list[ogr.Feature]:
    ds = ogr.Open(str(vector_path))
    layer = ds.GetLayer()
    return [feature for feature in layer]


def get_property_values_for_intersecting_features(geometry: ogr.Geometry, features: list[ogr.Feature]) -> bool:
    for feature in features:
        if feature.GetGeometryRef().Intersects(geometry):
            return True
    return False


def intersecting_feature_properties(
    geometry: ogr.Geometry, features: list[ogr.Feature], feature_property: str
) -> list[str]:
    property_values = []
    for feature in features:
        if feature.GetGeometryRef().Intersects(geometry):
            property_values.append(feature.GetField(feature_property))
    return property_values