<div dir="ltr"><br><div class="gmail_quote"><div dir="ltr"><div><div>Hi to eveyrone,<br><br></div>my name is Konstantinos and I am new both in the mailing list and postgis. My field is electrical engineering and currently  I m working on modelling of power networks. A sugnificant part of my work deals with spatial data. I have a question concerning an issue (explain it below) and would like to ask for your opinion.  The processing is done in 2 database from which I extract the following part of grids: <br><br><img src="cid:ii_iyx1xf8z2_15a1e212d032f917" style="margin-right:0px" width="255" height="147"><br><br></div><br><div>Both pictures illustrate the same object. The blue  is considered to be the valid one, while the green contains missing data. Therefore my aim is to compare the geometries of the lines of both grids and find similarities. For example the lines of the third image can be considers as identical and hence tha blue data can be copyed to the green. To mention that both databases have column with the geometries of the lines. I am running the following:<br><br><pre class="m_-2005620091781540490gmail-m_3961686253516418805gmail-lang-sql m_-2005620091781540490gmail-m_3961686253516418805gmail-prettyprint m_-2005620091781540490gmail-m_3961686253516418805gmail-prettyprinted"><font size="4"><code><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-kwd">create</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-kwd">index</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> geom_1_idx </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-kwd">ON</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> table1 </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-kwd">USING</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> GIST </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">(</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln">geometry_1</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">);<br></span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"></span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-kwd">select</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> table1_id</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">,</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> table2_id  </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-kwd">from</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> table1 t1</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">,</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> table2 t2
              </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-kwd">where</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">(</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-kwd">select</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> st_DWithin</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">(</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln">t1</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">.</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln">geometry_2 </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">,</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> t2</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">.</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln">geometry_2</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">,</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> units_of_srid</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">)</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> </span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">=</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pln"> true</span><span class="m_-2005620091781540490gmail-m_3961686253516418805gmail-pun">)</span></code></font></pre><br></div><div> I used spatial indexing according to your guidelines (<a href="http://revenant.ca/www/postgis/workshop/indexing.html" target="_blank">http://revenant.ca/www/<wbr>postgis/workshop/indexing.html</a><wbr>). As far I have understood by this code, I compare the bounding boxes of the lines and if the distance between the boxes is not higher than "units_of_srid", then the two geometries are considered identical. <br><br></div><div>My first question is if it s right what I wrote. How can I improve this in order not to miss "identical" lines?<br><br></div><div>Thank you in advance and for your understanding.<br></div><div>Looking forward for your comments and corrections.<br><br></div><div>Kind regards,<br><br></div><div>Konstantinos <br></div><div><br><br></div></div>
</div><br></div>