No subject
Wed Jul 15 09:27:57 PDT 2009
The rest of them have a Ring Self-intersection. These last dont through an
error exception, only the one with the "Self-intersection" error does.
I guess I should do propper geometry validation first.
Regards
On Wed, Aug 18, 2010 at 1:38 AM, Paragon Corporation <lr at pcorp.us> wrote:
> Ricardo,
> If you are looking in PgAdmin III, it won't show the geometry if it is too
> big.
>
> You should verify that your multipolygon is valid. Many of the GEOS
> functions won't work with invalid polygons.
>
> SELECT ST_IsValid(yourgeometryhere)
>
> Also let us know which version of GEOS you are using
>
> SELECT postgis_full_version();
>
> Lots of topology exception issues have been fixed in newer versions of GEOS
> so upgrading might fix your problem if you are running a GEOS version older
> than 3.2
>
> Hope that helps,
> Regina and Leo
> http://www.postgis.us
>
> ------------------------------
> *From:* postgis-users-bounces at postgis.refractions.net [mailto:
> postgis-users-bounces at postgis.refractions.net] *On Behalf Of *Ricardo
> Bayley
> *Sent:* Wednesday, August 18, 2010 12:06 AM
> *To:* postgis-users at postgis.refractions.net
> *Subject:* [postgis-users] GEOS pointonsurface() threw an error!
>
> Hi,
>
> I am calculation the PointOnSurface for a table of polygons. I get this
> error message
>
>
> NOTICE: TopologyException: found non-noded intersection between
> LINESTRING (857703 1.04022e+006, 857648 1.04019e+006) and LINESTRING (857536
> 1.04019e+006, 857870 1.04019e+006) at 857648 1.04019e+006
>
>
> ERROR: GEOS pointonsurface() threw an error!
>
> ********** Error **********
>
> ERROR: GEOS pointonsurface() threw an error!
> Estado SQL:XX000
>
> I checked the one of the rows that is getting this error with
> ST_GeometryType and I get "ST_Multipolygon", although I dont get the
> ST_AsText, or any output type.
> Any ideas why ?
>
>
> regards
>
>
> Ricardo
>
>
>
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
>
--0016364d2b6387c154048e11fa2c
Content-Type: text/html; charset=UTF-8
Content-Transfer-Encoding: quoted-printable
Hi Regina and Leo,<div><br></div><div>My version is:=C2=A0"POSTGIS=3D&=
quot;1.5.0" GEOS=3D"3.2.0-CAPI-1.6.0" PROJ=3D"Rel. 4.6.=
1, 21 August 2008" LIBXML=3D"2.7.6" USE_STATS"</div><di=
v>PostgreSQL 8.4 running on a WINXP</div>
<div><br></div><div>I was playing with the ST_IsValid() function and the ST=
_IsValidReason().</div><div><br></div><div>What I notice is the geometry, o=
f course, is not valid. Reason:=C2=A0"Self-intersection[-6.4825e+006 -=
3.15411e+006]"</div>
<div><br></div><div>I have other 96 invalid geometries in a table of 22k re=
cords.</div><div>From those 96, only 1 has the Self-intersection invalid re=
ason.=C2=A0</div><div>The rest of them have a Ring=C2=A0Self-intersection. =
These last dont through an error exception, only the one with the "Sel=
f-intersection" error does.</div>
<div><br></div><div>I guess I should do propper geometry validation first.<=
/div><div><br></div><div>Regards</div><div><br></div><div><br><br><div clas=
s=3D"gmail_quote">On Wed, Aug 18, 2010 at 1:38 AM, Paragon Corporation <spa=
n dir=3D"ltr"><<a href=3D"mailto:lr at pcorp.us">lr at pcorp.us</a>></span>=
wrote:<br>
<blockquote class=3D"gmail_quote" style=3D"margin:0 0 0 .8ex;border-left:1p=
x #ccc solid;padding-left:1ex;">
<div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">Ricardo,</font></span></div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">If you are looking in PgAdmin III, it won't show the geome=
try=20
if it is too big.</font></span></div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial"></font></span>=C2=A0</div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">You should verify that your multipolygon is valid.=C2=A0 Many=
=20
of the GEOS functions won't work with invalid polygons.</font></span></=
div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial"></font></span>=C2=A0</div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">SELECT ST_IsValid(yourgeometryhere)</font></span></div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial"></font></span>=C2=A0</div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">Also let us know which version of GEOS you are=20
using</font></span></div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial"></font></span>=C2=A0</div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">SELECT postgis_full_version();</font></span></div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial"></font></span>=C2=A0</div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">Lots of topology exception issues have been fixed in newer=20
versions of GEOS so upgrading might fix your problem if you are running a G=
EOS=20
version older than 3.2</font></span></div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial"></font></span>=C2=A0</div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">Hope that helps,</font></span></div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial">Regina and Leo</font></span></div>
<div dir=3D"ltr" align=3D"left"><span><font color=3D"#0000ff" size=3D"2" fa=
ce=3D"Arial"><a href=3D"http://www.postgis.us" target=3D"_blank">http://www=
.postgis.us</a></font></span></div><br>
<div dir=3D"ltr" lang=3D"en-us" align=3D"left">
<hr>
<font size=3D"2" face=3D"Tahoma"><b>From:</b>=20
<a href=3D"mailto:postgis-users-bounces at postgis.refractions.net" target=3D"=
_blank">postgis-users-bounces at postgis.refractions.net</a>=20
[mailto:<a href=3D"mailto:postgis-users-bounces at postgis.refractions.net" ta=
rget=3D"_blank">postgis-users-bounces at postgis.refractions.net</a>] <b>On Be=
half Of=20
</b>Ricardo Bayley<br><b>Sent:</b> Wednesday, August 18, 2010 12:06=20
AM<br><b>To:</b> <a href=3D"mailto:postgis-users at postgis.refractions.net" t=
arget=3D"_blank">postgis-users at postgis.refractions.net</a><br><b>Subject:</=
b>=20
[postgis-users] GEOS pointonsurface() threw an error!<br></font><br></div><=
div><div></div><div class=3D"h5">
<div></div>Hi,
<div><br></div>
<div>I am calculation the PointOnSurface for a table of polygons. I get thi=
s=20
error message</div>
<div><br></div>
<div><br></div>
<div>
<div>NOTICE: =C2=A0TopologyException: found non-noded intersection between=
=20
LINESTRING (857703 1.04022e+006, 857648 1.04019e+006) and LINESTRING (85753=
6=20
1.04019e+006, 857870 1.04019e+006) at 857648 1.04019e+006</div>
<div><br></div>
<div><br></div>
<div>ERROR: =C2=A0GEOS pointonsurface() threw an error!</div>
<div><br></div>
<div>********** Error **********</div>
<div><br></div>
<div>ERROR: GEOS pointonsurface() threw an error!</div>
<div>Estado SQL:XX000</div></div>
<div><br></div>
<div>I checked the one of the rows that is getting this error with=20
ST_GeometryType and I get "ST_Multipolygon", although I dont get =
the ST_AsText,=20
or any output type.</div>
<div>Any ideas why ?=C2=A0</div>
<div><br></div>
<div><br></div>
<div>regards</div>
<div><br></div>
<div><br></div>
<div>Ricardo</div>
<div><br></div>
<div><br></div></div></div></div>
<br>_______________________________________________<br>
postgis-users mailing list<br>
<a href=3D"mailto:postgis-users at postgis.refractions.net">postgis-users at post=
gis.refractions.net</a><br>
<a href=3D"http://postgis.refractions.net/mailman/listinfo/postgis-users" t=
arget=3D"_blank">http://postgis.refractions.net/mailman/listinfo/postgis-us=
ers</a><br>
<br></blockquote></div><br></div>
--0016364d2b6387c154048e11fa2c--