Hi Jose,<div><br></div><div>thanks for that. I&#39;ve used SetSpatialFilterRect(x,y,x+0.01,y+0.01)</div><div><br><br><div class="gmail_quote">On Wed, Dec 8, 2010 at 12:16 PM, Jose Gomez-Dans <span dir="ltr">&lt;<a href="mailto:jgomezdans@gmail.com">jgomezdans@gmail.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi,<br>
<div class="im"><br>
On 8 December 2010 09:38, Sebastian E. Ovide &lt;<a href="mailto:sebastian.ovide@gmail.com">sebastian.ovide@gmail.com</a>&gt; wrote:<br>
&gt; Using Python + Gdal/OGR is there any efficient way of finding all the<br>
&gt; polygons that contains the point X,Y ?<br>
<br>
</div>This is my try. Not the most efficient, but you could always use a<br>
spatial index and a spatial filter:<br>
&lt;<a href="http://gist.github.com/733207" target="_blank">http://gist.github.com/733207</a>&gt;<br>
<font color="#888888"><br>
Jose<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Sebastian E. Ovide<br><br><br><br><br>
</div>