[gdal-dev] spatial join using python ogr or shapely

[Giuseppe Amatulli]

HI,
I would like to make a spatial join of two shp using python ogr or shapely
Is this possible?
I do not have a common Item in the dbf so i can not make a table-join,
i have to do a spatial one.
I did not find any examples on the net.

so started ....

try:
    from osgeo import ogr
    import sys

shp00=ogr.Open("alabama.shp")
shp10=ogr.Open("tabblock2010_01_pophu.shp")

inPolyLshp00 = shp00.GetLayer(0)
inPolyLshp10 = shp10.GetLayer(0)

and then?

Any suggestion?
Thanks in advance
-- 
Giuseppe Amatulli
Web: www.spatial-ecology.net