G'day all,
I have some bad geometries (Self-intersections) that I need to fix:
when I run :
australia=# select gid, tailtag, area_name,
st_astext(st_centroid(the_geom)) as centroid from tailtagarea where
st_isValid(the_geom) is false order by st_centroid(the_geom);;
to make myself a nice list to work with I get :
NOTICE: Self-intersection at or near point 146.646 -41.1121
NOTICE: Self-intersection at or near point 147.787 -42.288
NOTICE: Self-intersection at or near point 115.999 -32.7807
NOTICE: Self-intersection at or near point 136.002 -33.3412
NOTICE: Self-intersection at or near point 151.385 -27.7379
NOTICE: Self-intersection at or near point 142.163 -26.2079
NOTICE: Self-intersection at or near point 151.385 -27.7379
NOTICE: Self-intersection at or near point 142.163 -26.2079
NOTICE: Self-intersection at or near point 115.999 -32.7807
NOTICE: Self-intersection at or near point 130.528 -12.7086
NOTICE: Self-intersection at or near point 115.556 -21.4113
NOTICE: Self-intersection at or near point 113.653 -24.8917
NOTICE: Self-intersection at or near point 115.142 -30.6844
then :
gid | tailtag | area_name | centroid
-----+---------+-----------------
+-------------------------------------------
276 | WCN | Carnarvon | POINT(114.382000002796
-24.3219058101861)
291 | WDN | Dandaragan | POINT(115.504407697781
-30.5410130360183)
380 | WWA | Waroona | POINT(115.91620760142
-32.8588263644047)
343 | WMY | Murray | POINT(115.966965752586
-32.6237045619179)
362 | WRO | Roebourne | POINT(117.03170469097
-21.0526004441269)
60 | TDG | Darwin And Gulf | POINT(133.897365573887
-14.2140172548334)
247 | S71 | S71 | POINT(136.913166015303
-32.9332738397356)
70 | QBC | Barcoo | POINT(142.269579925576
-25.1377197205728)
170 | QQP | Quilpie | POINT(143.889867295296
-26.354593975264)
498 | MLA | Latrobe | POINT(146.544098538635
-41.2417387155427)
510 | MRS | Ross | POINT(147.556140609695
-42.1118612308958)
147 | QMM | Millmerran | POINT(151.11832129932
-27.8551116855016)
169 | QPW | Pittsworth | POINT(151.538916320118
-27.7014266930915)
(13 rows)
It would be excellent if I could get the Notice to display in the
query results, so I could print it out and work with it to quickly
locate the problems (in QGIS)
... so can I do that? or am I best to just copy and paste the notices
Unfortunately as I write this "Safari could not open the page “http://postgis.refractions.net/documentation/
” because the server is not responding."
cheers
Ben
--
Ben Madin
REMOTE INFORMATION
t : +61 8 9192 5455
f : +61 8 9192 5535
m : 0448 887 220
Broome WA 6725
[email protected]
Out here, it pays to
know...
_______________________________________________
postgis-users mailing list
[email protected]
http://postgis.refractions.net/mailman/listinfo/postgis-users
