[geos-devel] 3.1.0rc1

Fri Jan 30 00:11:20 EST 2009

One more: Running the old aggregate union also fails,

aggtest=# select geometrytype(st_union_old(the_geom)) from tm_world;
NOTICE:  TopologyException: found non-noded intersection between
-89.142 17.8189, -89.142 17.8189 and -89.142 17.8189, -89.142 17.8189
-89.142 17.8189
ERROR:  GEOS union() threw an error!

But with a topology error in a different location. But, like the new
union, you can't re-produce the error by running standard union on the
pair of geometries at that location.

P.

On Thu, Jan 29, 2009 at 1:07 PM, Sean Gillies <sgillies at frii.com> wrote:
> Another data point: I found Dane's shapefile to contain 100 polygons and 145
> multipolygons. I was able to cascade union a) the first collection of 100
> polygons, and also b) the collection of 3238 polygons exploded from the 145
> multis. An attempt to cascade union everything (a + b) failed, as you're
> finding.
>
> http://lists.gispython.org/pipermail/community/2009-January/001963.html
>
> Sean
>
> On Jan 29, 2009, at 1:39 PM, Martin Davis wrote:
>
>> I would have expected so - so I'm puzzled.  No further ideas at the
>> moment, I'm afraid.
>>
>> Paul Ramsey wrote:
>>>
>>> I can do an individual union of Mozambique and Zimbabwe without
>>> problems. Should I be seeing a topology exception?
>>>
>>> P
>>>
>>> On Thu, Jan 29, 2009 at 12:09 PM, Martin Davis <mbdavis at refractions.net>
>>> wrote:
>>>
>>>> When I look at the area of the reported TopologyException in OpenJUMP,
>>>> I'm
>>>> seeing the following kinds of vertex values:
>>>>
>>>> Feature ID: 122
>>>> [   3:0:666] POINT (32.476230999999984 -21.321853999999973)
>>>> Feature ID: 222
>>>> [     0:144] POINT (32.476231 -21.321854)
>>>>
>>>> So the dataset is not actually 100% cleanly noded.  These kinds of small
>>>> discrepancies are exactly the kind of situation which used to cause
>>>> robustness errors in JTS.  This has pretty much been fixed by the
>>>> addition
>>>> of the Geometry Snapping concept.  I think this made it into GEOS - but
>>>> perhaps there's a bug in it?
>>>>
>>>> It should be possible to narrow down the problem by finding the two
>>>> geometries which containing the problematic points and just working with
>>>> them.  (Mozambique and Zimbabwe)
>>>>
>>>>
>>>> Paul Ramsey wrote:
>>>>
>>>>> Incidentally, the problem is not the invalid geometries. I first tried
>>>>> the union after "fixing" them with buffer(0), but I just now tried the
>>>>> union after deleting them, with the same result.
>>>>>
>>>>> aggtest=# delete from tm_world_2 where not st_isvalid(the_geom);
>>>>> NOTICE:  Ring Self-intersection at or near point -53.7564 48.5033
>>>>> NOTICE:  Ring Self-intersection at or near point -70.9172 -54.7086
>>>>> NOTICE:  Ring Self-intersection at or near point 5.33694 61.5928
>>>>> NOTICE:  Ring Self-intersection at or near point 143.662 49.3122
>>>>> DELETE 4
>>>>> aggtest=# select geometrytype(st_union(the_geom)) from tm_world_2;
>>>>> NOTICE:  TopologyException: found non-noded intersection between
>>>>> 32.4122 -21.2178, 32.4762 -21.3219 and 32.4762 -21.3219, 32.4122
>>>>> -21.2178 32.4756 -21.3208
>>>>> ERROR:  GEOS2POSTGIS returned an error
>>>>> aggtest=#
>>>>>
>>>>> P
>>>>>
>>>>> On Wed, Jan 28, 2009 at 1:49 PM, Paul Ramsey <pramsey at opengeo.org>
>>>>> wrote:
>>>>>
>>>>>
>>>>>> For anyone that needs a test geometry, here's Dane's input as a single
>>>>>> multipolygon:
>>>>>>
>>>>>> http://cleverelephant.dyndns.org/~pramsey/bigmulti.txt.bz2
>>>>>>
>>>>>> On Wed, Jan 28, 2009 at 7:44 AM, Hartmut Kaiser
>>>>>> <hartmut.kaiser at gmail.com> wrote:
>>>>>>
>>>>>>
>>>>>>>> I just did a bit further testing of the new cascaded union goodness
>>>>>>>> via Shapely. See the linked email below:
>>>>>>>>
>>>>>>>> http://lists.gispython.org/pipermail/community/2009-January/001960.html
>>>>>>>>
>>>>>>>> Seems that I reached where a point that 'should never be reached' ?
>>>>>>>>
>>>>>>>> http://trac.osgeo.org/geos/browser/trunk/source/index/strtree/AbstractS
>>>>>>>> TRtree.cpp#L358
>>>>>>>>
>>>>>>>> I am new to looking at the geos code, so any pointers to help me
>>>>>>>> test
>>>>>>>> more would be great.
>>>>>>>>
>>>>>>>>
>>>>>>> Could you provide a minimal test I could use to reproduce this
>>>>>>> problem?
>>>>>>> Regards Hartmut
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>> Thanks!
>>>>>>>>
>>>>>>>> Dane
>>>>>>>>
>>>>>>>>
>>>>>>>> On Jan 27, 2009, at 9:02 PM, Dane Springmeyer wrote:
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>> Paul,
>>>>>>>>>
>>>>>>>>> I tested on mac 10.5 and ubuntu 8.10 and all configure;make and
>>>>>>>>> make
>>>>>>>>> installs went perfectly.
>>>>>>>>>
>>>>>>>>> Just a few compiler warnings that I've pasted here:
>>>>>>>>>
>>>>>>>>>
>>>>>>>> http://dpaste.com/hold/113736/
>>>>>>>>
>>>>>>>>
>>>>>>>>> Dane
>>>>>>>>>
>>>>>>>>> On Jan 27, 2009, at 10:35 AM, Paul Ramsey wrote:
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>> This is a pre-warning to PSC members to download and test the
>>>>>>>>>> first
>>>>>>>>>> 3.1 rc tarball:
>>>>>>>>>>
>>>>>>>>>> http://download.osgeo.org/geos/geos-3.1.0rc1.tar.bz2
>>>>>>>>>>
>>>>>>>>>> If I don't hear any howls of dismay over the next 24 hours, I will
>>>>>>>>>> make a wider announcement for testing.
>>>>>>>>>>
>>>>>>>>>> Paul
>>>>>>>>>> _______________________________________________
>>>>>>>>>> geos-devel mailing list
>>>>>>>>>> geos-devel at lists.osgeo.org
>>>>>>>>>> http://lists.osgeo.org/mailman/listinfo/geos-devel
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> geos-devel mailing list
>>>>>>>>> geos-devel at lists.osgeo.org
>>>>>>>>> http://lists.osgeo.org/mailman/listinfo/geos-devel
>>>>>>>>>
>>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> geos-devel mailing list
>>>>>>>> geos-devel at lists.osgeo.org
>>>>>>>> http://lists.osgeo.org/mailman/listinfo/geos-devel
>>>>>>>>
>>>>>>>>
>>>>>>> _______________________________________________
>>>>>>> geos-devel mailing list
>>>>>>> geos-devel at lists.osgeo.org
>>>>>>> http://lists.osgeo.org/mailman/listinfo/geos-devel
>>>>>>>
>>>>>>>
>>>>>>>
>>>>> _______________________________________________
>>>>> geos-devel mailing list
>>>>> geos-devel at lists.osgeo.org
>>>>> http://lists.osgeo.org/mailman/listinfo/geos-devel
>>>>>
>>>>>
>>>>>
>>>> --
>>>> Martin Davis
>>>> Senior Technical Architect
>>>> Refractions Research, Inc.
>>>> (250) 383-3022
>>>>
>>>> _______________________________________________
>>>> geos-devel mailing list
>>>> geos-devel at lists.osgeo.org
>>>> http://lists.osgeo.org/mailman/listinfo/geos-devel
>>>>
>>>>
>>> _______________________________________________
>>> geos-devel mailing list
>>> geos-devel at lists.osgeo.org
>>> http://lists.osgeo.org/mailman/listinfo/geos-devel
>>>
>>>
>>
>> --
>> Martin Davis
>> Senior Technical Architect
>> Refractions Research, Inc.
>> (250) 383-3022
>>
>> _______________________________________________
>> geos-devel mailing list
>> geos-devel at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/geos-devel
>
> _______________________________________________
> geos-devel mailing list
> geos-devel at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/geos-devel
>