<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Thanks,<br><br>All things come to those who wait :-)<br><br>My data is 2D at present...<br><br>If I can't get thr PL/R solution working I think I'm better off waiting for GEOS 3.4/Postgis 2.1<br><br>Cheers,<br><br> Brent<br><br><br>--- On <b>Sat, 9/1/12, Paragon Corporation <i><lr@pcorp.us></i></b> wrote:<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>From: Paragon Corporation <lr@pcorp.us><br>Subject: Re: [postgis-users] Concave hull problem & delaunay triangulation -help needed please<br>To: "'PostGIS Users Discussion'" <postgis-users@postgis.refractions.net><br>Date: Saturday, September 1, 2012, 4:48 PM<br><br><div id="yiv1269145993">
<div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2">Brent,</font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2">Is your data 2D or 3D? You might want to try out the new
ST_DelaunayTriangles in 2.1</font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"><a rel="nofollow" target="_blank" href="http://www.postgis.org/documentation/manual-svn/ST_DelaunayTriangles.html">http://www.postgis.org/documentation/manual-svn/ST_DelaunayTriangles.html</a></font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2">Right now I beleive it's just 2D, but work is going on in GEOS
land to make it 3D.: <a rel="nofollow" target="_blank" href="http://trac.osgeo.org/geos/ticket/570">http://trac.osgeo.org/geos/ticket/570</a></font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"> It does require upgrading your system to 2.1 development
and Geos 3.4.0 dev so might not be an option for you
yet.</font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2">We also have windows build bot builds for those on
windows who want to try which gets built whenever there is a change in the 2.1
code base</font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"><a rel="nofollow" target="_blank" href="http://www.postgis.org/download/windows/experimental.php">http://www.postgis.org/download/windows/experimental.php</a>
(for 9.2 both 32-bit/64-bit and for 9.0,9.1 just 64-bit
windows)</font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"></font></span> </div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2">Hope that helps,</font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2">Regina</font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"><a rel="nofollow" target="_blank" href="http://www.postgis.us">http://www.postgis.us</a></font></span></div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"></span> </div>
<div dir="ltr" align="left"><span class="yiv1269145993949094204-01092012"><font color="#0000ff" face="Arial" size="2"></font></span> </div><br>
<div dir="ltr" class="yiv1269145993OutlookMessageHeader" align="left" lang="en-us">
<hr tabindex="-1">
<font face="Tahoma" size="2"><b>From:</b>
postgis-users-bounces@postgis.refractions.net
[mailto:postgis-users-bounces@postgis.refractions.net] <b>On Behalf Of
</b>pcreso@pcreso.com<br><b>Sent:</b> Friday, August 31, 2012 11:42
PM<br><b>To:</b> PostGIS Users Discussion<br><b>Subject:</b> [postgis-users]
Concave hull problem & delaunay triangulation -help needed
please<br></font><br></div>
<div></div>
<table border="0" cellpadding="0" cellspacing="0">
<tbody>
<tr>
<td valign="top">Hi,<br><br>I'm trying to come up with an approach using
Postgis to come up with a polygon representing the outer bounds of a set
of points, ie: a polygon passing through all the outer points.<br><br>Much
useful info (as usual - thanks guys!!) at spatialdbadvisor & bostongis
but my limited understanding is stalling progress for now.<br><br>I came
across the ST_Concavehull() function, which fails. Indeed, I was surprised
to see this actually expands the polygon extent representing the point
set, as shown at:<br><a rel="nofollow" target="_blank" href="http://www.bostongis.com/postgis_concavehull.snippet">http://www.bostongis.com/postgis_concavehull.snippet</a>
as this buffers the points in some cases, and often gives quite unexpected
(to me anyway) results. <br><br>Looking for something that fits better,
using existing capabilities, I figured I could pretty easily use Delaunay
triangulation to create the outer hull, then realised this bridges every
second outer point (pretty much by definition). That's easy, just remove
any triangles from the polygon where they share a line segment with the
polygon. I like Postgis!!<br><br>I have Simon's code working fine, at
least using his example invocation, from:<br><a rel="nofollow" target="_blank" href="http://www.spatialdbadvisor.com/postgis_tips_tricks/290/r-based-delaunay-triangulation-function-for-postgis-using-the-deldir-package">http://www.spatialdbadvisor.com/postgis_tips_tricks/290/r-based-delaunay-triangulation-function-for-postgis-using-the-deldir-package</a><br><br>But
I can't figure out how to make it work with my point data. I don't know R
or PL/R well enough to adapt it.<br><br>I have a set of sample points in a
table:<br><br><font face="courier,monaco,monospace,sans-serif"> Column
| Type
|
<br> --------+-----------------+-----------------<br>
id |
integer | not null
default<br> geom | geometry(Point) | </font><br><br>Can
anyone explain, ideally in words of one syllable or less, how Simon's
function from the link above needs to be adapted to read my points dataset
as input?<br><br>Once I get the triangles as an output set I can use
Postgis for the rest...<br><br>Thanks,<br><br> Brent
Wood<br><br><br></td></tr></tbody></table></div>
</div><br>-----Inline Attachment Follows-----<br><br><div class="plainMail">_______________________________________________<br>postgis-users mailing list<br><a ymailto="mailto:postgis-users@postgis.refractions.net" href="/mc/compose?to=postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br></div></blockquote></td></tr></table>