<div>Leticia,</div>
<div> </div>
<div>I note that it can be improved a little bit, the first loop starts in 0 instead of 1:</div>
<div> </div>
<div>CREATE OR REPLACE FUNCTION ExtendLine(eje_ geometry, bound_ geometry) RETURNS geometry AS $$<br>-- version: alfa , by Julio A. Galindo, April 17, 2007: <a href="mailto:juliogalindoq@gmail.com">juliogalindoq@gmail.com
</a><br> DECLARE<br> b_ geometry = boundary(bound_);<br> dist float;<br> max_dist float = 0;<br> n_points int;<br> pto_1 geometry;<br> pto_2 geometry;<br> first_pto geometry;<br> last_pto geometry;
<br> u_1 float;<br> u_2 float;<br> norm float;<br> result text = 'LINESTRING(';<br> BEGIN<br> IF GeometryType(eje_) NOT LIKE 'LINESTRING'<br> OR GeometryType(bound_) NOT LIKE 'POLYGON' THEN RETURN NULL; END IF;
<br> -- First Search how far is the boundary: (worst case)<br> pto_1 := StartPoint(eje_);<br> pto_2 := EndPoint(eje_);<br> FOR i IN 1..NumPoints(b_)-1 LOOP<br> dist := distance(PointN(b_,i),pto_1);<br> IF dist > max_dist THEN max_dist := dist; END IF;
<br> dist := distance(PointN(b_,i),pto_2);<br> IF dist > max_dist THEN max_dist := dist; END IF;<br> END LOOP;<br> -- Now extent the linestring:<br> pto_2 := PointN(eje_,2);<br> u_1 := X(pto_2)-X(pto_1);
<br> u_2 := Y(pto_2)-Y(pto_1);<br> norm := sqrt(u_1^2 + u_2^2);<br> first_pto := MakePoint(X(pto_1)-u_1/norm*dist,Y(pto_1)-u_2/norm*dist);<br> n_points := nPoints(eje_);<br> IF n_points > 2 THEN<br>
pto_1 := PointN(eje_,n_points-1);<br> pto_2 := PointN(eje_,n_points);<br> u_1 := X(pto_2)-X(pto_1);<br> u_2 := Y(pto_2)-Y(pto_1);<br> norm := sqrt(u_1^2 + u_2^2);<br> END IF;<br> last_pto := MakePoint(X(pto_2)+u_1/norm*dist,Y(pto_2)+u_2/norm*dist);
<br> result := result || X(first_pto) || ' ' || Y(first_pto) || ',';<br> FOR i IN 1..NumPoints(eje_) LOOP<br> result := result || X(PointN(eje_,i)) || ' ' || Y(PointN(eje_,i)) || ',';
<br> END LOOP;<br> result := result || X(last_pto) || ' ' || Y(last_pto) || ')';<br> -- Find the final Linestring:<br> b_ := intersection(GeomFromText(result,SRID(eje_)),bound_);<br> RETURN b_;
<br> END $$<br> LANGUAGE plpgsql<br> STABLE<br> RETURNS NULL ON NULL INPUT;<br><br><br> </div>
<div><span class="gmail_quote">On 4/17/07, <b class="gmail_sendername">Leticia</b> <<a href="mailto:lgomez@itba.edu.ar">lgomez@itba.edu.ar</a>> wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div lang="EN-US" vlink="blue" link="blue">
<div>
<p><font face="Arial" color="navy" size="2"><span style="FONT-SIZE: 10pt; COLOR: navy; FONT-FAMILY: Arial">Julio, THANKS!</span></font></p>
<p><font face="Arial" color="navy" size="2"><span style="FONT-SIZE: 10pt; COLOR: navy; FONT-FAMILY: Arial"> </span></font></p>
<p><font face="Arial" color="navy" size="2"><span style="FONT-SIZE: 10pt; COLOR: navy; FONT-FAMILY: Arial"> </span></font></p>
<p><font face="Arial" color="navy" size="2"><span style="FONT-SIZE: 10pt; COLOR: navy; FONT-FAMILY: Arial"> </span></font></p>
<div>
<div style="TEXT-ALIGN: center" align="center"><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">
<hr align="center" width="100%" size="3">
</span></font></div>
<p><b><font face="Tahoma" size="2"><span style="FONT-WEIGHT: bold; FONT-SIZE: 10pt; FONT-FAMILY: Tahoma">From:</span></font></b><font face="Tahoma" size="2"><span style="FONT-SIZE: 10pt; FONT-FAMILY: Tahoma"> <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">
postgis-users-bounces@postgis.refractions.net</a> [mailto:<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net
</a>] <b><span style="FONT-WEIGHT: bold">On Behalf Of </span></b>Julio Galindo<br><b><span style="FONT-WEIGHT: bold">Sent:</span></b> Tuesday, April 17, 2007 1:53 AM<br><b><span style="FONT-WEIGHT: bold">To:</span></b> PostGIS Users Discussion
<br><b><span style="FONT-WEIGHT: bold">Subject:</span></b> Re: [postgis-users] Extending Linestrings</span></font></p></div>
<div><span class="e" id="q_111fdf5eb17f49aa_1">
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">Hi Leticia,</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">I create the following function for your problem, please test it carefully, I only made two test (at the end):</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">CREATE OR REPLACE FUNCTION ExtendLine(eje_ geometry, bound_ geometry) RETURNS geometry AS $$<br>-- version: alfa , by Julio A. Galindo, April 17, 2007:
<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:juliogalindoq@gmail.com" target="_blank">juliogalindoq@gmail.com </a><br> DECLARE<br> b_ geometry = boundary(bound_);<br> dist float;<br> max_dist float = 0;
<br> n_points int;<br> pto_1 geometry;<br> pto_2 geometry;<br> first_pto geometry;<br> last_pto geometry; <br> u_1 float;<br> u_2 float;<br> norm float;<br> result text = 'LINESTRING(';
<br> BEGIN<br> IF GeometryType(eje_) NOT LIKE 'LINESTRING'<br> OR GeometryType(bound_) NOT LIKE 'POLYGON' THEN RETURN NULL; END IF; <br> -- First Search how far is the boundary: (worst case)<br> pto_1 := StartPoint(eje_);
<br> pto_2 := EndPoint(eje_);<br> FOR i IN 0..NumPoints(b_)-1 LOOP<br> dist := distance(PointN(b_,i),pto_1);<br> IF dist > max_dist THEN max_dist := dist; END IF; <br> dist := distance(PointN(b_,i),pto_2);<br>
IF dist > max_dist THEN max_dist := dist; END IF;<br> END LOOP;<br> -- Now extent the linestring:<br> pto_2 := PointN(eje_,2);<br> u_1 := X(pto_2)-X(pto_1); <br> u_2 := Y(pto_2)-Y(pto_1);<br> norm := sqrt(u_1^2 + u_2^2);
<br> first_pto := MakePoint(X(pto_1)-u_1/norm*dist,Y(pto_1)-u_2/norm*dist);<br> n_points := nPoints(eje_);<br> IF n_points > 2 THEN<br> pto_1 := PointN(eje_,n_points-1);<br> pto_2 := PointN(eje_,n_points);
<br> u_1 := X(pto_2)-X(pto_1);<br> u_2 := Y(pto_2)-Y(pto_1);<br> norm := sqrt(u_1^2 + u_2^2);<br> END IF;<br> last_pto := MakePoint(X(pto_2)+u_1/norm*dist,Y(pto_2)+u_2/norm*dist); <br> result := result || X(first_pto) || ' ' || Y(first_pto) || ',';
<br> FOR i IN 1..NumPoints(eje_) LOOP<br> result := result || X(PointN(eje_,i)) || ' ' || Y(PointN(eje_,i)) || ','; <br> END LOOP;<br> result := result || X(last_pto) || ' ' || Y(last_pto) || ')';
<br> -- Find the final Linestring:<br> b_ := intersection(GeomFromText(result,SRID(eje_)),bound_);<br> RETURN b_; <br> END $$<br> LANGUAGE plpgsql<br> STABLE<br> RETURNS NULL ON NULL INPUT;<br> </span></font></p>
</div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">Test 1:</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">select asText(extendLine(geomFromText('LINESTRING(5 5, 12 12)',-1),GeomFromText('POLYGON((0 0,15 0,15 20,0 20,0 0))',-1)));</span></font>
</p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> astext <br>---------------------------------<br> LINESTRING(0 0,5 5,12 12,15 15)<br> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">The solution is good, because the line is at 45 degress </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">And the same function can work with a general linestring (many points):</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">select asText(extendLine(geomFromText('LINESTRING(5 5, 12 12, 13 12)',-1),GeomFromText('POLYGON((0 0,15 0,15 20,0 20,0 0))',-1)));<br>
astext <br>--------------------------------------- <br> LINESTRING(0 0,5 5,12 12,13 12,15 12)</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">Let me know if you improve this simple function.</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">Regards,</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">Julio A. Galindo Q.</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">Topogen Ltda.</span></font></p></div>
<div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">La Paz</span></font>, Bolivia<br><br> </p></div>
<div>
<p><span><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt">On 4/16/07, <b><span style="FONT-WEIGHT: bold">Leticia</span></b> <<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:lgomez@itba.edu.ar" target="_blank">
lgomez@itba.edu.ar</a>> wrote:</span></font></span> </p>
<div vlink="purple" link="blue">
<div>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">Hi community,</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">I need to "extend" a linestring in another one. The new linestring must be obtained such as their new extremes are the intersection of it with a polygon that contains it. The original linestring is composed of only a linear segment.
</span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">Is there a simple method to do this?</span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">-------------</span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">| |</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> | | |</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> | | |</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> | |</span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">-------------</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">I need the other one:</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">-------|-----</span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">| | |</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> | | |</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> | | |</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> | | |</span></font></p>
<p style="TEXT-INDENT: 0.5in"><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">-------|-----</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">The original linestring is not always vertical u horizontal.</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">Thanks in advance,</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt">Leticia</span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Courier New" size="2"><span style="FONT-SIZE: 10pt"> </span></font></p>
<p><font face="Arial" size="2"><span style="FONT-SIZE: 10pt; FONT-FAMILY: Arial"> </span></font></p></div></div>
<p style="MARGIN-BOTTOM: 12pt"><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"><br>_______________________________________________<br>postgis-users mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:postgis-users@postgis.refractions.net" target="_blank">
postgis-users@postgis.refractions.net</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users
</a></span></font></p></div>
<p><font face="Times New Roman" size="3"><span style="FONT-SIZE: 12pt"> </span></font></p></span></div></div></div><br>_______________________________________________<br>postgis-users mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:postgis-users@postgis.refractions.net">
postgis-users@postgis.refractions.net</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users
</a><br><br></blockquote></div><br>