All,<br>
<br>
I am trying to calculate intersection of two polygon layers with<br>
polygons.  I searched threads, found following thread, and what I got<br>
was that OGR/GEOS provide me to intersection of two features, but I<br>
have to write something in order to perform this operation for two<br>
layers with multiple features.<br>
<br>
Open source vector geoprocessing libraries?<br>
<a href="http://groups.google.com/group/gdal/browse_thread/thread/7b5a4f461a4db319/" target="_blank">http://groups.google.com/group/gdal/browse_thread/thread/7b5a4f461a4db319/</a><br>
<br>
Intersect of two shapefiles<br>
<a href="http://groups.google.com/group/gdal/browse_thread/thread/4b178c776aecf425/" target="_blank">http://groups.google.com/group/gdal/browse_thread/thread/4b178c776aecf425/</a><br>
<br>
My questions are:<br>
<br>
1. Did I miss something?  Is there development of library which takes<br>
two polygon layers to calculate intersection or other basic geometry<br>
operations that i find in OGR?<br>
<br>
2. Assuming no to above, I got started writing code using python<br>
binding.  Not having too much experience, I codes to go through all<br>
possible pair (complete bipartite), test with intersect, then run<br>
intersection if it returns true.  I am wondering there is faster way<br>
to do this.  I am thinking something like get envelope polygon of each<br>
feature first, and then test intersection of them and then perform<br>
intersect operation only when first test passes (i assume that it is<br>
easier to test intersect if there are only handful of vertices as<br>
oppose to hundreds of vertices in my case).<br>
<br>
Right now I am trying with layer of hundreds of thousands of polygons<br>
vs. thousands of polygon, and it seems to be taking 10s of hours to<br>
finish in my environment.<br>
<br>
Thank you,<br>
<br>
Yosuke Kimura