<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">2014-05-21 21:02 GMT+02:00 Jeremy Palmer <span dir="ltr"><<a href="mailto:JPalmer@linz.govt.nz" target="_blank">JPalmer@linz.govt.nz</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class=""><br>
> As a side not<br>
> _ in your first use case, you could have forced the use of index by adding to the intersect :<br>
> AND admin.geom && other.geom,<br>
> or even computing this as a first step.<br>
<br>
</div>I thought ST_Intersect just expanded to "SELECT $1 && $2 AND _ST_Intersects($1,$2)"<br></blockquote><div><br>Yep,  but it seems that the bounding box was not used in your case<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">

<div class=""><br>
> _your admlinistrativ limit seems to be a partition of the space (no overlaping). In this case there exist very fast way to get the intersect answer, but it is not used in GIS field.<br>
<br>
</div>Interesting. Do you have a reference to this method? </blockquote><div><br></div><div>This is called a topology model, or a 2D Arrangment more formaly.<br></div><div>See for instance<br><a href="http://doc.cgal.org/latest/Arrangement_on_surface_2/group__PkgArrangement2.html">http://doc.cgal.org/latest/Arrangement_on_surface_2/group__PkgArrangement2.html</a><br>
</div><div>Or PostGIS Topology ( not performance-oriented).<br><br></div><div>Cheers,<br></div><div>Rémi-C<br></div></div></div></div>