[Qgis-user] Python code to select polygon at location

[Greg Troxel]

Zoltan Szecsei via QGIS-User <qgis-user at lists.osgeo.org> writes:

> Hi,
> I'm doing my "once a year need to do something in python with QGIS" :-/
> Using QGIS 3.38 and PyCharm.
>
> I have a SHP file with 28000 polygons.
> Another SHP file with 22000 Points
> Both in EPSG:4326
>
> Without using QGIS processing, what is the quickest way to, for each
> point 1 by 1, find the underlying polygon and read the attribute field
> values of that polygon?

Load it into postgis and write a select statement?